Considerable effort went into devising an expandable structure for MeshWorker, which can work with a lot of different data structures, and relieves you of all the hassle of integrating face terms on locally refined meshes, even with multigrid. I think it would be very advisable to use it.
Best, Guido On Wed, May 23, 2012 at 1:33 PM, Katharina Kormann <[email protected]> wrote: > Dear Saul, > > The MatrixFree class has its own loop structure that is not based on the > MeshWorker since it uses its own data structures and parallelization > strategy. These data structures are necessary since we need the precomputed > transformation data for each cell. > Face loops are not supported at the moment. I am using face terms for error > estimation in my research but for cartesian cells only. Including faces in a > general setting with the same efficiency as volume terms would require > considerable enhancement of the matrix free data structures. > > Best regards, > Katharina > > On 23 maj 2012, at 16.33, Saul Teukolsky wrote: > > Hi Guido, > > No, it has it's own cell_loop function that does volume terms. It is not > obvious to me how to do face terms. > > Saul > > Guido Kanschat wrote: > > I have not had a chance to check out the MatrixFree class. Does it use the > MeshWorker loop? In that case, it would be easy. If not, it should. > > Guido > > On May 23, 2012 8:14 AM, "Saul Teukolsky" <[email protected]> wrote: >> >> Hi, >> >> I am trying to solve a DG problem with the MatrixFree class. The cell_loop >> function provides a convenient way to compute volume terms. Is there a >> convenient way to loop over the faces of each cell to compute the flux >> terms? If so, a simple example of the syntax would be very helpful! >> >> Thanks, >> Saul >> _______________________________________________ >> dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii > > > > _______________________________________________ > dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii > -- Guido Kanschat, Department of Mathematics, Texas A&M University In accordance with the Texas Public Information Act all email sent to and from my Texas A&M email account can be obtained through a Public Information Request. If you do not want your correspondence public, please arrange a phone conversation with me. _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
