anchor.diff:  C     Partition the working space array W to provide the storage that is needed
   C     for the main calculation.
   C
  +      ITOTAL=N*(3*N+2*M+11)+4*M+6
         MPP=M+2
         ICON=1
         ISIM=ICON+MPP
  
positive.code:          PRINT *, '  VSIG = ', (VSIG(J),J=1,N), ' -- ', PARSIG
           PRINT *, '  VETA = ', (VETA(J),J=1,N), ' -- ', PARETA
           PRINT *, '  IBRNCH, IFLAG = ', IBRNCH, IFLAG
        END IF
  C
  C     If a new vertex is needed to improve acceptability, then decide which
  
negative.code:       IF (IPRINT .EQ. 3) THEN
           PRINT *, '  **DX = ', (DX(I),I=1,N)
           PRINT *, '  **IACT = ', (IACT(I),I=1,M+1)
           PRINT *, 'M,N,RHO,IFULL =', M, N, RHO, IFULL
           PRINT *, '  **CON = ', (CON(I),I=1,M)
           PRINT *, '  **A = ', ((A(I,K),I=1,N),K=1,MP)
           PRINT *, '  **W = ', (W(I),I=1,ITOTAL)
        END IF
        IF (IFULL .EQ. 0) THEN
            TEMP=0.0d0
  