diff --git a/tests/regression_tests/mgxs_library_mesh/inputs_true.dat b/tests/regression_tests/mgxs_library_mesh/inputs_true.dat index 4f138b24f6..f0f93d43c4 100644 --- a/tests/regression_tests/mgxs_library_mesh/inputs_true.dat +++ b/tests/regression_tests/mgxs_library_mesh/inputs_true.dat @@ -1,311 +1,35 @@ - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 1.26 1.26 - 17 17 - -10.71 -10.71 - -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1 -1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -1 1 2 1 1 2 1 1 2 1 1 2 1 1 2 1 1 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -1 1 1 2 1 1 1 1 1 1 1 1 1 2 1 1 1 -1 1 1 1 1 2 1 1 2 1 1 2 1 1 1 1 1 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 -1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 - - - 1.26 1.26 - 17 17 - -10.71 -10.71 - -3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 -3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 -3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3 -3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3 -3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 -3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3 -3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 -3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 -3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3 -3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 -3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 -3 3 4 3 3 4 3 3 4 3 3 4 3 3 4 3 3 -3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 -3 3 3 4 3 3 3 3 3 3 3 3 3 4 3 3 3 -3 3 3 3 3 4 3 3 4 3 3 4 3 3 3 3 3 -3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 -3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 - - - 21.42 21.42 - 21 21 - -224.91 -224.91 - -5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 -5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 -5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5 -5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 -5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 -5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 -5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 -5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 -5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 -5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 -5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 -5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 -5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 -5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 -5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 -5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 -5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 -5 5 5 5 5 6 6 6 6 6 6 6 6 6 6 6 5 5 5 5 5 -5 5 5 5 5 5 5 6 6 6 6 6 6 6 5 5 5 5 5 5 5 -5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 -5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 5 - - - 21.42 21.42 - 21 21 - -224.91 -224.91 - -7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 -7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 -7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7 -7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7 -7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 -7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 -7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 -7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 -7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 -7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 -7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 -7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 -7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 -7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 -7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 -7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 -7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 -7 7 7 7 7 8 8 8 8 8 8 8 8 8 8 8 7 7 7 7 7 -7 7 7 7 7 7 7 8 8 8 8 8 8 8 7 7 7 7 7 7 7 -7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 -7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 7 - - - - - - - - - - - - - - - - - - + + + + + + + + + + + + - - - - - - - + + + - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + eigenvalue - 100 - 10 - 5 - - - -160 -160 -183 160 160 183 - - + 1000 + 5 + 0 diff --git a/tests/regression_tests/mgxs_library_mesh/results_true.dat b/tests/regression_tests/mgxs_library_mesh/results_true.dat index 829edde9f4..f1ff29d6ce 100644 --- a/tests/regression_tests/mgxs_library_mesh/results_true.dat +++ b/tests/regression_tests/mgxs_library_mesh/results_true.dat @@ -1,310 +1,310 @@ mesh 1 group in nuclide mean std. dev. x y z -0 1 1 1 1 total 0.762544 0.085298 -2 1 2 1 1 total 0.644837 0.088457 -1 2 1 1 1 total 0.653375 0.153317 -3 2 2 1 1 total 0.676480 0.094215 +0 1 1 1 1 total 0.105390 0.006421 +2 1 2 1 1 total 0.105466 0.003175 +1 2 1 1 1 total 0.106221 0.004040 +3 2 2 1 1 total 0.102641 0.002129 mesh 1 group in nuclide mean std. dev. x y z -0 1 1 1 1 total 0.473988 0.088732 -2 1 2 1 1 total 0.399254 0.091318 -1 2 1 1 1 total 0.379821 0.167092 -3 2 2 1 1 total 0.424265 0.099551 +0 1 1 1 1 total 0.078603 0.006888 +2 1 2 1 1 total 0.075950 0.003755 +1 2 1 1 1 total 0.074519 0.004589 +3 2 2 1 1 total 0.072616 0.002838 mesh 1 group in nuclide mean std. dev. x y z -0 1 1 1 1 total 0.473988 0.088732 -2 1 2 1 1 total 0.399254 0.091318 -1 2 1 1 1 total 0.379821 0.167092 -3 2 2 1 1 total 0.424265 0.099551 +0 1 1 1 1 total 0.078605 0.006892 +2 1 2 1 1 total 0.075989 0.003746 +1 2 1 1 1 total 0.074571 0.004600 +3 2 2 1 1 total 0.072586 0.002824 mesh 1 group in nuclide mean std. dev. x y z -0 1 1 1 1 total 0.027288 0.005813 -2 1 2 1 1 total 0.020262 0.003701 -1 2 1 1 1 total 0.019449 0.004420 -3 2 2 1 1 total 0.021266 0.002869 +0 1 1 1 1 total 0.013600 0.000926 +2 1 2 1 1 total 0.013584 0.000551 +1 2 1 1 1 total 0.013692 0.000712 +3 2 2 1 1 total 0.013022 0.000430 mesh 1 group in nuclide mean std. dev. x y z -0 1 1 1 1 total 0.016037 0.006339 -2 1 2 1 1 total 0.013018 0.003521 -1 2 1 1 1 total 0.012153 0.003804 -3 2 2 1 1 total 0.012965 0.002454 +0 1 1 1 1 total 0.001333 0.001105 +2 1 2 1 1 total 0.001339 0.000693 +1 2 1 1 1 total 0.001330 0.000885 +3 2 2 1 1 total 0.001260 0.000534 mesh 1 group in nuclide mean std. dev. x y z -0 1 1 1 1 total 0.011251 0.003050 -2 1 2 1 1 total 0.007243 0.001219 -1 2 1 1 1 total 0.007296 0.001795 -3 2 2 1 1 total 0.008301 0.001066 +0 1 1 1 1 total 0.012266 0.000830 +2 1 2 1 1 total 0.012244 0.000486 +1 2 1 1 1 total 0.012361 0.000650 +3 2 2 1 1 total 0.011762 0.000376 mesh 1 group in nuclide mean std. dev. x y z -0 1 1 1 1 total 0.027498 0.007445 -2 1 2 1 1 total 0.017954 0.003077 -1 2 1 1 1 total 0.017912 0.004426 -3 2 2 1 1 total 0.020469 0.002617 +0 1 1 1 1 total 0.032001 0.002161 +2 1 2 1 1 total 0.031882 0.001271 +1 2 1 1 1 total 0.032193 0.001701 +3 2 2 1 1 total 0.030726 0.001000 mesh 1 group in nuclide mean std. dev. x y z -0 1 1 1 1 total 2.177345e+06 589804.301157 -2 1 2 1 1 total 1.404096e+06 236476.852674 -1 2 1 1 1 total 1.413154e+06 347806.623478 -3 2 2 1 1 total 1.608259e+06 206502.707123 +0 1 1 1 1 total 2.372379e+06 160440.303797 +2 1 2 1 1 total 2.368109e+06 93914.371991 +1 2 1 1 1 total 2.390701e+06 125743.883417 +3 2 2 1 1 total 2.274785e+06 72785.094827 mesh 1 group in nuclide mean std. dev. x y z -0 1 1 1 1 total 0.735256 0.080216 -2 1 2 1 1 total 0.624575 0.084974 -1 2 1 1 1 total 0.633925 0.149098 -3 2 2 1 1 total 0.655214 0.091422 +0 1 1 1 1 total 0.091790 0.005503 +2 1 2 1 1 total 0.091883 0.002653 +1 2 1 1 1 total 0.092530 0.003354 +3 2 2 1 1 total 0.089619 0.001721 mesh 1 group in nuclide mean std. dev. x y z -0 1 1 1 1 total 0.763779 0.070696 -2 1 2 1 1 total 0.628158 0.064356 -1 2 1 1 1 total 0.640809 0.158369 -3 2 2 1 1 total 0.645171 0.080467 +0 1 1 1 1 total 0.087817 0.005624 +2 1 2 1 1 total 0.090790 0.005246 +1 2 1 1 1 total 0.093736 0.005609 +3 2 2 1 1 total 0.092035 0.003633 mesh 1 group in group out legendre nuclide mean std. dev. x y z -0 1 1 1 1 1 P0 total 0.763779 0.070696 -1 1 1 1 1 1 P1 total 0.288556 0.024446 -2 1 1 1 1 1 P2 total 0.082441 0.011443 -3 1 1 1 1 1 P3 total -0.005627 0.012638 -8 1 2 1 1 1 P0 total 0.628158 0.064356 -9 1 2 1 1 1 P1 total 0.245583 0.022676 -10 1 2 1 1 1 P2 total 0.086370 0.007833 -11 1 2 1 1 1 P3 total 0.019590 0.005345 -4 2 1 1 1 1 P0 total 0.640809 0.158369 -5 2 1 1 1 1 P1 total 0.273553 0.066437 -6 2 1 1 1 1 P2 total 0.108446 0.024435 -7 2 1 1 1 1 P3 total 0.012229 0.003785 -12 2 2 1 1 1 P0 total 0.645171 0.080467 -13 2 2 1 1 1 P1 total 0.252215 0.032154 -14 2 2 1 1 1 P2 total 0.089251 0.009734 -15 2 2 1 1 1 P3 total 0.004748 0.002987 +0 1 1 1 1 1 P0 total 0.087684 0.005584 +1 1 1 1 1 1 P1 total 0.026787 0.002493 +2 1 1 1 1 1 P2 total 0.014937 0.001035 +3 1 1 1 1 1 P3 total 0.007893 0.001109 +8 1 2 1 1 1 P0 total 0.090687 0.005242 +9 1 2 1 1 1 P1 total 0.029516 0.002004 +10 1 2 1 1 1 P2 total 0.016952 0.001093 +11 1 2 1 1 1 P3 total 0.008019 0.001095 +4 2 1 1 1 1 P0 total 0.093670 0.005616 +5 2 1 1 1 1 P1 total 0.031703 0.002177 +6 2 1 1 1 1 P2 total 0.017922 0.001352 +7 2 1 1 1 1 P3 total 0.011171 0.001055 +12 2 2 1 1 1 P0 total 0.091808 0.003617 +13 2 2 1 1 1 P1 total 0.030025 0.001876 +14 2 2 1 1 1 P2 total 0.015181 0.002277 +15 2 2 1 1 1 P3 total 0.009550 0.001713 mesh 1 group in group out legendre nuclide mean std. dev. x y z -0 1 1 1 1 1 P0 total 0.763779 0.070696 -1 1 1 1 1 1 P1 total 0.288556 0.024446 -2 1 1 1 1 1 P2 total 0.082441 0.011443 -3 1 1 1 1 1 P3 total -0.005627 0.012638 -8 1 2 1 1 1 P0 total 0.628158 0.064356 -9 1 2 1 1 1 P1 total 0.245583 0.022676 -10 1 2 1 1 1 P2 total 0.086370 0.007833 -11 1 2 1 1 1 P3 total 0.019590 0.005345 -4 2 1 1 1 1 P0 total 0.640809 0.158369 -5 2 1 1 1 1 P1 total 0.273553 0.066437 -6 2 1 1 1 1 P2 total 0.108446 0.024435 -7 2 1 1 1 1 P3 total 0.012229 0.003785 -12 2 2 1 1 1 P0 total 0.645171 0.080467 -13 2 2 1 1 1 P1 total 0.252215 0.032154 -14 2 2 1 1 1 P2 total 0.089251 0.009734 -15 2 2 1 1 1 P3 total 0.004748 0.002987 - mesh 1 group in group out nuclide mean std. dev. - x y z -0 1 1 1 1 1 total 1.0 0.108337 -2 1 2 1 1 1 total 1.0 0.113128 -1 2 1 1 1 1 total 1.0 0.238517 -3 2 2 1 1 1 total 1.0 0.132597 +0 1 1 1 1 1 P0 total 0.087817 0.005624 +1 1 1 1 1 1 P1 total 0.026785 0.002504 +2 1 1 1 1 1 P2 total 0.014973 0.001041 +3 1 1 1 1 1 P3 total 0.007913 0.001144 +8 1 2 1 1 1 P0 total 0.090790 0.005246 +9 1 2 1 1 1 P1 total 0.029477 0.001987 +10 1 2 1 1 1 P2 total 0.016940 0.001094 +11 1 2 1 1 1 P3 total 0.008033 0.001104 +4 2 1 1 1 1 P0 total 0.093736 0.005609 +5 2 1 1 1 1 P1 total 0.031651 0.002201 +6 2 1 1 1 1 P2 total 0.017953 0.001364 +7 2 1 1 1 1 P3 total 0.011158 0.001044 +12 2 2 1 1 1 P0 total 0.092035 0.003633 +13 2 2 1 1 1 P1 total 0.030055 0.001856 +14 2 2 1 1 1 P2 total 0.015245 0.002274 +15 2 2 1 1 1 P3 total 0.009534 0.001700 mesh 1 group in group out nuclide mean std. dev. x y z -0 1 1 1 1 1 total 0.015584 0.003404 -2 1 2 1 1 1 total 0.017684 0.002499 -1 2 1 1 1 1 total 0.014200 0.003676 -3 2 2 1 1 1 total 0.022409 0.002481 +0 1 1 1 1 1 total 1.001515 0.075311 +2 1 2 1 1 1 total 1.001135 0.061671 +1 2 1 1 1 1 total 1.000704 0.055977 +3 2 2 1 1 1 total 1.002471 0.042246 + mesh 1 group in group out nuclide mean std. dev. + x y z +0 1 1 1 1 1 total 0.031246 0.001839 +2 1 2 1 1 1 total 0.032452 0.002365 +1 2 1 1 1 1 total 0.032568 0.002068 +3 2 2 1 1 1 total 0.031529 0.001639 mesh 1 group in group out nuclide mean std. dev. x y z -0 1 1 1 1 1 total 1.0 0.108337 -2 1 2 1 1 1 total 1.0 0.113128 -1 2 1 1 1 1 total 1.0 0.238517 -3 2 2 1 1 1 total 1.0 0.132597 +0 1 1 1 1 1 total 1.0 0.074891 +2 1 2 1 1 1 total 1.0 0.061618 +1 2 1 1 1 1 total 1.0 0.056068 +3 2 2 1 1 1 total 1.0 0.042067 mesh 1 group in group out legendre nuclide mean std. dev. x y z -0 1 1 1 1 1 P0 total 0.735256 0.113047 -1 1 1 1 1 1 P1 total 0.277780 0.041434 -2 1 1 1 1 1 P2 total 0.079362 0.014706 -3 1 1 1 1 1 P3 total -0.005417 0.012184 -8 1 2 1 1 1 P0 total 0.624575 0.110512 -9 1 2 1 1 1 P1 total 0.244182 0.041824 -10 1 2 1 1 1 P2 total 0.085877 0.014634 -11 1 2 1 1 1 P3 total 0.019478 0.006012 -4 2 1 1 1 1 P0 total 0.633925 0.212349 -5 2 1 1 1 1 P1 total 0.270615 0.089799 -6 2 1 1 1 1 P2 total 0.107281 0.034246 -7 2 1 1 1 1 P3 total 0.012098 0.004637 -12 2 2 1 1 1 P0 total 0.655214 0.126119 -13 2 2 1 1 1 P1 total 0.256141 0.049765 -14 2 2 1 1 1 P2 total 0.090641 0.016563 -15 2 2 1 1 1 P3 total 0.004822 0.003115 +0 1 1 1 1 1 P0 total 0.091790 0.008806 +1 1 1 1 1 1 P1 total 0.028042 0.003295 +2 1 1 1 1 1 P2 total 0.015636 0.001560 +3 1 1 1 1 1 P3 total 0.008263 0.001304 +8 1 2 1 1 1 P0 total 0.091883 0.006252 +9 1 2 1 1 1 P1 total 0.029905 0.002297 +10 1 2 1 1 1 P2 total 0.017175 0.001268 +11 1 2 1 1 1 P3 total 0.008124 0.001147 +4 2 1 1 1 1 P0 total 0.092530 0.006177 +5 2 1 1 1 1 P1 total 0.031317 0.002339 +6 2 1 1 1 1 P2 total 0.017704 0.001433 +7 2 1 1 1 1 P3 total 0.011035 0.001092 +12 2 2 1 1 1 P0 total 0.089619 0.004144 +13 2 2 1 1 1 P1 total 0.029309 0.001964 +14 2 2 1 1 1 P2 total 0.014820 0.002251 +15 2 2 1 1 1 P3 total 0.009322 0.001687 mesh 1 group in group out legendre nuclide mean std. dev. x y z -0 1 1 1 1 1 P0 total 0.735256 0.138292 -1 1 1 1 1 1 P1 total 0.277780 0.051210 -2 1 1 1 1 1 P2 total 0.079362 0.017035 -3 1 1 1 1 1 P3 total -0.005417 0.012198 -8 1 2 1 1 1 P0 total 0.624575 0.131169 -9 1 2 1 1 1 P1 total 0.244182 0.050123 -10 1 2 1 1 1 P2 total 0.085877 0.017565 -11 1 2 1 1 1 P3 total 0.019478 0.006403 -4 2 1 1 1 1 P0 total 0.633925 0.260681 -5 2 1 1 1 1 P1 total 0.270615 0.110590 -6 2 1 1 1 1 P2 total 0.107281 0.042750 -7 2 1 1 1 1 P3 total 0.012098 0.005462 -12 2 2 1 1 1 P0 total 0.655214 0.153147 -13 2 2 1 1 1 P1 total 0.256141 0.060250 -14 2 2 1 1 1 P2 total 0.090641 0.020464 -15 2 2 1 1 1 P3 total 0.004822 0.003180 +0 1 1 1 1 1 P0 total 0.091929 0.011205 +1 1 1 1 1 1 P1 total 0.028084 0.003918 +2 1 1 1 1 1 P2 total 0.015660 0.001956 +3 1 1 1 1 1 P3 total 0.008276 0.001447 +8 1 2 1 1 1 P0 total 0.091987 0.008443 +9 1 2 1 1 1 P1 total 0.029939 0.002948 +10 1 2 1 1 1 P2 total 0.017195 0.001653 +11 1 2 1 1 1 P3 total 0.008134 0.001253 +4 2 1 1 1 1 P0 total 0.092595 0.008065 +5 2 1 1 1 1 P1 total 0.031339 0.002924 +6 2 1 1 1 1 P2 total 0.017716 0.001743 +7 2 1 1 1 1 P3 total 0.011042 0.001255 +12 2 2 1 1 1 P0 total 0.089840 0.005621 +13 2 2 1 1 1 P1 total 0.029381 0.002326 +14 2 2 1 1 1 P2 total 0.014856 0.002342 +15 2 2 1 1 1 P3 total 0.009345 0.001737 mesh 1 group out nuclide mean std. dev. x y z -0 1 1 1 1 total 1.0 0.300047 -2 1 2 1 1 total 1.0 0.178169 -1 2 1 1 1 total 1.0 0.262180 -3 2 2 1 1 total 1.0 0.104797 +0 1 1 1 1 total 1.0 0.066520 +2 1 2 1 1 total 1.0 0.087934 +1 2 1 1 1 total 1.0 0.063390 +3 2 2 1 1 total 1.0 0.063791 mesh 1 group out nuclide mean std. dev. x y z -0 1 1 1 1 total 1.0 0.300047 -2 1 2 1 1 total 1.0 0.178169 -1 2 1 1 1 total 1.0 0.262180 -3 2 2 1 1 total 1.0 0.108931 +0 1 1 1 1 total 1.0 0.068463 +2 1 2 1 1 total 1.0 0.091776 +1 2 1 1 1 total 1.0 0.064705 +3 2 2 1 1 total 1.0 0.063003 mesh 1 group in nuclide mean std. dev. x y z -0 1 1 1 1 total 7.097008e-07 1.458546e-07 -2 1 2 1 1 total 4.407745e-07 7.903907e-08 -1 2 1 1 1 total 3.984535e-07 1.157576e-07 -3 2 2 1 1 total 4.750476e-07 6.207437e-08 +0 1 1 1 1 total 8.735713e-10 4.530341e-11 +2 1 2 1 1 total 8.821319e-10 3.206094e-11 +1 2 1 1 1 total 8.699208e-10 2.515246e-11 +3 2 2 1 1 total 8.738762e-10 1.734562e-11 mesh 1 group in nuclide mean std. dev. x y z -0 1 1 1 1 total 0.027311 0.007397 -2 1 2 1 1 total 0.017820 0.003054 -1 2 1 1 1 total 0.017783 0.004394 -3 2 2 1 1 total 0.020320 0.002598 +0 1 1 1 1 total 0.031799 0.002147 +2 1 2 1 1 total 0.031680 0.001263 +1 2 1 1 1 total 0.031989 0.001691 +3 2 2 1 1 total 0.030533 0.000994 mesh 1 group in group out nuclide mean std. dev. x y z -0 1 1 1 1 1 total 0.015584 0.003404 -2 1 2 1 1 1 total 0.017684 0.002499 -1 2 1 1 1 1 total 0.014200 0.003676 -3 2 2 1 1 1 total 0.022259 0.002508 +0 1 1 1 1 1 total 0.031056 0.001862 +2 1 2 1 1 1 total 0.032188 0.002420 +1 2 1 1 1 1 total 0.032304 0.002073 +3 2 2 1 1 1 total 0.031336 0.001614 mesh 1 delayedgroup group in nuclide mean std. dev. x y z -0 1 1 1 1 1 total 0.000006 1.689606e-06 -1 1 1 1 2 1 total 0.000033 8.718916e-06 -2 1 1 1 3 1 total 0.000032 8.323051e-06 -3 1 1 1 4 1 total 0.000072 1.866015e-05 -4 1 1 1 5 1 total 0.000031 7.654909e-06 -5 1 1 1 6 1 total 0.000013 3.206343e-06 -12 1 2 1 1 1 total 0.000004 6.723192e-07 -13 1 2 1 2 1 total 0.000022 3.706235e-06 -14 1 2 1 3 1 total 0.000022 3.674263e-06 -15 1 2 1 4 1 total 0.000052 8.774048e-06 -16 1 2 1 5 1 total 0.000024 4.168024e-06 -17 1 2 1 6 1 total 0.000010 1.726268e-06 -6 2 1 1 1 1 total 0.000004 1.003100e-06 -7 2 1 1 2 1 total 0.000022 5.425275e-06 -8 2 1 1 3 1 total 0.000021 5.324236e-06 -9 2 1 1 4 1 total 0.000050 1.251572e-05 -10 2 1 1 5 1 total 0.000022 5.762184e-06 -11 2 1 1 6 1 total 0.000009 2.391676e-06 -18 2 2 1 1 1 total 0.000005 5.962367e-07 -19 2 2 1 2 1 total 0.000025 3.200900e-06 -20 2 2 1 3 1 total 0.000025 3.127442e-06 -21 2 2 1 4 1 total 0.000058 7.296157e-06 -22 2 2 1 5 1 total 0.000026 3.298196e-06 -23 2 2 1 6 1 total 0.000011 1.370918e-06 +0 1 1 1 1 1 total 0.000007 4.734745e-07 +1 1 1 1 2 1 total 0.000036 2.443930e-06 +2 1 1 1 3 1 total 0.000035 2.333188e-06 +3 1 1 1 4 1 total 0.000078 5.231199e-06 +4 1 1 1 5 1 total 0.000032 2.144718e-06 +5 1 1 1 6 1 total 0.000013 8.984148e-07 +12 1 2 1 1 1 total 0.000007 2.770884e-07 +13 1 2 1 2 1 total 0.000036 1.430245e-06 +14 1 2 1 3 1 total 0.000035 1.365436e-06 +15 1 2 1 4 1 total 0.000078 3.061421e-06 +16 1 2 1 5 1 total 0.000032 1.255139e-06 +17 1 2 1 6 1 total 0.000013 5.257735e-07 +6 2 1 1 1 1 total 0.000007 3.731284e-07 +7 2 1 1 2 1 total 0.000037 1.925974e-06 +8 2 1 1 3 1 total 0.000035 1.838702e-06 +9 2 1 1 4 1 total 0.000079 4.122522e-06 +10 2 1 1 5 1 total 0.000032 1.690176e-06 +11 2 1 1 6 1 total 0.000014 7.080087e-07 +18 2 2 1 1 1 total 0.000007 2.050310e-07 +19 2 2 1 2 1 total 0.000035 1.058307e-06 +20 2 2 1 3 1 total 0.000033 1.010352e-06 +21 2 2 1 4 1 total 0.000075 2.265292e-06 +22 2 2 1 5 1 total 0.000031 9.287379e-07 +23 2 2 1 6 1 total 0.000013 3.890451e-07 mesh 1 delayedgroup group out nuclide mean std. dev. x y z 0 1 1 1 1 1 total 0.0 0.000000 -1 1 1 1 2 1 total 0.0 0.000000 -2 1 1 1 3 1 total 0.0 0.000000 -3 1 1 1 4 1 total 0.0 0.000000 -4 1 1 1 5 1 total 0.0 0.000000 +1 1 1 1 2 1 total 1.0 1.414214 +2 1 1 1 3 1 total 1.0 0.868831 +3 1 1 1 4 1 total 1.0 1.414214 +4 1 1 1 5 1 total 1.0 1.414214 5 1 1 1 6 1 total 0.0 0.000000 12 1 2 1 1 1 total 0.0 0.000000 -13 1 2 1 2 1 total 0.0 0.000000 +13 1 2 1 2 1 total 1.0 0.866827 14 1 2 1 3 1 total 0.0 0.000000 -15 1 2 1 4 1 total 0.0 0.000000 -16 1 2 1 5 1 total 0.0 0.000000 +15 1 2 1 4 1 total 1.0 0.455171 +16 1 2 1 5 1 total 1.0 0.868553 17 1 2 1 6 1 total 0.0 0.000000 6 2 1 1 1 1 total 0.0 0.000000 -7 2 1 1 2 1 total 0.0 0.000000 -8 2 1 1 3 1 total 0.0 0.000000 -9 2 1 1 4 1 total 0.0 0.000000 -10 2 1 1 5 1 total 0.0 0.000000 -11 2 1 1 6 1 total 0.0 0.000000 -18 2 2 1 1 1 total 0.0 0.000000 -19 2 2 1 2 1 total 0.0 0.000000 +7 2 1 1 2 1 total 1.0 1.414214 +8 2 1 1 3 1 total 1.0 1.414214 +9 2 1 1 4 1 total 1.0 0.674843 +10 2 1 1 5 1 total 1.0 1.414214 +11 2 1 1 6 1 total 1.0 0.866033 +18 2 2 1 1 1 total 1.0 1.414214 +19 2 2 1 2 1 total 1.0 1.414214 20 2 2 1 3 1 total 1.0 1.414214 -21 2 2 1 4 1 total 0.0 0.000000 +21 2 2 1 4 1 total 1.0 0.579059 22 2 2 1 5 1 total 0.0 0.000000 -23 2 2 1 6 1 total 0.0 0.000000 +23 2 2 1 6 1 total 1.0 1.414214 mesh 1 delayedgroup group in nuclide mean std. dev. x y z -0 1 1 1 1 1 total 0.000228 0.000084 -1 1 1 1 2 1 total 0.001195 0.000438 -2 1 1 1 3 1 total 0.001153 0.000420 -3 1 1 1 4 1 total 0.002629 0.000950 -4 1 1 1 5 1 total 0.001125 0.000398 -5 1 1 1 6 1 total 0.000470 0.000166 -12 1 2 1 1 1 total 0.000225 0.000044 -13 1 2 1 2 1 total 0.001232 0.000242 -14 1 2 1 3 1 total 0.001216 0.000239 -15 1 2 1 4 1 total 0.002882 0.000570 -16 1 2 1 5 1 total 0.001345 0.000270 -17 1 2 1 6 1 total 0.000558 0.000112 -6 2 1 1 1 1 total 0.000228 0.000057 -7 2 1 1 2 1 total 0.001222 0.000309 -8 2 1 1 3 1 total 0.001193 0.000304 -9 2 1 1 4 1 total 0.002780 0.000713 -10 2 1 1 5 1 total 0.001250 0.000328 -11 2 1 1 6 1 total 0.000520 0.000136 -18 2 2 1 1 1 total 0.000227 0.000027 -19 2 2 1 2 1 total 0.001225 0.000143 -20 2 2 1 3 1 total 0.001201 0.000140 -21 2 2 1 4 1 total 0.002815 0.000326 -22 2 2 1 5 1 total 0.001284 0.000147 -23 2 2 1 6 1 total 0.000533 0.000061 +0 1 1 1 1 1 total 0.000221 0.000019 +1 1 1 1 2 1 total 0.001139 0.000096 +2 1 1 1 3 1 total 0.001087 0.000092 +3 1 1 1 4 1 total 0.002437 0.000206 +4 1 1 1 5 1 total 0.000999 0.000084 +5 1 1 1 6 1 total 0.000419 0.000035 +12 1 2 1 1 1 total 0.000222 0.000012 +13 1 2 1 2 1 total 0.001144 0.000060 +14 1 2 1 3 1 total 0.001092 0.000057 +15 1 2 1 4 1 total 0.002448 0.000129 +16 1 2 1 5 1 total 0.001004 0.000053 +17 1 2 1 6 1 total 0.000420 0.000022 +6 2 1 1 1 1 total 0.000221 0.000015 +7 2 1 1 2 1 total 0.001143 0.000078 +8 2 1 1 3 1 total 0.001091 0.000075 +9 2 1 1 4 1 total 0.002446 0.000167 +10 2 1 1 5 1 total 0.001003 0.000069 +11 2 1 1 6 1 total 0.000420 0.000029 +18 2 2 1 1 1 total 0.000220 0.000009 +19 2 2 1 2 1 total 0.001136 0.000047 +20 2 2 1 3 1 total 0.001084 0.000045 +21 2 2 1 4 1 total 0.002431 0.000100 +22 2 2 1 5 1 total 0.000997 0.000041 +23 2 2 1 6 1 total 0.000417 0.000017 mesh 1 delayedgroup group in nuclide mean std. dev. x y z -0 1 1 1 1 1 total 0.013345 0.004923 -1 1 1 1 2 1 total 0.032674 0.011850 -2 1 1 1 3 1 total 0.120923 0.043307 -3 1 1 1 4 1 total 0.304289 0.106753 -4 1 1 1 5 1 total 0.855760 0.286466 -5 1 1 1 6 1 total 2.874120 0.965609 -12 1 2 1 1 1 total 0.013367 0.002548 -13 1 2 1 2 1 total 0.032520 0.006266 -14 1 2 1 3 1 total 0.121250 0.023544 -15 1 2 1 4 1 total 0.307552 0.060464 -16 1 2 1 5 1 total 0.867665 0.175131 -17 1 2 1 6 1 total 2.914635 0.587161 -6 2 1 1 1 1 total 0.013357 0.003345 -7 2 1 1 2 1 total 0.032590 0.008273 -8 2 1 1 3 1 total 0.121103 0.031074 -9 2 1 1 4 1 total 0.306111 0.080011 -10 2 1 1 5 1 total 0.862660 0.235694 -11 2 1 1 6 1 total 2.897534 0.788926 -18 2 2 1 1 1 total 0.013360 0.001587 -19 2 2 1 2 1 total 0.032564 0.003810 -20 2 2 1 3 1 total 0.121158 0.014038 -21 2 2 1 4 1 total 0.306653 0.035052 -22 2 2 1 5 1 total 0.864587 0.096680 -23 2 2 1 6 1 total 2.904111 0.325170 - mesh 1 delayedgroup group in group out nuclide mean std. dev. - x y z -0 1 1 1 1 1 1 total 0.00000 0.000000 -1 1 1 1 2 1 1 total 0.00000 0.000000 -2 1 1 1 3 1 1 total 0.00000 0.000000 -3 1 1 1 4 1 1 total 0.00000 0.000000 -4 1 1 1 5 1 1 total 0.00000 0.000000 -5 1 1 1 6 1 1 total 0.00000 0.000000 -12 1 2 1 1 1 1 total 0.00000 0.000000 -13 1 2 1 2 1 1 total 0.00000 0.000000 -14 1 2 1 3 1 1 total 0.00000 0.000000 -15 1 2 1 4 1 1 total 0.00000 0.000000 -16 1 2 1 5 1 1 total 0.00000 0.000000 -17 1 2 1 6 1 1 total 0.00000 0.000000 -6 2 1 1 1 1 1 total 0.00000 0.000000 -7 2 1 1 2 1 1 total 0.00000 0.000000 -8 2 1 1 3 1 1 total 0.00000 0.000000 -9 2 1 1 4 1 1 total 0.00000 0.000000 -10 2 1 1 5 1 1 total 0.00000 0.000000 -11 2 1 1 6 1 1 total 0.00000 0.000000 -18 2 2 1 1 1 1 total 0.00000 0.000000 -19 2 2 1 2 1 1 total 0.00000 0.000000 -20 2 2 1 3 1 1 total 0.00015 0.000151 -21 2 2 1 4 1 1 total 0.00000 0.000000 -22 2 2 1 5 1 1 total 0.00000 0.000000 -23 2 2 1 6 1 1 total 0.00000 0.000000 +0 1 1 1 1 1 total 0.013336 0.001120 +1 1 1 1 2 1 total 0.032739 0.002751 +2 1 1 1 3 1 total 0.120780 0.010147 +3 1 1 1 4 1 total 0.302780 0.025438 +4 1 1 1 5 1 total 0.849490 0.071370 +5 1 1 1 6 1 total 2.853000 0.239696 +12 1 2 1 1 1 total 0.013336 0.000695 +13 1 2 1 2 1 total 0.032739 0.001707 +14 1 2 1 3 1 total 0.120780 0.006296 +15 1 2 1 4 1 total 0.302780 0.015784 +16 1 2 1 5 1 total 0.849490 0.044285 +17 1 2 1 6 1 total 2.853000 0.148730 +6 2 1 1 1 1 total 0.013336 0.000906 +7 2 1 1 2 1 total 0.032739 0.002224 +8 2 1 1 3 1 total 0.120780 0.008205 +9 2 1 1 4 1 total 0.302780 0.020570 +10 2 1 1 5 1 total 0.849490 0.057711 +11 2 1 1 6 1 total 2.853000 0.193821 +18 2 2 1 1 1 total 0.013336 0.000528 +19 2 2 1 2 1 total 0.032739 0.001296 +20 2 2 1 3 1 total 0.120780 0.004781 +21 2 2 1 4 1 total 0.302780 0.011986 +22 2 2 1 5 1 total 0.849490 0.033628 +23 2 2 1 6 1 total 2.853000 0.112940 + mesh 1 delayedgroup group in group out nuclide mean std. dev. + x y z +0 1 1 1 1 1 1 total 0.000000 0.000000 +1 1 1 1 2 1 1 total 0.000055 0.000055 +2 1 1 1 3 1 1 total 0.000055 0.000034 +3 1 1 1 4 1 1 total 0.000052 0.000052 +4 1 1 1 5 1 1 total 0.000029 0.000029 +5 1 1 1 6 1 1 total 0.000000 0.000000 +12 1 2 1 1 1 1 total 0.000000 0.000000 +13 1 2 1 2 1 1 total 0.000058 0.000036 +14 1 2 1 3 1 1 total 0.000000 0.000000 +15 1 2 1 4 1 1 total 0.000149 0.000048 +16 1 2 1 5 1 1 total 0.000057 0.000035 +17 1 2 1 6 1 1 total 0.000000 0.000000 +6 2 1 1 1 1 1 total 0.000000 0.000000 +7 2 1 1 2 1 1 total 0.000027 0.000027 +8 2 1 1 3 1 1 total 0.000026 0.000026 +9 2 1 1 4 1 1 total 0.000105 0.000051 +10 2 1 1 5 1 1 total 0.000054 0.000054 +11 2 1 1 6 1 1 total 0.000051 0.000031 +18 2 2 1 1 1 1 total 0.000028 0.000028 +19 2 2 1 2 1 1 total 0.000025 0.000025 +20 2 2 1 3 1 1 total 0.000032 0.000032 +21 2 2 1 4 1 1 total 0.000080 0.000033 +22 2 2 1 5 1 1 total 0.000000 0.000000 +23 2 2 1 6 1 1 total 0.000027 0.000027 diff --git a/tests/regression_tests/mgxs_library_mesh/test.py b/tests/regression_tests/mgxs_library_mesh/test.py index b72fdf2f1e..3660d3eb74 100644 --- a/tests/regression_tests/mgxs_library_mesh/test.py +++ b/tests/regression_tests/mgxs_library_mesh/test.py @@ -2,42 +2,67 @@ import hashlib import openmc import openmc.mgxs +import pytest from tests.testing_harness import PyAPITestHarness +@pytest.fixture +def model(): + model = openmc.model.Model() + + fuel = openmc.Material() + fuel.set_density('g/cm3', 10.0) + fuel.add_nuclide('U235', 1.0) + zr = openmc.Material() + zr.set_density('g/cm3', 1.0) + zr.add_nuclide('Zr90', 1.0) + model.materials.extend([fuel, zr]) + + box1 = openmc.model.rectangular_prism(10.0, 10.0) + box2 = openmc.model.rectangular_prism(20.0, 20.0, boundary_type='reflective') + top = openmc.ZPlane(z0=10.0, boundary_type='vacuum') + bottom = openmc.ZPlane(z0=-10.0, boundary_type='vacuum') + cell1 = openmc.Cell(fill=fuel, region=box1 & +bottom & -top) + cell2 = openmc.Cell(fill=zr, region=~box1 & box2 & +bottom & -top) + model.geometry = openmc.Geometry([cell1, cell2]) + + model.settings.batches = 5 + model.settings.inactive = 0 + model.settings.particles = 1000 + + # Initialize a one-group structure + energy_groups = openmc.mgxs.EnergyGroups([0, 20.e6]) + + # Initialize MGXS Library for a few cross section types + # for one material-filled cell in the geometry + model.mgxs_lib = openmc.mgxs.Library(model.geometry) + model.mgxs_lib.by_nuclide = False + + # Test all MGXS types + model.mgxs_lib.mgxs_types = openmc.mgxs.MGXS_TYPES + openmc.mgxs.MDGXS_TYPES + model.mgxs_lib.energy_groups = energy_groups + model.mgxs_lib.num_delayed_groups = 6 + model.mgxs_lib.correction = None # Avoid warning about P0 correction + model.mgxs_lib.legendre_order = 3 + model.mgxs_lib.domain_type = 'mesh' + + # Instantiate a tally mesh + mesh = openmc.RegularMesh(mesh_id=1) + mesh.dimension = [2, 2] + mesh.lower_left = [-100., -100.] + mesh.width = [100., 100.] + + model.mgxs_lib.domains = [mesh] + model.mgxs_lib.build_library() + + # Add tallies + model.mgxs_lib.add_to_tallies_file(model.tallies, merge=False) + + return model + + class MGXSTestHarness(PyAPITestHarness): - def __init__(self, *args, **kwargs): - super().__init__(*args, **kwargs) - - # Initialize a one-group structure - energy_groups = openmc.mgxs.EnergyGroups(group_edges=[0, 20.e6]) - - # Initialize MGXS Library for a few cross section types - # for one material-filled cell in the geometry - self.mgxs_lib = openmc.mgxs.Library(self._model.geometry) - self.mgxs_lib.by_nuclide = False - - # Test all MGXS types - self.mgxs_lib.mgxs_types = openmc.mgxs.MGXS_TYPES + \ - openmc.mgxs.MDGXS_TYPES - self.mgxs_lib.energy_groups = energy_groups - self.mgxs_lib.num_delayed_groups = 6 - self.mgxs_lib.legendre_order = 3 - self.mgxs_lib.domain_type = 'mesh' - - # Instantiate a tally mesh - mesh = openmc.RegularMesh(mesh_id=1) - mesh.dimension = [2, 2] - mesh.lower_left = [-100., -100.] - mesh.width = [100., 100.] - - self.mgxs_lib.domains = [mesh] - self.mgxs_lib.build_library() - - # Add tallies - self.mgxs_lib.add_to_tallies_file(self._model.tallies, merge=False) - def _get_results(self, hash_output=False): """Digest info in the statepoint and return as a string.""" @@ -45,13 +70,14 @@ class MGXSTestHarness(PyAPITestHarness): sp = openmc.StatePoint(self._sp_name) # Load the MGXS library from the statepoint - self.mgxs_lib.load_from_statepoint(sp) + mgxs_lib = self._model.mgxs_lib + mgxs_lib.load_from_statepoint(sp) # Build a string from Pandas Dataframe for each 1-group MGXS outstr = '' - for domain in self.mgxs_lib.domains: - for mgxs_type in self.mgxs_lib.mgxs_types: - mgxs = self.mgxs_lib.get_mgxs(domain, mgxs_type) + for domain in mgxs_lib.domains: + for mgxs_type in mgxs_lib.mgxs_types: + mgxs = mgxs_lib.get_mgxs(domain, mgxs_type) df = mgxs.get_pandas_dataframe() outstr += df.to_string() + '\n' @@ -64,6 +90,6 @@ class MGXSTestHarness(PyAPITestHarness): return outstr -def test_mgxs_library_mesh(): - harness = MGXSTestHarness('statepoint.10.h5') +def test_mgxs_library_mesh(model): + harness = MGXSTestHarness('statepoint.5.h5', model) harness.main()