anchor.diff:        implicit none
         integer          n, m, lwsp, liwsp, itrace, iflag, iwsp(liwsp)
         double precision t, tol, anorm
  -      complex*16       v(n), w(n), wsp(lwsp)
  +      complex(8)       v(n), w(n), wsp(lwsp)
         external         matvec
   
   *-----Purpose----------------------------------------------------------|
  
positive.code:       implicit none
        double precision t
        integer          ideg, m, ldh, lwsp, iexph, ns, iflag, ipiv(m)
        complex*16       H(ldh,m), wsp(lwsp)
  
  *-----Purpose----------------------------------------------------------|
  *
  
negative.code:       iwsp(6) = ibrkflag
        iwsp(7) = mbrkdwn
  
        wsp(1)  = mycomplex( step_min, 0d0 )
        wsp(2)  = mycomplex( step_max, 0d0 )
        wsp(3)  = mycomplex( 0.0d0, 0d0 )
        wsp(4)  = mycomplex( 0.0d0, 0d0 )
        wsp(5)  = mycomplex( x_error, 0d0 )
        wsp(6)  = mycomplex( s_error, 0d0 )
        wsp(7)  = mycomplex( tbrkdwn, 0d0 )
        wsp(8)  = mycomplex( sgn*t_now, 0d0 )
        wsp(9)  = mycomplex( hump/vnorm, 0d0 )
        wsp(10) = mycomplex( beta/vnorm, 0d0 )
        END
  *----------------------------------------------------------------------|
  *----------------------------------------------------------------------|
  