MeshFaceCenters.hpp

//---- GPL ---------------------------------------------------------------------
//
// Copyright (C)  Stichting Deltares, 2011-2025.
//
// This program is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation version 3.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program.  If not, see <http://www.gnu.org/licenses/>.
//
// contact: delft3d.support@deltares.nl
// Stichting Deltares
// P.O. Box 177
// 2600 MH Delft, The Netherlands
//
// All indications and logos of, and references to, "Delft3D" and "Deltares"
// are registered trademarks of Stichting Deltares, and remain the property of
// Stichting Deltares. All rights reserved.
//
//------------------------------------------------------------------------------
 
#pragma once
 
#include <array>
#include <span>
#include <vector>
 
#include "MeshKernel/Definitions.hpp"
#include "MeshKernel/Mesh.hpp"
#include "MeshKernel/Point.hpp"
 
namespace meshkernel::algo
{
    std::vector<Point> ComputeFaceCircumcenters(const Mesh& mesh);
 
    void ComputeFaceCircumcenters(const Mesh& mesh, std::span<Point> faceCenters);
 
    Point CircumcenterOfTriangle(const Point& firstNode, const Point& secondNode, const Point& thirdNode, const Projection projection);
 
    Point ComputeFaceCircumenter(std::vector<Point>& polygon,
                                 const std::vector<UInt>& edgesNumFaces,
                                 const Projection projection);
 
    Point ComputeCircumCenter(const Point& centerOfMass,
                              const UInt pointCount,
                              const std::array<Point, constants::geometric::maximumNumberOfNodesPerFace>& middlePoints,
                              const std::array<Point, constants::geometric::maximumNumberOfNodesPerFace>& normals,
                              const Projection projection);
 
} // namespace meshkernel::algo