On 7/19/21 1:50 PM, Corbin Foucart wrote:
* It seems that since cell interior data (A inverse) is needed to
premultiply the B operator (which contains face terms), it makes more
sense to loop over each cell and handle the faces of the cell manually,
like in step-51, which performs a similar element-local static condensation.
* In the spirit of step-51, my idea is to use mesh_loop, but only use the
cell-assembly part and pass dummy functions to the face_worker and
boundary_worker arguments. Is that a sensible way to do it, or should I be
using a different MeshWorker function?
* Or alternatively, is there a way to use the face_worker function, but
preserve the data on the cell to the left and right? It seems expensive to
recompute the inverse mass matrix on every interface, for each adjacent
cell.
All of these seem reasonable approaches. Unless you have concrete evidence
that whatever performance you get, going with the easiest-to-implement
approach is always the right choice!
Best
Wolfgang
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
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/09bbfcf4-5885-1ee9-6a21-db3ffe855884%40colostate.edu.
