Cholesky decomposition: Difference between revisions
Content added Content deleted
(+scilab) |
|||
Line 2,920: | Line 2,920: | ||
assert(math.round( result * 100000 ) * 0.00001 == math.round( test * 100000 ) * 0.00001) |
assert(math.round( result * 100000 ) * 0.00001 == math.round( test * 100000 ) * 0.00001) |
||
}</lang> |
}</lang> |
||
=={{header|Scilab}}== |
|||
The Cholesky decomposition is builtin, and an upper triangular matrix is returned, such that <math>L^TL=A</math>. |
|||
<lang scilab>a = [25 15 -5; 15 18 0; -5 0 11]; |
|||
chol(a) |
|||
ans = |
|||
5. 3. -1. |
|||
0. 3. 1. |
|||
0. 0. 3. |
|||
a = [18 22 54 42; 22 70 86 62; |
|||
54 86 174 134; 42 62 134 106]; |
|||
chol(a) |
|||
ans = |
|||
4.2426407 5.1854497 12.727922 9.8994949 |
|||
0. 6.5659052 3.0460385 1.6245539 |
|||
0. 0. 1.6497422 1.849711 |
|||
0. 0. 0. 1.3926212</lang> |
|||
=={{header|Seed7}}== |
=={{header|Seed7}}== |