Hi Katharina,
Thanks for your response. I don't need the full efficiency for the face
terms for my application (which is to test various deal.ii features
against another code, basically). Would you be able to share with me
what you have for doing face terms so far, perhaps offline?
Thanks,
Saul
Katharina Kormann 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
