subroutine dlgams (x, dlgam, sgngam)
c july 1977 edition.  w. fullerton, c3, los alamos scientific lab.
c
c evaluate log abs (gamma(x)) and return the sign of gamma(x) in sgngam.
c sgngam is either +1.0 or -1.0.
c
      double precision x, dlgam, sgngam, dint, dlngam
      external  dlngam
c
      dlgam = dlngam(x)
      sgngam = 1.0d0
      if (x.gt.0.d0) return
c
      int = dmod (-dint(x), 2.0d0) + 0.1d0
      if (int.eq.0) sgngam = -1.0d0
c
      return
      end