Merge pull request #17 from wdhawkins/spherical

wdhawkins · web-flow · commit 4898221e5bea · 2025-06-18T22:55:47.000-05:00
Fixes for boundary conditions after product quadrature normalization changes
diff --git a/Six_1g_spherical_benchmarks/Problem_1.py b/Six_1g_spherical_benchmarks/Problem_1.py
@@ -77,7 +77,7 @@
     xs_mat.CreateSimpleOneGroup(sigma_t=sigt, c=c)
 
     # Boundary source
-    bsrc = [1.0]
+    bsrc = [4.0 * np.pi]
 
     # Angular quadrature
     pquad = GLCProductQuadrature3DXYZ(n_polar=2, n_azimuthal=4)
diff --git a/Six_1g_spherical_benchmarks/Problem_2.py b/Six_1g_spherical_benchmarks/Problem_2.py
@@ -77,7 +77,7 @@
     xs_mat.CreateSimpleOneGroup(sigma_t=sigt, c=c)
 
     # Boundary source
-    bsrc = [1.0]
+    bsrc = [4.0 * np.pi]
 
     # Angular quadrature
     pquad = GLCProductQuadrature3DXYZ(n_polar=2, n_azimuthal=4)
diff --git a/Six_1g_spherical_benchmarks/Problem_4.py b/Six_1g_spherical_benchmarks/Problem_4.py
@@ -40,7 +40,7 @@
     xs_mat.CreateSimpleOneGroup(sigma_t=1.0, c=0.0)
 
     # Boundary source
-    bsrc = [1.0]
+    bsrc = [4.0 * np.pi]
 
     # Angular quadrature
     pquad = GLCProductQuadrature3DXYZ(n_polar=2, n_azimuthal=4)
diff --git a/Six_1g_spherical_benchmarks/Problem_5.py b/Six_1g_spherical_benchmarks/Problem_5.py
@@ -40,7 +40,7 @@
     xs_mat.CreateSimpleOneGroup(sigma_t=1.0, c=0.0)
 
     # Boundary source
-    bsrc = [1.0]
+    bsrc = [4.0 * np.pi]
 
     # Angular quadrature
     pquad = GLCProductQuadrature3DXYZ(n_polar=2, n_azimuthal=4)
