Как получить значение числа ПИ с заданной точностью во время компиляции?
2 ответа
В основном взято отсюда Добавлено только ограничение точности.
#include <iostream>
#include <limits>
// Generic exponent compile-time calculations. Pow<2,3>::result == 2^3
template <unsigned long B, unsigned long E>
struct Pow
{
static constexpr double result = B * Pow<B, E - 1>::result;
};
template <unsigned long B>
struct Pow<B, 0>
{
static constexpr double result = 1;
};
// Bailey-Borwein-Plouffe formula for calculating pi
// http://en.wikipedia.org/wiki/Bailey-Borwein-Plouffe_formula
template <unsigned long N>
struct CalculatePi
{
static constexpr double pi =
(
1.0 / Pow<16, N>::result *
(
4.0 / (8 * N + 1.0) - 2.0 / (8 * N + 4.0) -
1.0 / (8 * N + 5.0) - 1.0 / (8 * N + 6.0)
)
) + CalculatePi<N - 1>::pi;
;
};
template <>
struct CalculatePi<-1>
{
static constexpr double pi = 0.0;
};
template <int IT = 10, int N = 10>
struct Pecision : public CalculatePi<IT>
{
static constexpr double value = (long long)(CalculatePi<IT>::pi * Pow <10, N>::result) / Pow <10, N>::result;
};
template <>
struct Pecision<0>
{
static constexpr double value = 0.0;
};
// main program. Print pi, calculated from 10 iterations of
// the BBP formula above
int main()
{
std::cout.precision(std::numeric_limits<double>::digits10);
std::cout << "pi: " << Pecision <10>::value << std::endl;
return 0;
}
-
-
@3.14 Набросал на скорую руку, надо было точность первым параметром поставить. поэтому надо писать так Pecision <10, 9> второй параметр точность 24 дек 2018 в 17:16
-
-
-
@3.14 Для double вот так будет более корректно Pecision<10, std::numeric_limits<double>::max_digits10>::value 24 дек 2018 в 17:22
Например с помощью серий Грэгори.
// Example program
#include <iostream>
#include <string>
template <class T>
T calculatePi(int depth)
{
T pi = 0;
bool sign = true; // true: +, false: -
for (int i = 1; i < depth; i += 2)
{
if (sign)
pi += 1/static_cast<T>(i);
else
pi -= 1/static_cast<T>(i);
sign = !sign;
}
return pi*4;
}
int main()
{
std::cout << calculatePi<double>(10000);
}
-
-
3@3.14 Поставьте перед функцией
constexpr
- будет во время компиляции. 24 дек 2018 в 15:38 -
Проблема в том что в теле constexpr функции может быть только return– user24515024 дек 2018 в 15:45
-
3Как по мне это бессмыслица.– user24515024 дек 2018 в 16:35
-
3@3.14 А почему вы компилятор готовый используете, так ведь любой может. Нужно свой. :Р 24 дек 2018 в 16:38
long double
, или может потребоваться произвольно большая точность?