diff --git a/test/localgeodesicfefunctiontest.cc b/test/localgeodesicfefunctiontest.cc
index aea15a5580e3e36ef2a85cf7c732dc9a6fd4671b..37f874f4a1e3aef70e84610db620b1cd68cfc7b8 100644
--- a/test/localgeodesicfefunctiontest.cc
+++ b/test/localgeodesicfefunctiontest.cc
@@ -342,39 +342,34 @@ int main()
 
     std::cout << std::setw(15) << std::setprecision(12);
     
-    GeometryType element;
-
     ////////////////////////////////////////////////////////////////
     //  Test functions on 1d elements
     ////////////////////////////////////////////////////////////////
-    element.makeSimplex(1);
     
-    test<RealTuple<double,1>,1>(element);
-    test<UnitVector<double,2>,1>(element);
-    test<UnitVector<double,3>,1>(element);
-    test<Rotation<double,3>,1>(element);
-    test<RigidBodyMotion<double,3>,1>(element);
+    test<RealTuple<double,1>,1>(GeometryTypes::simplex(1));
+    test<UnitVector<double,2>,1>(GeometryTypes::simplex(1));
+    test<UnitVector<double,3>,1>(GeometryTypes::simplex(1));
+    test<Rotation<double,3>,1>(GeometryTypes::simplex(1));
+    test<RigidBodyMotion<double,3>,1>(GeometryTypes::simplex(1));
     
     ////////////////////////////////////////////////////////////////
     //  Test functions on 2d simplex elements
     ////////////////////////////////////////////////////////////////
-    element.makeSimplex(2);
 
-    test<RealTuple<double,1>,2>(element);
-    test<UnitVector<double,2>,2>(element);
-    test<UnitVector<double,3>,2>(element);
-    test<Rotation<double,3>,2>(element);
-    test<RigidBodyMotion<double,3>,2>(element);
+    test<RealTuple<double,1>,2>(GeometryTypes::simplex(2));
+    test<UnitVector<double,2>,2>(GeometryTypes::simplex(2));
+    test<UnitVector<double,3>,2>(GeometryTypes::simplex(2));
+    test<Rotation<double,3>,2>(GeometryTypes::simplex(2));
+    test<RigidBodyMotion<double,3>,2>(GeometryTypes::simplex(2));
 
     ////////////////////////////////////////////////////////////////
     //  Test functions on 2d quadrilateral elements
     ////////////////////////////////////////////////////////////////
-    element.makeCube(2);
 
-    test<RealTuple<double,1>,2>(element);
-    test<UnitVector<double,2>,2>(element);
-    test<UnitVector<double,3>,2>(element);
-    test<Rotation<double,3>,2>(element);
-    test<RigidBodyMotion<double,3>,2>(element);
+    test<RealTuple<double,1>,2>(GeometryTypes::cube(2));
+    test<UnitVector<double,2>,2>(GeometryTypes::cube(2));
+    test<UnitVector<double,3>,2>(GeometryTypes::cube(2));
+    test<Rotation<double,3>,2>(GeometryTypes::cube(2));
+    test<RigidBodyMotion<double,3>,2>(GeometryTypes::cube(2));
 
 }