This commit is contained in:
GuySten 2026-07-20 16:30:47 +02:00 committed by GitHub
commit 379514eef7
No known key found for this signature in database
GPG key ID: B5690EEEBB952194
88 changed files with 2469 additions and 2347 deletions

View file

@ -112,11 +112,20 @@ private:
//! Determine if a particle is inside the cell for a complex cell.
//!
//! Uses the comobination of half-spaces and binary operators to determine
//! Uses the combination of half-spaces and binary operators to determine
//! if short circuiting can be used. Short cicuiting uses the relative and
//! absolute depth of parentheses in the expression.
bool contains_complex(Position r, Direction u, int32_t on_surface) const;
//! Find the oncoming boundary of this cell for a simple cell (only
//! intersection operators)
std::pair<double, int32_t> distance_simple(
Position r, Direction u, int32_t on_surface) const;
//! Find the oncoming boundary of this cell for a complex cell.
std::pair<double, int32_t> distance_complex(
Position r, Direction u, int32_t on_surface) const;
//! BoundingBox if the paritcle is in a simple cell.
BoundingBox bounding_box_simple() const;

View file

@ -945,6 +945,18 @@ std::string Region::str() const
std::pair<double, int32_t> Region::distance(
Position r, Direction u, int32_t on_surface) const
{
if (simple_) {
return distance_simple(r, u, on_surface);
} else {
return distance_complex(r, u, on_surface);
}
}
//==============================================================================
std::pair<double, int32_t> Region::distance_simple(
Position r, Direction u, int32_t on_surface) const
{
double min_dist {INFTY};
int32_t i_surf {std::numeric_limits<int32_t>::max()};
@ -973,6 +985,49 @@ std::pair<double, int32_t> Region::distance(
//==============================================================================
std::pair<double, int32_t> Region::distance_complex(
Position r, Direction u, int32_t on_surface) const
{
double min_dist;
int32_t i_surf;
double atleast {-1.0};
bool in_region = contains_complex(r, u, on_surface);
while (true) {
min_dist = INFTY;
i_surf = std::numeric_limits<int32_t>::max();
for (int32_t token : expression_) {
// Ignore this token if it corresponds to an operator rather than a
// region.
if (token >= OP_UNION)
continue;
// Calculate the distance to this surface.
// Note the off-by-one indexing
bool coincident {std::abs(token) == std::abs(on_surface)};
double d {model::surfaces[abs(token) - 1]->distance(r, u, coincident)};
// Check if this distance is the new minimum.
if ((d > atleast) && (d < min_dist)) {
if (min_dist - d >= FP_PRECISION * min_dist) {
min_dist = d;
i_surf = -token;
}
}
}
if (min_dist == INFTY)
break;
auto p = r + (min_dist + TINY_BIT) * u;
if (contains_complex(p, u, on_surface) != in_region)
break;
atleast = min_dist;
}
return {min_dist, i_surf};
}
//==============================================================================
bool Region::contains(Position r, Direction u, int32_t on_surface) const
{
if (simple_) {

View file

@ -0,0 +1,2 @@
k-combined:
6.085380E-02 2.048137E-03

View file

@ -0,0 +1,2 @@
k-combined:
6.085380E-02 2.048137E-03

View file

@ -0,0 +1,2 @@
k-combined:
6.085380E-02 2.048137E-03

View file

@ -0,0 +1,2 @@
k-combined:
6.085380E-02 2.048137E-03

View file

@ -0,0 +1,2 @@
k-combined:
6.085380E-02 2.048137E-03

View file

@ -0,0 +1,2 @@
k-combined:
6.085380E-02 2.048137E-03

View file

@ -0,0 +1,2 @@
k-combined:
6.085380E-02 2.048137E-03

View file

@ -0,0 +1,2 @@
k-combined:
6.085380E-02 2.048137E-03

View file

@ -1,11 +1,11 @@
k-combined:
2.603220E-01 1.429366E-03
2.589444E-01 9.409859E-04
tally 1:
2.624819E+00
1.378200E+00
2.730035E+00
1.492361E+00
1.013707E+00
2.055807E-01
1.123257E-01
2.530233E-03
2.578878E+00
1.331090E+00
2.702603E+00
1.462070E+00
9.847537E-01
1.942139E-01
1.084095E-01
2.358175E-03

View file

@ -1,53 +1,53 @@
k-combined:
9.035025E-02 2.654309E-03
9.946255E-02 2.184346E-03
tally 1:
5.994069E+00
3.594398E+00
6.089246E+00
3.712059E+00
tally 2:
4.559707E-04
2.080739E-08
4.646490E-04
2.160545E-08
tally 3:
5.994525E+00
3.594945E+00
6.089710E+00
3.712624E+00
tally 4:
0.000000E+00
0.000000E+00
tally 5:
1.473892E+00
2.187509E-01
1.556619E+00
2.452581E-01
tally 6:
5.223875E-05
2.992861E-10
0.000000E+00
0.000000E+00
tally 7:
1.473945E+00
2.187663E-01
1.556619E+00
2.452581E-01
tally 8:
1.885798E+01
3.558423E+01
1.001104E+02
1.003264E+03
tally 9:
7.467961E+00
5.580255E+00
7.645865E+00
5.850005E+00
tally 10:
5.082094E-04
2.584432E-08
4.646490E-04
2.160545E-08
tally 11:
7.468470E+00
5.581014E+00
7.646330E+00
5.850716E+00
tally 12:
1.885798E+01
3.558423E+01
1.001104E+02
1.003264E+03
tally 13:
0.000000E+00
0.000000E+00
tally 14:
2.739543E-04
7.600983E-09
3.197805E-04
1.030696E-08
tally 15:
2.739543E-04
7.600983E-09
3.197805E-04
1.030696E-08
tally 16:
7.881296E+01
6.221087E+02
0.000000E+00
0.000000E+00
tally 17:
1.051397E+02
1.106292E+03
1.077571E+02
1.162197E+03

View file

@ -1 +1 @@
e07ed2bc8893c69721abf61b123336f1f6128a3bee6ec63b84d1f549f31707a74a6ce885091ccc0eac6b7f16f7cab39ede4784584c08825829e108de878ea5fb
62452fff4ffd89deba11463d67bc60a72aa2f6d65f773b113e9ed34b3c4b4ecc7c189401c9752d0f7a60c57bd62ba1a3447890f39e7217e2de5957b2ea16dd65

View file

@ -1,244 +1,244 @@
k-combined:
7.729082E-01 3.775399E-02
8.023345E-01 6.979356E-02
tally 1:
5.296804E-02
5.661701E-04
8.356446E-02
1.412139E-03
5.041335E-02
5.143568E-04
1.299348E-01
3.467618E-03
3.929702E-01
3.147038E-02
1.379707E-01
3.888484E-03
1.405034E-01
4.473799E-03
3.785796E-01
2.940585E-02
1.422010E-01
4.113723E-03
5.647073E-02
6.735251E-04
7.911154E-02
1.329137E-03
5.160755E-02
5.361448E-04
6.669424E-02
9.090832E-04
1.008621E-01
2.134534E-03
6.808932E-02
9.355993E-04
1.873006E-01
7.135961E-03
6.221575E-01
7.819842E-02
1.856653E-01
6.954762E-03
2.014929E-01
8.327845E-03
5.853251E-01
6.945708E-02
1.709645E-01
5.917124E-03
7.214913E-02
1.058962E-03
1.027720E-01
2.138475E-03
6.099853E-02
7.493941E-04
6.892071E-02
9.630680E-04
1.035459E-01
2.173883E-03
6.973870E-02
9.904237E-04
2.125703E-01
9.112659E-03
9.012205E-01
2.163546E-01
2.066426E-01
8.617414E-03
2.258950E-01
1.039607E-02
9.476792E-01
2.350708E-01
2.225585E-01
1.017898E-02
7.111503E-02
1.036847E-03
1.117012E-01
2.530040E-03
6.870474E-02
9.551035E-04
5.738897E-02
6.699030E-04
9.522335E-02
1.835769E-03
6.570917E-02
8.656870E-04
1.945592E-01
7.593336E-03
5.514753E-01
6.122981E-02
2.144202E-01
9.421739E-03
1.971631E-01
7.944046E-03
6.088996E-01
7.442954E-02
1.965447E-01
7.765628E-03
7.005494E-02
1.012891E-03
1.010633E-01
2.084095E-03
6.145926E-02
7.694351E-04
4.999479E-02
5.129164E-04
7.238243E-02
1.062921E-03
4.902309E-02
4.852193E-04
1.324655E-01
3.642431E-03
3.305312E-01
2.265726E-02
1.332993E-01
3.728385E-03
1.547469E-01
4.894837E-03
3.625944E-01
2.747313E-02
1.435761E-01
4.334405E-03
5.789603E-02
7.065383E-04
7.589559E-02
1.205386E-03
5.210018E-02
5.790843E-04
5.745655E-02
6.712105E-04
7.867782E-02
1.248727E-03
5.868429E-02
7.131793E-04
1.438557E-01
4.308092E-03
3.787094E-01
2.891295E-02
1.374240E-01
3.869017E-03
1.347274E-01
3.834160E-03
3.987709E-01
3.261523E-02
1.334735E-01
3.626900E-03
4.778466E-02
4.648282E-04
7.840849E-02
1.297864E-03
4.916372E-02
4.947407E-04
6.995546E-02
9.876172E-04
1.087891E-01
2.387772E-03
6.400064E-02
8.393797E-04
1.945275E-01
7.631597E-03
5.867134E-01
7.069740E-02
1.792978E-01
6.602522E-03
1.955106E-01
7.732201E-03
6.226639E-01
7.780920E-02
1.796752E-01
6.469604E-03
6.679475E-02
9.049602E-04
1.008432E-01
2.075215E-03
5.586034E-02
6.531573E-04
6.137597E-02
7.755153E-04
1.147821E-01
2.767643E-03
6.809436E-02
9.285053E-04
2.104606E-01
9.134770E-03
9.269139E-01
2.245814E-01
2.021416E-01
8.202826E-03
2.176221E-01
9.854345E-03
9.336697E-01
2.326234E-01
2.095728E-01
8.873161E-03
7.246366E-02
1.052928E-03
1.083862E-01
2.417978E-03
6.342274E-02
8.587365E-04
6.418145E-02
8.295301E-04
1.066727E-01
2.322556E-03
6.830911E-02
9.380558E-04
1.832547E-01
6.752214E-03
5.428015E-01
5.938368E-02
1.951315E-01
7.663620E-03
1.572668E-01
5.032070E-03
5.625230E-01
6.367739E-02
1.918204E-01
7.408678E-03
6.451833E-02
8.482285E-04
9.426618E-02
1.787793E-03
6.648279E-02
9.126838E-04
5.317949E-02
5.712761E-04
8.572058E-02
1.489821E-03
5.734607E-02
6.866788E-04
1.365285E-01
3.892455E-03
3.670737E-01
2.795038E-02
1.429659E-01
4.492430E-03
1.222604E-01
3.045332E-03
3.261367E-01
2.184800E-02
1.255502E-01
3.263101E-03
4.699926E-02
4.464966E-04
7.180560E-02
1.057221E-03
5.028875E-02
5.175287E-04
tally 2:
4.597932E-02
4.246849E-04
8.494738E-02
1.467238E-03
5.155072E-02
5.373129E-04
1.479395E-01
4.468164E-03
3.931843E-01
3.141943E-02
1.264081E-01
3.243689E-03
1.229742E-01
3.276501E-03
3.768427E-01
2.927971E-02
1.574908E-01
5.079304E-03
5.621044E-02
6.877866E-04
8.490616E-02
1.510673E-03
4.709600E-02
4.576632E-04
5.382674E-02
5.890600E-04
1.116560E-01
2.599817E-03
6.648634E-02
8.862677E-04
1.996861E-01
8.137198E-03
6.218616E-01
7.805532E-02
1.672379E-01
5.667797E-03
1.841275E-01
6.936896E-03
5.887874E-01
7.031665E-02
1.872574E-01
7.086584E-03
7.574698E-02
1.160992E-03
1.100173E-01
2.453972E-03
5.427866E-02
5.955556E-04
5.644318E-02
6.503634E-04
1.058628E-01
2.254338E-03
8.012794E-02
1.297032E-03
2.292757E-01
1.064490E-02
9.001874E-01
2.153097E-01
1.921874E-01
7.537343E-03
2.058642E-01
8.537136E-03
9.442653E-01
2.359241E-01
2.419900E-01
1.202201E-02
7.563876E-02
1.169874E-03
1.165428E-01
2.750405E-03
5.646354E-02
6.408683E-04
4.963879E-02
4.975915E-04
1.005478E-01
2.046400E-03
7.041191E-02
1.002590E-03
1.959123E-01
7.727902E-03
5.625596E-01
6.366248E-02
1.987313E-01
8.108832E-03
1.922194E-01
7.602752E-03
6.062527E-01
7.384124E-02
2.039506E-01
8.379010E-03
7.019905E-02
1.034239E-03
1.081383E-01
2.344040E-03
5.125142E-02
5.430397E-04
4.230495E-02
3.666730E-04
7.722275E-02
1.208856E-03
4.998937E-02
5.049355E-04
1.461632E-01
4.456332E-03
3.315459E-01
2.288314E-02
1.238273E-01
3.205844E-03
1.355317E-01
3.754923E-03
3.646520E-01
2.766067E-02
1.543013E-01
5.044487E-03
6.082173E-02
7.636985E-04
8.490569E-02
1.506309E-03
4.046126E-02
3.495832E-04
4.543475E-02
4.334802E-04
8.327402E-02
1.412708E-03
5.876305E-02
7.155117E-04
1.614394E-01
5.411560E-03
3.717165E-01
2.787574E-02
1.260209E-01
3.218342E-03
1.253684E-01
3.332746E-03
4.032393E-01
3.350955E-02
1.448782E-01
4.379454E-03
5.071496E-02
5.305658E-04
7.987421E-02
1.317997E-03
4.219165E-02
3.699167E-04
5.746706E-02
6.694157E-04
1.140169E-01
2.603241E-03
6.643597E-02
9.208619E-04
2.189665E-01
9.709850E-03
5.827286E-01
6.955166E-02
1.650015E-01
5.490500E-03
1.803242E-01
6.573690E-03
6.254933E-01
7.861927E-02
1.931243E-01
7.508718E-03
6.560024E-02
8.649472E-04
1.026985E-01
2.141614E-03
5.105660E-02
5.415526E-04
5.906447E-02
7.318418E-04
1.148513E-01
2.734426E-03
6.999733E-02
9.823103E-04
2.242530E-01
1.049196E-02
9.228789E-01
2.232918E-01
1.852616E-01
6.872969E-03
1.980949E-01
8.177349E-03
9.355465E-01
2.355627E-01
2.277687E-01
1.049830E-02
7.978124E-02
1.279972E-03
1.150264E-01
2.686454E-03
5.347570E-02
6.120057E-04
5.589535E-02
6.332204E-04
1.122848E-01
2.544997E-03
7.096065E-02
1.015080E-03
2.048604E-01
8.420923E-03
5.431051E-01
5.950426E-02
1.871604E-01
7.046694E-03
1.407324E-01
4.003323E-03
5.578866E-01
6.250709E-02
2.002467E-01
8.086142E-03
6.597658E-02
8.956961E-04
1.026623E-01
2.120442E-03
5.494589E-02
6.199369E-04
4.691398E-02
4.413619E-04
9.072176E-02
1.677400E-03
5.541886E-02
6.672463E-04
1.495874E-01
4.606567E-03
3.649631E-01
2.770976E-02
1.248362E-01
3.405376E-03
1.171804E-01
2.836268E-03
3.258318E-01
2.188203E-02
1.371604E-01
3.846064E-03
4.834961E-02
4.795272E-04
7.603041E-02
1.176251E-03
4.112691E-02
3.447493E-04

View file

@ -1,342 +1,342 @@
k-combined:
7.729082E-01 3.775399E-02
8.023345E-01 6.979356E-02
tally 1:
5.296804E-02
5.661701E-04
8.356446E-02
1.412139E-03
5.041335E-02
5.143568E-04
1.299348E-01
3.467618E-03
3.929702E-01
3.147038E-02
1.379707E-01
3.888484E-03
1.405034E-01
4.473799E-03
3.785796E-01
2.940585E-02
1.422010E-01
4.113723E-03
5.647073E-02
6.735251E-04
7.911154E-02
1.329137E-03
5.160755E-02
5.361448E-04
6.669424E-02
9.090832E-04
1.008621E-01
2.134534E-03
6.808932E-02
9.355993E-04
1.873006E-01
7.135961E-03
6.221575E-01
7.819842E-02
1.856653E-01
6.954762E-03
2.014929E-01
8.327845E-03
5.853251E-01
6.945708E-02
1.709645E-01
5.917124E-03
7.214913E-02
1.058962E-03
1.027720E-01
2.138475E-03
6.099853E-02
7.493941E-04
6.892071E-02
9.630680E-04
1.035459E-01
2.173883E-03
6.973870E-02
9.904237E-04
2.125703E-01
9.112659E-03
9.012205E-01
2.163546E-01
2.066426E-01
8.617414E-03
2.258950E-01
1.039607E-02
9.476792E-01
2.350708E-01
2.225585E-01
1.017898E-02
7.111503E-02
1.036847E-03
1.117012E-01
2.530040E-03
6.870474E-02
9.551035E-04
5.738897E-02
6.699030E-04
9.522335E-02
1.835769E-03
6.570917E-02
8.656870E-04
1.945592E-01
7.593336E-03
5.514753E-01
6.122981E-02
2.144202E-01
9.421739E-03
1.971631E-01
7.944046E-03
6.088996E-01
7.442954E-02
1.965447E-01
7.765628E-03
7.005494E-02
1.012891E-03
1.010633E-01
2.084095E-03
6.145926E-02
7.694351E-04
4.999479E-02
5.129164E-04
7.238243E-02
1.062921E-03
4.902309E-02
4.852193E-04
1.324655E-01
3.642431E-03
3.305312E-01
2.265726E-02
1.332993E-01
3.728385E-03
1.547469E-01
4.894837E-03
3.625944E-01
2.747313E-02
1.435761E-01
4.334405E-03
5.789603E-02
7.065383E-04
7.589559E-02
1.205386E-03
5.210018E-02
5.790843E-04
5.745655E-02
6.712105E-04
7.867782E-02
1.248727E-03
5.868429E-02
7.131793E-04
1.438557E-01
4.308092E-03
3.787094E-01
2.891295E-02
1.374240E-01
3.869017E-03
1.347274E-01
3.834160E-03
3.987709E-01
3.261523E-02
1.334735E-01
3.626900E-03
4.778466E-02
4.648282E-04
7.840849E-02
1.297864E-03
4.916372E-02
4.947407E-04
6.995546E-02
9.876172E-04
1.087891E-01
2.387772E-03
6.400064E-02
8.393797E-04
1.945275E-01
7.631597E-03
5.867134E-01
7.069740E-02
1.792978E-01
6.602522E-03
1.955106E-01
7.732201E-03
6.226639E-01
7.780920E-02
1.796752E-01
6.469604E-03
6.679475E-02
9.049602E-04
1.008432E-01
2.075215E-03
5.586034E-02
6.531573E-04
6.137597E-02
7.755153E-04
1.147821E-01
2.767643E-03
6.809436E-02
9.285053E-04
2.104606E-01
9.134770E-03
9.269139E-01
2.245814E-01
2.021416E-01
8.202826E-03
2.176221E-01
9.854345E-03
9.336697E-01
2.326234E-01
2.095728E-01
8.873161E-03
7.246366E-02
1.052928E-03
1.083862E-01
2.417978E-03
6.342274E-02
8.587365E-04
6.418145E-02
8.295301E-04
1.066727E-01
2.322556E-03
6.830911E-02
9.380558E-04
1.832547E-01
6.752214E-03
5.428015E-01
5.938368E-02
1.951315E-01
7.663620E-03
1.572668E-01
5.032070E-03
5.625230E-01
6.367739E-02
1.918204E-01
7.408678E-03
6.451833E-02
8.482285E-04
9.426618E-02
1.787793E-03
6.648279E-02
9.126838E-04
5.317949E-02
5.712761E-04
8.572058E-02
1.489821E-03
5.734607E-02
6.866788E-04
1.365285E-01
3.892455E-03
3.670737E-01
2.795038E-02
1.429659E-01
4.492430E-03
1.222604E-01
3.045332E-03
3.261367E-01
2.184800E-02
1.255502E-01
3.263101E-03
4.699926E-02
4.464966E-04
7.180560E-02
1.057221E-03
5.028875E-02
5.175287E-04
tally 2:
2.572693E-01
1.338921E-02
2.567576E-01
1.333558E-02
6.127899E-01
7.781914E-02
6.163075E-01
7.600273E-02
2.566545E-01
1.372553E-02
2.509744E-01
1.286275E-02
3.366720E-01
2.282560E-02
3.216175E-01
2.103222E-02
1.118187E+00
2.697988E-01
1.050003E+00
2.391098E-01
3.407360E-01
2.332214E-02
3.150094E-01
2.006040E-02
3.130223E-01
2.001074E-02
3.238651E-01
2.101351E-02
1.061186E+00
2.372828E-01
1.099461E+00
2.597134E-01
3.445524E-01
2.400393E-02
3.429055E-01
2.372392E-02
2.299776E-01
1.064851E-02
2.208901E-01
9.829109E-03
6.015657E-01
7.399276E-02
5.851229E-01
7.165716E-02
2.583609E-01
1.364046E-02
2.456380E-01
1.257735E-02
2.670690E-01
1.455424E-02
2.447172E-01
1.234981E-02
6.320025E-01
8.177928E-02
6.127544E-01
7.527645E-02
2.397989E-01
1.166167E-02
2.369379E-01
1.189279E-02
3.579293E-01
2.610406E-02
3.210301E-01
2.076744E-02
1.097889E+00
2.664104E-01
1.043586E+00
2.338589E-01
3.326713E-01
2.225475E-02
3.338383E-01
2.271154E-02
3.220172E-01
2.104639E-02
3.416243E-01
2.354218E-02
9.868052E-01
2.128044E-01
1.091708E+00
2.578225E-01
3.417944E-01
2.354436E-02
3.053713E-01
1.901696E-02
2.589156E-01
1.364023E-02
2.619494E-01
1.399573E-02
5.532464E-01
6.269623E-02
5.468628E-01
6.098899E-02
2.270733E-01
1.040737E-02
2.361659E-01
1.129872E-02
tally 3:
5.296804E-02
5.661701E-04
8.356446E-02
1.412139E-03
5.041335E-02
5.143568E-04
1.299348E-01
3.467618E-03
3.929702E-01
3.147038E-02
1.379707E-01
3.888484E-03
1.405034E-01
4.473799E-03
3.785796E-01
2.940585E-02
1.422010E-01
4.113723E-03
5.647073E-02
6.735251E-04
7.911154E-02
1.329137E-03
5.160755E-02
5.361448E-04
6.669424E-02
9.090832E-04
1.008621E-01
2.134534E-03
6.808932E-02
9.355993E-04
1.873006E-01
7.135961E-03
6.221575E-01
7.819842E-02
1.856653E-01
6.954762E-03
2.014929E-01
8.327845E-03
5.853251E-01
6.945708E-02
1.709645E-01
5.917124E-03
7.214913E-02
1.058962E-03
1.027720E-01
2.138475E-03
6.099853E-02
7.493941E-04
6.892071E-02
9.630680E-04
1.035459E-01
2.173883E-03
6.973870E-02
9.904237E-04
2.125703E-01
9.112659E-03
9.012205E-01
2.163546E-01
2.066426E-01
8.617414E-03
2.258950E-01
1.039607E-02
9.476792E-01
2.350708E-01
2.225585E-01
1.017898E-02
7.111503E-02
1.036847E-03
1.117012E-01
2.530040E-03
6.870474E-02
9.551035E-04
5.738897E-02
6.699030E-04
9.522335E-02
1.835769E-03
6.570917E-02
8.656870E-04
1.945592E-01
7.593336E-03
5.514753E-01
6.122981E-02
2.144202E-01
9.421739E-03
1.971631E-01
7.944046E-03
6.088996E-01
7.442954E-02
1.965447E-01
7.765628E-03
7.005494E-02
1.012891E-03
1.010633E-01
2.084095E-03
6.145926E-02
7.694351E-04
4.999479E-02
5.129164E-04
7.238243E-02
1.062921E-03
4.902309E-02
4.852193E-04
1.324655E-01
3.642431E-03
3.305312E-01
2.265726E-02
1.332993E-01
3.728385E-03
1.547469E-01
4.894837E-03
3.625944E-01
2.747313E-02
1.435761E-01
4.334405E-03
5.789603E-02
7.065383E-04
7.589559E-02
1.205386E-03
5.210018E-02
5.790843E-04
5.745655E-02
6.712105E-04
7.867782E-02
1.248727E-03
5.868429E-02
7.131793E-04
1.438557E-01
4.308092E-03
3.787094E-01
2.891295E-02
1.374240E-01
3.869017E-03
1.347274E-01
3.834160E-03
3.987709E-01
3.261523E-02
1.334735E-01
3.626900E-03
4.778466E-02
4.648282E-04
7.840849E-02
1.297864E-03
4.916372E-02
4.947407E-04
6.995546E-02
9.876172E-04
1.087891E-01
2.387772E-03
6.400064E-02
8.393797E-04
1.945275E-01
7.631597E-03
5.867134E-01
7.069740E-02
1.792978E-01
6.602522E-03
1.955106E-01
7.732201E-03
6.226639E-01
7.780920E-02
1.796752E-01
6.469604E-03
6.679475E-02
9.049602E-04
1.008432E-01
2.075215E-03
5.586034E-02
6.531573E-04
6.137597E-02
7.755153E-04
1.147821E-01
2.767643E-03
6.809436E-02
9.285053E-04
2.104606E-01
9.134770E-03
9.269139E-01
2.245814E-01
2.021416E-01
8.202826E-03
2.176221E-01
9.854345E-03
9.336697E-01
2.326234E-01
2.095728E-01
8.873161E-03
7.246366E-02
1.052928E-03
1.083862E-01
2.417978E-03
6.342274E-02
8.587365E-04
6.418145E-02
8.295301E-04
1.066727E-01
2.322556E-03
6.830911E-02
9.380558E-04
1.832547E-01
6.752214E-03
5.428015E-01
5.938368E-02
1.951315E-01
7.663620E-03
1.572668E-01
5.032070E-03
5.625230E-01
6.367739E-02
1.918204E-01
7.408678E-03
6.451833E-02
8.482285E-04
9.426618E-02
1.787793E-03
6.648279E-02
9.126838E-04
5.317949E-02
5.712761E-04
8.572058E-02
1.489821E-03
5.734607E-02
6.866788E-04
1.365285E-01
3.892455E-03
3.670737E-01
2.795038E-02
1.429659E-01
4.492430E-03
1.222604E-01
3.045332E-03
3.261367E-01
2.184800E-02
1.255502E-01
3.263101E-03
4.699926E-02
4.464966E-04
7.180560E-02
1.057221E-03
5.028875E-02
5.175287E-04
tally 4:
2.572693E-01
1.338921E-02
2.567576E-01
1.333558E-02
6.127899E-01
7.781914E-02
6.163075E-01
7.600273E-02
2.566545E-01
1.372553E-02
2.509744E-01
1.286275E-02
3.366720E-01
2.282560E-02
3.216175E-01
2.103222E-02
1.118187E+00
2.697988E-01
1.050003E+00
2.391098E-01
3.407360E-01
2.332214E-02
3.150094E-01
2.006040E-02
3.130223E-01
2.001074E-02
3.238651E-01
2.101351E-02
1.061186E+00
2.372828E-01
1.099461E+00
2.597134E-01
3.445524E-01
2.400393E-02
3.429055E-01
2.372392E-02
2.299776E-01
1.064851E-02
2.208901E-01
9.829109E-03
6.015657E-01
7.399276E-02
5.851229E-01
7.165716E-02
2.583609E-01
1.364046E-02
2.456380E-01
1.257735E-02
2.670690E-01
1.455424E-02
2.447172E-01
1.234981E-02
6.320025E-01
8.177928E-02
6.127544E-01
7.527645E-02
2.397989E-01
1.166167E-02
2.369379E-01
1.189279E-02
3.579293E-01
2.610406E-02
3.210301E-01
2.076744E-02
1.097889E+00
2.664104E-01
1.043586E+00
2.338589E-01
3.326713E-01
2.225475E-02
3.338383E-01
2.271154E-02
3.220172E-01
2.104639E-02
3.416243E-01
2.354218E-02
9.868052E-01
2.128044E-01
1.091708E+00
2.578225E-01
3.417944E-01
2.354436E-02
3.053713E-01
1.901696E-02
2.589156E-01
1.364023E-02
2.619494E-01
1.399573E-02
5.532464E-01
6.269623E-02
5.468628E-01
6.098899E-02
2.270733E-01
1.040737E-02
2.361659E-01
1.129872E-02

View file

@ -49,10 +49,10 @@ tally 1:
9.089898E-06
7.058968E-03
4.982903E-05
9.013188E-04
8.123755E-07
7.501461E-04
5.627191E-07
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
5.134898E-03
@ -69,74 +69,74 @@ tally 1:
1.611949E-05
2.698215E-03
7.280363E-06
1.751215E-02
1.073504E-04
1.356908E-02
9.589157E-05
1.171344E-02
7.372535E-05
5.400405E-03
2.916437E-05
5.451466E-03
2.971848E-05
3.165676E-04
1.002151E-07
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
6.400753E-04
4.096964E-07
0.000000E+00
0.000000E+00
3.724938E-03
1.387516E-05
2.812863E-03
7.278837E-06
1.000700E+01
1.001402E+01
2.345904E-02
2.127641E-04
2.451653E-02
1.750725E-04
8.695904E-03
5.553981E-05
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.000609E+01
1.001220E+01
1.991905E-02
1.933708E-04
2.154398E-02
1.662365E-04
1.370500E-03
1.878271E-06
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.051159E-03
1.641189E-05
3.590490E-03
1.289162E-05
0.000000E+00
0.000000E+00
4.947565E-03
2.447840E-05
1.725721E-02
8.878803E-05
1.242587E-02
6.544613E-05
5.656444E+01
3.199538E+02
2.159961E-02
1.034364E-04
3.091063E-02
5.543317E-04
4.232209E-03
1.791159E-05
9.266696E-03
8.587166E-05
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.716170E-03
2.224226E-05
1.321855E-02
5.288265E-05
1.104038E-03
1.218900E-06
4.682079E-03
2.192187E-05
8.704920E-03
3.815239E-05
7.790475E-03
2.341864E-05
5.309576E-02
7.984865E-04
5.656179E+01
3.199238E+02
2.399311E-02
2.936670E-04
4.680753E-02
9.526509E-04
2.093139E-02
2.830759E-04
1.074699E-02
1.154978E-04
0.000000E+00
@ -145,20 +145,20 @@ tally 1:
0.000000E+00
4.526425E-03
2.048853E-05
3.983003E-03
9.553596E-06
1.726789E-02
9.579848E-05
4.522349E-03
1.296396E-05
1.996841E-02
1.373021E-04
3.355605E-02
3.331438E-04
3.578422E-02
3.560242E-04
5.649209E+01
3.191358E+02
2.591495E-02
4.169753E-04
9.484744E-03
8.065715E-05
1.772102E-02
1.673611E-04
5.653183E-04
2.713258E-07
0.000000E+00
0.000000E+00
4.232587E-03

View file

@ -1,362 +1,362 @@
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 0.103374 0.004981
2 1 2 1 total 0.103852 0.004752
1 2 1 1 total 0.103322 0.003613
3 2 2 1 total 0.102889 0.003387
0 1 1 1 total 0.103392 0.006724
2 1 2 1 total 0.104216 0.003920
1 2 1 1 total 0.102878 0.004544
3 2 2 1 total 0.104559 0.005436
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 0.072760 0.005235
2 1 2 1 total 0.074856 0.004972
1 2 1 1 total 0.074583 0.004000
3 2 2 1 total 0.072925 0.003485
0 1 1 1 total 0.073701 0.006894
2 1 2 1 total 0.072983 0.004236
1 2 1 1 total 0.071789 0.004763
3 2 2 1 total 0.074547 0.005794
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 0.072768 0.005236
2 1 2 1 total 0.074864 0.004972
1 2 1 1 total 0.074603 0.004002
3 2 2 1 total 0.072881 0.003491
0 1 1 1 total 0.073688 0.006895
2 1 2 1 total 0.072983 0.004236
1 2 1 1 total 0.071778 0.004762
3 2 2 1 total 0.074483 0.005805
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 0.013309 0.000729
2 1 2 1 total 0.013223 0.000707
1 2 1 1 total 0.013222 0.000596
3 2 2 1 total 0.013282 0.000564
0 1 1 1 total 0.013218 0.000943
2 1 2 1 total 0.013285 0.000674
1 2 1 1 total 0.013050 0.000695
3 2 2 1 total 0.013457 0.000810
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 0.013241 0.000728
2 1 2 1 total 0.013142 0.000705
1 2 1 1 total 0.013137 0.000596
3 2 2 1 total 0.013221 0.000564
0 1 1 1 total 0.013155 0.000942
2 1 2 1 total 0.013227 0.000674
1 2 1 1 total 0.012988 0.000694
3 2 2 1 total 0.013397 0.000809
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 0.001281 0.000874
2 1 2 1 total 0.001297 0.000731
1 2 1 1 total 0.001281 0.000744
3 2 2 1 total 0.001288 0.000719
0 1 1 1 total 0.001277 0.001050
2 1 2 1 total 0.001274 0.000849
1 2 1 1 total 0.001252 0.000839
3 2 2 1 total 0.001301 0.000901
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 0.012029 0.000664
2 1 2 1 total 0.011926 0.000639
1 2 1 1 total 0.011941 0.000545
3 2 2 1 total 0.011994 0.000515
0 1 1 1 total 0.011941 0.000847
2 1 2 1 total 0.012011 0.000598
1 2 1 1 total 0.011798 0.000634
3 2 2 1 total 0.012156 0.000736
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 0.031361 0.001735
2 1 2 1 total 0.031091 0.001675
1 2 1 1 total 0.031169 0.001418
3 2 2 1 total 0.031255 0.001366
0 1 1 1 total 0.031158 0.002189
2 1 2 1 total 0.031301 0.001533
1 2 1 1 total 0.030704 0.001668
3 2 2 1 total 0.031644 0.001953
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 2.326447e+06 128504.886535
2 1 2 1 total 2.306518e+06 123638.096130
1 2 1 1 total 2.309435e+06 105395.059055
3 2 2 1 total 2.319667e+06 99692.620001
0 1 1 1 total 2.309421e+06 163809.243792
2 1 2 1 total 2.322988e+06 115559.889881
1 2 1 1 total 2.281847e+06 122611.604752
3 2 2 1 total 2.351060e+06 142256.311927
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 0.090065 0.004262
2 1 2 1 total 0.090629 0.004086
1 2 1 1 total 0.090100 0.003045
3 2 2 1 total 0.089607 0.002828
0 1 1 1 total 0.090175 0.005785
2 1 2 1 total 0.090931 0.003266
1 2 1 1 total 0.089828 0.003877
3 2 2 1 total 0.091102 0.004654
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 0.090188 0.005151
2 1 2 1 total 0.094339 0.004349
1 2 1 1 total 0.088294 0.003953
3 2 2 1 total 0.089633 0.002592
0 1 1 1 total 0.089509 0.005190
2 1 2 1 total 0.089214 0.003829
1 2 1 1 total 0.090843 0.004007
3 2 2 1 total 0.089166 0.004904
mesh 1 group in group out legendre nuclide mean std. dev.
x y
0 1 1 1 1 P0 total 0.090122 0.005149
1 1 1 1 1 P1 total 0.030614 0.001612
2 1 1 1 1 P2 total 0.016490 0.000897
3 1 1 1 1 P3 total 0.010163 0.000829
8 1 2 1 1 P0 total 0.094307 0.004337
9 1 2 1 1 P1 total 0.028996 0.001462
10 1 2 1 1 P2 total 0.017486 0.000808
11 1 2 1 1 P3 total 0.009753 0.000742
4 2 1 1 1 P0 total 0.088198 0.003959
5 2 1 1 1 P1 total 0.028739 0.001715
6 2 1 1 1 P2 total 0.015608 0.000680
7 2 1 1 1 P3 total 0.008232 0.000363
12 2 2 1 1 P0 total 0.089536 0.002551
13 2 2 1 1 P1 total 0.029964 0.000821
14 2 2 1 1 P2 total 0.015742 0.001260
15 2 2 1 1 P3 total 0.008944 0.000385
0 1 1 1 1 P0 total 0.089477 0.005189
1 1 1 1 1 P1 total 0.029692 0.001522
2 1 1 1 1 P2 total 0.016160 0.001232
3 1 1 1 1 P3 total 0.009278 0.001253
8 1 2 1 1 P0 total 0.089214 0.003829
9 1 2 1 1 P1 total 0.031233 0.001604
10 1 2 1 1 P2 total 0.017305 0.001152
11 1 2 1 1 P3 total 0.010357 0.000432
4 2 1 1 1 P0 total 0.090810 0.004019
5 2 1 1 1 P1 total 0.031089 0.001428
6 2 1 1 1 P2 total 0.016755 0.000688
7 2 1 1 1 P3 total 0.008554 0.000830
12 2 2 1 1 P0 total 0.089068 0.004878
13 2 2 1 1 P1 total 0.030012 0.002004
14 2 2 1 1 P2 total 0.016990 0.001267
15 2 2 1 1 P3 total 0.009947 0.000865
mesh 1 group in group out legendre nuclide mean std. dev.
x y
0 1 1 1 1 P0 total 0.090188 0.005151
1 1 1 1 1 P1 total 0.030606 0.001614
2 1 1 1 1 P2 total 0.016462 0.000898
3 1 1 1 1 P3 total 0.010173 0.000824
8 1 2 1 1 P0 total 0.094339 0.004349
9 1 2 1 1 P1 total 0.028988 0.001462
10 1 2 1 1 P2 total 0.017473 0.000811
11 1 2 1 1 P3 total 0.009763 0.000749
4 2 1 1 1 P0 total 0.088294 0.003953
5 2 1 1 1 P1 total 0.028719 0.001720
6 2 1 1 1 P2 total 0.015571 0.000691
7 2 1 1 1 P3 total 0.008255 0.000361
12 2 2 1 1 P0 total 0.089633 0.002592
13 2 2 1 1 P1 total 0.030008 0.000847
14 2 2 1 1 P2 total 0.015748 0.001284
15 2 2 1 1 P3 total 0.008951 0.000387
0 1 1 1 1 P0 total 0.089509 0.005190
1 1 1 1 1 P1 total 0.029704 0.001529
2 1 1 1 1 P2 total 0.016151 0.001230
3 1 1 1 1 P3 total 0.009264 0.001253
8 1 2 1 1 P0 total 0.089214 0.003829
9 1 2 1 1 P1 total 0.031233 0.001604
10 1 2 1 1 P2 total 0.017305 0.001152
11 1 2 1 1 P3 total 0.010357 0.000432
4 2 1 1 1 P0 total 0.090843 0.004007
5 2 1 1 1 P1 total 0.031100 0.001425
6 2 1 1 1 P2 total 0.016745 0.000688
7 2 1 1 1 P3 total 0.008540 0.000824
12 2 2 1 1 P0 total 0.089166 0.004904
13 2 2 1 1 P1 total 0.030076 0.002035
14 2 2 1 1 P2 total 0.017008 0.001271
15 2 2 1 1 P3 total 0.009936 0.000860
mesh 1 group in group out nuclide mean std. dev.
x y
0 1 1 1 1 total 1.000736 0.060260
2 1 2 1 1 total 1.000346 0.042656
1 2 1 1 1 total 1.001091 0.052736
3 2 2 1 1 total 1.001078 0.038212
0 1 1 1 1 total 1.000358 0.064142
2 1 2 1 1 total 1.000000 0.052450
1 2 1 1 1 total 1.000358 0.045966
3 2 2 1 1 total 1.001104 0.058382
mesh 1 group in group out nuclide mean std. dev.
x y
0 1 1 1 1 total 0.031890 0.002523
2 1 2 1 1 total 0.032714 0.001502
1 2 1 1 1 total 0.030414 0.002369
3 2 2 1 1 total 0.031051 0.001521
0 1 1 1 1 total 0.031429 0.002151
2 1 2 1 1 total 0.030726 0.002050
1 2 1 1 1 total 0.031474 0.002368
3 2 2 1 1 total 0.030461 0.001891
mesh 1 group in group out nuclide mean std. dev.
x y
0 1 1 1 1 total 1.0 0.060234
2 1 2 1 1 total 1.0 0.042523
1 2 1 1 1 total 1.0 0.052777
3 2 2 1 1 total 1.0 0.037842
0 1 1 1 1 total 1.0 0.064127
2 1 2 1 1 total 1.0 0.052450
1 2 1 1 1 total 1.0 0.046096
3 2 2 1 1 total 1.0 0.058105
mesh 1 group in group out legendre nuclide mean std. dev.
x y
0 1 1 1 1 P0 total 0.090065 0.006899
1 1 1 1 1 P1 total 0.030595 0.002243
2 1 1 1 1 P2 total 0.016480 0.001229
3 1 1 1 1 P3 total 0.010157 0.000977
8 1 2 1 1 P0 total 0.090629 0.005616
9 1 2 1 1 P1 total 0.027865 0.001820
10 1 2 1 1 P2 total 0.016804 0.001044
11 1 2 1 1 P3 total 0.009373 0.000813
4 2 1 1 1 P0 total 0.090100 0.005647
5 2 1 1 1 P1 total 0.029359 0.002172
6 2 1 1 1 P2 total 0.015944 0.000984
7 2 1 1 1 P3 total 0.008410 0.000523
12 2 2 1 1 P0 total 0.089607 0.004415
13 2 2 1 1 P1 total 0.029988 0.001459
14 2 2 1 1 P2 total 0.015755 0.001411
15 2 2 1 1 P3 total 0.008951 0.000528
0 1 1 1 1 P0 total 0.090175 0.008180
1 1 1 1 1 P1 total 0.029923 0.002590
2 1 1 1 1 P2 total 0.016286 0.001683
3 1 1 1 1 P3 total 0.009351 0.001421
8 1 2 1 1 P0 total 0.090931 0.005780
9 1 2 1 1 P1 total 0.031834 0.002214
10 1 2 1 1 P2 total 0.017639 0.001436
11 1 2 1 1 P3 total 0.010556 0.000663
4 2 1 1 1 P0 total 0.089828 0.005672
5 2 1 1 1 P1 total 0.030753 0.001978
6 2 1 1 1 P2 total 0.016574 0.001010
7 2 1 1 1 P3 total 0.008461 0.000905
12 2 2 1 1 P0 total 0.091102 0.007049
13 2 2 1 1 P1 total 0.030698 0.002649
14 2 2 1 1 P2 total 0.017378 0.001607
15 2 2 1 1 P3 total 0.010174 0.001045
mesh 1 group in group out legendre nuclide mean std. dev.
x y
0 1 1 1 1 P0 total 0.090131 0.008782
1 1 1 1 1 P1 total 0.030617 0.002905
2 1 1 1 1 P2 total 0.016492 0.001581
3 1 1 1 1 P3 total 0.010164 0.001154
8 1 2 1 1 P0 total 0.090661 0.006820
9 1 2 1 1 P1 total 0.027875 0.002174
10 1 2 1 1 P2 total 0.016810 0.001267
11 1 2 1 1 P3 total 0.009376 0.000906
4 2 1 1 1 P0 total 0.090199 0.007385
5 2 1 1 1 P1 total 0.029391 0.002669
6 2 1 1 1 P2 total 0.015962 0.001295
7 2 1 1 1 P3 total 0.008419 0.000686
12 2 2 1 1 P0 total 0.089704 0.005591
13 2 2 1 1 P1 total 0.030020 0.001856
14 2 2 1 1 P2 total 0.015772 0.001536
15 2 2 1 1 P3 total 0.008961 0.000629
0 1 1 1 1 P0 total 0.090207 0.010021
1 1 1 1 1 P1 total 0.029934 0.003225
2 1 1 1 1 P2 total 0.016292 0.001981
3 1 1 1 1 P3 total 0.009354 0.001543
8 1 2 1 1 P0 total 0.090931 0.007494
9 1 2 1 1 P1 total 0.031834 0.002773
10 1 2 1 1 P2 total 0.017639 0.001709
11 1 2 1 1 P3 total 0.010556 0.000863
4 2 1 1 1 P0 total 0.089860 0.007018
5 2 1 1 1 P1 total 0.030764 0.002432
6 2 1 1 1 P2 total 0.016580 0.001266
7 2 1 1 1 P3 total 0.008464 0.000985
12 2 2 1 1 P0 total 0.091203 0.008836
13 2 2 1 1 P1 total 0.030732 0.003200
14 2 2 1 1 P2 total 0.017397 0.001902
15 2 2 1 1 P3 total 0.010186 0.001203
mesh 1 group out nuclide mean std. dev.
x y
0 1 1 1 total 1.0 0.098093
2 1 2 1 total 1.0 0.042346
1 2 1 1 total 1.0 0.104352
3 2 2 1 total 1.0 0.067889
0 1 1 1 total 1.0 0.082200
2 1 2 1 total 1.0 0.089260
1 2 1 1 total 1.0 0.097625
3 2 2 1 total 1.0 0.071328
mesh 1 group out nuclide mean std. dev.
x y
0 1 1 1 total 1.0 0.099401
2 1 2 1 total 1.0 0.042085
1 2 1 1 total 1.0 0.104135
3 2 2 1 total 1.0 0.067307
0 1 1 1 total 1.0 0.086574
2 1 2 1 total 1.0 0.087769
1 2 1 1 total 1.0 0.101064
3 2 2 1 total 1.0 0.071547
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 8.547713e-10 3.477664e-11
2 1 2 1 total 9.057083e-10 4.171858e-11
1 2 1 1 total 8.666731e-10 2.371072e-11
3 2 2 1 total 8.532016e-10 1.638929e-11
0 1 1 1 total 8.656883e-10 5.405903e-11
2 1 2 1 total 8.635887e-10 2.325887e-11
1 2 1 1 total 8.656521e-10 3.594341e-11
3 2 2 1 total 8.517154e-10 4.271120e-11
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 0.031163 0.001724
2 1 2 1 total 0.030895 0.001664
1 2 1 1 total 0.030973 0.001409
3 2 2 1 total 0.031058 0.001358
0 1 1 1 total 0.030961 0.002175
2 1 2 1 total 0.031103 0.001523
1 2 1 1 total 0.030510 0.001657
3 2 2 1 total 0.031444 0.001941
mesh 1 group in group out nuclide mean std. dev.
x y
0 1 1 1 1 total 0.031781 0.002541
2 1 2 1 1 total 0.032490 0.001487
1 2 1 1 1 total 0.030274 0.002354
3 2 2 1 1 total 0.030882 0.001500
0 1 1 1 1 total 0.031222 0.002220
2 1 2 1 1 total 0.030464 0.002002
1 2 1 1 1 total 0.031317 0.002427
3 2 2 1 1 total 0.030326 0.001887
mesh 1 group in nuclide mean std. dev.
x y surf
3 1 1 x-max in 1 total 0.1888 0.008218
2 1 1 x-max out 1 total 0.1828 0.008212
3 1 1 x-max in 1 total 0.1798 0.013283
2 1 1 x-max out 1 total 0.1912 0.016344
1 1 1 x-min in 1 total 0.0000 0.000000
0 1 1 x-min out 1 total 0.0000 0.000000
7 1 1 y-max in 1 total 0.1820 0.012534
6 1 1 y-max out 1 total 0.1794 0.015302
7 1 1 y-max in 1 total 0.1764 0.012444
6 1 1 y-max out 1 total 0.1820 0.015843
5 1 1 y-min in 1 total 0.0000 0.000000
4 1 1 y-min out 1 total 0.0000 0.000000
19 1 2 x-max in 1 total 0.1870 0.011696
18 1 2 x-max out 1 total 0.1850 0.012919
19 1 2 x-max in 1 total 0.1848 0.011608
18 1 2 x-max out 1 total 0.1772 0.012741
17 1 2 x-min in 1 total 0.0000 0.000000
16 1 2 x-min out 1 total 0.0000 0.000000
23 1 2 y-max in 1 total 0.0000 0.000000
22 1 2 y-max out 1 total 0.0000 0.000000
21 1 2 y-min in 1 total 0.1794 0.015302
20 1 2 y-min out 1 total 0.1820 0.012534
21 1 2 y-min in 1 total 0.1820 0.015843
20 1 2 y-min out 1 total 0.1764 0.012444
11 2 1 x-max in 1 total 0.0000 0.000000
10 2 1 x-max out 1 total 0.0000 0.000000
9 2 1 x-min in 1 total 0.1828 0.008212
8 2 1 x-min out 1 total 0.1888 0.008218
15 2 1 y-max in 1 total 0.1850 0.011256
14 2 1 y-max out 1 total 0.1870 0.009597
9 2 1 x-min in 1 total 0.1912 0.016344
8 2 1 x-min out 1 total 0.1798 0.013283
15 2 1 y-max in 1 total 0.1714 0.009020
14 2 1 y-max out 1 total 0.1866 0.008812
13 2 1 y-min in 1 total 0.0000 0.000000
12 2 1 y-min out 1 total 0.0000 0.000000
27 2 2 x-max in 1 total 0.0000 0.000000
26 2 2 x-max out 1 total 0.0000 0.000000
25 2 2 x-min in 1 total 0.1850 0.012919
24 2 2 x-min out 1 total 0.1870 0.011696
25 2 2 x-min in 1 total 0.1772 0.012741
24 2 2 x-min out 1 total 0.1848 0.011608
31 2 2 y-max in 1 total 0.0000 0.000000
30 2 2 y-max out 1 total 0.0000 0.000000
29 2 2 y-min in 1 total 0.1870 0.009597
28 2 2 y-min out 1 total 0.1850 0.011256
29 2 2 y-min in 1 total 0.1866 0.008812
28 2 2 y-min out 1 total 0.1714 0.009020
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 4.581283 0.329623
2 1 2 1 total 4.452968 0.295762
1 2 1 1 total 4.469295 0.239674
3 2 2 1 total 4.570894 0.218416
0 1 1 1 total 4.522798 0.423053
2 1 2 1 total 4.567288 0.265089
1 2 1 1 total 4.643242 0.308048
3 2 2 1 total 4.471436 0.347510
mesh 1 group in nuclide mean std. dev.
x y
0 1 1 1 total 4.580746 0.329583
2 1 2 1 total 4.452519 0.295716
1 2 1 1 total 4.468066 0.239680
3 2 2 1 total 4.573657 0.219069
0 1 1 1 total 4.523555 0.423292
2 1 2 1 total 4.567288 0.265089
1 2 1 1 total 4.643970 0.308091
3 2 2 1 total 4.475279 0.348767
mesh 1 delayedgroup group in nuclide mean std. dev.
x y
0 1 1 1 1 total 0.000007 3.812309e-07
1 1 1 2 1 total 0.000036 1.967797e-06
2 1 1 3 1 total 0.000034 1.878630e-06
3 1 1 4 1 total 0.000077 4.212043e-06
4 1 1 5 1 total 0.000031 1.726878e-06
5 1 1 6 1 total 0.000013 7.233832e-07
12 1 2 1 1 total 0.000007 3.663985e-07
13 1 2 2 1 total 0.000035 1.891236e-06
14 1 2 3 1 total 0.000034 1.805539e-06
15 1 2 4 1 total 0.000076 4.048166e-06
16 1 2 5 1 total 0.000031 1.659691e-06
17 1 2 6 1 total 0.000013 6.952388e-07
6 2 1 1 1 total 0.000007 3.123173e-07
7 2 1 2 1 total 0.000035 1.612086e-06
8 2 1 3 1 total 0.000034 1.539037e-06
9 2 1 4 1 total 0.000076 3.450649e-06
10 2 1 5 1 total 0.000031 1.414717e-06
11 2 1 6 1 total 0.000013 5.926201e-07
18 2 2 1 1 total 0.000007 2.946716e-07
19 2 2 2 1 total 0.000036 1.521004e-06
20 2 2 3 1 total 0.000034 1.452083e-06
21 2 2 4 1 total 0.000076 3.255689e-06
22 2 2 5 1 total 0.000031 1.334787e-06
23 2 2 6 1 total 0.000013 5.591374e-07
0 1 1 1 1 total 0.000007 4.871913e-07
1 1 1 2 1 total 0.000035 2.514731e-06
2 1 1 3 1 total 0.000034 2.400781e-06
3 1 1 4 1 total 0.000076 5.382749e-06
4 1 1 5 1 total 0.000031 2.206851e-06
5 1 1 6 1 total 0.000013 9.244423e-07
12 1 2 1 1 total 0.000007 3.452278e-07
13 1 2 2 1 total 0.000036 1.781960e-06
14 1 2 3 1 total 0.000034 1.701213e-06
15 1 2 4 1 total 0.000077 3.814261e-06
16 1 2 5 1 total 0.000031 1.563793e-06
17 1 2 6 1 total 0.000013 6.550675e-07
6 2 1 1 1 total 0.000007 3.621903e-07
7 2 1 2 1 total 0.000035 1.869515e-06
8 2 1 3 1 total 0.000034 1.784801e-06
9 2 1 4 1 total 0.000075 4.001672e-06
10 2 1 5 1 total 0.000031 1.640629e-06
11 2 1 6 1 total 0.000013 6.872538e-07
18 2 2 1 1 total 0.000007 4.170403e-07
19 2 2 2 1 total 0.000036 2.152634e-06
20 2 2 3 1 total 0.000035 2.055091e-06
21 2 2 4 1 total 0.000078 4.607684e-06
22 2 2 5 1 total 0.000032 1.889085e-06
23 2 2 6 1 total 0.000013 7.913312e-07
mesh 1 delayedgroup group out nuclide mean std. dev.
x y
0 1 1 1 1 total 1.0 1.414214
1 1 1 2 1 total 0.0 0.000000
2 1 1 3 1 total 0.0 0.000000
3 1 1 4 1 total 1.0 0.578922
2 1 1 3 1 total 1.0 1.414214
3 1 1 4 1 total 1.0 0.677547
4 1 1 5 1 total 0.0 0.000000
5 1 1 6 1 total 0.0 0.000000
12 1 2 1 1 total 0.0 0.000000
13 1 2 2 1 total 1.0 0.578922
14 1 2 3 1 total 1.0 1.414214
15 1 2 4 1 total 1.0 1.414214
16 1 2 5 1 total 1.0 0.875472
13 1 2 2 1 total 0.0 0.000000
14 1 2 3 1 total 1.0 0.877657
15 1 2 4 1 total 1.0 0.806180
16 1 2 5 1 total 1.0 1.414214
17 1 2 6 1 total 1.0 1.414214
6 2 1 1 1 total 0.0 0.000000
7 2 1 2 1 total 0.0 0.000000
8 2 1 3 1 total 1.0 1.414214
9 2 1 4 1 total 1.0 0.579392
6 2 1 1 1 total 1.0 1.414214
7 2 1 2 1 total 1.0 0.866887
8 2 1 3 1 total 0.0 0.000000
9 2 1 4 1 total 1.0 1.414214
10 2 1 5 1 total 1.0 1.414214
11 2 1 6 1 total 0.0 0.000000
18 2 2 1 1 total 1.0 0.868163
18 2 2 1 1 total 0.0 0.000000
19 2 2 2 1 total 1.0 1.414214
20 2 2 3 1 total 0.0 0.000000
21 2 2 4 1 total 1.0 0.868969
22 2 2 5 1 total 1.0 1.414214
20 2 2 3 1 total 1.0 1.414214
21 2 2 4 1 total 1.0 1.414214
22 2 2 5 1 total 1.0 0.875616
23 2 2 6 1 total 0.0 0.000000
mesh 1 delayedgroup group in nuclide mean std. dev.
x y
0 1 1 1 1 total 0.000221 0.000015
1 1 1 2 1 total 0.001140 0.000079
2 1 1 3 1 total 0.001088 0.000075
3 1 1 4 1 total 0.002440 0.000169
4 1 1 5 1 total 0.001001 0.000069
5 1 1 6 1 total 0.000419 0.000029
12 1 2 1 1 total 0.000221 0.000013
13 1 2 2 1 total 0.001139 0.000066
14 1 2 3 1 total 0.001088 0.000063
15 1 2 4 1 total 0.002439 0.000142
16 1 2 5 1 total 0.001000 0.000058
17 1 2 6 1 total 0.000419 0.000024
6 2 1 1 1 total 0.000220 0.000013
7 2 1 2 1 total 0.001138 0.000067
8 2 1 3 1 total 0.001086 0.000064
9 2 1 4 1 total 0.002435 0.000144
10 2 1 5 1 total 0.000998 0.000059
11 2 1 6 1 total 0.000418 0.000025
18 2 2 1 1 total 0.000221 0.000013
19 2 2 2 1 total 0.001141 0.000066
20 2 2 3 1 total 0.001090 0.000063
21 2 2 4 1 total 0.002443 0.000141
22 2 2 5 1 total 0.001002 0.000058
23 2 2 6 1 total 0.000420 0.000024
0 1 1 1 1 total 0.000221 0.000018
1 1 1 2 1 total 0.001139 0.000094
2 1 1 3 1 total 0.001087 0.000089
3 1 1 4 1 total 0.002438 0.000201
4 1 1 5 1 total 0.001000 0.000082
5 1 1 6 1 total 0.000419 0.000034
12 1 2 1 1 total 0.000221 0.000015
13 1 2 2 1 total 0.001143 0.000075
14 1 2 3 1 total 0.001091 0.000072
15 1 2 4 1 total 0.002447 0.000161
16 1 2 5 1 total 0.001003 0.000066
17 1 2 6 1 total 0.000420 0.000028
6 2 1 1 1 total 0.000222 0.000015
7 2 1 2 1 total 0.001145 0.000078
8 2 1 3 1 total 0.001093 0.000074
9 2 1 4 1 total 0.002451 0.000166
10 2 1 5 1 total 0.001005 0.000068
11 2 1 6 1 total 0.000421 0.000029
18 2 2 1 1 total 0.000222 0.000016
19 2 2 2 1 total 0.001144 0.000081
20 2 2 3 1 total 0.001093 0.000077
21 2 2 4 1 total 0.002450 0.000173
22 2 2 5 1 total 0.001004 0.000071
23 2 2 6 1 total 0.000421 0.000030
mesh 1 delayedgroup nuclide mean std. dev.
x y
0 1 1 1 total 0.013336 0.000919
1 1 1 2 total 0.032739 0.002257
2 1 1 3 total 0.120780 0.008325
3 1 1 4 total 0.302780 0.020871
4 1 1 5 total 0.849490 0.058555
5 1 1 6 total 2.853000 0.196656
12 1 2 1 total 0.013336 0.000770
13 1 2 2 total 0.032739 0.001891
14 1 2 3 total 0.120780 0.006977
15 1 2 4 total 0.302780 0.017490
16 1 2 5 total 0.849490 0.049071
17 1 2 6 total 2.853000 0.164804
6 2 1 1 total 0.013336 0.000790
7 2 1 2 total 0.032739 0.001940
8 2 1 3 total 0.120780 0.007157
9 2 1 4 total 0.302780 0.017942
10 2 1 5 total 0.849490 0.050338
11 2 1 6 total 2.853000 0.169061
18 2 2 1 total 0.013336 0.000757
19 2 2 2 total 0.032739 0.001858
20 2 2 3 total 0.120780 0.006855
21 2 2 4 total 0.302780 0.017186
22 2 2 5 total 0.849490 0.048217
23 2 2 6 total 2.853000 0.161935
0 1 1 1 total 0.013336 0.001104
1 1 1 2 total 0.032739 0.002710
2 1 1 3 total 0.120780 0.009996
3 1 1 4 total 0.302780 0.025059
4 1 1 5 total 0.849490 0.070306
5 1 1 6 total 2.853000 0.236122
12 1 2 1 total 0.013336 0.000884
13 1 2 2 total 0.032739 0.002170
14 1 2 3 total 0.120780 0.008006
15 1 2 4 total 0.302780 0.020070
16 1 2 5 total 0.849490 0.056310
17 1 2 6 total 2.853000 0.189116
6 2 1 1 total 0.013336 0.000893
7 2 1 2 total 0.032739 0.002193
8 2 1 3 total 0.120780 0.008091
9 2 1 4 total 0.302780 0.020283
10 2 1 5 total 0.849490 0.056907
11 2 1 6 total 2.853000 0.191120
18 2 2 1 total 0.013336 0.000915
19 2 2 2 total 0.032739 0.002246
20 2 2 3 total 0.120780 0.008286
21 2 2 4 total 0.302780 0.020773
22 2 2 5 total 0.849490 0.058281
23 2 2 6 total 2.853000 0.195737
mesh 1 delayedgroup group in group out nuclide mean std. dev.
x y
0 1 1 1 1 1 total 0.000026 0.000026
0 1 1 1 1 1 total 0.000025 0.000025
1 1 1 2 1 1 total 0.000000 0.000000
2 1 1 3 1 1 total 0.000000 0.000000
3 1 1 4 1 1 total 0.000083 0.000034
2 1 1 3 1 1 total 0.000025 0.000025
3 1 1 4 1 1 total 0.000157 0.000075
4 1 1 5 1 1 total 0.000000 0.000000
5 1 1 6 1 1 total 0.000000 0.000000
12 1 2 1 1 1 total 0.000000 0.000000
13 1 2 2 1 1 total 0.000081 0.000033
14 1 2 3 1 1 total 0.000026 0.000026
15 1 2 4 1 1 total 0.000026 0.000026
16 1 2 5 1 1 total 0.000059 0.000036
13 1 2 2 1 1 total 0.000000 0.000000
14 1 2 3 1 1 total 0.000058 0.000036
15 1 2 4 1 1 total 0.000139 0.000079
16 1 2 5 1 1 total 0.000033 0.000033
17 1 2 6 1 1 total 0.000033 0.000033
6 2 1 1 1 1 total 0.000000 0.000000
7 2 1 2 1 1 total 0.000000 0.000000
8 2 1 3 1 1 total 0.000032 0.000032
9 2 1 4 1 1 total 0.000080 0.000033
10 2 1 5 1 1 total 0.000028 0.000028
6 2 1 1 1 1 total 0.000025 0.000025
7 2 1 2 1 1 total 0.000056 0.000034
8 2 1 3 1 1 total 0.000000 0.000000
9 2 1 4 1 1 total 0.000050 0.000050
10 2 1 5 1 1 total 0.000025 0.000025
11 2 1 6 1 1 total 0.000000 0.000000
18 2 2 1 1 1 total 0.000054 0.000033
19 2 2 2 1 1 total 0.000029 0.000029
20 2 2 3 1 1 total 0.000000 0.000000
21 2 2 4 1 1 total 0.000054 0.000033
22 2 2 5 1 1 total 0.000032 0.000032
18 2 2 1 1 1 total 0.000000 0.000000
19 2 2 2 1 1 total 0.000025 0.000025
20 2 2 3 1 1 total 0.000025 0.000025
21 2 2 4 1 1 total 0.000026 0.000026
22 2 2 5 1 1 total 0.000059 0.000036
23 2 2 6 1 1 total 0.000000 0.000000

File diff suppressed because it is too large Load diff

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
4.929000E-02 8.212396E-03
4.784265E-02 7.584515E-03

View file

@ -1,2 +1,2 @@
k-combined:
4.929000E-02 8.212396E-03
4.784265E-02 7.584515E-03

View file

@ -1,2 +1,2 @@
k-combined:
4.929000E-02 8.212396E-03
4.784265E-02 7.584515E-03

View file

@ -1,2 +1,2 @@
k-combined:
4.929000E-02 8.212396E-03
4.784265E-02 7.584515E-03

View file

@ -1,2 +1,2 @@
k-combined:
1.403371E+00 1.456192E-02
1.941773E+00 3.249992E-01

View file

@ -1,2 +1,2 @@
k-combined:
1.403371E+00 1.456192E-02
1.941773E+00 3.249992E-01

View file

@ -1,2 +1,2 @@
k-combined:
1.403371E+00 1.456192E-02
1.941773E+00 3.249992E-01

View file

@ -1,2 +1,2 @@
k-combined:
1.403371E+00 1.456192E-02
1.941773E+00 3.249992E-01

View file

@ -1,2 +1,2 @@
k-combined:
1.266853E+00 4.552028E-02
1.518808E+00 1.406725E-02

View file

@ -1,2 +1,2 @@
k-combined:
1.266853E+00 4.552028E-02
1.518808E+00 1.406725E-02

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
9.947197E-01 3.711779E-02
1.015340E+00 5.680758E-02

View file

@ -1,2 +1,2 @@
k-combined:
9.947197E-01 3.711779E-02
1.015340E+00 5.680758E-02

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
5.642735E-02 1.494035E-02
6.085380E-02 2.048137E-03

View file

@ -1,2 +1,2 @@
k-combined:
4.929000E-02 8.212396E-03
4.784265E-02 7.584515E-03

View file

@ -1,6 +1,7 @@
"""Tests for surface flux tallying via flux score + SurfaceFilter."""
import math
import numpy as np
import pytest
import openmc
@ -162,3 +163,44 @@ def test_cellfrom_filter_flux_directional(two_cell_model, run_in_tmpdir):
assert mean_from1 == pytest.approx(1.0)
# No particles cross xmid from cell2 → flux = 0
assert mean_from2 == pytest.approx(0.0)
def test_surface_filter_do_not_tally_virtual_surface_crossing(run_in_tmpdir):
openmc.reset_auto_ids()
model = openmc.Model()
zmin = openmc.ZPlane(z0 = -1.0, surface_id=1, name="plane 1")
zmax = openmc.ZPlane(z0 = 1.0, surface_id=2, name="plane 2")
ymin = openmc.YPlane(y0 = -1.0, surface_id=3, name="plane 3")
ymax = openmc.YPlane(y0 = 1.0, surface_id=4, name="plane 4")
xmin = openmc.XPlane(x0 = -1.0, surface_id=5, name="plane 5")
xmax = openmc.XPlane(x0 = 1.0, surface_id=6, name="plane 6")
sph = openmc.Sphere(r=100.0, boundary_type='vacuum')
cube = (+zmin & -zmax & +ymin & -ymax & +xmin & -xmax )
sphere = -sph&~(cube)
cube_cell = openmc.Cell(region=cube)
sphere_cell = openmc.Cell(region=sphere)
univ = openmc.Universe(cells=[cube_cell, sphere_cell])
model.geometry = openmc.Geometry(univ)
src = openmc.IndependentSource()
src.space = openmc.stats.Point((0.0, 0.0, 0.0))
model.settings.run_mode = 'fixed source'
model.settings.batches = 1
model.settings.particles = 100
model.settings.source = src
surface_filter = openmc.SurfaceFilter([1,2,3,4,5,6])
tally = openmc.Tally()
tally.scores = ["current"]
tally.filters = [surface_filter]
model.tallies = [tally]
model.run(apply_tally_results=True)
current_mean = np.abs(tally.mean.flatten())
# Every particle crosses exit the cube with weight 1, so current = 1.0
assert current_mean.sum() == pytest.approx(1.0, rel=1e-8)