@@ -1404,11 +1404,11 @@
 /* INDENT ON */
 
 /* gamma(x+i) for 0 <= x < 1  */
 static struct Double
 gam_n(int i, double x) {
-        struct Double rr, yy;
+        struct Double rr = {0.0L, 0.0L}, yy;
         double r1, r2, t2, z, xh, xl, yh, yl, zh, z1, z2, zl, x5, wh, wl;
 
         /* compute yy = gamma(x+1) */
         if (x > 0.2845) {
                 if (x > 0.6374) {

@@ -1556,11 +1556,11 @@
                 ww = large_gam(x, &m);
                 w = ww.h + ww.l;
                 __HI(w) += m << 20;
                 return (w);
         }
-        if (hx > 0) {           /* x from 0 to 8 */
+        if (hx > 0) {           /* 0 < x < 8 */
                 i = (int) x;
                 ww = gam_n(i, x - (double) i);
                 return (ww.h + ww.l);
         }