A linear equation with tridiagonal matrix structure is solved here. The main diagonal is stored on bu, the upper diagonal on au, and the lower diagonal on cu, the right hand side is stored on du. The method used here is the simplified Gauss elimination, also called Thomas algorithm.
IMPLICIT NONEINPUT PARAMETERS:
integer, intent(in) :: N,fi,ltOUTPUT PARAMETERS:
REALTYPE :: value(0:N)REVISION HISTORY:
Original author(s): Hans Burchard & Karsten Bolding