Hi Praveen, Correction, I meant to say face worker instead of boundary worker. Step-74 <https://dealii.org/developer/doxygen/deal.II/step_74.html#Theassemble_systemfunction> is a DG code that uses mesh_loop() with a face worker, and also makes use of the nice FEInterfaceValues <https://dealii.org/developer/doxygen/deal.II/classFEInterfaceValues.html> class that has also very recently been extended to use FEValuesExtractors (see item 17 in the specific improvements section <https://dealii.org/developer/doxygen/deal.II/changes_between_9_2_0_and_9_3_0.html> for the upcoming release).
Best, Jean-Paul > On 31. May 2021, at 19:57, Jean-Paul Pelteret <jppelte...@gmail.com> wrote: > > 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 >> <mailto: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/ >> <http://www.dealii.org/> >> For mailing list/forum options, see >> https://groups.google.com/d/forum/dealii?hl=en >> <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 >> <mailto: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 >> >> <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/2D3471C3-6094-4BEF-AA32-7F1B7F97CDF2%40gmail.com.
