Function SetBernoulliNumbersSum

extern/ttmath/ttmath.h:2233–2272

	*/
	template<class ValueType>
	ValueType SetBernoulliNumbersSum(CGamma<ValueType> & cgamma, const ValueType & n_, uint m,
									  const volatile StopCalculating * stop = 0)
	{
	ValueType k_, temp, temp2, temp3, sum;

		sum.SetZero();
		
		for(uint k=0 ; k<m ; ++k)			// k<m means k<=m-1
		{
			if( stop && (k & 15)==0 )		// means: k % 16 == 0
				if( stop->WasStopSignal() )
					return ValueType();		// NaN

			if( k>1 && (k & 1) == 1 )		// for that k the Bernoulli number is zero
				continue;

			k_ = k;

			temp = n_;				// n_ is equal 2
			temp.Pow(k_);
			// temp = 2^k

			temp2 = cgamma.fact[m];
			temp3 = cgamma.fact[k];
			temp3.Mul(cgamma.fact[m-k]);
			temp2.Div(temp3);
			// temp2 = (m k) = m! / ( k! * (m-k)! )

			temp.Mul(temp2);
			temp.Mul(cgamma.bern[k]);

			sum.Add(temp);
			// sum += 2^k * (m k) * B(k)

			if( sum.IsNan() )
				break;
		}

	return sum;
	}


	/*!