Prevent segfault in distance to boundary calculation (#2659) #3185

Job	Run time
Python 3.10 (omp=n, mpi=n, dagmc=n, ncrystal=n, libmesh=n, event=n vectfit=n)	1h 0m 27s
Python 3.10 (omp=y, mpi=n, dagmc=n, ncrystal=n, libmesh=n, event=n vectfit=n)	2h 43m 56s
Python 3.10 (omp=n, mpi=y, dagmc=n, ncrystal=n, libmesh=n, event=n vectfit=n)	45m 59s
Python 3.10 (omp=y, mpi=y, dagmc=n, ncrystal=n, libmesh=n, event=n vectfit=n)	2h 26m 59s
Python 3.7 (omp=n, mpi=n, dagmc=, ncrystal=, libmesh=, event= vectfit=)	1h 30m 13s
Python 3.8 (omp=n, mpi=n, dagmc=, ncrystal=, libmesh=, event= vectfit=)	1h 0m 22s
Python 3.9 (omp=n, mpi=n, dagmc=, ncrystal=, libmesh=, event= vectfit=)	1h 0m 7s
Python 3.11 (omp=n, mpi=n, dagmc=, ncrystal=, libmesh=, event= vectfit=)	1h 25m 21s
Python 3.10 (omp=y, mpi=y, dagmc=y, ncrystal=, libmesh=, event= vectfit=)	3h 0m 28s
Python 3.10 (omp=n, mpi=n, dagmc=, ncrystal=y, libmesh=, event= vectfit=)	1h 18m 39s
Python 3.10 (omp=y, mpi=y, dagmc=, ncrystal=, libmesh=y, event= vectfit=)	2h 20m 4s
Python 3.10 (omp=y, mpi=n, dagmc=, ncrystal=, libmesh=y, event= vectfit=)	3h 3m 24s
Python 3.10 (omp=y, mpi=n, dagmc=, ncrystal=, libmesh=, event=y vectfit=)	3h 7m 4s
Python 3.10 (omp=n, mpi=y, dagmc=, ncrystal=, libmesh=, event= vectfit=y)	1h 3m 23s
finish	5s
	1d 1h 46m 31s

Provide feedback