diff --git a/MathLib/Point3d.h b/MathLib/Point3d.h
index ebc86d6fbe1a611179a05030629bf6bd9648cdde..d59fc2ab998c537e6808514599d15bb5c3736945 100644
--- a/MathLib/Point3d.h
+++ b/MathLib/Point3d.h
@@ -53,6 +53,7 @@ inline MathLib::Point3d operator*(MATRIX const& mat, MathLib::Point3d const& p)
 {
     MathLib::Point3d new_p;
     for (std::size_t i(0); i<3; ++i) {
+        new_p[i] = 0.0;
         for (std::size_t j(0); j<3; ++j) {
             new_p[i] += mat(i,j)*p[j];
         }