id,summary,reporter,owner,description,type,status,milestone,component,version,severity,resolution,keywords,cc
10104,program aborts with log1p(-2.0Q),Charles Karney <charles@…>,John Maddock,"The math functions, e.g., log1p, for float128 do not handle parameters
outside the normal domain.  Thus sqrt(-2.0Q) produces a nan as does
log1p(-2.0), but log1p(-2.0Q) crashes.  (Similar remarks apply to hypot
and probably other similarly defined functions.)


{{{
#include <iostream>
#include <cmath>
#include <boost/multiprecision/float128.hpp>
#include <boost/math/special_functions/log1p.hpp>

int main() {
  double x = -2;
  boost::multiprecision::float128 y = -2;
  
  std::cout << std::sqrt(x) << ""\n"" << std::log1p(x) << ""\n"";
  std::cout << sqrt(y) << ""\n"";
  std::cout << boost::math::log1p(y) << ""\n"";
}
}}}
produces
{{{
-nan
-nan
-nan
terminate called after throwing an instance of 'boost::exception_detail::clone_impl<boost::exception_detail::error_info_injector<std::domain_error> >'
  what():  Error in function boost::math::log1p<N5boost14multiprecision6numberINS0_8backends16float128_backendELNS0_26expression_template_optionE0EEE>(N5boost14multiprecision6numberINS0_8backends16float128_backendELNS0_26expression_template_optionE0EEE): log1p(x) requires x > -1, but got x = -2.
Aborted (core dumped)
}}}
",Bugs,closed,To Be Determined,multiprecision,Boost 1.54.0,Problem,invalid,,