share/colnew/fortran/appsln.f

   1 C
   2 C----------------------------------------------------------------------
   3 C                             p a r t  4
   4 C               polynomial and service routines
   5 C----------------------------------------------------------------------
   6 C
   7       SUBROUTINE APPSLN (X, Z, FSPACE, ISPACE)
   8 C
   9 C*****************************************************************
  10 C
  11 C     purpose
  12 C
  13 C           set up a standard call to  approx  to evaluate the
  14 C           approximate solution  z = z( u(x) )  at a point x
  15 C           (it has been computed by a call to  colnew ).
  16 C           the parameters needed for  approx  are retrieved
  17 C           from the work arrays  ispace  and  fspace .
  18 C
  19 C*****************************************************************
  20 C
  21       IMPLICIT REAL*8 (A-H,O-Z)
  22       DIMENSION Z(1), FSPACE(1), ISPACE(1), A(28), DUMMY(1)
  23       IS6 = ISPACE(6)
  24       IS5 = ISPACE(1) + 2
  25       IS4 = IS5 + ISPACE(4) * (ISPACE(1) + 1)
  26       I = 1
  27       CALL APPROX (I, X, Z, A, FSPACE(IS6), FSPACE(1), ISPACE(1),
  28      1             FSPACE(IS5), FSPACE(IS4), ISPACE(2), ISPACE(3),
  29      2             ISPACE(5), ISPACE(8), ISPACE(4), 2, DUMMY, 0)
  30       RETURN
  31       END