| 113 | |
| 114 | |
| 115 | qcomp getSum(qvector vec) { |
| 116 | |
| 117 | qcomp sum = qcomp(0,0); |
| 118 | qcomp y, t, c=sum; |
| 119 | |
| 120 | // complex Kahan summation |
| 121 | for (auto& x : vec) { |
| 122 | y = x - c; |
| 123 | t = sum + y; |
| 124 | c = ( t - sum ) - y; |
| 125 | sum = t; |
| 126 | } |
| 127 | |
| 128 | return sum; |
| 129 | } |
| 130 | |
| 131 | |
| 132 | qreal getSum(vector<qreal> vec) { |
no test coverage detected