Hello everyone, I've been reading through some of the matrix-free tutorials---step-67 in particular. I'd like to write something similar, an RKDG solver for hyperbolic problems for which I can use explicit RK time integration.
- Does it make sense to use the MatrixFree functionality with regular Vector<Number> operands? The tutorials often make use of LinearAlgebra::distributed::Vector<Number> objects to store the vectors. - I don't want to use the latter because I'm driving the code from Python, and any use of MPI or distributed memory could be a headache in communicating data back and forth between the two languages. Additionally, I'm solving many small problems as opposed to a single large one. - Ultimately, I'm looking for a fast way to apply the elemental inverse mass matrices and face integrations, so MatrixFree functionality seemed the logical choice. Best, Corbin -- 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/09075ee8-93e7-412a-9b86-8a220785ea2fn%40googlegroups.com.
