Contacts.hpp

//---- GPL ---------------------------------------------------------------------
//
// Copyright (C)  Stichting Deltares, 2011-2020.
//
// 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 <MeshKernel/Mesh1D.hpp>
#include <MeshKernel/Mesh2D.hpp>
 
namespace meshkernel
{
    class Contacts
    {
    public:
        Contacts(Mesh1D& mesh1d, Mesh2D& mesh2d);
 
        void ComputeSingleContacts(const std::vector<bool>& oneDNodeMask, const Polygons& polygons, double projectionFactor);
 
        void ComputeMultipleContacts(const std::vector<bool>& oneDNodeMask);
 
        void ComputeContactsWithPolygons(const std::vector<bool>& oneDNodeMask, const Polygons& polygons);
 
        void ComputeContactsWithPoints(const std::vector<bool>& oneDNodeMask, const std::vector<Point>& points);
 
        void ComputeBoundaryContacts(const std::vector<bool>& oneDNodeMask,
                                     const Polygons& polygons,
                                     double searchRadius);
 
        std::vector<UInt> const& Mesh1dIndices() const { return m_mesh1dIndices; }
 
        std::vector<UInt> const& Mesh2dIndices() const { return m_mesh2dIndices; }
 
        void SetIndices(const std::vector<meshkernel::UInt>& mesh1dIndices,
                        const std::vector<meshkernel::UInt>& mesh2dIndices);
 
    private:
        [[nodiscard]] bool IsContactIntersectingMesh1d(UInt node, UInt face) const;
 
        [[nodiscard]] bool IsContactIntersectingContact(UInt node, UInt face) const;
 
        void Connect1dNodesWithCrossingFaces(UInt node,
                                             double projectionFactor);
 
        void Validate() const;
 
        Mesh1D& m_mesh1d;                        
        Mesh2D& m_mesh2d;                        
        std::vector<UInt> m_mesh1dIndices;       
        std::vector<UInt> m_mesh2dIndices;       
        std::vector<Point> m_facesCircumcenters; 
    };
} // namespace meshkernel