The numerical treatment of stiff ODE systems is carried out by using implicit methods. A long sequence of non-linear systems has to be treated when implicit methods are used. The Newton iterative method is often used in the solution of these systems. This leads to the calculation of Jacobian...

It is hard to achieve high efficiency when large sparse problems are solved on parallel computers. There are three major reasons for this:• one must use indirect addressing in all sparse matrix techniques,• the data involved in computations are irregularly distributed, which means that there...

The horizontal advection is one of the most important physical processes in an air pollution model. While it is clear how to describe mathematically this process, the computer treatment of the arising first-order partial differential equation (PDE) causes great difficulties. It is assumed that...