Made corner angle computation independent of differentiable structure

Author: rjc8237

src/optimization/cone_metric.cpp

@@ -190,7 +190,6 @@ void PennerConeMetric::make_discrete_metric()
     SolveStats<Scalar> solve_stats;
     bool use_ptolemy = true;
     ConformalIdealDelaunay<Scalar>::MakeDelaunay(*this, u, del_stats, solve_stats, use_ptolemy);
-    m_flip_seq = del_stats.flip_seq; // TODO Append
     m_is_discrete_metric = true;
     return;
 }

src/optimization/constraint.cpp

@@ -5,7 +5,7 @@
 
 namespace CurvatureMetric {
 
-bool satisfies_triangle_inequality(const DifferentiableConeMetric& cone_metric)
+bool satisfies_triangle_inequality(const Mesh<Scalar>& cone_metric)
 {
     int num_halfedges = cone_metric.n_halfedges();
 
@@ -28,7 +28,7 @@ bool satisfies_triangle_inequality(const DifferentiableConeMetric& cone_metric)
     return true;
 }
 
-void corner_angles(const DifferentiableConeMetric& cone_metric, VectorX& he2angle, VectorX& he2cot)
+void corner_angles(const Mesh<Scalar>& cone_metric, VectorX& he2angle, VectorX& he2cot)
 {
     int num_halfedges = cone_metric.n_halfedges();
     he2angle.setZero(num_halfedges);

src/optimization/constraint.hh

@@ -8,17 +8,17 @@ namespace CurvatureMetric {
 /// Check the triangle inequality for every triangle in the mesh with respect to the
 /// halfedge metric coordinates
 ///
-/// @param[in] cone_metric: mesh with differentiable metric
+/// @param[in] cone_metric: mesh with metric
 /// @return: true iff each triangle in the mesh satisfies the triangle inequality
-bool satisfies_triangle_inequality(const DifferentiableConeMetric& cone_metric);
+bool satisfies_triangle_inequality(const Mesh<Scalar>& cone_metric);
 
 /// Compute the triangle angles and cotangent angles of a Delaunay (possibly
 /// symmetric) mesh m with metric. Angles are indexed by their opposing halfedge.
 ///
-/// @param[in] cone_metric: mesh with differentiable metric
+/// @param[in] cone_metric: mesh with metric
 /// @param[out] he2angle: map from halfedges to opposing angle
 /// @param[out] he2cot: map from halfedges to cotan of opposing angle
-void corner_angles(const DifferentiableConeMetric& cone_metric, VectorX& he2angle, VectorX& he2cot);
+void corner_angles(const Mesh<Scalar>& cone_metric, VectorX& he2angle, VectorX& he2cot);
 
 /// Compute the vertex angles of a Delaunay (possibly symmetric) mesh with metric
 /// and the Jacobian with respect to the reduced coordinates if needed.