Descenso del bloque de vértices

Presentamos el descenso por bloques de vértices, una solución de descenso por coordenadas de bloques para la forma variacional de Euler implícito mediante iteraciones de Gauss-Seidel a nivel de vértices. Funciona con actualizaciones locales de la posición de los vértices que logran reducciones en la energía variacional global con un paralelismo maximizado. Esto conforma un solucionador físico capaz de alcanzar la convergencia numérica con estabilidad incondicional y un rendimiento computacional excepcional. También puede adaptarse a un presupuesto de cálculo dado simplemente limitando el número de iteraciones, al tiempo que mantiene su estabilidad y su tasa de convergencia superior. 

Presentamos y evaluamos nuestro método en el contexto de la dinámica de cuerpos elásticos, proporcionando detalles de todos los componentes esenciales y demostrando que supera a las técnicas alternativas. Además, discutimos y mostramos ejemplos de cómo nuestro método puede utilizarse para otros sistemas de simulación, incluyendo simulaciones basadas en partículas y cuerpos rígidos.