anchor.diff:          endif
           call stop_error("hankel1: zbesh error")
       endif
  -    res = cmplx(resReal(1), resImag(1), kind=dp)
  +
  +    ! take care of negative orders if needed:
  +    res = Hrotate(1, order)*(resReal(1) + i_*resImag(1))
   end function hankel1v_complex
   
   
   function hankel2n_real(order, x) result(res)
  -    implicit none
       integer, intent(in) :: order
       real(dp), intent(in) :: x
       complex(dp) :: res
  
positive.code:     res = resReal + i_*resImag
  end function dairybi
  
  end module
  
negative.code: 
  
  function besselkv_complex(order, z) result(res)
      real(dp), intent(in) :: order
      complex(dp), intent(in) :: z
      complex(dp) :: res
  
      integer, parameter :: scaling=1, length=1
      real(dp) :: resReal(length), resImag(length), callorder
      integer :: underflows, ierr
  
      ! K_{-v} = K_v; NIST Handbook (10.27.3)
      callorder = abs(order)
  
      call zbesk(real(z), aimag(z), callorder, scaling, length, resReal, resImag, underflows, ierr)
      if(ierr > 0) then
          if(ierr == 1) then
              print *, "IERR=1, INPUT ERROR - NO COMPUTATION"
  