share/hompack/fortran/otputp.f

   1       SUBROUTINE OTPUTP(N,NUMPAT,CL,FACV,CLX,X,XNP1)
   2 C
   3 C OTPUTP  POSTPROCESSES THE ENDPOINTS OF THE PATHS, UNTRANSFORMING
   4 C AND UNSCALING THEM.
   5 C
   6 C ON INPUT:
   7 C
   8 C N  IS THE NUMBER OF EQUATIONS AND VARIABLES.
   9 C
  10 C NUMPAT  IS THE CURRENT PATH NUMBER.
  11 C
  12 C CL  IS THE ARRAY THAT DEFINES THE PROJECTIVE TRANSFORMATION.
  13 C
  14 C FACV  CONTAINS THE VARIABLE SCALING FACTORS.
  15 C
  16 C X  IS THE ENDPOINT OF THE PATH, POSSIBLY TRANSFORMED AND/OR SCALED
  17 C   DEPENDING ON THE  POLSYS  INPUT FLAG  IFLG1.
  18 C
  19 C CLX  IS WORKSPACE.
  20 C
  21 C ON OUTPUT:
  22 C
  23 C N, NUMPAT, CL, AND  FACV  ARE UNCHANGED.
  24 C
  25 C X  IS THE UNTRANSFORMED AND UNSCALED VERSION OF X.
  26 C
  27 C XNP1  IS THE PROJECTIVE COORDINATE "X(N+1)".  XNP1  EQUALS UNITY IF
  28 C   THE PROJECTIVE TRANSFORMATION IS NOT SPECIFIED.
  29 C
  30 C DECLARATIONS OF INPUT, WORKSPACE, AND OUTPUT:
  31       INTEGER N,NUMPAT
  32       DOUBLE PRECISION CL,FACV,CLX,X,XNP1
  33       DIMENSION CL(2,N+1),FACV(N),CLX(2,N),X(2,N),XNP1(2)
  34 C
  35 C DECLARATION OF VARIABLES
  36       INTEGER I,IERR,J,NP1
  37       DOUBLE PRECISION D1MACH,FAC,TEMP
  38       DIMENSION TEMP(2)
  39 C
  40       NP1=N+1
  41 C COMPUTE XNP1
  42       DO 1 J=1,N
  43         CALL MULP(CL(1,J),X(1,J),CLX(1,J))
  44  1    CONTINUE
  45       DO 2  I=1,2
  46         XNP1(I)=CL(I,NP1)
  47         DO 2  J=1,N
  48           XNP1(I) = XNP1(I) + CLX(I,J)
  49   2   CONTINUE
  50 C UNTRANSFORM VARIABLES
  51       DO 10  J=1,N
  52         CALL DIVP(X(1,J),XNP1,TEMP,IERR)
  53         X(1,J)=TEMP(1)
  54         X(2,J)=TEMP(2)
  55   10  CONTINUE
  56 C UNSCALE VARIABLES
  57       TEMP(1)=D1MACH(2)
  58       DO 30 J=1,N
  59         FAC=10.**FACV(J)
  60         DO 30 I=1,2
  61           IF( (ABS(X(I,J))/TEMP(1))*FAC .LT. 1.0 ) X(I,J)=FAC*X(I,J)
  62   30  CONTINUE
  63       RETURN
  64       END