Re: [Getfem-users] ILU problem

[Yves Renard]

Dear Marco,

I tested your code with the matrix you attached and I have no error on the ilu 
computation. (This was on a 64 bit dual core opteron with Ubuntu and g++ 
4.4.3).


Yves.



On mardi 14 juin 2011, Marco Zecchi wrote:
> 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.


-- 

  Yves Renard (address@hidden)       tel : (33) 04.72.43.87.08
  Pole de Mathematiques, INSA-Lyon             fax : (33) 04.72.43.85.29
  20, rue Albert Einstein
  69621 Villeurbanne Cedex, FRANCE
  http://math.univ-lyon1.fr/~renard

---------