fix a minor bug in hypgeom_m (nan result only setting real part to nan)

fredrik-johansson · fredrik-johansson · commit c0965c3733ee · 2022-01-24T16:16:55.000+01:00
diff --git a/acb_hypgeom/m.c b/acb_hypgeom/m.c
@@ -56,7 +56,11 @@ acb_hypgeom_m_asymp(acb_t res, const acb_t a, const acb_t b, const acb_t z, int
     if (!regularized)
     {
         acb_gamma(v, b, prec);
-        acb_mul(t, t, v, prec);
+
+        if (acb_is_finite(v))
+            acb_mul(t, t, v, prec);
+        else
+            acb_indeterminate(t);
     }
 
     if (acb_is_real(a) && acb_is_real(b) && acb_is_real(z))

Original file line number	Diff line number	Diff line change
`@@ -56,7 +56,11 @@ acb_hypgeom_m_asymp(acb_t res, const acb_t a, const acb_t b, const acb_t z, int`
`56`	`56`	`if (!regularized)`
`57`	`57`	`{`
`58`	`58`	`acb_gamma(v, b, prec);`
`59`		`- acb_mul(t, t, v, prec);`
	`59`	`+`
	`60`	`+ if (acb_is_finite(v))`
	`61`	`+ acb_mul(t, t, v, prec);`
	`62`	`+ else`
	`63`	`+ acb_indeterminate(t);`
`60`	`64`	`}`
`61`	`65`
`62`	`66`	`if (acb_is_real(a) && acb_is_real(b) && acb_is_real(z))`