boundary_lib.h

// ---------------------------------------------------------------------
//
// Copyright (C) 2014 - 2018 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE at
// the top level of the deal.II distribution.
//
// ---------------------------------------------------------------------


#ifndef dealii_occ_boundary_lib_h
#define dealii_occ_boundary_lib_h

#include <deal.II/base/config.h>

#ifdef DEAL_II_WITH_OPENCASCADE

#include <deal.II/opencascade/utilities.h>
#include <deal.II/grid/manifold.h>

// opencascade needs "HAVE_CONFIG_H" to be exported...
#define HAVE_CONFIG_H
#include <Adaptor3d_Curve.hxx>
#include <Adaptor3d_HCurve.hxx>
#include <BRepAdaptor_Curve.hxx>
#undef HAVE_CONFIG_H

DEAL_II_NAMESPACE_OPEN

namespace OpenCASCADE
{
  template <int dim, int spacedim>
  class NormalProjectionBoundary : public FlatManifold<dim,spacedim>
  {
  public:

    NormalProjectionBoundary(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;


  private:
    const TopoDS_Shape sh;

    const double tolerance;
  };

  template <int dim, int spacedim>
  class DirectionalProjectionBoundary : public FlatManifold<dim,spacedim>
  {
  public:
    DirectionalProjectionBoundary(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;

  private:
    const TopoDS_Shape sh;

    const Tensor<1, spacedim> direction;

    const double tolerance;
  };


  template <int dim, int spacedim>
  class NormalToMeshProjectionBoundary : public FlatManifold<dim,spacedim>
  {
  public:
    NormalToMeshProjectionBoundary(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;

  private:
    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;

  private:
    const TopoDS_Shape sh;

    Handle_Adaptor3d_HCurve curve;

    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;

  private:
    std::tuple<double, double, double, double>
    get_uv_bounds() const;

    TopoDS_Face face;

    double tolerance;
  };

}

DEAL_II_NAMESPACE_CLOSE


#endif // DEAL_II_WITH_OPENCASCADE

/*------------------------------ occ_boundary_lib.h ------------------------------*/
#endif
/*------------------------------ occ_boundary_lib.h ------------------------------*/