From d1982bfca9bf5fe44cb681fd9c86cc3290c82cfd Mon Sep 17 00:00:00 2001
From: ckormanyos <ckormanyos@yahoo.com>
Date: Fri, 29 Aug 2025 14:07:10 +0200
Subject: [PATCH] Handle review comments

---
 .../multiprecision/detail/default_ops.hpp     | 34 ++-----------------
 .../multiprecision/detail/functions/pow.hpp   |  8 ++---
 2 files changed, 6 insertions(+), 36 deletions(-)

diff --git a/include/boost/multiprecision/detail/default_ops.hpp b/include/boost/multiprecision/detail/default_ops.hpp
index 6fd4b65e..87ada6e4 100644
--- a/include/boost/multiprecision/detail/default_ops.hpp
+++ b/include/boost/multiprecision/detail/default_ops.hpp
@@ -2619,13 +2619,6 @@ inline BOOST_MP_CXX14_CONSTEXPR typename multiprecision::detail::expression<tag,
 template <class Backend, multiprecision::expression_template_option ExpressionTemplates>
 inline BOOST_MP_CXX14_CONSTEXPR multiprecision::number<Backend, ExpressionTemplates> erf BOOST_PREVENT_MACRO_SUBSTITUTION(const multiprecision::number<Backend, ExpressionTemplates>& arg)
 {
-   const int fpc_arg { (boost::multiprecision::fpclassify)(arg) };
-
-   if (fpc_arg == FP_ZERO)
-   {
-      return arg;
-   }
-
    detail::scoped_default_precision<multiprecision::number<Backend, ExpressionTemplates> > precision_guard(arg);
    return boost::math::erf(arg, c99_error_policy());
 }
@@ -2634,14 +2627,6 @@ inline BOOST_MP_CXX14_CONSTEXPR typename multiprecision::detail::expression<tag,
 {
    using value_type = typename multiprecision::detail::expression<tag, A1, A2, A3, A4>::result_type;
    detail::scoped_default_precision<value_type> precision_guard(arg);
-
-   const int fpc_arg { (boost::multiprecision::fpclassify)(arg) };
-
-   if (fpc_arg == FP_ZERO)
-   {
-      return arg;
-   }
-
    return erf(value_type(arg));
 }
 template <class Backend, multiprecision::expression_template_option ExpressionTemplates>
@@ -2693,27 +2678,12 @@ inline BOOST_MP_CXX14_CONSTEXPR typename multiprecision::detail::expression<tag,
 template <class Backend, multiprecision::expression_template_option ExpressionTemplates>
 inline BOOST_MP_CXX14_CONSTEXPR multiprecision::number<Backend, ExpressionTemplates> tgamma BOOST_PREVENT_MACRO_SUBSTITUTION(const multiprecision::number<Backend, ExpressionTemplates>& arg)
 {
-   using local_number_type = multiprecision::number<Backend, ExpressionTemplates>;
-
-   detail::scoped_default_precision<local_number_type> precision_guard(arg);
-
-   const int fpc_arg { (boost::multiprecision::fpclassify)(arg) };
-
-   if ((fpc_arg == FP_ZERO) && std::numeric_limits<local_number_type>::has_infinity)
+   detail::scoped_default_precision<multiprecision::number<Backend, ExpressionTemplates> > precision_guard(arg);
+   if ((arg == 0) && std::numeric_limits<multiprecision::number<Backend, ExpressionTemplates> >::has_infinity)
    {
       errno = ERANGE;
       return 1 / arg;
    }
-   else if ((fpc_arg == FP_NAN) && std::numeric_limits<local_number_type>::has_quiet_NaN)
-   {
-      return arg;
-   }
-   else if ((fpc_arg == FP_INFINITE) && std::numeric_limits<local_number_type>::has_infinity && std::numeric_limits<local_number_type>::has_quiet_NaN)
-   {
-      const bool is_neg { arg.compare(local_number_type(0)) < 0 };
-
-      return ((!is_neg) ? std::numeric_limits<local_number_type>::infinity() : std::numeric_limits<local_number_type>::quiet_NaN());
-   }
    return boost::math::tgamma(arg, c99_error_policy());
 }
 template <class tag, class A1, class A2, class A3, class A4>
diff --git a/include/boost/multiprecision/detail/functions/pow.hpp b/include/boost/multiprecision/detail/functions/pow.hpp
index e79aacff..9cf265fb 100644
--- a/include/boost/multiprecision/detail/functions/pow.hpp
+++ b/include/boost/multiprecision/detail/functions/pow.hpp
@@ -509,15 +509,15 @@ inline void eval_pow(T& result, const T& x, const T& a)
       return;
    }
 
-   const int fpc_x { eval_fpclassify(x) };
-   const int fpc_a { eval_fpclassify(a) };
-
-   if (fpc_a == FP_ZERO)
+   if (eval_is_zero(a))
    {
       result = si_type(1);
       return;
    }
 
+   const int fpc_x { eval_fpclassify(x) };
+   const int fpc_a { eval_fpclassify(a) };
+
    switch (fpc_x)
    {
    case FP_ZERO: