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