In solving linear systems of billions of equations resulting from reservoir simulator pressure equations for giant reservoirs organized into fine-grid models, a multi-level methodology acquires a fast and coarser grid solution to be used as a more accurate initial estimate for the fine grid problem. Solver iterations performed on the coarser grid in reduced processing time generate a better initial estimate, resulting in fewer of the more expensive iterations at the fine-grid level.