Impala Documentation

template<class C>

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

Definition at line 543 of file Matrix.h.

Referenced by CxMatrixTem< C >::i().

00544 {
00545     int ii = -1;
00546 
00547         int i;
00548     for (i=0;i<n;i++) {
00549         int ip=indx[i];
00550         double sum=b[ip];
00551         b[ip]=b[i];
00552         if (ii>=0)
00553             for (int j=ii;j<=i-1;j++)
00554                 sum -= a[i*n+j]*b[j];
00555         else if (sum) ii=i;
00556         b[i]=sum;
00557     }
00558     for (i=n-1;i>=0;i--) {
00559         double sum=b[i];
00560         for (int j=i+1;j<n;j++)
00561             sum -= a[i*n+j]*b[j];
00562         b[i]=sum/a[i*n+i];
00563     }
00564 }