Updated test results, made fission sampling in MG mode a smidge faster

src/macroxs_header.F90

@@ -628,8 +628,8 @@ contains
     real(8) :: prob             ! Running probability
 
     xi = prn()
-    prob = ZERO
-    gout = 0
+    gout = 1
+    prob = this % chi(gout,gin)
 
     do while (prob < xi)
       gout = gout + 1
@@ -650,8 +650,8 @@ contains
     call find_angle(this % polar, this % azimuthal, uvw, iazi, ipol)
 
     xi = prn()
-    prob = ZERO
-    gout = 0
+    gout = 1
+    prob = this % chi(gout,gin,iazi,ipol)
 
     do while (prob < xi)
       gout = gout + 1

src/scattdata_header.F90

@@ -308,6 +308,8 @@ contains
           ! Now make sure integral norms to zero
           norm = this % dist(gin) % data(order,gout)
           if (norm > ZERO) then
+            this % fmu(gin) % data(:,gout) = &
+                 this % fmu(gin) % data(:,gout) / norm
             this % dist(gin) % data(:,gout) = &
                  this % dist(gin) % data(:,gout) / norm
           end if

tests/test_mg_basic/results_true.dat

@@ -1,2 +1,2 @@
 k-combined:
-1.045320E+00 5.851680E-02
+1.033731E+00 4.974463E-02

tests/test_mg_max_order/inputs_true.dat

@@ -1 +1 @@
-abe20c626d613e73ccb1a3f8468ad1b9aecca528afa9e8131a411d754eb86b8ab64a6fb1fdc9c0b8b8158ff7c82f548de5912041bf035aa5a2d4532cfe0c9510
+322483933c38fe6ecfa41d632c7214b5cd35af4a56415872585914d9c775dc99171e918eebf3221ab6292689c37269b8c3ce5ff85b3633b5c05ee481bf1b212a

tests/test_mg_max_order/results_true.dat

@@ -1,2 +1,2 @@
 k-combined:
-1.083030E+00 1.855038E-02
+1.055274E+00 1.715904E-02

tests/test_mg_max_order/test_mg_max_order.py

@@ -76,9 +76,10 @@ class MGMaxOrderTestHarness(PyAPITestHarness):
         self._input_set = MGNuclideInputSet()
 
     def _build_inputs(self):
-        super(MGMaxOrderTestHarness, self)._build_inputs()
         # Set P1 scattering
         self._input_set.settings.max_order = 1
+        # Call standard input build
+        super(MGMaxOrderTestHarness, self)._build_inputs()
 
 if __name__ == '__main__':
     harness = MGMaxOrderTestHarness('statepoint.10.*', False, mg=True)

tests/test_mg_nuclide/results_true.dat

@@ -1,2 +1,2 @@
 k-combined:
-1.380785E-01 5.556526E-03
+1.317412E-01 5.926047E-03