<https://scicomp.stackexchange.com/posts/44208/timeline> 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 to compute the action of the operator in a vector. What is the easiest way to implement a matrix free solver with GMRES (non-symmetric matrix)? I have considered the MatrixFree <https://www.dealii.org/current/doxygen/deal.II/classMatrixFree.html> class and also the Trilinos Wrapper <https://www.dealii.org/current/doxygen/deal.II/group__TrilinosWrappers.html>. However, the MatrixFree class seems to be implemented for parallel computing the elements of the matrix on the fly (I am happy with a look-up table), and Trilinos is a bit scary for me (I can try if no other option). Am I forced to use parallel computing to get access to Matrix Free algorithms in deal ii? -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/58416b47-4e76-426e-81f5-a37dc4dae150n%40googlegroups.com.
