I am aware that both FEniCS and deal.ii are capable of solving problems with Discontinuous Galerkin (DG) method. I would like to specifically know if any of these two softwares can cater these requirements. Other software suggestions are also welcome; I am aware of these two because they are actively developed. I am specifically interested in solving hyperbolic (wave-dominated) problems.

  1. Local DG implementation
  2. User defined numerical flux function.
  3. Access and modify nodal/modal basis function coefficients. This is required because I want to implement limiters (for shock capturing).
  4. Support for both structured and unstructured meshes.

I request FEniCS/deal.ii users to kindly answer.

  • $\begingroup$ There is the DUNE / DUNE PDELAB ecosystem, where some DG methods are already implemented, but I can't really recommend them for lack of Documentation. If I where you, I would go with whatever package has the most active users and the best documentation. That might influence your overall productivity more than minor differences in the implementations. (my 2 cents) $\endgroup$ – MPIchael Jun 25 '19 at 10:14
  • $\begingroup$ I would also take into account the type of element you want to use. Dealii supports only quads/hexes while fenics focuses on simplex elements ( triangles and tetrahedron). I have found dealii to be quite easy to use once you get past the initial confusion that always occurs with large codes. $\endgroup$ – BlaB Jun 27 '19 at 10:54

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Browse other questions tagged or ask your own question.