2.5: Remove OOPS code from the outliner space, as discussed
[blender-staging.git] / extern / ode / dist / ode / fbuild / ldlt.m
1 function [L,d] = ldlt(A)
2
3 n=length(A);
4 d=zeros(n,1);
5
6 d(1) = 1/A(1,1);
7 for i=2:n
8   for j=2:i-1
9     A(i,j) = A(i,j) - A(j,1:j-1) * A(i,1:j-1)';
10   end
11   sum = 0;
12   for j=1:i-1
13     q1 = A(i,j);
14     q2 = q1 * d(j);
15     A(i,j) = q2;
16     sum = sum + q1*q2;
17   end
18   d(i) = 1/(A(i,i) - sum);
19 end
20
21 L=A;
22 for i=1:n
23   L(i,i:n)=zeros(1,n+1-i);
24   L(i,i)=1;
25 end
26 d = d .\ 1;