template<class C>
Definition at line 551 of file MatrixTem.h. Referenced by Impala::Core::Matrix::MatrixTem< C >::i(). 00552 { 00553 int ii = -1; 00554 00555 int i; 00556 for (i=0;i<n;i++) { 00557 int ip=indx[i]; 00558 double sum=b[ip]; 00559 b[ip]=b[i]; 00560 if (ii>=0) 00561 for (int j=ii;j<=i-1;j++) 00562 sum -= a[i*n+j]*b[j]; 00563 else if (sum) ii=i; 00564 b[i]=sum; 00565 } 00566 for (i=n-1;i>=0;i--) { 00567 double sum=b[i]; 00568 for (int j=i+1;j<n;j++) 00569 sum -= a[i*n+j]*b[j]; 00570 b[i]=sum/a[i*n+i]; 00571 } 00572 }
|