[geom][cylsurface] reintroduce area computation

frackit/geometry/cylindersurface.hh

@@ -106,6 +106,8 @@ public:
     ctype height() const { return height_; }
     //! Return the radius of the cylinder
     ctype radius() const { return bottom_.radius(); }
+    //! Return the area of this surface
+    ctype area() { return 2.0*M_PI*radius()*height(); }
 
     //! Return the circle describing the top boundary of this surface
     const Circle& upperBoundingCircle() const { return top_; }

frackit/magnitude/area.hh

@@ -42,12 +42,6 @@ template<class Geometry>
 typename Geometry::ctype computeArea(const Geometry& geometry)
 { return geometry.area(); }
 
-//! \todo TODO doc me.
-template<class ctype>
-ctype computeArea(const CylinderSurface<ctype>& cylSurface)
-{ return 2.0*M_PI*cylSurface.radius()*cylSurface.height(); }
-
-//! \todo TODO doc me.
 template<class ctype = double>
 ctype computeArea(const TopoDS_Face& face,
                   ctype eps = Precision<ctype>::confusion(),