High-Performance computational magnetohydrodynamics with Python

We present the AGATE simulation code, a Python-based framework developed primarily for solving the magnetohydrodynamics (MHD) equations while maintaining adaptability to other equation sets. The code employs a modular, object-oriented architecture that separates interface specifications from numerical implementations, allowing users to customize numerical methods and physics models. Built on a Godunov-type finite-volume scheme, AGATE currently supports the ideal, Hall, and Chew-Goldberger-Low (CGL) MHD equations, with multiple acceleration options ranging from Numpy to GPU-enabled computation via NVIDIA CUDA. Performance testing demonstrates that our GPU implementations achieve > 40x speedups over CPU versions. Comprehensive verification through established benchmarks confirms accurate reproduction of both linear and nonlinear phenomena across different MHD regimes. This combination of modularity, performance, and extensibility makes AGATE suitable for multiple applications: from rapid prototyping to production simulations, and from numerical algorithm development to physics education.