diff --git a/src/GaussianEMRestraint.cpp b/src/GaussianEMRestraint.cpp
index b79788c2f139c67cceecea520b09e6281aae1e0f..732359e73daf355480ff96856b93010e17abd038 100755
--- a/src/GaussianEMRestraint.cpp
+++ b/src/GaussianEMRestraint.cpp
@@ -97,7 +97,7 @@ void GaussianEMRestraint::compute_initial_scores() {
        iter != map_score_dd_.end(); ++iter) {
     cached_score_term_ += std::log(iter->second);
   }
-  cached_score_term_ += (3.0 + dsize_) / 2.0 * std::log(2);
+  cached_score_term_ += (3.0 + dsize_) / 2.0 * std::log(2.0);
   cached_score_term_ += 0.5 * std::log(IMP::PI);
   cached_score_term_ -= boost::math::lgamma(dsize_ / 2.0);