stumpf.f

Go to the documentation of this file.

      SUBROUTINE stumpf(IPAIR,Z)
*
*       Modified Stumpff functions.
*       ---------------------------
*
*       Developed by Seppo Mikkola 10/1997.
*       ns is the max order included in z^n.
*
      include 'common6.h'
      parameter(ns=12)
      real*8 c(12)
*
*
*       Generate Stumpff functions for argument Z (SCOEFF set by START).
      c(ns)=1-z*scoeff(12)
      c(ns-1)=1-z*scoeff(11)
      do i=ns-2,3,-1
          c(i)=1-z*c(i+2)*scoeff(i)
      end do
      c(1) = 1.0
      c(2) = 1.0
*
*       Copy the five first coefficients (for predictor & corrector).
      DO 2 i = 1,5
          sf(i,ipair) = c(i)
*         SF(I,IPAIR) = 1.0
    2 CONTINUE
*
*       Generate Stumpff functions for argument 4*Z (only down to c_5).
      z4 = 4.0d0*z
      c(ns)=1-z4*scoeff(12)
      c(ns-1)=1-z4*scoeff(11)
      do i=ns-2,5,-1
          c(i)=1-z4*c(i+2)*scoeff(i)
      end do
*
*       Copy c_5 & c_6 to sixth & seventh location (for time integration).
      DO 3 i = 6,7
          sf(i,ipair) = c(i-1)
*         SF(I,IPAIR) = 1.0
    3 CONTINUE
*
      RETURN
*
      END