diff --git a/fdaPDE/pde/pde.h b/fdaPDE/pde/pde.h
index b90b6658..579d9859 100644
--- a/fdaPDE/pde/pde.h
+++ b/fdaPDE/pde/pde.h
@@ -75,13 +75,14 @@ class PDE : public PDEBase {
     PDE(const D& domain) : domain_(domain) { }
     PDE(const D& domain, const DVector<double>& time) : domain_(domain), time_(time) { };
     PDE(const D& domain, E diff_op) : domain_(domain), diff_op_(diff_op) { };
-    PDE(const D& domain, E diff_op, const DVector<double>& time) : domain_(domain), diff_op_(diff_op), time_(time) { };
+    PDE(const D& domain, const DVector<double>& time, E diff_op) : domain_(domain), time_(time), diff_op_(diff_op) { };
     void set_forcing(const F& forcing_data) { forcing_data_ = forcing_data; }
     void set_differential_operator(E diff_op) { diff_op_ = diff_op; }
     // full constructors
     PDE(const D& domain, E diff_op, const F& forcing_data) :
         domain_(domain), diff_op_(diff_op), forcing_data_(forcing_data) { }
-
+    PDE(const D& domain, const DVector<double>& time, E diff_op, const F& forcing_data) :
+        domain_(domain), time_(time), diff_op_(diff_op), forcing_data_(forcing_data) { }
     // setters
     virtual void set_dirichlet_bc(const DMatrix<double>& data) { boundary_data_ = data; }
     virtual void set_initial_condition(const DVector<double>& data) { initial_condition_ = data; };
diff --git a/fdaPDE/utils/type_erasure.h b/fdaPDE/utils/type_erasure.h
index 6438a52f..62160833 100644
--- a/fdaPDE/utils/type_erasure.h
+++ b/fdaPDE/utils/type_erasure.h
@@ -16,7 +16,7 @@
 
 #ifndef __TYPE_ERASURE_H__
 #define __TYPE_ERASURE_H__
-
+#include<utility>
 #include "assert.h"
 #include "traits.h"