Code used to benchmark the 3D diffusion code by calculating the criticality condition.