Beispiel #1
0
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;
}
Beispiel #2
0
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;
}
 bool solve()
 {
   internal_solve();
   return true;
 }
Beispiel #4
0
void MUMPSLU::solve(sparse_matrix const& A, vector& x, vector const& b)
{
    this->allocate_coordinate_format_storage(A);
    internal_solve(A, x, b);
}