kahansum.h

Go to the documentation of this file.

 
#ifndef KAHANSUM_H
#define KAHANSUM_H
 
#include <complex.h>
 
static inline void kahaninit(double * const sum, double * const compensation) {
    *sum = 0;
    *compensation = 0;
}
 
static inline void kahanadd(double *sum, double *compensation, double input) {
    double compensated_input = input - *compensation;
    double nsum = *sum + compensated_input;
    *compensation = (nsum - *sum) - compensated_input;
    *sum = nsum;
}
 
 
static inline void ckahaninit(complex double * const sum, complex double * const compensation) {
    *sum = 0;
    *compensation = 0;
}
 
static inline void ckahanadd(complex double *sum, complex double *compensation, complex double input) {
    complex double compensated_input = input - *compensation;
    complex double nsum = *sum + compensated_input;
    *compensation = (nsum - *sum) - compensated_input;
    *sum = nsum;
}
 
#endif //KAHANSUM_H