Avoided NaNs in keff and fission analog tallies when micro_xs(i)%absorption was zero as it is for He-4.
Fixed nu-scatter (w/ energyout filter) survival biasing issue.
Added Real Spherical Harmonics function, calc_rn to math.F90; orders up to N=10 are provided; these specific harmonics are chosen to be of the Ferrer form, consistent with chapter 8 of the nuclear engineering handbook. These have been verified to produce the same results as in the NE handbook (for up to l=3) as well as those produced by SymPy for other orders (which also match the NE Handbook for L<=3).
Conformed the code in math.F90 to OMC standards (code was generated by another code).
Added ability to read the new harmonic tally types in input_xml.
Renamed tally % scatt_order to tally % moment_order, updated revision statepoint number, added new moment tally types to output.F90.
added last_uvw attribute to particle type to store the pre-collision direction for accurate analog flux-moment weighted tallies.
And corrected a bug i so easily introduced in previous commit; was storing last_uvw in wrong place.
Added tallying of nu-scatter-n and nu-scatter-pn
Fixed input_xmls setting of size for -yn moments (size is (N+1)^2, not N)
Added tallies for all types, now to debug.
Seem to have all the bugs fixed for the harmonic tallying; next up will be documentation.
Modified tracklength tallies of harmonics to use current uvw, not last_uvw. whoops.
Simplified printing of -pn or -yn tallies in write_tallies; also had to extend the score name size, but tallies.out will still not exceed 80 characters.
Updated documentation for the new score types.;
Merge remote-tracking branch 'origin/develop' into survival
Fixed erroneous tallying of fission and kappa-fission in cases without survival biasing.
Fixed flux tallying for case of survival biasing.
Merge remote-tracking branch 'upstream/develop' into harmonics
Merge remote-tracking branch 'upstream/develop' into survival