anchor.diff:  implicit none
   
   ! Arguments
  -integer, intent(in)             :: n, np
  -double precision, intent(inout) :: d(np), v(np,np)
  +integer          :: n, np
  +double precision :: d(np), v(np,np)
   
   ! local variables
   integer :: i, j, k
  
positive.code: implicit none
  
  ! Arguments
  integer, intent(in)          :: nvalues
  integer, intent(out)         :: ierr
  double precision, intent(in) :: values(nvalues)
  double precision             :: mean
  
  ! Local variables
  integer :: i
  
negative.code: 
  subroutine cross_product(vec1,vec2,cross)
  implicit none
  double precision :: vec1(3), vec2(3)
  double precision :: cross(3)
  cross(1) = vec1(2)*vec2(3) - vec1(3)*vec2(2)
  cross(2) = vec1(3)*vec2(1) - vec1(1)*vec2(3)
  cross(3) = vec1(1)*vec2(2) - vec1(2)*vec2(1)
  