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; }