[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Getfem-users] ILU problem
From: |
Marco Zecchi |
Subject: |
[Getfem-users] ILU problem |
Date: |
Mon, 13 Jun 2011 20:52:06 -0400 |
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
- [Getfem-users] ILU problem,
Marco Zecchi <=