anchor.diff:        KK=IACT(ICON)
         DO 90 I=1,N
      90 DXNEW(I)=A(I,KK)
  -      TOT=0.0d0
  +      TOT=0.0D0
         K=N
  +      print *, '   k, nact, dxnew =', k, nact, (dxnew(i),i=1,n)
     100 IF (K .GT. NACT) THEN
             SP=0.0d0
             SPABS=0.0d0
  
positive.code:   110     SPABS=SPABS+DABS(TEMP)
            ACCA=SPABS+0.1d0*DABS(SP)
            ACCB=SPABS+0.2d0*DABS(SP)
            IF (SPABS .GE. ACCA .OR. ACCA .GE. ACCB) SP=0.0d0
            IF (TOT .EQ. 0.0d0) THEN
                TOT=SP
            ELSE
                KP=K+1
                TEMP=DSQRT(SP*SP+TOT*TOT)
  
negative.code:       SUBROUTINE TRSTLP (N,M,A,B,RHO,DX,IFULL,IACT,Z,ZDOTA,VMULTC,
       1  SDIRN,DXNEW,VMULTD,IPRINT)
        IMPLICIT DOUBLE PRECISION (A-H,O-Z)
        DOUBLE PRECISION TEMP
        DIMENSION A(N,*),B(*),DX(*),IACT(*),Z(N,*),ZDOTA(*),
       1  VMULTC(*),SDIRN(*),DXNEW(*),VMULTD(*)
  C
  