3D discontinuous Galerkin method

Dear all,
Can we write 3D discontinuous Galerkin method FreeFem++ code?

Yes, see from the source:
LapDG3d.edp (1.7 KB)

Thank you very much.