ADI for the Convection Equation

Program adiprcon.f ssolves the two-dimensional convection equation using a Peaceman-Rachford ADI scheme. You also need the adiprcon.in input file, the u_exa.f exact solution, the tridia.f tridiagonal solver, the subroutine out.f taking care of output, the subroutine in.f taking care of input, and the openw.f file open subroutine.

It should not be difficult to convert the program to one for the heat equation. Actually, the heat equation is easier, since it has boundary conditions at all four boundaries, not just two of them.


Return