[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: [Getfem-users] ILU problem
From: |
Yves Renard |
Subject: |
Re: [Getfem-users] ILU problem |
Date: |
Wed, 15 Jun 2011 13:18:44 +0200 |
User-agent: |
KMail/1.13.5 (Linux/2.6.32-32-server; KDE/4.4.5; x86_64; ; ) |
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
---------