Mixed-precision arithmetic in the ENDGame dynamical core of the Unified Model, a numerical weather prediction and climate model code

The Met Office's weather and climate simulation code the Unified Model is used for both operational Numerical Weather Prediction and Climate modelling. The computational performance of the model running on parallel supercomputers is a key consideration. A Krylov sub-space solver is employed to solve the equations of the dynamical core of the model, known as ENDGame. These describe the evolution of the Earth's atmosphere. Typically, 64-bit precision is used throughout weather and climate applications. This work presents a mixed-precision implementation of the solver, the beneficial effect on run-time and the impact on solver convergence. The complex interplay of errors arising from accumulated round-off in floating-point arithmetic and other numerical effects is discussed. A careful analysis is required. The mixed-precision solver is now employed in the operational forecast to satisfy run-time constraints without compromising the accuracy of the solution.

[1]  Henk A. van der Vorst,et al.  Bi-CGSTAB: A Fast and Smoothly Converging Variant of Bi-CG for the Solution of Nonsymmetric Linear Systems , 1992, SIAM J. Sci. Comput..

[2]  T. Allen.  Comparison of iterative pressure solvers for turbulent flow over hills , 2007 .

[3]  Peter D. Düben,et al.  Single Precision in Weather Forecasting Models: An Evaluation with the IFS , 2017 .

[4]  Peter D. Düben,et al.  The use of imprecise processing to improve accuracy in weather & climate prediction , 2014, J. Comput. Phys..

[5]  Wayne Luk,et al.  On the use of programmable hardware and reduced numerical precision in earth‐system modeling , 2015, Journal of advances in modeling earth systems.

[6]  Peter R. Graves-Morris.  The Breakdowns of BiCGStab , 2004, Numerical Algorithms.

[7]  Peter D. Düben,et al.  Benchmark Tests for Numerical Weather Forecasts on Inexact Hardware , 2014 .

[8]  Hirofumi Tomita,et al.  Single Precision in the Dynamical Core of a Nonhydrostatic Global Atmospheric Model: Evaluation Using a Baroclinic Wave Test Case , 2018 .

[9]  Richard Barrett,et al.  Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods , 1994, Other Titles in Applied Mathematics.

[10]  Julien Langou,et al.  Accelerating scientific computations with mixed precision algorithms , 2008, Comput. Phys. Commun..

[11]  John Thuburn,et al.  Conservation and Linear Rossby-Mode Dispersion on the Spherical C Grid , 2004 .

[12]  Brian Golding,et al.  Unified Modeling and Prediction of Weather and Climate: A 25-Year Journey , 2012 .

[13]  Rajiv Gupta,et al.  Fence Scoping , 2014, SC14: International Conference for High Performance Computing, Networking, Storage and Analysis.

[14]  Nicholas J. Higham,et al.  A New Analysis of Iterative Refinement and Its Application to Accurate Solution of Ill-Conditioned Sparse Linear Systems , 2017, SIAM J. Sci. Comput..

[15]  Matthew Chantry,et al.  Scale-Selective Precision for Weather and Climate Forecasting , 2019, Monthly Weather Review.

[16]  M. Diamantakis,et al.  An inherently mass‐conserving semi‐implicit semi‐Lagrangian discretization of the deep‐atmosphere global non‐hydrostatic equations , 2014 .