On Tue, 9 May 2017, Salazar De Troya, Miguel wrote: > There is a function System::project_vector() to project the vectors > defined in one mesh onto a new refined mesh, but what about > matrices? The same projection operators could be used on that > matrix. I can imagine that we could recycle the preconditioner onto > the new mesh without having to rebuild it again. Is this possible in > libmesh?
We've been futzing with the new GenericProjector code to allow users to extract projection matrices. I assume the recycled preconditioner you'd want, given an original preconditioner P and a projection matrix M, would just be M P M^T? If so, then you might want to help us get the projection matrix code working. I've got incomplete code in the roystgnr/projection_matrix branch; the matrix assignment code is written but I didn't do preallocation yet. Paul Bauman and one of his students may or may not have taken that further in their multigrid work. --- Roy ------------------------------------------------------------------------------ Check out the vibrant tech community on one of the world's most engaging tech sites, Slashdot.org! http://sdm.link/slashdot _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
