/* umfpack-simple.c
 *
 * Demo file extracted from UMFPACK's documentation.
 *
 * Compule with:
 *    gcc -Wall -ansi -pedantic -O2 umfpack-simple.c -lumfpack -lamd -lm
*/
#include <stdio.h>
#include "umfpack/umfpack.h"

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    n=5;

int main(void)
{
    void *Symbolic, *Numeric ;
    int i ;

    umfpack_di_symbolic(n, n, Ap, Ai, Ax, &Symbolic, NULL, NULL);
    umfpack_di_numeric(Ap, Ai, Ax, Symbolic, &Numeric, NULL, NULL);
    umfpack_di_free_symbolic(&Symbolic);
    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;
}