diff --git a/Periodic_3_triangulation_3/include/CGAL/Periodic_3_Delaunay_triangulation_3.h b/Periodic_3_triangulation_3/include/CGAL/Periodic_3_Delaunay_triangulation_3.h
index 2bd69020be5..d058a0a5af8 100644
--- a/Periodic_3_triangulation_3/include/CGAL/Periodic_3_Delaunay_triangulation_3.h
+++ b/Periodic_3_triangulation_3/include/CGAL/Periodic_3_Delaunay_triangulation_3.h
@@ -696,33 +696,7 @@ public:
   // and thus cannot be done without changing the Traits concept.
 
   FT dual_volume(Vertex_handle v) const {
-    std::list<Edge> edges;
-    incident_edges(v, std::back_inserter(edges));
-
-    FT vol(0);
-    for (typename std::list<Edge>::iterator eit = edges.begin() ; 
-	 eit != edges.end() ; ++eit) {
-
-      // compute the dual of the edge *eit but handle the translations
-      // with respect to the dual of v. That is why we cannot use one
-      // of the existing dual functions here.
-      Facet_circulator fstart = incident_facets(*eit);
-      Facet_circulator fcit = fstart;
-      std::vector<Point> pts;
-      do {
-	// TODO: possible speed-up by caching the circumcenters
-	Point dual_orig = periodic_circumcenter(fcit->first).first;
-	int idx = fcit->first->index(v);
-	Offset off = periodic_point(fcit->first,idx).second;
-	pts.push_back(point(std::make_pair(dual_orig,-off)));
-	++fcit;
-      } while (fcit != fstart);
-
-      Point orig(0,0,0);
-      for (unsigned int i=1 ; i<pts.size()-1 ; i++) 
-	vol += Tetrahedron(orig,pts[0],pts[i],pts[i+1]).volume();
-    }
-    return vol;
+    return Base::dual_volume(v, geom_traits().construct_circumcenter_3_object());
   }
 
   /// Centroid computations
diff --git a/Periodic_3_triangulation_3/include/CGAL/Periodic_3_regular_triangulation_3.h b/Periodic_3_triangulation_3/include/CGAL/Periodic_3_regular_triangulation_3.h
index b414a171a13..10a6c72298d 100644
--- a/Periodic_3_triangulation_3/include/CGAL/Periodic_3_regular_triangulation_3.h
+++ b/Periodic_3_triangulation_3/include/CGAL/Periodic_3_regular_triangulation_3.h
@@ -6093,6 +6093,15 @@ public:
   Stream& draw_dual(Stream& os) const {
     return Base::draw_dual(os, geom_traits().construct_weighted_circumcenter_3_object());
   }
+
+  /// Volume computations
+
+  // Note: Polygon area computation requires to evaluate square roots
+  // and thus cannot be done without changing the Traits concept.
+
+  FT dual_volume(Vertex_handle v) const {
+    return Base::dual_volume(v, geom_traits().construct_weighted_circumcenter_3_object());
+  }
 //@}
 };
 
diff --git a/Periodic_3_triangulation_3/include/CGAL/Periodic_3_triangulation_3.h b/Periodic_3_triangulation_3/include/CGAL/Periodic_3_triangulation_3.h
index d1734f6cec1..746613b6a16 100644
--- a/Periodic_3_triangulation_3/include/CGAL/Periodic_3_triangulation_3.h
+++ b/Periodic_3_triangulation_3/include/CGAL/Periodic_3_triangulation_3.h
@@ -1651,6 +1651,42 @@ protected:
     CGAL_triangulation_assertion( i == number_of_facets() );
     return os;
   }
+
+  /// Volume computations
+
+  // Note: Polygon area computation requires to evaluate square roots
+  // and thus cannot be done without changing the Traits concept.
+
+  template <class ConstructCircumcenter>
+  FT dual_volume(Vertex_handle v, ConstructCircumcenter construct_circumcenter) const {
+    std::list<Edge> edges;
+    incident_edges(v, std::back_inserter(edges));
+
+    FT vol(0);
+    for (typename std::list<Edge>::iterator eit = edges.begin() ;
+   eit != edges.end() ; ++eit) {
+
+      // compute the dual of the edge *eit but handle the translations
+      // with respect to the dual of v. That is why we cannot use one
+      // of the existing dual functions here.
+      Facet_circulator fstart = incident_facets(*eit);
+      Facet_circulator fcit = fstart;
+      std::vector<Point> pts;
+      do {
+  // TODO: possible speed-up by caching the circumcenters
+  Point dual_orig = periodic_circumcenter(fcit->first, construct_circumcenter).first;
+  int idx = fcit->first->index(v);
+  Offset off = periodic_point(fcit->first,idx).second;
+  pts.push_back(point(std::make_pair(dual_orig,-off)));
+  ++fcit;
+      } while (fcit != fstart);
+
+      Point orig(0,0,0);
+      for (unsigned int i=1 ; i<pts.size()-1 ; i++)
+  vol += Tetrahedron(orig,pts[0],pts[i],pts[i+1]).volume();
+    }
+    return vol;
+  }
 };
 
 template < class GT, class TDS >