Как вы печатаете двойной поток, чтобы при чтении в нем вы не теряли точность?
Я пробовал:
std::stringstream ss;
double v = 0.1 * 0.1;
ss << std::setprecision(std::numeric_limits<T>::digits10) << v << " ";
double u;
ss >> u;
std::cout << "precision " << ((u == v) ? "retained" : "lost") << std::endl;
Это не сработало, как я ожидал.
Но я могу увеличить точность (что меня удивило, так как я считал, что цифры 10 являются максимальными).
ss << std::setprecision(std::numeric_limits<T>::digits10 + 2) << v << " ";
// ^^^^^^ +2
Он имеет отношение к числу значащих цифр, а первые два не учитываются в (0,01).
Итак, кто-нибудь посмотрел на то, что точно представляет числа с плавающей запятой? Какое точное магическое заклинание в потоке мне нужно сделать?
После некоторых экспериментов:
Проблема была с моей оригинальной версией. В десятичной строке были записаны незначительные цифры после десятичной точки, которые влияли на точность.
Итак, чтобы компенсировать это, мы можем использовать научную нотацию для компенсации:
ss << std::scientific
<< std::setprecision(std::numeric_limits<double>::digits10 + 1)
<< v;
Это все еще не объясняет необходимость в +1.
Также, если я распечатаю номер с большей точностью, я получаю более точную распечатку!
std::cout << std::scientific << std::setprecision(std::numeric_limits<double>::digits10) << v << "\n";
std::cout << std::scientific << std::setprecision(std::numeric_limits<double>::digits10 + 1) << v << "\n";
std::cout << std::scientific << std::setprecision(std::numeric_limits<double>::digits) << v << "\n";
Это приводит к:
1.000000000000000e-02
1.0000000000000002e-02
1.00000000000000019428902930940239457413554200000000000e-02
Основываясь на ответе @Stephen Canon ниже:
Мы можем распечатать именно с помощью форматирования printf(), "% a" или "% A". Для достижения этого на С++ нам нужно использовать фиксированные и научные манипуляторы (см. N3225: 22.4.2.2.2p5, таблица 88).
std::cout.flags(std::ios_base::fixed | std::ios_base::scientific);
std::cout << v;
Теперь я определил:
template<typename T>
std::ostream& precise(std::ostream& stream)
{
std::cout.flags(std::ios_base::fixed | std::ios_base::scientific);
return stream;
}
std::ostream& preciselngd(std::ostream& stream){ return precise<long double>(stream);}
std::ostream& precisedbl(std::ostream& stream) { return precise<double>(stream);}
std::ostream& preciseflt(std::ostream& stream) { return precise<float>(stream);}
Далее: Как мы обрабатываем NaN/Inf?