Hi Corbin, Vector<Number> and LiearAlgebra::Vector<Number> should work work with MatrixFree. If not let us know what (compilation) errors you get.
Alternatively, you can use LinearAlgebra::distributed::Vector<Number> and run the program with a single process. Hope that helps, Peter On Monday, 23 May 2022 at 10:29:03 UTC+2 corbin....@gmail.com wrote: > 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/fd28184a-4ec2-47f7-9911-f0008c0eb803n%40googlegroups.com.
