Update doc

Author: andrjohns

stan/math/prim/fun/hypergeometric_3F2.hpp

@@ -12,7 +12,6 @@
 #include <stan/math/prim/fun/sum.hpp>
 #include <stan/math/prim/fun/sign.hpp>
 #include <stan/math/prim/fun/value_of_rec.hpp>
-#include <stan/math/prim/fun/log_sum_exp_signed.hpp>
 
 namespace stan {
 namespace math {
@@ -31,7 +30,7 @@ T_return hypergeometric_3F2_infsum(const Ta& a, const Tb& b, const Tz& z,
   check_3F2_converges("hypergeometric_3F2", a_array[0], a_array[1], a_array[2],
                       b_array[0], b_array[1], z);
 
-  T_return t_acc = 0.0;
+  T_return t_acc = 1.0;
   T_return log_t = 0.0;
   T_return log_z = log(fabs(z));
   Eigen::ArrayXi a_signs = sign(value_of_rec(a_array));
@@ -40,7 +39,7 @@ T_return hypergeometric_3F2_infsum(const Ta& a, const Tb& b, const Tz& z,
   plain_type_t<decltype(b_array)> bpk = b_array;
   int z_sign = sign(value_of_rec(z));
   int t_sign = z_sign * a_signs.prod() * b_signs.prod();
-  int acc_sign = 1;
+
   int k = 0;
   while (k <= max_steps && log_t >= log(precision)) {
     // Replace zero values with 1 prior to taking the log so that we accumulate
@@ -53,8 +52,7 @@ T_return hypergeometric_3F2_infsum(const Ta& a, const Tb& b, const Tz& z,
     }
 
     log_t += p + log_z;
-    std::forward_as_tuple(t_acc, acc_sign)
-        = log_sum_exp_signed(t_acc, acc_sign, log_t, t_sign);
+    t_acc += t_sign * exp(log_t);
 
     if (is_inf(t_acc)) {
       throw_domain_error("hypergeometric_3F2", "sum (output)", t_acc,
@@ -72,7 +70,7 @@ T_return hypergeometric_3F2_infsum(const Ta& a, const Tb& b, const Tz& z,
                        "exceeded  iterations, hypergeometric function did not ",
                        "converge.");
   }
-  return acc_sign * exp(t_acc);
+  return t_acc;
 }
 }  // namespace internal