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

20  20 
dvariable tmp=0.0; 
21  21 
for (int k=m.slicemin();k<=m.slicemax();k++) 
22  22 
{ 
23 
tmp+=norm2(m(k));


23 
tmp+=norm2(m(k));


24  24 
} 
25  25 
tmp=sqrt(tmp); 
26  26 
RETURN_ARRAYS_DECREMENT(); 
...  ...  
37  37 
dvariable tmp=0.0; 
38  38 
for (int k=m.slicemin();k<=m.slicemax();k++) 
39  39 
{ 
40 
tmp+=norm2(m(k));


40 
tmp+=norm2(m(k));


41  41 
} 
42  42 
RETURN_ARRAYS_DECREMENT(); 
43  43 
return tmp; 
trunk/src/linad99/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) 
trunk/src/linad99/dmat33.cpp (revision 1657)  

32  32 
{ 
33  33 
if (j<=i) 
34  34 
{ 
35 
j1=j;


36 
i1=i;


35 
j1=j;


36 
i1=i;


37  37 
} 
38  38 
else 
39  39 
{ 
40 
j1=i;


41 
i1=j;


40 
j1=i;


41 
i1=j;


42  42 
} 
43  43 
if ( (i1j1) < bw) 
44  44 
ofs << setw(w) << S(i1,j1) << " "; 
trunk/src/linad99/dmat29.cpp (revision 1657)  

59  59 
{ 
60  60 
for (int j=1;j<=m;j++) 
61  61 
{ 
62 
tmp(i,j)=v(ii++);


62 
tmp(i,j)=v(ii++);


63  63 
} 
64  64 
} 
Also available in: Unified diff