anchor.diff:      midx = 0.5_kreal*(self%xmin + self%xmax)
       midy = 0.5_kreal*(self%ymin + self%ymax)
       midz = 0.5_kreal*(self%zmin + self%zmax)
  -    
  +
       call bisectAll(1)%init(self%xmin, midx, self%ymin, midy, self%zmin, midz)
       call bisectAll(2)%init(midx, self%xmax, self%ymin, midy, self%zmin, midz)
       call bisectAll(3)%init(self%xmin, midx, midy, self%ymax, self%zmin, midz)
  
positive.code:     class(Box3d), intent(in) :: self
      type(Logger), intent(inout) :: aLog
      character(len=256) :: infostring
      
      call StartSection(aLog)
      write(infostring,'(6(A,G10.4),A)') 'x in [', self%xmin, ', ', self%xmax, &
          ']; y in [', self%ymin, ', ', self%ymax, ']; z in [', self%zmin, ', ', self%zmax, ']'
  
negative.code:   real(kreal), dimension(3, 6) :: face_locs
    integer(kint) :: i
    real(kreal) :: test_dist
  
    face_locs = self%faceCentroids()
    calc_min_radius = huge(1.0_kreal)
    do i=1, 6
  