LSST Applications g02d81e74bb+86cf3d8bc9,g180d380827+7a4e862ed4,g2079a07aa2+86d27d4dc4,g2305ad1205+e1ca1c66fa,g29320951ab+012e1474a1,g295015adf3+341ea1ce94,g2bbee38e9b+0e5473021a,g337abbeb29+0e5473021a,g33d1c0ed96+0e5473021a,g3a166c0a6a+0e5473021a,g3ddfee87b4+c429d67c83,g48712c4677+f88676dd22,g487adcacf7+27e1e21933,g50ff169b8f+96c6868917,g52b1c1532d+585e252eca,g591dd9f2cf+b41db86c35,g5a732f18d5+53520f316c,g64a986408d+86cf3d8bc9,g858d7b2824+86cf3d8bc9,g8a8a8dda67+585e252eca,g99cad8db69+84912a7fdc,g9ddcbc5298+9a081db1e4,ga1e77700b3+15fc3df1f7,ga8c6da7877+a2b54eae19,gb0e22166c9+60f28cb32d,gba4ed39666+c2a2e4ac27,gbb8dafda3b+6681f309db,gc120e1dc64+f0fcc2f6d8,gc28159a63d+0e5473021a,gcf0d15dbbd+c429d67c83,gdaeeff99f8+f9a426f77a,ge6526c86ff+0433e6603d,ge79ae78c31+0e5473021a,gee10cc3b42+585e252eca,gff1a9f87cc+86cf3d8bc9,w.2024.17
LSST Data Management Base Package
|
A numerically stable summation algorithm for floating-point numbers. More...
#include <SafeSum.h>
Public Member Functions | |
SafeSum (T initial=static_cast< T >(0)) noexcept | |
SafeSum & | operator= (T value) noexcept |
SafeSum & | operator+= (T value) noexcept |
SafeSum & | operator-= (T value) noexcept |
operator T () const noexcept | |
A numerically stable summation algorithm for floating-point numbers.
SafeSum implements += and -= operators that can be used to accumulate floating point numbers with very different magnitudes, with accuracy limited only by the usual floating-point inaccuracy in represented the final sum.
SafeSum is explicitly convertible to and from its underlying floating-point type and only supports in-place addition and subtraction, in order to avoid cases where a combination of implicit conversion and multiple overloaded operators could lead to accidental use of regular floating-point operations.
SafeSum uses the Kahan-Neumaier algorithm (though this should be considered an implementation detail by callers), which maintains a lower-order-bit correction that compensates for the loss of precision in the main sum. Particularly aggressive compiler optimizations that do not preserve IEEE floating point behavior (e.g. gcc's -fassociative-math
) may optimize away the correction and reduce SafeSum's behavior to a standard unsafe sum.
|
inlineexplicitnoexcept |
|
inlineexplicitnoexcept |
|
inlinenoexcept |
Definition at line 76 of file SafeSum.h.
|
inlinenoexcept |
|
inlinenoexcept |