anchor.diff:  
       tmpVector(:) = 0.
   
  +    call OMP_SET_NUM_THREADS(omp_threads)
  +
  +    !$OMP PARALLEL DO &
  +    !$OMP DEFAULT(SHARED) PRIVATE(tmpVector,ind,mask,VSa)
       do i=1, numnodesi
          ind = minloc(neighbors(i,:),dim=1) - 1
          allocate (mask(ind))
  
positive.code:           end if
         end do
      end do
    end subroutine filter_neighbors
  
  end module modphi
  
negative.code:        DEALLOCATE(VSa)
         DEALLOCATE(mask)
      end do
      !$OMP END PARALLEL DO
    end subroutine normal_vector_stroke
  
    double precision function dpair(i,j)
  