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 both the C++ and Python sides. 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?
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 · 2026-04-20T20:16:40Z

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

Would @MicahGale, @pshriwise, or @magnoxemo like to review this (at least as a first pass)?

MicahGale · 2026-04-20T20:18:19Z

I can do a "conceptual review" for sure, but no guarantees on competence of the review.

magnoxemo · 2026-04-21T18:39:46Z

@tjlaboss Just one quick question, are you planning to implement multi-D support in a following PR? or this is gonna be 1D for now?

tjlaboss · 2026-04-21T22:00:16Z

No, 1D. Multiple dimensions are left as an exercise to the reader.

MicahGale · 2026-04-22T01:14:26Z

No, 1D. Multiple dimensions are left as an exercise to the reader.

You need to write the M&C paper first to say that.

MicahGale

Just a few comments. I don't understand the tally infrastructure well enough to have strong opinions. Though I couldn't find where you actually calculate the Fourier series.

MicahGale · 2026-04-22T12:52:38Z

+    @order.setter
+    def order(self, order):


Should this make a call to check_value?

Not in openmc.lib I think.

tjlaboss · 2026-06-11T21:38:41Z

Went ahead and made SpatialExpansion on the C++ side too.

Also fixed some C API stuff, doc references, etc. that I missed previously.

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

MicahGale reviewed Apr 22, 2026

View reviewed changes

tjlaboss added 4 commits June 11, 2026 14:50

Register spatialfourier in _FILTER_TYPE_MAP

57673b5

Include SpatialFourier in docs/source/pythonapi

05f5d22

Implement SpatialExpansionFilter in C++ too

c198b23

_dll get and set order

e940f7e

whitespace cleanup

0dcc91a

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 21 commits into
openmc-dev:developfrom
tjlaboss:fourier-expansion

tjlaboss commented Mar 25, 2026 •

edited

Loading

Uh oh!

MicahGale commented Mar 25, 2026

Uh oh!

tjlaboss commented Apr 20, 2026

Uh oh!

MicahGale commented Apr 20, 2026

Uh oh!

magnoxemo commented Apr 21, 2026 •

edited

Loading

Uh oh!

tjlaboss commented Apr 21, 2026

Uh oh!

MicahGale commented Apr 22, 2026

Uh oh!

MicahGale left a comment

Uh oh!

MicahGale Apr 22, 2026

Uh oh!

tjlaboss Jun 11, 2026

Uh oh!

Uh oh!

tjlaboss commented Jun 11, 2026

Uh oh!

Reviewers

Assignees

Labels

Projects

Milestone

Development

Uh oh!

3 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!

tjlaboss commented Apr 20, 2026

Uh oh!

MicahGale commented Apr 20, 2026

Uh oh!

magnoxemo commented Apr 21, 2026 • edited Loading Uh oh! There was an error while loading. Please reload this page.

Uh oh!

Uh oh!

tjlaboss commented Apr 21, 2026

Uh oh!

MicahGale commented Apr 22, 2026

Uh oh!

MicahGale left a comment

Choose a reason for hiding this comment

Uh oh!

MicahGale Apr 22, 2026

Choose a reason for hiding this comment

Uh oh!

tjlaboss Jun 11, 2026

Choose a reason for hiding this comment

Uh oh!

Uh oh!

tjlaboss commented Jun 11, 2026

Uh oh!

Reviewers

Assignees

Labels

Projects

Milestone

Development

Uh oh!

3 participants

tjlaboss commented Mar 25, 2026 •

edited

Loading

magnoxemo commented Apr 21, 2026 •

edited

Loading