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:           DO 80 I=1,N
     80     OPTNEW=OPTNEW-DX(I)*A(I,MCON)
        END IF
        IF (ICOUNT .EQ. 0 .OR. OPTNEW .LT. OPTOLD) THEN
            OPTOLD=OPTNEW
            NACTX=NACT
  
negative.code:               TEMP=ALPHA*Z(I,K)+BETA*Z(I,KP)
                Z(I,KP)=ALPHA*Z(I,KP)-BETA*Z(I,K)
    120         Z(I,K)=TEMP
                print *, '   k, alpha, beta, tot', k, alpha, beta, tot
            END IF
            K=K-1
            GOTO 100
  