# HG changeset patch
# User alpar
# Date 1207224408 0
# Node ID b9905565d18539111ad8a1a0b0fbf6fe4928371d
# Parent  71f4bd3a9ae8f965d265b88f11e6904edb7f2b47
Serious bugfix in Random::gamma()
Backport from hg b6bede534255

diff -r 71f4bd3a9ae8 -r b9905565d185 lemon/random.h
--- a/lemon/random.h	Tue Mar 25 16:28:06 2008 +0000
+++ b/lemon/random.h	Thu Apr 03 12:06:48 2008 +0000
@@ -802,7 +802,7 @@
 	    nu=V0*std::exp(-xi);
 	  }
       } while(nu>std::pow(xi,delta-1.0)*std::exp(-xi));
-      return theta*(xi-gamma(int(std::floor(k))));
+      return theta*(xi+gamma(int(std::floor(k))));
     }