message: Begin generalizing types
diff: @@ -53,8 +53,8 @@ namespace primal
  * \param [in] npts the number of Gaussian quadrature nodes for each component
  * \return the value of the integral
  */
  -template <class Lambda>
  -double evaluate_line_integral(const axom::Array<primal::BezierCurve<double, 2>>& cs,
  +template <class Lambda, typename T, int NDIMS>
  +double evaluate_line_integral(const axom::Array<primal::BezierCurve<T, NDIMS>>& cs,
  Lambda&& integrand,
  int npts)
  {
  @@ -88,8 +88,8 @@ double evaluate_line_integral(const axom::Array<primal::BezierCurve<double, 2>>&
  * \param [in] npts_P the number of quadrature points to evaluate the antiderivative
  * \return the value of the integral
  */
  -template <class Lambda>
  -double evaluate_line_integral(const primal::CurvedPolygon<double, 2> cpoly,
  +template <class Lambda, typename T, int NDIMS>
  +double evaluate_line_integral(const primal::CurvedPolygon<T, NDIMS> cpoly,
  Lambda&& integrand,
  int npts)
  {
  @@ -124,8 +124,8 @@ double evaluate_line_integral(const primal::CurvedPolygon<double, 2> cpoly,
  * \param [in] npts the number of quadrature nodes
  * \return the value of the integral
  */
  -template <class Lambda>
  -double evaluate_line_integral(const primal::BezierCurve<double, 2>& c,
  +template <class Lambda, typename T, int NDIMS>
  +double evaluate_line_integral(const primal::BezierCurve<T, NDIMS>& c,
  Lambda&& integrand,
  int npts)
  {
  @@ -154,8 +154,8 @@ double evaluate_line_integral(const primal::BezierCurve<double, 2>& c,
  * \param [in] npts_P the number of quadrature points to evaluate the antiderivative
  * \return the value of the integral
  */
  -template <class Lambda>
  -double evaluate_area_integral(const axom::Array<primal::BezierCurve<double, 2>>& cs,
  +template <class Lambda, typename T, int NDIMS>
  +double evaluate_area_integral(const axom::Array<primal::BezierCurve<T, NDIMS>>& cs,
  Lambda&& integrand,
  int npts_Q,
  int npts_P = 0)
  