manifold_lib.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 2014 - 2023 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Part of the source code is dual licensed under Apache-2.0 WITH
// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
// governing the source code and code contributions can be found in
// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
//
// ------------------------------------------------------------------------
 
 
#ifndef dealii_occ_manifold_lib_h
#define dealii_occ_manifold_lib_h
 
#include <deal.II/base/config.h>
 
#ifdef DEAL_II_WITH_OPENCASCADE
 
#  include <deal.II/grid/manifold.h>
 
#  include <deal.II/opencascade/utilities.h>
 
// opencascade needs "HAVE_CONFIG_H" to be exported...
#  define HAVE_CONFIG_H
#  include <Adaptor3d_Curve.hxx>
#  if !DEAL_II_OPENCASCADE_VERSION_GTE(7, 6, 0)
#    include <Adaptor3d_HCurve.hxx>
#  endif
#  include <BRepAdaptor_Curve.hxx>
#  undef HAVE_CONFIG_H
 
DEAL_II_NAMESPACE_OPEN
 
namespace OpenCASCADE
{
  template <int dim, int spacedim>
  class NormalProjectionManifold : public FlatManifold<dim, spacedim>
  {
  public:
    NormalProjectionManifold(const TopoDS_Shape &sh,
                             const double        tolerance = 1e-7);
 
    virtual std::unique_ptr<Manifold<dim, spacedim>>
    clone() const override;
 
    virtual Point<spacedim>
    project_to_manifold(
      const ArrayView<const Point<spacedim>> &surrounding_points,
      const Point<spacedim>                  &candidate) const override;
 
 
  protected:
    const TopoDS_Shape sh;
 
    const double tolerance;
  };
 
  template <int dim, int spacedim>
  class DirectionalProjectionManifold : public FlatManifold<dim, spacedim>
  {
  public:
    DirectionalProjectionManifold(const TopoDS_Shape        &sh,
                                  const Tensor<1, spacedim> &direction,
                                  const double               tolerance = 1e-7);
 
    virtual std::unique_ptr<Manifold<dim, spacedim>>
    clone() const override;
 
    virtual Point<spacedim>
    project_to_manifold(
      const ArrayView<const Point<spacedim>> &surrounding_points,
      const Point<spacedim>                  &candidate) const override;
 
  protected:
    const TopoDS_Shape sh;
 
    const Tensor<1, spacedim> direction;
 
    const double tolerance;
  };
 
 
  template <int dim, int spacedim>
  class NormalToMeshProjectionManifold : public FlatManifold<dim, spacedim>
  {
  public:
    NormalToMeshProjectionManifold(const TopoDS_Shape &sh,
                                   const double        tolerance = 1e-7);
 
    virtual std::unique_ptr<Manifold<dim, spacedim>>
    clone() const override;
 
    virtual Point<spacedim>
    project_to_manifold(
      const ArrayView<const Point<spacedim>> &surrounding_points,
      const Point<spacedim>                  &candidate) const override;
 
  protected:
    const TopoDS_Shape sh;
 
    const double tolerance;
  };
 
  template <int dim, int spacedim>
  class ArclengthProjectionLineManifold : public ChartManifold<dim, spacedim, 1>
  {
  public:
    ArclengthProjectionLineManifold(const TopoDS_Shape &sh,
                                    const double        tolerance = 1e-7);
 
    virtual std::unique_ptr<Manifold<dim, spacedim>>
    clone() const override;
 
    virtual Point<1>
    pull_back(const Point<spacedim> &space_point) const override;
 
    virtual Point<spacedim>
    push_forward(const Point<1> &chart_point) const override;
 
  protected:
    const TopoDS_Shape sh;
 
#  if DEAL_II_OPENCASCADE_VERSION_GTE(7, 6, 0)
    Handle_Adaptor3d_Curve curve;
#  else
    Handle_Adaptor3d_HCurve curve;
#  endif
 
    const double tolerance;
 
    const double length;
  };
 
  template <int dim, int spacedim>
  class NURBSPatchManifold : public ChartManifold<dim, spacedim, 2>
  {
  public:
    NURBSPatchManifold(const TopoDS_Face &face, const double tolerance = 1e-7);
 
    virtual std::unique_ptr<Manifold<dim, spacedim>>
    clone() const override;
 
    virtual Point<2>
    pull_back(const Point<spacedim> &space_point) const override;
 
    virtual Point<spacedim>
    push_forward(const Point<2> &chart_point) const override;
 
    virtual DerivativeForm<1, 2, spacedim>
    push_forward_gradient(const Point<2> &chart_point) const override;
 
  protected:
    std::tuple<double, double, double, double>
    get_uv_bounds() const;
 
    TopoDS_Face face;
 
    double tolerance;
  };
 
} // namespace OpenCASCADE
 
DEAL_II_NAMESPACE_CLOSE
 
 
#endif // DEAL_II_WITH_OPENCASCADE
#endif // dealii_occ_manifold_lib_h