def test_trajPhotonFlooding(self):
options = FakeOptionsRK()
options.pbfile = "mixed.prb"
options.configfilename = "standardRKmulti.conf"
runTrajPrbAdd(options)
outdata = numpy.loadtxt(options.outtrajfile, comments = "#")
self.assertAlmostEqual(outdata[:, 3].sum(), 16200. / 1000, delta = 100.0 / 1000)
def test_trajClippingPhotonFlooding(self):
options = FakeOptionsRK()
options.pbfile = "mixed.prb"
options.configfilename = "standardRKclip.conf"
runTrajPrbAdd(options)
outdata = numpy.loadtxt(options.outtrajfile, comments = "#")
self.assertEqual(len(outdata[:, 3]), 3760)
