Changeset 116:b6bede534255 in lemon-main

Timestamp:

04/03/08 14:00:18 (17 years ago)

Author:

Alpar Juttner <alpar@…>

Branch:

default

Children:

117:7b0ce9fb1169, 123:8899d1891a3c

Phase:

public

Message:

Serious buxfixes in Random::gamma() and Random::pareto()

File:

• lemon/random.h (modified) (2 diffs)

lemon/random.h

@@ -796 +796 @@
           }
       } 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))));
     }
     
@@ -820 +820 @@
     double pareto(double k,double x_min)
     {
-      return exponential(gamma(k,1.0/x_min));
+      return exponential(gamma(k,1.0/x_min))+x_min;
     }