Hi Praveen, > Is it possible to use WorkStream to loop over faces ?
It is, in the sense that you can use it to multithread the loop over all cells, and then loop over all all faces on a cell yourself. This has been automated in the mesh_loop() <https://dealii.org/developer/doxygen/deal.II/group__MeshWorker.html#ga76ec61fbd188fb320fe8ca166a79b322> function (in fact, mesh_loop() is used as a part of the backend to MeshWorker). There you can specify a boundary worker that will be called either before all after the cells — you get to choose using the MeshWorker::AssembleFlags <https://dealii.org/developer/doxygen/deal.II/namespaceMeshWorker.html#ac7a9db8b34d398d7d398d1e8809874aaa44a76e905b1d4cd80af387b5fac4d8aa> . It’s quite a convenient assembly device, and its being put to use in most of the new tutorials. There are a few tutorials that use mesh_loop() already. Just look in the keywords of the tutorial list <https://dealii.org/developer/doxygen/deal.II/Tutorial.html> to find them. I’m sure that there must also be a DG example in the test suite somewhere, if you need it. > I know MeshWorker can be used for this, but it seems to be not used by many > people and I suspect it may not be maintained well. Yes, it seems not to be used terribly often (at least, the questions on the mailing list are infrequent, and the number of people that can answer those questions seems somewhat limited). I advocate using mesh_loop() in preference to MeshWorker. Best, Jean-Paul > On 31. May 2021, at 06:21, Praveen C <cprav...@gmail.com> wrote: > > Dear all > > I have a DG-type scheme which is totally quadrature-free and matrix-free. I > need to first loop over cells and then loop over faces. I know MeshWorker can > be used for this, but it seems to be not used by many people and I suspect it > may not be maintained well. > > Is it possible to use WorkStream to loop over faces ? > > Thanks > praveen > > -- > 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/98B64908-C102-4C4F-B18E-99066E32C1A4%40gmail.com. -- 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/A5AA339A-E94A-44AD-B098-6035CC37033A%40gmail.com.
