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
