Impala Documentation

template<class C>

static void Impala::Core::Matrix::lubksb	(	C *	a,
		int	n,
		short *	indx,
		double *	b
	)			`[static]`

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 }