diff --git a/dynare++/integ/cc/quasi_mcarlo.cc b/dynare++/integ/cc/quasi_mcarlo.cc
index 2073c13af..8df7be8f8 100644
--- a/dynare++/integ/cc/quasi_mcarlo.cc
+++ b/dynare++/integ/cc/quasi_mcarlo.cc
@@ -112,7 +112,7 @@ HaltonSequence::HaltonSequence(int n, int mxn, int dim, const PermutationScheme
   // todo: raise if |dim > num_primes|
   // todo: raise if |n > mxn|
   for (int i = 0; i < dim; i++)
-    ri.push_back(RadicalInverse(num, primes[i], maxn));
+    ri.emplace_back(num, primes[i], maxn);
   eval();
 }
 
@@ -123,7 +123,7 @@ HaltonSequence::operator=(const HaltonSequence &hs)
   maxn = hs.maxn;
   ri.clear();
   for (const auto & i : hs.ri)
-    ri.push_back(RadicalInverse(i));
+    ri.emplace_back(i);
   pt = hs.pt;
   return *this;
 }
diff --git a/dynare++/tl/cc/pyramid_prod.cc b/dynare++/tl/cc/pyramid_prod.cc
index a7ba4ba14..30eb37faf 100644
--- a/dynare++/tl/cc/pyramid_prod.cc
+++ b/dynare++/tl/cc/pyramid_prod.cc
@@ -68,7 +68,7 @@ USubTensor::addKronColumn(int i, const vector<const FGSTensor *> &ts,
       IntSequence ind(pindex, lastdim, lastdim+t->dimen());
       lastdim += t->dimen();
       index in(t, ind);
-      tmpcols.push_back(ConstVector(*t, *in));
+      tmpcols.emplace_back(*t, *in);
     }
 
   URSingleTensor kronmult(tmpcols);
diff --git a/dynare++/utils/cc/pascal_triangle.hh b/dynare++/utils/cc/pascal_triangle.hh
index a2c021fda..488e14571 100644
--- a/dynare++/utils/cc/pascal_triangle.hh
+++ b/dynare++/utils/cc/pascal_triangle.hh
@@ -33,7 +33,7 @@ namespace ogu
   public:
     PascalTriangle()
     {
-      tr.push_back(PascalRow());
+      tr.emplace_back();
     }
     PascalTriangle(const PascalTriangle &triang)
       : tr(triang.tr)