Fix missing mul

Author: andrewjradcliffe

src/e_lgammaf_r.jl

@@ -57,7 +57,7 @@ function _lgammaf_r(x::Float32)
             y = 1.0f0 + c - x
             z = y * y
             p1 = evalpoly(z, (7.7215664089f-2, 6.7352302372f-2, 7.3855509982f-3, 1.1927076848f-3, 2.2086278477f-4, 2.5214456400f-5))
-            p2 = evalpoly(z, (3.2246702909f-1, 2.0580807701f-2, 2.8905137442f-3, 5.1006977446f-4, 1.0801156895f-4, 4.4864096708f-5))
+            p2 = z * evalpoly(z, (3.2246702909f-1, 2.0580807701f-2, 2.8905137442f-3, 5.1006977446f-4, 1.0801156895f-4, 4.4864096708f-5))
             p = muladd(p1, y, p2)
             r += muladd(y, -0.5f0, p)
         elseif f ≥ 0.23163998f0 # or, the lb? 0.2316322f0