776 	const int	nCnt)  in PVRTMatrixLinearEqSolveF()  argument
785 	_RPT1(_CRT_WARN, "LinearEqSolve(%d)\n", nCnt);  in PVRTMatrixLinearEqSolveF()
786 	for(i = 0; i < nCnt; ++i)  in PVRTMatrixLinearEqSolveF()
789 		for(j = 1; j <= nCnt; ++j)  in PVRTMatrixLinearEqSolveF()
795 	if(nCnt == 1)  in PVRTMatrixLinearEqSolveF()
803 	i = nCnt;  in PVRTMatrixLinearEqSolveF()
808 		if(pSrc[i][nCnt] != 0)  in PVRTMatrixLinearEqSolveF()
811 			if(i != (nCnt-1))  in PVRTMatrixLinearEqSolveF()
813 				for(j = 0; j <= nCnt; ++j)  in PVRTMatrixLinearEqSolveF()
816 					f = pSrc[nCnt-1][j];  in PVRTMatrixLinearEqSolveF()
817 					pSrc[nCnt-1][j] = pSrc[i][j];  in PVRTMatrixLinearEqSolveF()
823 			for(j = 0; j < (nCnt-1); ++j)  in PVRTMatrixLinearEqSolveF()
825 				_ASSERT(pSrc[nCnt-1][nCnt] != 0);  in PVRTMatrixLinearEqSolveF()
826 				f = pSrc[j][nCnt] / pSrc[nCnt-1][nCnt];  in PVRTMatrixLinearEqSolveF()
829 				for(k = 0; k < nCnt; ++k)  in PVRTMatrixLinearEqSolveF()
831 					pSrc[j][k] -= f * pSrc[nCnt-1][k];  in PVRTMatrixLinearEqSolveF()
840 	PVRTMatrixLinearEqSolveF(pRes, pSrc, nCnt - 1);  in PVRTMatrixLinearEqSolveF()
843 	f = pSrc[nCnt-1][0];  in PVRTMatrixLinearEqSolveF()
844 	for(k = 1; k < nCnt; ++k)  in PVRTMatrixLinearEqSolveF()
846 		f -= pSrc[nCnt-1][k] * pRes[k-1];  in PVRTMatrixLinearEqSolveF()
848 	_ASSERT(pSrc[nCnt-1][nCnt] != 0);  in PVRTMatrixLinearEqSolveF()
849 	f /= pSrc[nCnt-1][nCnt];  in PVRTMatrixLinearEqSolveF()
850 	pRes[nCnt-1] = f;  in PVRTMatrixLinearEqSolveF()
860 		for(i = 1; i <= nCnt; ++i)  in PVRTMatrixLinearEqSolveF()
861 			fCnt += pSrc[nCnt-1][i] * pRes[i-1];  in PVRTMatrixLinearEqSolveF()
863 		_ASSERT(abs(fCnt - pSrc[nCnt-1][0]) < 1e-3);  in PVRTMatrixLinearEqSolveF()