diff --git a/test/localgeodesicfefunctiontest.cc b/test/localgeodesicfefunctiontest.cc
index 30ec726bdc110a316f584f7184faf4af18b20bb4..a2fa30ba11d9fd40d18180f39a36141ee3fef364 100644
--- a/test/localgeodesicfefunctiontest.cc
+++ b/test/localgeodesicfefunctiontest.cc
@@ -76,8 +76,7 @@ void testPermutationInvariance(const std::vector<TargetSpace>& corners)
     PQkLocalFiniteElementCache<double,double,domainDim,1> feCache;
     typedef typename PQkLocalFiniteElementCache<double,double,domainDim,1>::FiniteElementType LocalFiniteElement;
     
-    GeometryType simplex;
-    simplex.makeSimplex(domainDim);
+    GeometryType simplex = GeometryTypes::simplex(domainDim);
 
     //
     std::vector<TargetSpace> cornersRotated1(domainDim+1);
diff --git a/test/localprojectedfefunctiontest.cc b/test/localprojectedfefunctiontest.cc
index 7d266d6240c26a01f6ebd843f65ac544936f8420..d2f6829398d13b3c60ba4052988e38ad8f701a7d 100644
--- a/test/localprojectedfefunctiontest.cc
+++ b/test/localprojectedfefunctiontest.cc
@@ -118,8 +118,7 @@ void testPermutationInvariance(const std::vector<TargetSpace>& corners)
     PQkLocalFiniteElementCache<double,double,domainDim,1> feCache;
     typedef typename PQkLocalFiniteElementCache<double,double,domainDim,1>::FiniteElementType LocalFiniteElement;
 
-    GeometryType simplex;
-    simplex.makeSimplex(domainDim);
+    GeometryType simplex = GeometryTypes::simplex(domainDim);
 
     //
     std::vector<TargetSpace> cornersRotated1(domainDim+1);