[Getfem-users] GMM++ Preconditioner

[Danesh Daroui]

Dear all,

I would like to use GMM++ iterative solvers. The problem is that I can
not use for example incomplete LU factorization preconditioner. Let say
that I have the equation Ax=b. The I define the preconditioner as:

gmm::dense_matrix<double> A(100, 100);

// filling the matrix

gmm::ilu_precond<double> P1(A);

but I get compilation errors. It just compiled successfullt when I define:

gmm::ilu_precond<double> P1;

but there is no hint and neither sample code to show which function should I 
call to create the preconditioner.

Can anybody help?

Regards,

D.