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

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