Mercurial > octave-nkf
diff liboctave/UMFPACK/UMFPACK/Demo/umfpack_simple.c @ 5164:57077d0ddc8e
[project @ 2005-02-25 19:55:24 by jwe]
author | jwe |
---|---|
date | Fri, 25 Feb 2005 19:55:28 +0000 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/liboctave/UMFPACK/UMFPACK/Demo/umfpack_simple.c Fri Feb 25 19:55:28 2005 +0000 @@ -0,0 +1,30 @@ +/* -------------------------------------------------------------------------- */ +/* UMFPACK Version 4.4, Copyright (c) 2005 by Timothy A. Davis. CISE Dept, */ +/* Univ. of Florida. All Rights Reserved. See ../Doc/License for License. */ +/* web: http://www.cise.ufl.edu/research/sparse/umfpack */ +/* -------------------------------------------------------------------------- */ + +#include <stdio.h> +#include "umfpack.h" + +int n = 5 ; +int Ap [ ] = {0, 2, 5, 9, 10, 12} ; +int Ai [ ] = { 0, 1, 0, 2, 4, 1, 2, 3, 4, 2, 1, 4} ; +double Ax [ ] = {2., 3., 3., -1., 4., 4., -3., 1., 2., 2., 6., 1.} ; +double b [ ] = {8., 45., -3., 3., 19.} ; +double x [5] ; + +int main (void) +{ + double *null = (double *) NULL ; + int i ; + void *Symbolic, *Numeric ; + (void) umfpack_di_symbolic (n, n, Ap, Ai, Ax, &Symbolic, null, null) ; + (void) umfpack_di_numeric (Ap, Ai, Ax, Symbolic, &Numeric, null, null) ; + umfpack_di_free_symbolic (&Symbolic) ; + (void) umfpack_di_solve (UMFPACK_A, Ap, Ai, Ax, x, b, Numeric, null, null) ; + umfpack_di_free_numeric (&Numeric) ; + for (i = 0 ; i < n ; i++) printf ("x [%d] = %g\n", i, x [i]) ; + return (0) ; +} +