AMReX_CompensatedSum.H

Go to the documentation of this file.

#ifndef AMREX_COMPENSATEDSUM_H_
#define AMREX_COMPENSATEDSUM_H_
#include <AMReX_Config.H>
 
#include <AMReX_Extension.H>
#include <AMReX_GpuQualifiers.H>
 
namespace amrex
{
 
template <typename T, typename U>
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE
void compensatedAdd (T& sum, T& compensation, U value) noexcept
{
    T volatile y = static_cast<T>(value) - compensation;
    T volatile t = sum + y;
    T volatile c = t - sum;
    compensation = c - y;
    sum = t;
}
 
}
 
#endif