anchor.diff:      call OMP_SET_NUM_THREADS(4)
   
       ptr_jacobian => fjacobian2
  -    ptr_f => f
       !$OMP PARALLEL DO &
       !$OMP DEFAULT(SHARED) PRIVATE(nt)
       do i=1,numnodes
          nt = OMP_GET_THREAD_NUM()
          call integrate(nstroke_coordinates(i,:),node_coordinates(i,:),i,centres,weights,nt)
  -       intphi_under(i) = folding(i,dim_quad,nt)
  +       intphi_under(i) = folding(i,f,dim_quad,nt)
       end do
       !$OMP END PARALLEL DO
     end subroutine calcomp2
  
positive.code:   subroutine setup_calcomp ()
      use dbsym
      ptr_jacobian => fjacobian
      ptr_f => f
    end subroutine setup_calcomp
  
    subroutine calcomp()
  
negative.code:     call OMP_SET_NUM_THREADS(4)
  
      ptr_jacobian => fjacobian
      !$OMP PARALLEL DO &
      !$OMP DEFAULT(SHARED) PRIVATE(nt)
      do i=1,numnodes
         nt = OMP_GET_THREAD_NUM()
         call integrate(nstroke_coordinates(i,:),node_coordinates(i,:),i,centres,weights,nt)
         gauss(i,6) = folding(i,f2,dim_quad,nt)
      end do
      !$OMP END PARALLEL DO
      gauss(1,6) = 0.0
  