blob: ce8a7f30f4dc0cc57e2b1262df7c8913fc6dedb7 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
|
/* -*- Mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*- */
/*
* This file is part of the LibreOffice project.
*
* This Source Code Form is subject to the terms of the Mozilla Public
* License, v. 2.0. If a copy of the MPL was not distributed with this
* file, You can obtain one at http://mozilla.org/MPL/2.0/.
*
*/
#pragma once
#include <math.h>
namespace sc::op
{
/**
* Performs one step of the Neumanier sum between doubles
* Overwrites the summand and error
* @parma sum
* @param err
* @param value
*/
inline void sumNeumanierNormal(double& sum, double& err, const double& value)
{
double t = sum + value;
if (fabs(sum) >= fabs(value))
err += (sum - t) + value;
else
err += (value - t) + sum;
sum = t;
}
} // end namespace sc::op
/* vim:set shiftwidth=4 softtabstop=4 expandtab: */
|