When par->traceRayAlgorithm==1 the software makes use of an algorithm which finds the chain of Delaunay cells traversed by each ray then interpolates quantities of interest inside the cells to solve the RTE. I note with some chagrin however that I neglected to make this procedure thread-safe, and at least one user has had a crash during raytracing in parallel when this option was chosen.
I'll put a note in the documentation for now.
When par->traceRayAlgorithm==1 the software makes use of an algorithm which finds the chain of Delaunay cells traversed by each ray then interpolates quantities of interest inside the cells to solve the RTE. I note with some chagrin however that I neglected to make this procedure thread-safe, and at least one user has had a crash during raytracing in parallel when this option was chosen.
I'll put a note in the documentation for now.