Dear all, I am developing a quasi Newton solver in deal.ii. During each iteration, I need to solve a linear system B * x = p, where the matrix B has the following form: B = B_0 + W * M * W^T. B_0 is an n by n sparse matrix (n is large > 100k); M is a m by m *full matrix* (m is small m = 20 for example); W is a n by m full matrix; Since W * M * W^T will be an n by n fully dense matrix, I cannot afford to explicitly store it.
Does the deal.ii community have any suggestion what is the most efficient way to solve this linear system? I am thinking to use LinearOperator to define and operate on W * M * W^T, but I am not sure whether it supports full matrix and how efficient it would be. Best, Tao -- 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 dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/b7aed966-09b5-4486-ba68-3d74183fc203n%40googlegroups.com.