minter1d_rat.F File Reference

#include "my_real.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	minter1d_rat (x0, x1, x2, x3, y0, y1, y2, y3, x, y, yp, i, n)

Function/Subroutine Documentation

minter1d_rat()

subroutine minter1d_rat	(		x0,
			x1,
			x2,
			x3,
			y0,
			y1,
			y2,
			y3,
			x,
			y,
			yp,
		integer	i,
		integer	n )

Definition at line 32 of file minter1d_rat.F.

C     INTERPOLATION (RATIONAL FUNCTION)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
        USE constant_mod , ONLY : zero, em20, one
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
        IMPLICIT NONE
#include      "my_real.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER I, N
        my_real x0, x1, x2, x3, y0, y1, y2, y3, x, y, yp
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        my_real q, d, r, s, sp, c2, dm, sm, c1, c6, c4, c3, c3d, c5
C-----------------------------------------------
        q = x-x1
        d = x2-x1
        r = d-q
        s = (y2-y1) / d
        sp = (y3-y2) / (x3-x2)
        c2 = (sp-s) / (x3-x1)
        dm = x1-x0
        dm = sign(max(em20,abs(dm)),dm)
        sm = (y1-y0) / dm
        c1 = (s-sm) / (d+dm)
        c6 = zero
        IF(i == 1)THEN
          IF(s*(s-d*c2) <= zero) c2 = s/d
          c4 = c2
        ELSE IF(i == n-1)THEN
          c4 = c1
        ELSE
          IF(i == 2 .AND. sm*(sm-dm*c1) <= zero) c1=(s-sm-sm) / d
          c3 = abs(c2*r)
          c3d = c3+abs(c1*q)
          c5 = zero
          IF(c3d > zero) THEN
            c3 = c3/c3d
            c5 = c3*(c1-c2)
          ENDIF
          c4 = c2+c5
          c6 = d*c5*(one-c3)
        ENDIF
        y = y1 + q*(s-r*c4)
        yp = s+(q-r)*c4+c6
        RETURN