anchor.diff:              t_step = t_out-t_now
               goto 300
            endif
  -         wsp(ih+(j-1)*mh+j) = COMPLEX( hj1j, 0d0 )
  +         wsp(ih+(j-1)*mh+j) = mycomplex( hj1j, 0d0 )
            call ZDSCAL( n, 1.0d0/hj1j, wsp(j1v),1 )
            j1v = j1v + n
    200  continue
  
positive.code:      .                 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
        double precision DZNRM2
  *
  *---  check restrictions on input parameters ...
  
negative.code:       endif
  *
        mx = mbrkdwn + MAX( 0,k1-2 )
        hjj = mycomplex( beta, 0d0 )
        print *,'bef zgemv 4'
        call ZGEMV( 'n', n,mx,hjj,wsp(iv),n,wsp(iphih),1,ONE,w,1 )
        print *,'aft zgemv 4'
  