Supongamos que tengo dos archivos VTK, ambos en formato de cuadrícula estructurada. Las cuadrículas estructuradas son las mismas (tienen la misma lista de puntos, en el mismo orden), y hay un campo, llamado "Phi", en cada archivo VTK. Quiero crear un tercer archivo VTK, nuevamente con la misma cuadrícula estructurada, y trazar un campo que sea la diferencia entre Phi en el primer archivo VTK y Phi en el segundo archivo VTK.
Sé cómo hacer esto manualmente; Puedo analizar el texto sin formato en los dos archivos VTK, copiar los datos en matrices, restar una matriz de la otra y luego volcar los datos en el formato correcto en un nuevo archivo. ¿Hay una mejor manera de calcular esta diferencia y exportarla a VTK? Una solución en Python, o en software de visualización como VisIt o Paraview sería preferible a usar un lenguaje compilado como C ++.
El propósito de calcular esta diferencia es comparar diferentes métodos numéricos para calcular la solución de un PDE; Como estoy usando el mismo software para generar las soluciones, puedo garantizar que todos los datos, excepto el campo Phi, serán los mismos en cada archivo que genere.