message: Bugfix for XL for multiply-defined template functions that are fully specialized
  Need to either inline or move to cpp file. I chose the former here.
diff: @@ -776,7 +776,7 @@ constexpr typename Delaunay<DIM>::IndexType Delaunay<DIM>::INVALID_INDEX;
  // 2D specialization for generateInitialMesh(...)
  template <>
  -void Delaunay<2>::generateInitialMesh(std::vector<DataType>& points,
  +inline void Delaunay<2>::generateInitialMesh(std::vector<DataType>& points,
  std::vector<IndexType>& elem,
  const BoundingBox& bb)
  {
  @@ -801,7 +801,7 @@ void Delaunay<2>::generateInitialMesh(std::vector<DataType>& points,
  // 3D specialization for generateInitialMesh(...)
  template <>
  -void Delaunay<3>::generateInitialMesh(std::vector<DataType>& points,
  +inline void Delaunay<3>::generateInitialMesh(std::vector<DataType>& points,
  std::vector<IndexType>& elem,
  const BoundingBox& bb)
  {
  @@ -833,7 +833,8 @@ void Delaunay<3>::generateInitialMesh(std::vector<DataType>& points,
  // 2D specialization for getBaryCoords(...)
  template <>
  -Delaunay<2>::BaryCoordType Delaunay<2>::getBaryCoords(IndexType element_idx,
  +inline Delaunay<2>::BaryCoordType Delaunay<2>::getBaryCoords(
  + IndexType element_idx,
  const PointType& query_pt) const
  {
  const auto verts = m_mesh.boundaryVertices(element_idx);
  @@ -846,7 +847,8 @@ Delaunay<2>::BaryCoordType Delaunay<2>::getBaryCoords(IndexType element_idx,
  // 3D specialization for getBaryCoords(...)
  template <>
  -Delaunay<3>::BaryCoordType Delaunay<3>::getBaryCoords(IndexType element_idx,
  +inline Delaunay<3>::BaryCoordType Delaunay<3>::getBaryCoords(
  + IndexType element_idx,
  const PointType& query_pt) const
  {
  const auto verts = m_mesh.boundaryVertices(element_idx);
  @@ -860,7 +862,7 @@ Delaunay<3>::BaryCoordType Delaunay<3>::getBaryCoords(IndexType element_idx,
  // 2D specialization for isPointInSphere(...)
  template <>
  -bool Delaunay<2>::InsertionHelper::isPointInSphere(const PointType& query_pt,
  +inline bool Delaunay<2>::InsertionHelper::isPointInSphere(const PointType& query_pt,
  IndexType element_idx) const
  {
  const auto verts = m_mesh.boundaryVertices(element_idx);
  @@ -872,7 +874,7 @@ bool Delaunay<2>::InsertionHelper::isPointInSphere(const PointType& query_pt,
  // 3D specialization for isPointInSphere(...)
  template <>
  -bool Delaunay<3>::InsertionHelper::isPointInSphere(const PointType& query_pt,
  +inline bool Delaunay<3>::InsertionHelper::isPointInSphere(const PointType& query_pt,
  IndexType element_idx) const
  {
  const auto verts = m_mesh.boundaryVertices(element_idx);
  