Main Page
Related Pages
Modules
Namespaces
Classes
Files
File List
File Members
external_packages
csparse
cs_lsolve.c
Go to the documentation of this file.
1
#include "
cs.h
"
2
/* solve Lx=b where x and b are dense. x=b on input, solution on output. */
3
int
cs_lsolve
(
const
cs
*
L
,
double
*x)
4
{
5
int
p, j, n, *Lp, *Li ;
6
double
*Lx ;
7
if
(!
CS_CSC
(L) || !x)
return
(0) ;
/* check inputs */
8
n = L->
n
; Lp = L->
p
; Li = L->
i
; Lx = L->
x
;
9
for
(j = 0 ; j < n ; j++)
10
{
11
x [j] /= Lx [Lp [j]] ;
12
for
(p = Lp [j]+1 ; p < Lp [j+1] ; p++)
13
{
14
x [Li [p]] -= Lx [p] * x [j] ;
15
}
16
}
17
return
(1) ;
18
}
cs_sparse::n
int n
Definition:
cs.h:20
cs.h
cs_lsolve
int cs_lsolve(const cs *L, double *x)
Definition:
cs_lsolve.c:3
CS_CSC
#define CS_CSC(A)
Definition:
cs.h:140
cs_sparse::p
int * p
Definition:
cs.h:21
cs_sparse
Definition:
cs.h:16
cs_sparse::i
int * i
Definition:
cs.h:22
L
#define L
cs_sparse::x
double * x
Definition:
cs.h:23
acado
Author(s): Milan Vukov, Rien Quirynen
autogenerated on Mon Jun 10 2019 12:34:31