anchor.diff:      do iz=1,dim_quad
          do ik=1,Nk
             tmp = farr(nthread,iz,ik)
  -          gtmp = gtmp + real(jacobian(nthread,iz,ik))*tmp
  +          gtmp = gtmp + (jacobian(nthread,iz,ik))*tmp
          end do
       end do
       foldingarr = gtmp
  
positive.code:     do i=1,numnodes
         nt = OMP_GET_THREAD_NUM()
         call singrate(node_coordinates(i,:),node_coordinates(i,:),i,k_wave,centres,weights,nt)
         gauss(i,4) = sum(jacobian(nt + 1,:,:))
      end do
      !$OMP END PARALLEL DO
    end function calcsing
  
negative.code: 
      f2 = cdexp((0,1)*k_wave*norm(y))*phi(x,i,hval**2)/norm(y)
    end function f2
  
    double complex function one(x,i,j)
      integer, intent(in) :: i,j
      double precision, intent(in), dimension(:) :: x
      one = 1
    end function one
  
    double complex function f(x,i,j)
      integer, intent(in) :: i,j
  