Collide two photon packets with each other, and see what happens. Two photon groups contain packets using a top hat functions. One of them has the lower value at zero, the other is at nonzero. The third photon group packets are Gaussians. Ideally, they should just pass through each other. But that is not what moment based methods do with radiation. To run with GEAR-RT, compile SWIFT with --with-rt=GEAR_2 --with-rt-riemann-solver=GLF --with-hydro-dimension=1 --with-hydro=gizmo-mfv \ --with-riemann-solver=hllc --with-stars=GEAR --with-feedback=none --with-grackle=$GRACKLE_ROOT