1
<?xml version="1.0" encoding="ISO-8859-1" standalone="no"?>
2
<!DOCTYPE MAN SYSTEM "../../manrev.dtd">
4
<LANGUAGE>eng</LANGUAGE>
6
<TYPE>Scilab Function</TYPE>
7
<DATE>April 1993</DATE>
8
<SHORT_DESCRIPTION name="spchol"> sparse cholesky factorization</SHORT_DESCRIPTION>
10
<CALLING_SEQUENCE_ITEM>[R,P] = spchol(X) </CALLING_SEQUENCE_ITEM>
15
<PARAM_NAME>X</PARAM_NAME>
17
<SP>: symmetric positive definite real sparse matrix</SP>
21
<PARAM_NAME>P</PARAM_NAME>
23
<SP>: permutation matrix</SP>
27
<PARAM_NAME>R</PARAM_NAME>
29
<SP>: cholesky factor</SP>
35
<P><VERB>[R,P] = spchol(X)</VERB> produces a
36
lower triangular matrix <VERB>R</VERB> such that <VERB>P*R*R'*P' = X</VERB>.</P>
41
3., 0., 0., 2., 0., 0., 2., 0., 2., 0., 0. ;
42
0., 5., 4., 0., 0., 0., 0., 0., 0., 0., 0. ;
43
0., 4., 5., 0., 0., 0., 0., 0., 0., 0., 0. ;
44
2., 0., 0., 3., 0., 0., 2., 0., 2., 0., 0. ;
45
0., 0., 0., 0. , 5., 0., 0., 0., 0., 0., 4. ;
46
0., 0., 0., 0., 0., 4., 0., 3., 0., 3., 0. ;
47
2., 0., 0., 2., 0., 0., 3., 0., 2., 0., 0. ;
48
0., 0., 0., 0., 0., 3., 0., 4., 0., 3., 0. ;
49
2., 0., 0., 2., 0., 0., 2., 0., 3., 0., 0. ;
50
0., 0., 0., 0., 0., 3., 0., 3., 0., 4., 0. ;
51
0., 0., 0., 0., 4., 0., 0., 0., 0., 0., 5.];
52
X=sparse(X);[R,P] = spchol(X);