source: PSPA/parmelaPSPA/trunk/lodgau.f @ 488

      subroutine lodgau(nbase,n,x,y,f,d,j)
c-----------------------------------------------------------------------
c     
c     load a gaussian distribution using the hammersley's sequence.
c     (nbase=0 ==> random sampling !)
c     
c-----------------------------------------------------------------------
      dimension       x(*), y(*), f(*), d(*), j(*)
      data            xbig/5./
c-----------------------------------------------------------------------
c     
c     cumulative distribution of gaussian
c     
      delx = 2.*xbig/(n-1)
      xk = -xbig-0.5*delx
      x(1) = 0.0
      do 10 i=2,n
         xk = xk+delx
         x(i) = x(i-1) + exp(-0.5*xk*xk)
 10   continue
      zcons = 1./x(n)
      call sscal(n,zcons,x,1)
c     
c     uniform sampling
c     
      call loduni(nbase,n,y)
c     
c     inverse linear interpolation
c     
      call srchfd(n,y,n,x,j,d,f)
      do 200 i=1,n
         x1 = -xbig+j(i)*delx
         x0 = x1-delx
 200  x(i) = x0*(1.-f(i)) + x1*f(i)
c     
      return
      end
c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++*