Spatial Fourier Functional Expansion Tally by tjlaboss · Pull Request #3899 · openmc-dev/openmc

tjlaboss · 2026-03-25T00:09:12Z

Description

Addresses Part 1 of #3815 . Implements a new functional expansion tally basis in space: spatial Fourier. This borrows heavily from the implementation of spatial Legendre.

Adds new class SpatialFourierFilter.
Created a new base class SpatialExpansionFilter as a parent of SpatialFourierFilter and SpatialLegendreFilter on the Python side. Moved most of SpatialLegendreFilter to its parent.
Added unit and regression tests.

Checklist

I have performed a self-review of my own code
I have run clang-format (version 18) on any C++ source files (if applicable)
I have followed the style guidelines for Python source files (if applicable)
~~I have made corresponding changes to the documentation (if applicable)~~
I have added tests that prove my fix is effective or that my feature works (if applicable)

Discussion

Should this filter have a flag for only odd (sin) or even (cos) modes?
Is there a better convention for the bin labels?
Should the SpatialExpansionFilter base class be implemented on the C++ side as well?
This does filter does not normalize by the factor of sqrt(2). I believe that the convention for FETs in OpenMC is to normalize in postprocessing (though I'm not sure if this is explicitly established: see Return actual functional representation of funcational expansion tallies in Python API #3419).

Results

Details

Unhashed _results_true.dat_:

k-combined:
9.509528E-01 1.260257E-02
tally 1:
4.544192E+01
4.134369E+02
4.587609E+01
4.210518E+02
4.670908E+01
4.367621E+02
4.622590E+01
4.297902E+02
4.633576E+01
4.295793E+02
tally 2:
4.575713E+01
4.189118E+02
4.607291E+01
4.250499E+02
4.648406E+01
4.324258E+02
4.559699E+01
4.168418E+02
4.602948E+01
4.240500E+02
tally 3:
1.047472E+01
2.217769E+01
1.311403E+01
3.476624E+01
1.190020E+01
2.846511E+01
9.817675E+00
1.944668E+01
9.640395E+00
1.866903E+01
1.441632E+01
4.189765E+01
1.152881E+01
2.668579E+01
1.022064E+01
2.094964E+01
9.345487E+00
1.747941E+01
1.455081E+01
4.307116E+01
1.206247E+01
2.931906E+01
1.065944E+01
2.300125E+01
9.887436E+00
1.984094E+01
1.348406E+01
3.690150E+01
1.247311E+01
3.144075E+01
1.021770E+01
2.123148E+01
9.985555E+00
1.999493E+01
1.390450E+01
3.900624E+01
1.232130E+01
3.050514E+01
1.000565E+01
2.022633E+01
tally 4:
0.000000E+00
0.000000E+00
6.791923E+01
1.020235E+03
1.961029E+00
3.845636E+00
8.337637E+00
6.951618E+01
tally 5:
9.510932E-04
1.812975E-07
1.268378E-05
3.224356E-11
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.031713E-03
2.132866E-07
1.378257E-05
3.806312E-11
4.909250E-03
4.830322E-06
1.607239E-04
5.177343E-09
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
5.561199E-03
6.196522E-06
1.811530E-04
6.575150E-09
4.686797E-03
4.402486E-06
5.660713E-04
6.422262E-08
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
5.444020E-03
5.937938E-06
6.594664E-04
8.713247E-08
1.050818E-02
2.213101E-05
3.181668E-03
2.028876E-06
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.273369E-02
3.248536E-05
3.902110E-03
3.050526E-06
4.308207E-03
3.719966E-06
3.659779E-03
2.684451E-06
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
5.775627E-03
6.682975E-06
4.989193E-03
4.986948E-06
1.804693E-03
6.527584E-07
5.148790E-03
5.313198E-06
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.400355E-03
1.154307E-06
6.964737E-03
9.718121E-06
tally 6:
4.561158E+01
4.188739E+02
4.571003E+01
4.183680E+02
5.346481E+01
5.721449E+02
9.615221E+00
1.855790E+01
tally 7:
4.524000E+01
4.118471E+02
4.559250E+01
4.160547E+02
5.197000E+01
5.407097E+02
6.352500E+00
8.133144E+00
tally 8:
4.202750E+01
3.557575E+02
0.000000E+00
0.000000E+00
8.000000E-02
1.487500E-03
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.057925E+00
2.242740E-01
0.000000E+00
0.000000E+00
2.276201E+00
1.044393E+00
3.210000E+00
2.062475E+00
0.000000E+00
0.000000E+00
4.096750E+01
3.359720E+02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.455522E-01
1.261673E-02
0.000000E+00
0.000000E+00
5.425338E-01
6.028348E-02
2.500000E-03
6.250000E-06
0.000000E+00
0.000000E+00
4.545000E+00
4.132288E+00
0.000000E+00
0.000000E+00
4.880750E+01
4.768958E+02
3.215509E-02
2.456414E-04
0.000000E+00
0.000000E+00
7.083345E-02
1.082307E-03
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.162500E+00
2.005156E+00
1.477679E-01
4.541066E-03
6.352500E+00
8.133144E+00
2.695133E-01
1.513303E-02
tally 9:
2.947536E+01
1.738278E+02
7.044852E+00
9.933902E+00
5.937075E+01
7.093035E+02
5.204605E+01
5.422926E+02
tally 10:
1.489500E+01
4.440685E+01
1.489500E+01
4.440685E+01
1.754000E+01
6.169969E+01
1.754250E+01
6.171787E+01
3.721000E+01
2.772030E+02
3.721000E+01
2.772030E+02
7.951000E+01
1.265991E+03
7.951000E+01
1.265991E+03
tally 11:
3.110000E+00
1.945800E+00
3.110000E+00
1.945800E+00
4.395000E+00
3.868588E+00
4.395000E+00
3.868588E+00
4.070000E+00
3.348888E+00
4.070000E+00
3.348888E+00
3.200000E+00
2.056688E+00
3.200000E+00
2.056688E+00
3.737500E+00
2.817369E+00
3.737500E+00
2.817369E+00
5.297500E+00
5.621269E+00
5.297500E+00
5.621269E+00
4.612500E+00
4.313169E+00
4.612500E+00
4.313169E+00
3.747500E+00
2.825719E+00
3.750000E+00
2.829813E+00
7.965000E+00
1.278139E+01
7.965000E+00
1.278139E+01
1.139000E+01
2.607964E+01
1.139000E+01
2.607964E+01
9.760000E+00
1.919485E+01
9.760000E+00
1.919485E+01
7.872500E+00
1.247288E+01
7.872500E+00
1.247288E+01
1.704500E+01
5.842291E+01
1.704500E+01
5.842291E+01
2.361750E+01
1.125088E+02
2.361750E+01
1.125088E+02
2.114000E+01
8.950206E+01
2.114000E+01
8.950206E+01
1.737000E+01
6.055176E+01
1.737000E+01
6.055176E+01
tally 12:
2.276723E+01
1.039193E+02
5.854642E+01
6.862265E+02
7.129457E+01
1.018744E+03
5.751055E+01
6.616700E+02
2.046997E+01
8.402181E+01
tally 13:
2.219625E+01
9.860830E+01
5.908079E+01
6.994068E+02
7.086900E+01
1.005632E+03
5.677892E+01
6.453639E+02
2.101561E+01
8.844643E+01
tally 14:
4.879514E+00
4.816471E+00
6.813531E+00
9.427081E+00
6.064033E+00
7.435738E+00
4.956988E+00
5.074194E+00
1.254725E+01
3.161254E+01
1.729706E+01
6.049359E+01
1.537356E+01
4.752235E+01
1.314206E+01
3.458236E+01
1.559803E+01
4.879859E+01
2.137329E+01
9.282614E+01
1.880012E+01
7.094898E+01
1.535340E+01
4.729341E+01
1.205983E+01
2.923369E+01
1.756018E+01
6.243520E+01
1.485848E+01
4.439850E+01
1.292792E+01
3.360042E+01
4.248974E+00
3.677725E+00
6.425662E+00
8.394859E+00
5.189692E+00
5.423136E+00
4.540746E+00
4.176273E+00
tally 15:
1.491550E+02
4.453432E+03
1.491575E+02
4.453578E+03
5.730081E+01
6.578885E+02
5.730055E+01
6.578828E+02
2.052647E+01
8.446977E+01
2.052526E+01
8.446050E+01
2.179001E+00
9.833680E-01
2.179384E+00
9.836111E-01
-2.033562E+00
8.651558E-01
-2.032725E+00
8.643578E-01
tally 16:
2.297048E+02
1.056002E+04
8.436346E+00
1.964815E+01
-3.585165E+01
2.600872E+02
-6.324897E+00
8.966235E+00
-8.771532E+00
1.568559E+01
tally 17:
2.297048E+02
1.056002E+04
-4.932157E+01
4.946025E+02
-1.516594E+01
5.579878E+01
-2.988884E+01
1.885191E+02
-3.250576E+00
1.280698E+01
-9.311595E+00
5.136714E+01
7.672533E+00
1.796054E+01
-1.005586E+01
2.222811E+01
-5.632241E+00
1.070039E+01
-1.288737E+01
3.724165E+01
1.635532E+01
6.213608E+01
tally 18:
1.491550E+02
4.453432E+03
1.491575E+02
4.453578E+03
2.299406E+02
1.058178E+04
1.541575E+02
4.756893E+03
-3.241616E-01
1.043689E-01
-3.225711E-01
1.043060E-01
-2.534571E-01
1.590360E-01
-3.142774E-01
1.051304E-01
2.910549E-01
8.081484E-02
2.906772E-01
8.074167E-02
2.125905E+00
1.431960E+00
3.476996E-01
7.756024E-02
3.583372E-01
1.155094E-01
3.564458E-01
1.155993E-01
-3.051202E-01
5.241652E-02
3.480452E-01
8.739425E-02
-1.193142E-01
2.331107E-02
-1.213984E-01
2.362402E-02
-1.046787E-01
3.203498E-01
-6.440058E-02
1.344148E-02
2.603361E-01
5.538613E-02
2.599198E-01
5.533392E-02
1.036723E+00
3.735011E-01
1.962834E-01
4.857429E-02
-1.250932E-02
4.577743E-02
-1.367369E-02
4.569515E-02
-2.957782E-01
3.174211E-01
-1.012426E-01
6.825028E-02
-3.236805E-01
9.291578E-02
-3.231854E-01
9.270339E-02
-5.504895E-01
3.940803E-01
-3.749395E-01
9.086805E-02
2.082093E-01
3.535544E-02
2.085722E-01
3.538224E-02
9.334556E-01
4.896500E-01
2.583522E-01
6.244720E-02
-9.270389E-02
6.126815E-02
-9.105364E-02
6.165724E-02
-1.560886E-01
1.181402E-01
-1.326451E-01
6.452274E-02
-3.626055E-02
8.257324E-02
-3.555628E-02
8.256519E-02
5.412889E-01
1.399070E-01
-2.325259E-02
7.555563E-02
-1.590580E-03
2.873992E-02
-2.453385E-03
2.892563E-02
2.575346E-01
3.828925E-01
7.377688E-02
4.442002E-02
4.158547E-01
1.301465E-01
4.163998E-01
1.302149E-01
-1.125162E-01
4.917395E-01
4.288947E-01
1.239845E-01
-4.843101E-02
1.084864E-02
-4.740500E-02
1.084966E-02
4.226017E-01
9.888255E-02
-6.166814E-02
1.468814E-02
-1.844615E-02
3.110642E-02
-1.856876E-02
3.109377E-02
-5.131375E-01
9.662145E-02
-5.256893E-02
2.558826E-02
-3.523569E-01
6.730738E-02
-3.513970E-01
6.709665E-02
-7.074196E-01
2.913525E-01
-3.787501E-01
6.876670E-02
1.722060E-01
5.172088E-02
1.716094E-01
5.152236E-02
1.058503E+00
4.523600E-01
1.716184E-01
4.144209E-02
5.450104E-03
3.601203E-02
4.790328E-03
3.581970E-02
-1.065336E-01
2.681613E-01
-1.697541E-02
3.869111E-02
8.231673E-02
7.019014E-02
8.344706E-02
7.000945E-02
3.296209E-01
1.988737E-01
9.360732E-02
7.338492E-02
-1.259676E-01
3.897886E-03
-1.254279E-01
3.874917E-03
-6.599053E-02
1.032379E-01
-1.499346E-01
6.353538E-03
-3.227770E-01
4.186141E-02
-3.220479E-01
4.189754E-02
-7.540375E-01
3.582144E-01
-3.917373E-01
5.240141E-02
1.525495E-01
3.145077E-02
1.519077E-01
3.145980E-02
-1.412429E-01
4.056226E-02
1.924999E-01
4.230709E-02
-8.737337E-02
8.943151E-02
-8.757017E-02
8.936687E-02
1.829284E-01
3.681819E-01
-1.233220E-01
9.668452E-02
-1.386761E-02
2.774899E-02
-1.425138E-02
2.767738E-02
2.251663E-01
1.480266E-01
5.488943E-03
2.674471E-02
-1.200588E-01
1.002133E-01
-1.217202E-01
1.005802E-01
1.232881E-01
3.141823E-01
-1.247928E-01
1.058144E-01
tally 19:
2.299406E+02
1.058178E+04
1.541575E+02
4.756893E+03
-2.534571E-01
1.590360E-01
-3.142774E-01
1.051304E-01
2.125905E+00
1.431960E+00
3.476996E-01
7.756024E-02
-3.051202E-01
5.241652E-02
3.480452E-01
8.739425E-02
-1.046787E-01
3.203498E-01
-6.440058E-02
1.344148E-02
1.036723E+00
3.735011E-01
1.962834E-01
4.857429E-02
-2.957782E-01
3.174211E-01
-1.012426E-01
6.825028E-02
-5.504895E-01
3.940803E-01
-3.749395E-01
9.086805E-02
9.334556E-01
4.896500E-01
2.583522E-01
6.244720E-02
-1.560886E-01
1.181402E-01
-1.326451E-01
6.452274E-02
5.412889E-01
1.399070E-01
-2.325259E-02
7.555563E-02
2.575346E-01
3.828925E-01
7.377688E-02
4.442002E-02
-1.125162E-01
4.917395E-01
4.288947E-01
1.239845E-01
4.226017E-01
9.888255E-02
-6.166814E-02
1.468814E-02
-5.131375E-01
9.662145E-02
-5.256893E-02
2.558826E-02
-7.074196E-01
2.913525E-01
-3.787501E-01
6.876670E-02
1.058503E+00
4.523600E-01
1.716184E-01
4.144209E-02
-1.065336E-01
2.681613E-01
-1.697541E-02
3.869111E-02
3.296209E-01
1.988737E-01
9.360732E-02
7.338492E-02
-6.599053E-02
1.032379E-01
-1.499346E-01
6.353538E-03
-7.540375E-01
3.582144E-01
-3.917373E-01
5.240141E-02
-1.412429E-01
4.056226E-02
1.924999E-01
4.230709E-02
1.829284E-01
3.681819E-01
-1.233220E-01
9.668452E-02
2.251663E-01
1.480266E-01
5.488943E-03
2.674471E-02
1.232881E-01
3.141823E-01
-1.247928E-01
1.058144E-01
tally 20:
2.305888E+02
1.064078E+04
1.544016E+02
4.772042E+03
-8.591898E-01
1.131788E+00
-6.321030E-01
2.997205E-01
2.283040E+00
1.413501E+00
1.470581E-01
1.217435E-01
-8.974386E-01
1.557385E+00
3.727917E-01
6.344496E-01
3.077991E-01
6.968899E-01
2.208815E-01
8.099539E-02
1.657098E+00
1.586609E+00
2.205466E-01
2.066504E-01
-5.585804E-01
4.495641E-01
5.524431E-03
1.357578E-01
-6.504464E-01
5.854695E-01
-3.003911E-01
8.713589E-02
9.437127E-01
9.564931E-01
5.038755E-01
2.068738E-01
-4.424961E-01
3.335991E-01
-1.053527E-01
7.997562E-02
8.268687E-01
5.389782E-01
1.220661E-01
9.141750E-02
8.781942E-01
1.614827E-01
2.210256E-01
5.369655E-02
9.174355E-01
1.341251E+00
7.857567E-01
4.179630E-01
7.181864E-01
4.933596E-01
1.473982E-01
1.419231E-01
-2.226391E-01
1.358904E-01
2.444261E-02
2.167447E-02
-8.464324E-01
2.648335E-01
-5.440846E-01
8.063111E-02
7.320304E-01
6.539587E-01
8.560126E-02
5.926878E-02
-1.180227E-01
5.316560E-01
-4.047059E-02
1.113558E-01
2.235887E-01
1.334216E-01
3.093397E-01
8.666226E-02
3.037116E-03
2.720869E-01
-1.437483E-01
3.049882E-02
-6.593625E-01
3.724096E-01
-4.359114E-01
9.932295E-02
-4.589440E-01
1.829684E-01
3.351449E-02
4.341308E-02
1.509591E-01
6.479608E-01
-6.245819E-02
1.962339E-01
5.873502E-01
1.761062E-01
8.244269E-02
3.807625E-02
-4.164914E-02
6.309154E-01
1.296661E-01
1.456159E-01
tally 21:
6.585262E+01
8.688513E+02
8.033743E+00
1.301436E+01
6.130500E+01
7.520187E+02
7.040875E+00
9.932881E+00
7.388636E+01
1.094201E+03
6.834587E+01
9.346781E+02
tally 22:
4.967463E+00
4.938885E+00
4.768175E+00
4.549953E+00
4.846100E+00
4.701389E+00
4.593704E+00
4.223051E+00
4.723002E+00
4.462099E+00
tally 23:
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.863118E+00
6.959612E-01
1.705476E-02
5.829204E-05
5.535250E+01
6.146207E+02
9.572182E-01
1.838063E-01
1.467875E-05
4.328234E-11
1.854356E+08
6.897901E+15
3.015022E-03
3.180709E-06
3.364815E-01
2.270778E-02
9.027048E-01
1.636200E-01
2.357652E+00
1.114919E+00
1.278833E+01
3.277354E+01
1.199151E+01
2.881594E+01
1.465145E+01
4.301318E+01
2.340596E+00
1.098853E+00
1.726778E+08
5.981585E+15
1.852040E+08
6.880779E+15
8.647179E-03
1.497891E-05
1.944671E-02
7.868354E-05
0.000000E+00
0.000000E+00
9.107500E+01
1.663868E+03
0.000000E+00
0.000000E+00
5.322652E-06
5.694836E-12
0.000000E+00
0.000000E+00
5.869544E-05
2.519352E-09
2.593983E-02
1.353744E-04
1.941787E-02
7.845807E-05
0.000000E+00
0.000000E+00
3.064185E+00
1.885279E+00
2.975045E+00
1.776995E+00
3.083631E+00
1.909304E+00
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.794079E-01
4.618258E-02
0.000000E+00
0.000000E+00
1.748075E+02
6.126322E+03
0.000000E+00
0.000000E+00
2.779303E-05
1.552490E-10
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.797311E-05
6.148531E-10
1.003831E-01
2.025238E-03
0.000000E+00
0.000000E+00
4.763813E+01
4.546982E+02
4.763687E+01
4.546740E+02
4.811754E+01
4.638933E+02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.447311E+00
1.200639E+00
1.705476E-02
5.829204E-05
3.596400E+02
2.593705E+04
9.572182E-01
1.838063E-01
5.333993E-05
5.721669E-10
1.854356E+08
6.897901E+15
3.078283E-03
3.361020E-06
3.654666E-01
2.679249E-02
1.042153E+00
2.179058E-01
2.357652E+00
1.114919E+00
7.143905E+01
1.022953E+03
7.054186E+01
9.974400E+02
7.388636E+01
1.094201E+03
2.340596E+00
1.098853E+00
1.726778E+08
5.981585E+15
1.852040E+08
6.880779E+15
8.647179E-03
1.497891E-05
tally 24:
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.072418E+00
2.308414E-01
1.407702E-02
3.978985E-05
5.535250E+01
6.146207E+02
8.876160E-01
1.581985E-01
1.467875E-05
4.328234E-11
1.716697E+08
5.917518E+15
5.341501E-05
7.668405E-10
6.344301E-05
8.082762E-10
1.848020E-01
6.844659E-03
2.164241E+00
9.405018E-01
1.593040E-01
5.083825E-03
1.456985E-01
4.252227E-03
1.231722E+00
3.042957E-01
2.150164E+00
9.283068E-01
1.602628E+08
5.157247E+15
1.716554E+08
5.916536E+15
6.595811E-03
8.735488E-06
3.311743E-03
2.462032E-06
0.000000E+00
0.000000E+00
5.535250E+01
6.146207E+02
0.000000E+00
0.000000E+00
1.467875E-05
4.328234E-11
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.564516E-05
2.420531E-10
1.169918E-04
2.739639E-09
0.000000E+00
0.000000E+00
5.425380E+00
5.895416E+00
5.424012E+00
5.892397E+00
5.428692E+00
5.902658E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.863118E+00
6.959612E-01
1.705476E-02
5.829204E-05
5.535250E+01
6.146207E+02
9.572182E-01
1.838063E-01
1.467875E-05
4.328234E-11
1.854356E+08
6.897901E+15
3.015022E-03
3.180709E-06
3.364815E-01
2.270778E-02
9.027048E-01
1.636200E-01
2.357652E+00
1.114919E+00
1.278833E+01
3.277354E+01
1.199151E+01
2.881594E+01
1.465145E+01
4.301318E+01
2.340596E+00
1.098853E+00
1.726778E+08
5.981585E+15
1.852040E+08
6.880779E+15
8.647179E-03
1.497891E-05
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
9.107500E+01
1.663868E+03
0.000000E+00
0.000000E+00
5.322652E-06
5.694836E-12
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
9.107500E+01
1.663868E+03
0.000000E+00
0.000000E+00
5.322652E-06
5.694836E-12
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.944671E-02
7.868354E-05
0.000000E+00
0.000000E+00
9.107500E+01
1.663868E+03
0.000000E+00
0.000000E+00
5.322652E-06
5.694836E-12
0.000000E+00
0.000000E+00
5.869544E-05
2.519352E-09
2.593983E-02
1.353744E-04
1.941787E-02
7.845807E-05
0.000000E+00
0.000000E+00
3.064185E+00
1.885279E+00
2.975045E+00
1.776995E+00
3.083631E+00
1.909304E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.748075E+02
6.126322E+03
0.000000E+00
0.000000E+00
2.779303E-05
1.552490E-10
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.745603E-03
1.802573E-06
0.000000E+00
0.000000E+00
1.748075E+02
6.126322E+03
0.000000E+00
0.000000E+00
2.779303E-05
1.552490E-10
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.240472E-05
2.241945E-10
1.026468E-04
2.110010E-09
0.000000E+00
0.000000E+00
4.760492E+00
4.541066E+00
4.759269E+00
4.538675E+00
4.763238E+00
4.546372E+00
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.794079E-01
4.618258E-02
0.000000E+00
0.000000E+00
1.748075E+02
6.126322E+03
0.000000E+00
0.000000E+00
2.779303E-05
1.552490E-10
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.797311E-05
6.148531E-10
1.003831E-01
2.025238E-03
0.000000E+00
0.000000E+00
4.763813E+01
4.546982E+02
4.763687E+01
4.546740E+02
4.811754E+01
4.638933E+02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.072418E+00
2.308414E-01
1.407702E-02
3.978985E-05
3.596400E+02
2.593705E+04
8.876160E-01
1.581985E-01
5.333993E-05
5.721669E-10
1.716697E+08
5.917518E+15
5.341501E-05
7.668405E-10
6.344301E-05
8.082762E-10
1.848020E-01
6.844659E-03
2.164241E+00
9.405018E-01
1.593040E-01
5.083825E-03
1.456985E-01
4.252227E-03
1.231722E+00
3.042957E-01
2.150164E+00
9.283068E-01
1.602628E+08
5.157247E+15
1.716554E+08
5.916536E+15
6.595811E-03
8.735488E-06
6.478686E-03
9.795832E-06
0.000000E+00
0.000000E+00
3.596400E+02
2.593705E+04
0.000000E+00
0.000000E+00
5.333993E-05
5.721669E-10
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
5.347096E-05
1.195492E-09
2.355669E-04
1.111058E-08
0.000000E+00
0.000000E+00
1.092124E+01
2.389726E+01
1.091835E+01
2.388439E+01
1.092772E+01
2.392590E+01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.447311E+00
1.200639E+00
1.705476E-02
5.829204E-05
3.596400E+02
2.593705E+04
9.572182E-01
1.838063E-01
5.333993E-05
5.721669E-10
1.854356E+08
6.897901E+15
3.078283E-03
3.361020E-06
3.654666E-01
2.679249E-02
1.042153E+00
2.179058E-01
2.357652E+00
1.114919E+00
7.143905E+01
1.022953E+03
7.054186E+01
9.974400E+02
7.388636E+01
1.094201E+03
2.340596E+00
1.098853E+00
1.726778E+08
5.981585E+15
1.852040E+08
6.880779E+15
8.647179E-03
1.497891E-05
tally 25:
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.872500E+00
7.043938E-01
1.396649E-02
7.598320E-05
1.443500E+01
4.188168E+01
9.738738E-01
1.909539E-01
1.456937E-05
4.256948E-11
1.886870E+08
7.168253E+15
2.500000E-03
6.250000E-06
3.075000E-01
2.003125E-02
0.000000E+00
0.000000E+00
2.353102E+00
1.116236E+00
1.256250E+01
3.173376E+01
1.181750E+01
2.806099E+01
1.443500E+01
4.188168E+01
2.339136E+00
1.103322E+00
1.756374E+08
6.211106E+15
1.883952E+08
7.146298E+15
5.888970E-03
1.301949E-05
1.000000E-02
3.750000E-05
0.000000E+00
0.000000E+00
3.042500E+00
1.870344E+00
0.000000E+00
0.000000E+00
5.377597E-06
5.793845E-12
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.250000E-02
2.437500E-04
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.032500E+00
1.857869E+00
2.942500E+00
1.749119E+00
3.042500E+00
1.870344E+00
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.500000E-01
4.235000E-02
0.000000E+00
0.000000E+00
4.809000E+01
4.635862E+02
0.000000E+00
0.000000E+00
3.241167E-05
2.111966E-10
0.000000E+00
0.000000E+00
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.764000E+01
4.549033E+02
4.764000E+01
4.549033E+02
4.809000E+01
4.635862E+02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.410000E+00
1.164662E+00
1.396649E-02
7.598320E-05
7.362500E+01
1.086573E+03
9.738738E-01
1.909539E-01
5.893474E-05
6.969745E-10
1.886870E+08
7.168253E+15
2.500000E-03
6.250000E-06
3.475000E-01
2.546875E-02
0.000000E+00
0.000000E+00
2.353102E+00
1.116236E+00
7.121500E+01
1.016622E+03
7.037000E+01
9.926181E+02
7.362500E+01
1.086573E+03
2.339136E+00
1.103322E+00
1.756374E+08
6.211106E+15
1.883952E+08
7.146298E+15
5.888970E-03
1.301949E-05
tally 26:
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.082500E+00
2.355437E-01
4.768447E-03
1.139585E-05
1.232500E+00
3.052313E-01
8.975008E-01
1.621256E-01
4.660895E-06
4.397492E-12
1.735815E+08
6.064414E+15
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.168735E+00
9.473268E-01
1.500000E-01
4.737500E-03
1.300000E-01
3.662500E-03
1.232500E+00
3.052313E-01
2.163966E+00
9.432460E-01
1.620580E+08
5.285958E+15
1.735793E+08
6.064269E+15
3.322020E-04
9.208899E-08
2.500000E-03
6.250000E-06
0.000000E+00
0.000000E+00
5.337500E+00
5.718394E+00
0.000000E+00
0.000000E+00
4.162328E-06
3.477362E-12
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
5.335000E+00
5.713188E+00
5.335000E+00
5.713188E+00
5.337500E+00
5.718394E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.872500E+00
7.043938E-01
1.396649E-02
7.598320E-05
1.443500E+01
4.188168E+01
9.738738E-01
1.909539E-01
1.456937E-05
4.256948E-11
1.886870E+08
7.168253E+15
2.500000E-03
6.250000E-06
3.075000E-01
2.003125E-02
0.000000E+00
0.000000E+00
2.353102E+00
1.116236E+00
1.256250E+01
3.173376E+01
1.181750E+01
2.806099E+01
1.443500E+01
4.188168E+01
2.339136E+00
1.103322E+00
1.756374E+08
6.211106E+15
1.883952E+08
7.146298E+15
5.888970E-03
1.301949E-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
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.000000E-02
3.750000E-05
0.000000E+00
0.000000E+00
3.042500E+00
1.870344E+00
0.000000E+00
0.000000E+00
5.377597E-06
5.793845E-12
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.250000E-02
2.437500E-04
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.032500E+00
1.857869E+00
2.942500E+00
1.749119E+00
3.042500E+00
1.870344E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.500000E-03
6.250000E-06
0.000000E+00
0.000000E+00
4.900000E+00
4.810900E+00
0.000000E+00
0.000000E+00
1.475862E-06
4.417388E-13
0.000000E+00
0.000000E+00
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.897500E+00
4.805869E+00
4.897500E+00
4.805869E+00
4.900000E+00
4.810900E+00
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.500000E-01
4.235000E-02
0.000000E+00
0.000000E+00
4.809000E+01
4.635862E+02
0.000000E+00
0.000000E+00
3.241167E-05
2.111966E-10
0.000000E+00
0.000000E+00
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.764000E+01
4.549033E+02
4.764000E+01
4.549033E+02
4.809000E+01
4.635862E+02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.082500E+00
2.355437E-01
4.768447E-03
1.139585E-05
1.232500E+00
3.052313E-01
8.975008E-01
1.621256E-01
4.660895E-06
4.397492E-12
1.735815E+08
6.064414E+15
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.168735E+00
9.473268E-01
1.500000E-01
4.737500E-03
1.300000E-01
3.662500E-03
1.232500E+00
3.052313E-01
2.163966E+00
9.432460E-01
1.620580E+08
5.285958E+15
1.735793E+08
6.064269E+15
3.322020E-04
9.208899E-08
5.000000E-03
2.500000E-05
0.000000E+00
0.000000E+00
1.100000E+01
2.426909E+01
0.000000E+00
0.000000E+00
5.879616E-06
6.950007E-12
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.099500E+01
2.424706E+01
1.099500E+01
2.424706E+01
1.100000E+01
2.426909E+01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.410000E+00
1.164662E+00
1.396649E-02
7.598320E-05
7.362500E+01
1.086573E+03
9.738738E-01
1.909539E-01
5.893474E-05
6.969745E-10
1.886870E+08
7.168253E+15
2.500000E-03
6.250000E-06
3.475000E-01
2.546875E-02
0.000000E+00
0.000000E+00
2.353102E+00
1.116236E+00
7.121500E+01
1.016622E+03
7.037000E+01
9.926181E+02
7.362500E+01
1.086573E+03
2.339136E+00
1.103322E+00
1.756374E+08
6.211106E+15
1.883952E+08
7.146298E+15
5.888970E-03
1.301949E-05
tally 27:
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.849376E+00
6.865029E-01
1.675823E-02
5.637034E-05
1.443500E+01
4.188168E+01
9.419399E-01
1.779426E-01
1.456937E-05
4.256948E-11
1.824754E+08
6.677996E+15
4.669911E-03
6.353478E-06
3.285294E-01
2.182240E-02
9.037721E-01
1.644322E-01
2.321331E+00
1.080765E+00
1.258562E+01
3.185085E+01
1.180336E+01
2.800309E+01
1.443500E+01
4.188168E+01
2.304572E+00
1.065211E+00
1.699151E+08
5.790227E+15
1.822399E+08
6.660724E+15
8.492170E-03
1.448733E-05
1.758680E-02
6.526252E-05
0.000000E+00
0.000000E+00
3.042500E+00
1.870344E+00
0.000000E+00
0.000000E+00
5.377597E-06
5.793845E-12
0.000000E+00
0.000000E+00
2.208312E-04
4.535270E-08
2.765573E-02
1.577589E-04
1.754734E-02
6.495904E-05
0.000000E+00
0.000000E+00
3.024913E+00
1.848625E+00
2.927116E+00
1.731112E+00
3.042500E+00
1.870344E+00
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.779751E-01
4.599007E-02
0.000000E+00
0.000000E+00
4.809000E+01
4.635862E+02
0.000000E+00
0.000000E+00
3.241167E-05
2.111966E-10
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
5.202339E-05
6.721590E-10
1.000818E-01
2.016681E-03
0.000000E+00
0.000000E+00
4.761202E+01
4.544142E+02
4.761082E+01
4.543909E+02
4.809000E+01
4.635862E+02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.427789E+00
1.180817E+00
1.675823E-02
5.637034E-05
7.362500E+01
1.086573E+03
9.419399E-01
1.779426E-01
5.893474E-05
6.969745E-10
1.824754E+08
6.677996E+15
4.890742E-03
7.253871E-06
3.596282E-01
2.609803E-02
1.040856E+00
2.176627E-01
2.321331E+00
1.080765E+00
7.119721E+01
1.016165E+03
7.030311E+01
9.908064E+02
7.362500E+01
1.086573E+03
2.304572E+00
1.065211E+00
1.699151E+08
5.790227E+15
1.822399E+08
6.660724E+15
8.492170E-03
1.448733E-05
tally 28:
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.060216E+00
2.253377E-01
1.385339E-02
3.848021E-05
1.443500E+01
4.188168E+01
8.735389E-01
1.529994E-01
1.456937E-05
4.256948E-11
1.689471E+08
5.723043E+15
7.046785E-05
1.240304E-09
6.269959E-05
7.965202E-10
1.866769E-01
6.985523E-03
2.129960E+00
9.096368E-01
1.562417E-01
4.907671E-03
1.428934E-01
4.102873E-03
1.216458E+00
2.966425E-01
2.116107E+00
8.978426E-01
1.577208E+08
4.987735E+15
1.689327E+08
5.722065E+15
6.491030E-03
8.447968E-06
3.777489E-03
3.255196E-06
0.000000E+00
0.000000E+00
1.443500E+01
4.188168E+01
0.000000E+00
0.000000E+00
1.456937E-05
4.256948E-11
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
4.104266E-05
4.889770E-10
1.139025E-04
2.611032E-09
0.000000E+00
0.000000E+00
5.304440E+00
5.651143E+00
5.302200E+00
5.646237E+00
5.308217E+00
5.659334E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.849376E+00
6.865029E-01
1.675823E-02
5.637034E-05
1.443500E+01
4.188168E+01
9.419399E-01
1.779426E-01
1.456937E-05
4.256948E-11
1.824754E+08
6.677996E+15
4.669911E-03
6.353478E-06
3.285294E-01
2.182240E-02
9.037721E-01
1.644322E-01
2.321331E+00
1.080765E+00
1.258562E+01
3.185085E+01
1.180336E+01
2.800309E+01
1.443500E+01
4.188168E+01
2.304572E+00
1.065211E+00
1.699151E+08
5.790227E+15
1.822399E+08
6.660724E+15
8.492170E-03
1.448733E-05
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.042500E+00
1.870344E+00
0.000000E+00
0.000000E+00
5.377597E-06
5.793845E-12
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
3.042500E+00
1.870344E+00
0.000000E+00
0.000000E+00
5.377597E-06
5.793845E-12
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.758680E-02
6.526252E-05
0.000000E+00
0.000000E+00
3.042500E+00
1.870344E+00
0.000000E+00
0.000000E+00
5.377597E-06
5.793845E-12
0.000000E+00
0.000000E+00
2.208312E-04
4.535270E-08
2.765573E-02
1.577589E-04
1.754734E-02
6.495904E-05
0.000000E+00
0.000000E+00
3.024913E+00
1.848625E+00
2.927116E+00
1.731112E+00
3.042500E+00
1.870344E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
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.809000E+01
4.635862E+02
0.000000E+00
0.000000E+00
3.241167E-05
2.111966E-10
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.745161E-03
1.702574E-06
0.000000E+00
0.000000E+00
4.809000E+01
4.635862E+02
0.000000E+00
0.000000E+00
3.241167E-05
2.111966E-10
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.558491E-05
2.364490E-10
1.028539E-04
2.120089E-09
0.000000E+00
0.000000E+00
4.755015E+00
4.532855E+00
4.753850E+00
4.530567E+00
4.757760E+00
4.538149E+00
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.779751E-01
4.599007E-02
0.000000E+00
0.000000E+00
4.809000E+01
4.635862E+02
0.000000E+00
0.000000E+00
3.241167E-05
2.111966E-10
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
5.202339E-05
6.721590E-10
1.000818E-01
2.016681E-03
0.000000E+00
0.000000E+00
4.761202E+01
4.544142E+02
4.761082E+01
4.543909E+02
4.809000E+01
4.635862E+02
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
1.060216E+00
2.253377E-01
1.385339E-02
3.848021E-05
7.362500E+01
1.086573E+03
8.735389E-01
1.529994E-01
5.893474E-05
6.969745E-10
1.689471E+08
5.723043E+15
7.046785E-05
1.240304E-09
6.269959E-05
7.965202E-10
1.866769E-01
6.985523E-03
2.129960E+00
9.096368E-01
1.562417E-01
4.907671E-03
1.428934E-01
4.102873E-03
1.216458E+00
2.966425E-01
2.116107E+00
8.978426E-01
1.577208E+08
4.987735E+15
1.689327E+08
5.722065E+15
6.491030E-03
8.447968E-06
6.880377E-03
1.045360E-05
0.000000E+00
0.000000E+00
7.362500E+01
1.086573E+03
0.000000E+00
0.000000E+00
5.893474E-05
6.969745E-10
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
6.946156E-05
1.302423E-09
2.323582E-04
1.082628E-08
0.000000E+00
0.000000E+00
1.079564E+01
2.337660E+01
1.079216E+01
2.336117E+01
1.080252E+01
2.340680E+01
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
0.000000E+00
2.427789E+00
1.180817E+00
1.675823E-02
5.637034E-05
7.362500E+01
1.086573E+03
9.419399E-01
1.779426E-01
5.893474E-05
6.969745E-10
1.824754E+08
6.677996E+15
4.890742E-03
7.253871E-06
3.596282E-01
2.609803E-02
1.040856E+00
2.176627E-01
2.321331E+00
1.080765E+00
7.119721E+01
1.016165E+03
7.030311E+01
9.908064E+02
7.362500E+01
1.086573E+03
2.304572E+00
1.065211E+00
1.699151E+08
5.790227E+15
1.822399E+08
6.660724E+15
8.492170E-03
1.448733E-05
tally 29:
3.335068E+01
2.229495E+02
1.128844E+01
2.558948E+01
5.471015E+01
6.004217E+02
3.352822E+01
2.250658E+02
1.158784E+01
2.687498E+01
5.594127E+01
6.263579E+02
1.090994E+02
2.387594E+03
tally 30:
3.275696E+01
2.160444E+02
1.121655E+01
2.537659E+01
5.459568E+01
5.978631E+02
3.336676E+01
2.230181E+02
1.123814E+01
2.527618E+01
5.600406E+01
6.275331E+02
1.084967E+02
2.362981E+03
tally 31:
3.275696E+01
2.160444E+02
1.121655E+01
2.537659E+01
5.459568E+01
5.978631E+02
3.336676E+01
2.230181E+02
1.123814E+01
2.527618E+01
5.600406E+01
6.275331E+02
1.084967E+02
2.362981E+03
tally 32:
1.363871E-04
4.518856E-09
2.581266E-01
1.339951E-02
2.512393E-05
1.271096E-10
2.540779E-21
5.645239E-42
9.830895E-01
1.942917E-01
3.240205E+08
2.103581E+16
3.198869E+05
2.049459E+10

Details

Running on the example problem from example-filters.ipynb:

spatiallegendre	nuclide	score	mean	std. dev.
P0	total	flux	3.63e+01	8.06e-02
P1	total	flux	-2.60e-02	4.23e-02
P2	total	flux	-4.35e+00	2.69e-02
P3	total	flux	-1.36e-02	2.02e-02
P4	total	flux	-2.81e-01	1.40e-02
P5	total	flux	4.92e-03	1.15e-02
P6	total	flux	1.13e-01	9.69e-03
P7	total	flux	5.20e-03	8.94e-03
P8	total	flux	-7.84e-02	7.61e-03

spatialfourier	nuclide	score	mean	std. dev.
a0 (constant)	total	flux	3.63e+01	8.06e-02
a1 (cos)	total	flux	-6.45e+00	4.60e-02
b1 (sin)	total	flux	8.13e-03	6.10e-02
a2 (cos)	total	flux	-2.32e+00	2.96e-02
b2 (sin)	total	flux	2.65e-02	3.14e-02
a3 (cos)	total	flux	-6.38e-01	2.47e-02
b3 (sin)	total	flux	2.61e-02	2.51e-02
a4 (cos)	total	flux	-6.40e-01	2.42e-02
b4 (sin)	total	flux	2.72e-02	2.35e-02

Copy-paste SpatialLegendreFilter. No other changes needed. Could be combined to a SpatialExpansionFilter header?

MicahGale · 2026-03-25T20:13:20Z

I'll give this a look as well. @pshriwise and @magnoxemo might be interested as well.

tjlaboss added 16 commits March 24, 2026 17:44

Add SpatialExpansionFilter, SpatialLegendreFilter to Python API

1441cf6

SpatialFourierFilter header

1c56e0a

Copy-paste SpatialLegendreFilter. No other changes needed. Could be combined to a SpatialExpansionFilter header?

openmc.lib.SpatialFourierFilter boilerplate

cd23c73

Start SpatialFourier C++

34c04a6

Implement SpatialFourierFilter

0fc321c

clang-format

27d824c

Fourier filter bin labels

cbd7f78

Unit test SpatialFourierFilter

3b8d77b

Need to inherit from SpatialExpansionFilter

857f8ba

Add SpatialLegendre and SpatialFourier to regression test

b2f4df5

Clean up spacing

2f21b7f

Register "spatialfourier"

c99df6e

Add to a few more locations

4ba558d

Add spatial fourier to CMakeLists.txt

de2cbb7

Unbotch inputs_true.dat

0633eea

pytest update

df6ec06

tjlaboss marked this pull request as ready for review March 25, 2026 19:41

tjlaboss requested a review from paulromano as a code owner March 25, 2026 19:41

tjlaboss changed the title ~~Fourier expansion~~ Spatial Fourier Functional Expansion Tally Mar 26, 2026

Provide feedback

Saved searches

Use saved searches to filter your results more quickly

Spatial Fourier Functional Expansion Tally#3899

Spatial Fourier Functional Expansion Tally#3899
tjlaboss wants to merge 16 commits intoopenmc-dev:developfrom
tjlaboss:fourier-expansion

tjlaboss commented Mar 25, 2026 •

edited

Loading

Uh oh!

MicahGale commented Mar 25, 2026

Uh oh!

Reviewers

Assignees

Labels

Projects

Milestone

Development

Uh oh!

2 participants

Conversation

tjlaboss commented Mar 25, 2026 • edited Loading Uh oh! There was an error while loading. Please reload this page.

Uh oh!

Description

Checklist

Discussion

Results

Uh oh!

MicahGale commented Mar 25, 2026

Uh oh!

Reviewers

Assignees

Labels

Projects

Milestone

Development

Uh oh!

2 participants

tjlaboss commented Mar 25, 2026 •

edited

Loading