anchor.diff:      implicit none
   
       double precision, dimension(:,:), intent(in) :: q
  -    double precision, dimension(:), intent(in) :: n
  +    integer, dimension(:), intent(in) :: n
       double precision, dimension(:), intent(in) :: y
       double precision, intent(in) :: lambda
       double precision, intent(in) :: sigma
  
positive.code:     implicit none
  
      double precision, dimension(:,:), intent(in) :: D2
      double precision, dimension(:), intent(in) :: N2
      double precision, dimension(:), intent(in) :: alpha
      double precision, intent(in) :: dgamma
  
  
negative.code:     implicit none
  
      double precision, dimension(:,:), intent(in) :: d
      integer, dimension(:), intent(in) :: n
      double precision, dimension(:), intent(in) :: y
      double precision, intent(in) :: dgamma
      double precision, intent(in) :: lambda
  