He was formerly Professor of Numerical Analysis at Oxford. Our group are now working on parallel Newton's method on distributed memory architectures, whose kernel is based on parallel sparse Gaussian elimination with partial pivoting. We would like to release the software package to public in the near future. Now we are looking for some large testing examples from real applications. Your help is appreciated.
Any suggestions are welcome.
Each reservoir is modeled by a single ODE in time that accounts for the fluxes of carbon into and out of the reservoir, and the accumulation of carbon within the reservoir. Thus, the model is a set of seven first-order, initial value ODEs. This input is expressed as an exponential function of time. The model does not have a climate component. However, the output of the model, e.g., the carbon concentrations in the lower and upper atmospheres, could be the input to a climate model.