Fast spin up of Ocean biogeochemical models using matrix-free Newton-Krylov

Ocean Modelling

Volumn 23, Issue 3-4, 2008, Pages 121-129

  Khatiwala, Samar ^a  

^a Department of Earth and Environmental Sciences   (United States)

Author keywords

[No Author keywords available]

Indexed keywords

BIOGEOCHEMISTRY; MATRIX ALGEBRA; SOLUTIONS;

BIOGEOCHEMICAL MODELING; CLASS OF METHODS; COMPUTATIONAL APPROACH; EFFICIENT COMPUTATION; EQUILIBRIUM SOLUTIONS; GLOBAL MODELING; JACOBIAN; KRYLOV SUBSPACE METHODS; LARGE SYSTEMS; NEWTON KRYLOV; NEWTON-TYPE METHODS; NONLINEAR ALGEBRAIC EQUATIONS; SPIN-UP;

NONLINEAR EQUATIONS;

ALGORITHM; BIOGEOCHEMICAL CYCLE; EFFICIENCY MEASUREMENT; MARINE ENVIRONMENT; MATRIX; NONLINEARITY; OCEANOGRAPHY; PERFORMANCE ASSESSMENT;

EID: 46849115990     PISSN: 14635003     EISSN: None     Source Type: Journal    
DOI: 10.1016/j.ocemod.2008.05.002     Document Type: Article

References (43)

1. Amestoy P.R., Duff I.S., and LÉxcellent J.-Y. Multifrontal parallel distributed symmetric and unsymmetric solvers. Comput. Method. Appl. Mech. Eng. 184 (2000) 501-520
2. Aprille T., and Trick T. A computer algorithm to determine the steady-state response of nonlinear oscillators. IEEE Trans. Circuits Syst. 19 (1972) 354-360
3. Barrett R., Berry M., Chan T.F., Demmel J., Donato J., Dongarra J., Eijkhout V., Pozo R., Romine C., and der Vorst H.V. Templates for the Solution of Linear Systems: Building Blocks for Iterative Methods. second ed. (1994), Society for Industrial and Applied Mathematics, Philadelphia, PA
4. Bernsen E., Dijkstra H.A., and Wubs F.W. A method to reduce the spin-up time of ocean models. Ocean Modell. 20 (2008) 380-392
5. Bischof C.H., Roh L., and Mauer-Oats A.J. ADIC: an extensible automatic differentiation tool for ANSI-C. Softw. Pract. Exper. 27 12 (1997) 1427-1456
6. Brown P.N. A local convergence theory for combined inexact-Newton/finite-difference projection methods. SIAM J. Numer. Anal. 24 (1987) 407-434
7. Brown P.N., and Hindmarsh A.C. Matrix-free methods for stiff systems of ODE's. SIAM J. Numer. Anal. 23 (1986) 610-638
8. Brown P.N., and Saad Y. Hybrid Krylov methods for nonlinear systems of equations. SIAM J. Sci. Stat. Comput. 11 (1990) 450-481
9. Chan T.F., and Jackson K.R. Nonlinearly preconditioned Krylov subspace methods for discrete Newton algorithms. SIAM J. Sci. Stat. Comput. 5 (1984) 533-542
10. Coleman T.F., and Moré J.J. Estimation of sparse jacobian matrices and graph coloring problems. SIAM J. Numer. Anal. 20 (1983) 187-209
11. Curtis A.R., Powell M.J.D., and Reid J.K. On the estimation of sparse jacobian matrices. J. Inst. Math. Appl. 13 (1974) 117-119
12. Davidson B.D. Large-scale continuation and numerical bifurcation for partial differential equations. SIAM J. Numer. Anal. 34 (1997) 2008-2027
13. Dembo R.S., Eisenstat S.C., and Steihaug T. Inexact Newton methods. SIAM J. Numer. Anal. 19 (1982) 400-408
14. Dembo R.S., and Steihaug T. Truncated-Newton algorithms for large scale unconstrained optimization. Math. Prog. 26 (1983) 190-212
15. Demmel J.W., Eisenstat S.C., Gilbert J.R., Li X.S., and Liu J.W.H. A supernodal approach to sparse partial pivoting. SIAM J. Matrix Anal. Appl. 20 3 (1999) 720-755
16. Dennis J.E., and Schnabel R.B. Numerical Methods for Unconstrained Optimization and Nonlinear Equations. Classics in Applied Mathematics vol. 16 (1996), Society for Industrial and Applied Mathematics
17. Deuflhard P. Newton Methods for Nonlinear Problems: Affine Invariance and Adaptive Algorithms. Springer Series in Computational Mathematics vol. 35 (2004), Springer Verlag
18. Dutkiewicz, S., Follows, M., Parekh, P., 2005. Interactions of the iron and phosphorus cycles: A three-dimensional model study. Glob. Biogeochem. Cy. 19 (doi:10.1029/2004GB002342).
19. Dutkiewicz, S., Follows, M.J., Heimbach, P., Marshall, J., 2006. Controls on ocean productivity and air-sea carbon flux: An adjoint model sensitivity study. Geophys. Res. Lett. 33. doi:10.1029/2005GL024987.
20. Eisenstat S.C., and Walker H.F. Choosing the forcing terms in an inexact Newton method. SIAM J. Sci. Comput. 17 1 (1996) 16-32
21. Gear C.W., and Saad Y. Iterative solution of linear equations in ODE codes. SIAM J. Sci. Stat. Comput. 4 (1983) 583-601
22. Gebremedhin A., Manne F., and Pothen A. What color is your Jacobian? Graph coloring for computing derivatives. SIAM Rev. 47 (2005) 629-705
23. Giering R., and Kaminski T. Recipes for adjoint code generation. ACM Trans. Math. Softw. 24 (1998) 437-474
24. Hindmarsh A.C., Brown P.N., Grant K.E., Lee S.L., Serban R., Shumaker D.E., and Woodward C.S. SUNDIALS: suite of nonlinear and differential/algebraic equation solvers. ACM Trans. Math. Softw. 31 (2005) 363-396
25. Kelley C.T. Iterative Methods for Linear and Nonlinear Equations. Frontiers in applied mathematics vol. 16 (1995), Society for Industrial and Applied Mathematics
26. Khatiwala S., Visbeck M., and Cane M. Accelerated simulation of passive tracers in ocean circulation models. Ocean Modell. 32 (2005) 51-69
27. Knoll D.A., and Keyes D.E. Jacobian-free Newton-Krylov methods: a survey of approaches and applications. J. Comput. Phys. 193 (2004) 357-397
28. Levitus et al., 1998. World Ocean Database 1998. NOAA Atlas NESDIS 18, US Government Printing Office, Washington, DC.
29. Li, X., Primeau, F., 2008. A fast Newton-Krylov solver for seasonally varying global ocean biogeochemistry models. Ocean Modell., in press.
30. Li X.S., and Demmel J.W. SuperLU_DIST: a scalable distributed-memory sparse direct solver for unsymmetric linear systems. ACM Trans. Math. Softw. 29 2 (2003) 110-140 June
31. Mahowald N., Lou C., del Corral J., and Zender C. Interannual variability in atmospheric mineral aerosols from a 22-year model simulation and observational data. J. Geophys. Res. 108 (2003) 10.1029/2002JD002821
32. Merlis T., and Khatiwala S. Fast dynamical spin-up of ocean general circulation models using Newton-Krylov methods. Ocean Modell. 21 (2008) 97-105
33. Mousseau V.A., Knoll D.A., and Rider W.J. Physics-based preconditioning and the Newton-Krylov method for non-equilibrium radiation diffusion. J. Comput. Phys. 160 (2000) 743-765
34. Orr, J.C., Najjar, R.G., Sabine, C.L., Joos, F., 1999. Abiotic-HOWTO: Internal OCMIP Report. LSCE/CEA Saclay, Gif-sur-Yvette, France, 25 pp.
35. Parekh P., Follows M.J., and Boyle E.A. Decoupling of iron and phosphate in the global ocean. Global Biogeochem. Cy. 19 (2005) 10.1029/2004GB002280
36. Reisner J., Mousseau V.A., Wyszogrodzki A.A., and Knoll D.A. An efficient physics-based preconditioner for the fully implicit solution of small-scale thermally driven atmospheric flows. J. Comput. Phys. 189 (2003) 30-44
37. Saad Y. Iterative Methods for Sparse Linear Systems. Classics in Applied Mathematics vol. 16 (2003), Society for Industrial and Applied Mathematics
38. Shroff G.M., and Keller H.B. Stabilization of unstable procedures: the recursive projection method. SIAM J. Numer. Anal. 30 (1993) 1099-1120
39. Skelboe S. Computation of the periodic steady-state response of nonlinear networks by extrapolation methods. IEEE Trans. Circuits Syst. 27 (1980) 161-175
40. Telichevesky, R., Kundert, K.S., White, J.K., 1995. Efficient steady-state analysis based on matrix-free Krylov-subspace methods. In: 32nd ACM/IEEE Design Automation Conference.
41. Tylavsky D.J., and Sohie G.R.L. Generalization of the matrix inversion lemma. Proc. IEEE 74 (1986) 1050-1052
42. Wunsch C., and Heimbach P. Practical global oceanic state estimation. Physica D 230 (2007) 197-208
43. Wunsch C., and Heimbach P. How long to oceanic tracer and proxy equilibrium?. Quat. Sci. Rev. 27 (2008) 637-651