Stephen On Wednesday, July 10, 2019 at 3:56:02 PM UTC-4, Stephen wrote: > > I have an error estimator which I need to calculate between edges and I'd > like to loop over all edges in the triangulation. I know that I can do this > the "traditional" way via looping over all cells and then looping over all > faces in the cell but, if I do this, I end up calculating the same thing > twice since I then run over non-boundary faces once from one cell and then > again from it's neighbour; this is not acceptable since efficiency is > really of paramount importance in the code I'm currently writing. Is it > possible to instead loop over all active edges in the triangulation and > then find the cell IDs of the cells which neighbour this face? This way I > should be able to calculate the terms directly on the edge instead of > trying to go via the cell first. Thanks! > What you can do is to check that the face is not on the boundary and in that case only compute your error on the left and the bottom edges (for example). I think you could also use this function <https://www.dealii.org/current/doxygen/deal.II/group__MeshWorker.html#gab75fe399736977323bb319959ab66feb> to do all the heavy lifting for you.
Best, Bruno -- 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/a9f98552-5ca3-4bf6-a54e-01fed8f46ade%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.
