Fortran Subroutine for the Solution of Periodic Block-Tridiagonal Systems