Hello everyone, I am currently developing a matrix-free code which has a non-analytical spatially varying coefficient. As in step-37 I am using a Table<2, VectorizedArray<number> <https://www.dealii.org/current/doxygen/deal.II/classTable.html>> to store evaluations of the coefficient at all the quadrature points for a given cell batch.
I read my coefficient in from a CSV file corresponding to the values at DOF in a structured mesh which I am using to generate a FEFieldFunction object which i call for all quadrature points to initialize the values in the matrix-free coefficient table. This approach is fairly slow but works okay for the finest level (where the parallel distributed triangulation is the same as the structured mesh). However for coarser multigrid levels I run into issues where the cell within which I am trying to evaluate points is no longer owned by the current processor. I was wondering if anyone has any suggestions as to how to get around this, or perhaps has any ideas for a faster approach? Many thanks, Tom -- 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 visit https://groups.google.com/d/msgid/dealii/3f0dd574-5f3a-4300-8803-9a80dd2f318bn%40googlegroups.com.
