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

220  220 
{ 
221  221 
if (j<=i) 
222  222 
{ 
223 
j1=j;


224 
i1=i;


223 
j1=j;


224 
i1=i;


225  225 
} 
226  226 
else 
227  227 
{ 
228 
j1=i;


229 
i1=j;


228 
j1=i;


229 
i1=j;


230  230 
} 
231  231 
if ( (i1j1) < bw) 
232  232 
(*this)(i,j)=S(i1,j1); 
...  ...  
332  332 
for (j=ibw+1;j<=i1;j++) 
333  333 
{ 
334  334 
if (j>1) 
335 
{


335 
{ 

336  336 
tmp=M(i,j); 
337  337 
for (k=ibw+1;k<=j1;k++) 
338  338 
{ 
339 
if (k>0 && k>jbw)


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


340  340 
tmp=L(i,k)*L(j,k); 
341  341 
} 
342  342 
L(i,j)=tmp/L(j,j); 
...  ...  
345  345 
tmp=M(i,i); 
346  346 
for (k=ibw+1;k<=i1;k++) 
347  347 
{ 
348 
if (k>0)


348 
if (k>0) 

349  349 
tmp=L(i,k)*L(i,k); 
350  350 
} 
351  351 
if (tmp<=0) 
Also available in: Unified diff