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 ( (i1-j1) < bw)
232 232
        (*this)(i,j)=S(i1,j1);
......
332 332
    for (j=i-bw+1;j<=i-1;j++)
333 333
    {
334 334
      if (j>1)
335
      {	
335
      {
336 336
        tmp=M(i,j);
337 337
        for (k=i-bw+1;k<=j-1;k++)
338 338
        {
339
	  if (k>0 && k>j-bw)
339
          if (k>0 && k>j-bw)
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=i-bw+1;k<=i-1;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