|
5164
|
1 /* -------------------------------------------------------------------------- */ |
|
|
2 /* UMFPACK Version 4.4, Copyright (c) 2005 by Timothy A. Davis. CISE Dept, */ |
|
|
3 /* Univ. of Florida. All Rights Reserved. See ../Doc/License for License. */ |
|
|
4 /* web: http://www.cise.ufl.edu/research/sparse/umfpack */ |
|
|
5 /* -------------------------------------------------------------------------- */ |
|
|
6 |
|
|
7 #include <stdio.h> |
|
|
8 #include "umfpack.h" |
|
|
9 |
|
|
10 int n = 5 ; |
|
|
11 int Ap [ ] = {0, 2, 5, 9, 10, 12} ; |
|
|
12 int Ai [ ] = { 0, 1, 0, 2, 4, 1, 2, 3, 4, 2, 1, 4} ; |
|
|
13 double Ax [ ] = {2., 3., 3., -1., 4., 4., -3., 1., 2., 2., 6., 1.} ; |
|
|
14 double b [ ] = {8., 45., -3., 3., 19.} ; |
|
|
15 double x [5] ; |
|
|
16 |
|
|
17 int main (void) |
|
|
18 { |
|
|
19 double *null = (double *) NULL ; |
|
|
20 int i ; |
|
|
21 void *Symbolic, *Numeric ; |
|
|
22 (void) umfpack_di_symbolic (n, n, Ap, Ai, Ax, &Symbolic, null, null) ; |
|
|
23 (void) umfpack_di_numeric (Ap, Ai, Ax, Symbolic, &Numeric, null, null) ; |
|
|
24 umfpack_di_free_symbolic (&Symbolic) ; |
|
|
25 (void) umfpack_di_solve (UMFPACK_A, Ap, Ai, Ax, x, b, Numeric, null, null) ; |
|
|
26 umfpack_di_free_numeric (&Numeric) ; |
|
|
27 for (i = 0 ; i < n ; i++) printf ("x [%d] = %g\n", i, x [i]) ; |
|
|
28 return (0) ; |
|
|
29 } |
|
|
30 |
