//============================================================================
//
// Listing 16.5
//
// Mathematics for 3D Game Programming and Computer Graphics, 3rd ed.
// By Eric Lengyel
//
// The code in this file may be freely used in any software. It is provided
// as-is, with no warranty of any kind.
//
//============================================================================


void LURefineSolution(int n, const float *m, const float *d,
		const unsigned short *index, const float *r, float *x)
{
	float *t = new float[n];
	
	for (int i = 0; i < n; i++)
	{
		double q = -r[i];
		for (int k = 0; k < n; k++) q += m[k * n + i] * x[k];
		t[i] = (float) q;
	}
	
	LUBacksubstitute(n, d, index, t, t);
	for (int i = 0; i < n; i++) x[i] -= t[i];
	
	delete[] t;
}