Pronto implementaré el método nodal discontinuo de Galerkin, y habiendo hecho esto antes de conocer las matrices de indexación básicas que tendré que calcular, dada una malla y datos polinomiales.
El problema que encontré en el código anterior fueron errores sutiles que cometí al calcular cosas como la indexación de trazas interior / exterior. Los problemas que no surgieron en casos de prueba más simples surgirían en mallas más grandes, y generalmente esto produce un esquema inestable ya que las condiciones de contorno no se imponen correctamente (por lo que no hay posibilidad de simplemente ver la simulación cada 10 pasos más o menos y ver un problema localizado) .
Espero que algunas personas con más experiencia conozcan buenas pruebas para ejecutar en las matrices de índices para tener la confianza de que tienen razón. Las cuadraturas, derivados y similares son muy fáciles de probar, pero no puedo entender otras cosas.
Algunas pruebas que he hecho en el pasado están agregando normales interiores a normales exteriores, lo que debería producir 0 o + -2. Ser útil poder ver rápidamente el resultado de algún cambio de código es útil, pero no puedo pensar en una manera significativa de hacerlo con la indexación.
También debo mencionar que estos están yendo para quads y hexes, con el potencial de elementos curvilíneos. Ya no hay mucho código existente para comparar con una biblioteca en funcionamiento.
Puntos de bonificación si hay buenas pruebas unitarias que puedo escribir que no se basarían en una respuesta correcta existente para comparar. Me conformaré con muchas buenas heurísticas.