def test_05_finite_anisotropic(self):
ur1, ut1, uru, utu = iad.rt(1.0, 1.0, 0.8, 1.0, 0.8)
self.assertAlmostEqual(ur1, 0.03041, delta=0.0001)
self.assertAlmostEqual(ut1, 0.76388, delta=0.0001)
self.assertAlmostEqual(uru, 0.08416, delta=0.0001)
self.assertAlmostEqual(utu, 0.61111, delta=0.0001)
def test_04_finite(self):
ur1, ut1, uru, utu = iad.rt(1.0, 1.0, 0.8, 1.0, 0.0)
self.assertAlmostEqual(ur1, 0.21085, delta=0.0001)
self.assertAlmostEqual(ut1, 0.54140, delta=0.0001)
self.assertAlmostEqual(uru, 0.28015, delta=0.0001)
self.assertAlmostEqual(utu, 0.41624, delta=0.0001)
def test_02_thick(self):
ur1, ut1, uru, utu = iad.rt(1.0, 1.0, 0.8, 100000.0, 0.0)
self.assertAlmostEqual(ur1, 0.28525, delta=0.0001)
self.assertAlmostEqual(ut1, 0.00000, delta=0.0001)
self.assertAlmostEqual(uru, 0.34187, delta=0.0001)
self.assertAlmostEqual(utu, 0.00000, delta=0.0001)
def test_03_thick_non_absorbing(self):
ur1, ut1, uru, utu = iad.rt(1.0, 1.0, 1.0, 100000.0, 0.0)
self.assertAlmostEqual(ur1, 1.0000, delta=0.0001)
self.assertAlmostEqual(ut1, 0.0000, delta=0.0001)
self.assertAlmostEqual(uru, 1.0000, delta=0.0001)
self.assertAlmostEqual(utu, 0.0000, delta=0.0001)
def test_01_thick_non_scattering(self):
ur1, ut1, uru, utu = iad.rt(1.0, 1.0, 0.0, 100000.0, 0.0)
self.assertAlmostEqual(ur1, 0.00000, delta=0.0001)
self.assertAlmostEqual(ut1, 0.00000, delta=0.0001)
self.assertAlmostEqual(uru, 0.00000, delta=0.0001)
self.assertAlmostEqual(utu, 0.00000, delta=0.0001)