double mxm_solve(double *x, const double *A, const double *b, const int N)
{
mxm_local_block(a2, double, N);
mxm_set(a2, A, N);
mxv_set(x, b, N);
double det = internal_solve(a2, x, N);
mxm_free_local(a2);
return det;
}
real mxm_solve(real *x, const real *A, const real *b, const int N)
{
mxm_local_block(a2, real, N);
mxm_set(a2, A, N);
mxv_set(x, b, N);
real det = internal_solve(a2, x, N);
mxm_free_local(a2);
return det;
}
