• Introduces an efficient iterative solver (GMRES with incomplete LU preconditioning) for 3D finite-difference gravitational field modeling, overcoming memory and time limitations of direct solvers.
• Employs central finite-difference discretization on right rectangular prismatic grids and high-degree Lagrange interpolation for accurate partial derivatives of gravitational potential.
• Validates the algorithm using three density models, demonstrating high accuracy, reliability, and flexibility for 3D forward modeling.
• Provides a practical solution for large-scale 3D gravity data inversion and geological interpretation.