MAXDGTD is a discontinuous Galerkin code for solving Maxwell's equations in the time-domain.

The discontinuous Galerkin code uses a combination of high-order spatial and various temporal discretizations methods

  1. Spatial discretizations
    1. Centered numerical fluxes
    2. High-order nodal polynomial interpolation
    3. High-order geometrical mapping for curvilinear domains
    4. Efficient cubature formulas for line, surface and volume integrals

  2. Temporal discretizations
    1. Arbitrarily high-order explicit leap-frog integrators
    2. Second-order implicit Crank-Nicolson scheme
    3. Second-order hybrid explicit/implicit scheme
    4. Exponential, one-step, explicit, polynomial integrators

MAXDGTD can handle

  • Non-conforming simplicial and hexahedral meshes with arbitrary level hanging nodes
  • Non-conforming hybrid tetrahedral-hexahedral meshes
  • Locally varying polynomial degrees
  • Complex geometries with curved boundaries
  • Strongly heterogeneous media, conductive materials, source terms

The code is written in Fortran 77/90 (around 50K SLOCs).

Publications related to the development of this code: [J1] - [J2] - [J3] - [J4] - [J6] - [J7] - [J8] - [C2] - [C3], see also my PhD thesis.

If you are interested in this code for your own research, please contact me at

Hassan FAHS
LAMA Laboratory - UMR 8050
Université Paris-Est - Marne-la-Vallée
5, boulevard Descartes
Cité Descartes - Champs-sur-Marne
77454 Marne-la-Vallée cedex 2 - France

or by email