Function opensn::ComputePointwisePsiChange

Function Documentation

double opensn::ComputePointwisePsiChange(DiscreteOrdinatesProblem &do_problem, const LBSGroupset &groupset, const std::vector<double> &psi_new, const std::vector<double> &psi_old)

Maximum pointwise relative change between the new and old delayed angular unknowns of a groupset (the lagged incoming angular fluxes of opposing reflecting boundaries and broken cyclic dependencies, and their CSDA energy slopes), laid out as AngleAggregation::GetNewDelayedAngularDOFsAsSTLVector() returns them, reduced over all ranks.

Each entry’s change is relative to the larger of its own magnitude and the magnitude of the new or old scalar flux at the face node of the cell that reads it, for its group. An entry whose scale is below the smallest normal double contributes its absolute change. Collective.