anchor.diff:       .                 s_error, x_error, t_now, t_new, t_step, t_old,
        .                 xm, beta, break_tol, p1, p2, p3, eps, rndoff,
        .                 vnorm, avnorm, hj1j, hump, SQR1
  -      complex*16 hij
  -
  -      intrinsic AINT,ABS,COMPLEX,DBLE,INT,LOG10,MAX,MIN,NINT,SIGN,SQRT
  -      !complex*16 ZDOTC
  +      complex(8) hij
  +
  +      intrinsic AINT,ABS,DBLE,INT,LOG10,MAX,MIN,NINT,SIGN,SQRT
  +      complex(8) mycomplex
         double precision DZNRM2
   *
   *---  check restrictions on input parameters ...
  
positive.code:       iwsp(6) = ibrkflag
        iwsp(7) = mbrkdwn
  
        wsp(1)  = COMPLEX( step_min, 0d0 )
        wsp(2)  = COMPLEX( step_max, 0d0 )
        wsp(3)  = COMPLEX( 0.0d0, 0d0 )
        wsp(4)  = COMPLEX( 0.0d0, 0d0 )
        wsp(5)  = COMPLEX( x_error, 0d0 )
        wsp(6)  = COMPLEX( s_error, 0d0 )
        wsp(7)  = COMPLEX( tbrkdwn, 0d0 )
        wsp(8)  = COMPLEX( sgn*t_now, 0d0 )
        END
  *----------------------------------------------------------------------|
  *----------------------------------------------------------------------|
  
negative.code:             t_step = t_out-t_now
              goto 300
           endif
           wsp(ih+(j-1)*mh+j) = mycomplex( hj1j, 0d0 )
           wsp(ih+j*mh+j-1) = mycomplex( hj1j, 0d0 )
           call ZDSCAL( n, 1.0d0/hj1j, wsp(j1v),1 )
           j1v = j1v + n
   200  continue
  