[Getfem-users] ILU problem

[Marco Zecchi]

Hello,

I'm having problem with the gmm++ ilu preconditioner.
I'm attaching the matrix for which I'm trying to get the ilu factorization.

This is the code:

int main()
{
   
    gmm::col_matrix< gmm::wsvector<double> > M;
    MatrixMarket_load("./M.txt", M);
   
    gmm::csr_matrix<double> A;
    gmm::copy(M,A);
    cout << A.ncols() << "\t" << A.nrows() << "\t" << gmm::nnz(M) << endl;
    gmm::ilu_precond<gmm::csr_matrix<double>> P(A);

return 0;
}

I get a "vector subscript out of range" error...

However, in Matlab I am able to fine the ILU(0) factorization with...

load 'M.txt'
S = spconvert(M);
setup.type = 'nofill';
P = ilu(S,setup)

Is this a bug, or I am doing something wrong?

Thank you very much for your help!

Best regards

Marco.

M.zip
Description: Zip archive