class LinearSolver {
public:
	void Init(UINT n);
	void LoadDiagTerms();
	void CGradSolve();
	void BiCGradSolve(double _BMag = 0.0);
	__forceinline void PushElement(UINT Row, UINT Col, double Entry)
	{
		M.PushElement(Row, Col, Entry);
	}

	Vector<double> p, p2, q, q2, r, r2, z, z2, Temp, x, b, DiagTerms;
	SparseMatrix M;
};