A very simple computational model of pollution in the air over the U.S. might try to compute the concentration of three pollutants at each point in a 3000 x 1500 x 25 grid.

From physics and chemistry a mathematical model of the interaction of these three pollutants can be derived. This mathematical model is a system of 3 time-dependent, nonlinear, partial differential equations.

A **numerical method** for this problem would compute values of each of the three unknowns so that the
mathematical equations are approximately satisfied at each of the grid points.

So the total number of unknowns is

The classic algorithm requires ~~O~~ (*n*^{3}) flops, where *n* is the number of unknowns. Hence,
we need more than 10^{25} flops to solve the problem. On a teraflop machine, this would take at least

Unfortunately, since the problem is time-dependent, we actually need to solve one of the above problems for each of hundreds or thousands of time steps!

**flops**: floating point operations

CS1104 Main Page

Last Updated 04/16/2000

© L.Heath, 2000