Skip to content

Commit

Permalink
remove unused functions
Browse files Browse the repository at this point in the history
  • Loading branch information
Alexander Voigt authored and Alexander Voigt committed Aug 16, 2016
1 parent 3aa01a8 commit ec35ffd
Show file tree
Hide file tree
Showing 2 changed files with 0 additions and 24 deletions.
16 changes: 0 additions & 16 deletions src/numerics.cpp
Expand Up @@ -19,22 +19,6 @@
using namespace softsusy;
using namespace Eigen;

double lnLPoisson(unsigned k, double lambda) {
if (lambda < 6.66e-66) throw("In lnPoisson: lambda not positive\n");
double lnL = -lambda + double(k) * log(lambda);
for (unsigned i=2; i<=k; i++) lnL -= log(double(i));

return lnL;
}

double LPoisson(unsigned k, double lambda) {
if (lambda < 6.66e-66) throw("In lnPoisson: lambda not positive\n");
double L = exp(-lambda) * pow(lambda,int(k));
for (unsigned i=2; i<=k; i++) L /= double(i);

return L;
}

double calcDerivative(double (*func)(double), double x, double h, double
*err){
const double CON = 1.4, CON2 = CON * CON, BIG = 1.0e30,
Expand Down
8 changes: 0 additions & 8 deletions src/numerics.h
Expand Up @@ -26,14 +26,6 @@
/// calculates root(1+x), where x<<1 accurately
double accurateSqrt1Plusx(double x);

/// Calculates log likelihood of a Poisson with k observed events, expecting
/// lambda>0.
double lnLPoisson(unsigned k, double lambda);

/// Calculates likelihood of a Poisson with k observed events, expecting
/// lambda>0.
double LPoisson(unsigned k, double lambda);

/// func is user-supplied, h is an estimate of what step-size to start with
/// and err returns error flags
double calcDerivative(double (*func)(double),
Expand Down

0 comments on commit ec35ffd

Please sign in to comment.