anchor.diff:  pure subroutine initSpecific(self, xmin, xmax, ymin, ymax, zmin, zmax)
     class(Box3d), intent(inout) :: self
     real(kreal), intent(in) :: xmin, xmax, ymin, ymax, zmin, zmax
  -  
  +
     self%xmin = xmin
     self%xmax = xmax
     self%ymin = ymin
  
positive.code:     real(kreal), dimension(3,8) :: corners
      integer(kint) :: i
      real(kreal) :: test_dist
      
      corners = self%corners()
      calc_max_radius = 0.0_kreal
      do i=1,8
  
negative.code:     real(kreal), dimension(3), intent(in) :: queryPt
      containsPoint = (self%xmin <= queryPt(1) .and. queryPt(1) <= self%xmax) .and. &
                      (self%ymin <= queryPt(2) .and. queryPt(2) <= self%ymax) .and. &
                      (self%zmin <= queryPt(3) .and. queryPt(3) <= self%zmax)
  end function
  
  pure function bisectAll(self)
      type(Box3d), dimension(8) :: bisectAll
  