real_sq_solve()

void linalg::real_sq_solve	(	double *	ad,
		double *	bd,
		double *	xd,
		int	dim
	)

Solve a linear equation ax=b.

Purpose

Solve a linear equation ax=b for a given real square matrx ( \(a\)) and a real vector ( \(b\)).

Parameters

ad	A real square matrix given as a pointer.
bd	A real vector given as a pointer.
xd	A real vector into which the solution is set. Given as a pointer.
dim	Size of the matrix and the vector

Definition at line 106 of file linalg_eigen.cpp.

References get_dvector(), set_dmatrix(), and set_dvector().

{
   dmatrix A_eigen;
   dvector b_eigen,X_eigen;   
   set_dmatrix(n,a,A_eigen);
   set_dvector(n,b,b_eigen);
   Eigen::FullPivLU<dmatrix> lu(A_eigen);
   X_eigen=lu.solve(b_eigen);
   get_dvector(n,x,X_eigen);
}

Here is the call graph for this function: