Re: [deal.II] Matrix Free alternatives in dealii

Dear Felipe,Gmres (and all other iterative methods) in deal.Ii is already matrix free. It never accesses the matrix , only its matrix vector method.  FMM has already been implemented in deal.II and it is available in the repository pi-BEM. Take a look there for inspiration. If you know how to compu

[deal.II] Matrix Free alternatives in dealii

I am implementing a Fast Multipole Method (FMM) in deal ii. I do not want to store a dense matrix, but lower rank matrices and to use matrix free methods. By now, I store the elements of the low-rank matrices, and I want to retrieve them