Mercurial > octave-nkf

comparison liboctave/UMFPACK/UMFPACK/Demo/umfpack_simple.c @ 5164:57077d0ddc8e

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 2005-02-25 19:55:24 by jwe]
author jwe
date Fri, 25 Feb 2005 19:55:28 +0000
parents
children
comparison
equal deleted inserted replaced
5163:9f3299378193 5164:57077d0ddc8e
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