LSSTApplications
18.0.0+106,18.0.0+50,19.0.0,19.0.0+1,19.0.0+10,19.0.0+11,19.0.0+13,19.0.0+17,19.0.0+2,19.0.0-1-g20d9b18+6,19.0.0-1-g425ff20,19.0.0-1-g5549ca4,19.0.0-1-g580fafe+6,19.0.0-1-g6fe20d0+1,19.0.0-1-g7011481+9,19.0.0-1-g8c57eb9+6,19.0.0-1-gb5175dc+11,19.0.0-1-gdc0e4a7+9,19.0.0-1-ge272bc4+6,19.0.0-1-ge3aa853,19.0.0-10-g448f008b,19.0.0-12-g6990b2c,19.0.0-2-g0d9f9cd+11,19.0.0-2-g3d9e4fb2+11,19.0.0-2-g5037de4,19.0.0-2-gb96a1c4+3,19.0.0-2-gd955cfd+15,19.0.0-3-g2d13df8,19.0.0-3-g6f3c7dc,19.0.0-4-g725f80e+11,19.0.0-4-ga671dab3b+1,19.0.0-4-gad373c5+3,19.0.0-5-ga2acb9c+2,19.0.0-5-gfe96e6c+2,w.2020.01
LSSTDataManagementBasePackage
|
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 |
|
inlinenoexcept |
Definition at line 89 of file SafeSum.h.
|
inlinenoexcept |