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