Revision 1657 trunk/src/linad99/dmat43.cpp
dmat43.cpp (revision 1657)  

82  82 
for (j=ibw+1;j<=i1;j++) 
83  83 
{ 
84  84 
if (j>1) 
85 
{


85 
{ 

86  86 
tmp=M(i,j); 
87  87 
for (k=ibw+1;k<=j1;k++) 
88  88 
{ 
89 
if (k>0 && k>jbw)


89 
if (k>0 && k>jbw)


90  90 
tmp=L(i,k)*L(j,k); 
91  91 
} 
92  92 
L(i,j)=tmp/L(j,j); 
...  ...  
95  95 
tmp=M(i,i); 
96  96 
for (k=ibw+1;k<=i1;k++) 
97  97 
{ 
98 
if (k>0)


98 
if (k>0) 

99  99 
tmp=L(i,k)*L(i,k); 
100  100 
} 
101  101 
if (tmp<=0) 
Also available in: Unified diff