Sometimes a value gets negative zero causing a NaN on the sqrt.

Closes #61

1 files changed, 1 insertions, 0 deletions

diff --git a/src/mathfunc.h b/src/mathfunc.h
index 5b1e2ec..e786e87 100644
--- a/src/mathfunc.h
+++ b/src/mathfunc.h
@@ -46,6 +46,7 @@ bool has_inf(const gsl_matrix *m);
 bool is_integer(const std::string & s);
 
 double safe_log(const double d);
+double safe_sqrt(const double d);
 
 double VectorVar(const gsl_vector *v);
 void CenterMatrix(gsl_matrix *G);