AVANCES EN EL DESARROLLO DE GEOJIP-3D: UN SOFTWARE DE ALTO RENDIMIENTO COMPUTACIONAL PARA MODELADO E INVERSIÓN CONJUNTA MULTIMODAL
Este trabajo presenta los avances en el desarrollo de GeoJIP-3D, una plataforma computacional orientada al modelado e inversión conjunta de datos gravimétricos, magnéticos y (pronto) magnetotelúricos. La resolución de este problema inverso multifísico exige estabilidad matemática y eficiencia computacional. En este marco, el software optimiza la convergencia mediante el método de Gradientes Conjugados asistido por el cálculo analítico del tamaño de paso óptimo a través de derivadas de Gâteaux, trazando trayectorias eficientes y directas hacia el mínimo global.
Entre los desarrollos recientes destaca la formulación de un nuevo operador directo MT en 3D, el cual complementa los módulos de campos potenciales existentes. Con el fin de mejorar la adaptabilidad geométrica a la topografía y el costo computacional, el dominio se discretiza mediante mallas estructuradas no uniformes con refinamiento adaptativo (tipo Octree), preservando la robustez numérica del esquema de diferencias finitas en los límites del dominio. Por su parte, el acoplamiento multifísico se rige por una restricción de similitud estructural, optimizada mediante la reducción analítica del determinante de la Matriz de Gram, permitiendo vincular simultáneamente los modelos de densidad, susceptibilidad y resistividad con bajo costo computacional.
Originalmente prototipado en Fortran y acelerado con CUDA (NVIDIA), el desarrollo de la versión beta ha migrado a un núcleo modular en Python apoyado en librerías de cómputo científico de alto rendimiento. Esta transición permite generalizar el código hacia una arquitectura agnóstica de hardware, habilitando paralelización masiva tanto en GPUs NVIDIA como AMD (ROCm). Así, GeoJIP-3D pretende consolidarse como una herramienta versátil, escalable y de alto desempeño para la caracterización del subsuelo en problemas geofísicos variados.