diff --git a/exputil/ExpDeproj.cc b/exputil/ExpDeproj.cc
index 3982bf1c0..3a8994696 100644
--- a/exputil/ExpDeproj.cc
+++ b/exputil/ExpDeproj.cc
@@ -48,11 +48,6 @@ double ExpDeproj::density(double R)
 
 double ExpDeproj::mass(double R)
 {
-  // Initialize precomputed arrays?
-  if (mv.size() == 0) {
-    initialize();
-  }
-
   if (R < 0) {
     throw std::invalid_argument("R must be non-negative");
   }
diff --git a/include/ExpDeproj.H b/include/ExpDeproj.H
index 1c9da8a02..a0ff35659 100644
--- a/include/ExpDeproj.H
+++ b/include/ExpDeproj.H
@@ -15,7 +15,7 @@ class ExpDeproj
 
 public:
   //! Constructor
-  ExpDeproj(int n=4000) : ngrid(n) {}
+  ExpDeproj(int n=4000) : ngrid(n) { initialize(); }
 
   //! Destructor
   virtual ~ExpDeproj() {}