MortarSegmentHelper.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html

#pragma once

#include "MortarSegmentInfo.h"

#include "libmesh/point.h"
#include "libmesh/int_range.h"

#include <vector>

using libMesh::Point;
using libMesh::Real;

class MortarSegmentHelper
{
public:
  MortarSegmentHelper(const std::vector<Point> secondary_nodes,
                      const Point & center,
                      const Point & normal,
                      const MortarSegmentTriangulationMode triangulation_mode,
                      const bool triangulate_triangles);

  Point getIntersection(
      const Point & p1, const Point & p2, const Point & q1, const Point & q2, Real & s) const;

  bool isInsideSecondary(const Point & pt) const;

  bool isDisjoint(const std::vector<Point> & poly) const;

  std::vector<Point> projectPrimaryPoly(const std::vector<Point> & primary_nodes) const;

  std::vector<Point> clipPoly(const std::vector<Point> & primary_nodes) const;

  void triangulatePoly(std::vector<Point> & poly_nodes,
                       std::vector<std::vector<unsigned int>> & tri_map) const;

  void getMortarSegments(const std::vector<Point> & primary_nodes,
                         std::vector<Point> & nodes,
                         std::vector<std::vector<unsigned int>> & elem_to_nodes);

  Real area(const std::vector<Point> & nodes) const;

  const Point & center() const { return _center; }

  Real remainder() const { return _remaining_area_fraction; }

  Point point(unsigned int i) const
  {
    return (_secondary_poly[i](0) * _u) + (_secondary_poly[i](1) * _v) + _center;
  }

private:
  Point _center;

  Point _normal;

  Point _u, _v;

  Real _secondary_area;

  Real _remaining_area_fraction;

  bool _debug;

  Real _tolerance = 1e-8;

  const MortarSegmentTriangulationMode _triangulation_mode;

  const bool _triangulate_triangles;

  Real _area_tol;

  Real _length_tol;

  std::vector<Point> _secondary_poly;
};