docs: def gaussian_pdf(x, g):
          """ Gaussian probability density function at ``x`` for |GVar| ``g``. """
          return (
              numpy.exp(-(x - g.mean) ** 2 / 2. /g.var) /
              numpy.sqrt(g.var * 2 * numpy.pi)
              )