Implementación del método LDG para mallas no estructuradas en 3D

En este artículo se describe una implementación del método "Local Discontinuous Galerkin" (LDG) aplicado a problemas elípticos en 3D. Se discute la implementación de los principales operadores. En particular el uso de aproximaciones de alto orden y de mallas no estructuradas. Estructuras de datos eficientes que permiten un rápido ensamblado del sistema lineal en su formulación mixta son descritas en detalle.