부스트 muliprecision cpp_dec_float를 원하는 정밀도와 비교하십시오.

Question

10 진수 부동 유형에 boost :: multiprecision 라이브러리를 사용하고 있고 두 개의 부동 소수점을 지정된 정밀도와 비교하려고합니다.

그러나 cpp_dec_float 지정된 정밀도로하지 수를 비교하는 것뿐만 아니라, 가드 자리 포함

#include <iostream> 

#include <boost/multiprecision/cpp_dec_float.hpp> 
//#include <boost/math/special_functions.hpp> 

typedef boost::multiprecision::number<boost::multiprecision::cpp_dec_float<50> > flp_type; 

int main(int argc, char* argv[]) 
{ 
    // 50 decimal digits 
    flp_type sqrt2("1.4142135623730950488016887242096980785696718753769"); 
    // Contains calculated guard digits 
    flp_type result(boost::multiprecision::sqrt(flp_type("2"))); 

    // The 50 digits of precision actually ompare equal 
    std::cout << std::setprecision(50) << sqrt2 << std::endl; 
    std::cout << std::setprecision(50) << result << std::endl; 
    // I want this to compare to the specified precision of the type, not the guard digits 
    std::cout << (result==sqrt2) << std::endl; 

    return 0; 
} 

출력 : 예상

1.4142135623730950488016887242096980785696718753769 
1.4142135623730950488016887242096980785696718753769 
0 

:

1.4142135623730950488016887242096980785696718753769 
1.4142135623730950488016887242096980785696718753769 
1 

을 See on Coliru

precision()을 사용하여 "잘라내려고"시도했지만 아무 소용이 없습니다. 엡실론 비교없이 두 숫자를 비교하는 방법이 있습니까?

Answer 1

가드 비트를 제거하면 유형의 충실도가 효과적으로 상실됩니다.

확실히 확실한 방법은 (비) 직렬화를 사용하는 것입니다.

그래서 나는 엡실론이 1e-49이

인쇄입니다

Live On Coliru

// Either 
std::cout << std::numeric_limits<flp_type>::epsilon() << "\n"; 
std::cout << (abs(result-sqrt2) < std::numeric_limits<flp_type>::epsilon()) << std::endl; 

// Or 
result = flp_type { result.str(49, std::ios::fixed) }; 
std::cout << (result==sqrt2) << std::endl; 

주를 제안

1.4142135623730950488016887242096980785696718753769 
1.4142135623730950488016887242096980785696718753769 
1e-49 
1 
1 

분명히 epsilon() 기반 비교가 더 효율적으로 보일 것입니다.