Go to the documentation of this file.00001 #include "cs.h"
00002
00003 int cs_lusol (int order, const cs *A, double *b, double tol)
00004 {
00005 double *x ;
00006 css *S ;
00007 csn *N ;
00008 int n, ok ;
00009 if (!CS_CSC (A) || !b) return (0) ;
00010 n = A->n ;
00011 S = cs_sqr (order, A, 0) ;
00012 N = cs_lu (A, S, tol) ;
00013 x = cs_malloc (n, sizeof (double)) ;
00014 ok = (S && N && x) ;
00015 if (ok)
00016 {
00017 cs_ipvec (N->pinv, b, x, n) ;
00018 cs_lsolve (N->L, x) ;
00019 cs_usolve (N->U, x) ;
00020 cs_ipvec (S->q, x, b, n) ;
00021 }
00022 cs_free (x) ;
00023 cs_sfree (S) ;
00024 cs_nfree (N) ;
00025 return (ok) ;
00026 }