普通的实数求和函数:
template< class data_type >
data_type sum( const data_type a[], int n )
{
data_type sum = 0.0;
for (int i=0; i<n; i++)
sum += a[i];
return sum;
}
修正一下中间计算产生的误差:
template< class data_type >
data_type kahansum( const data_type a[], int n )
{
data_type sum=0, c=0, y, t;
for (int i=0; i<n; i++)
{