geometryreference.hh

Go to the documentation of this file.

// SPDX-FileCopyrightText: Copyright © DUNE Project contributors, see file LICENSE.md in module root
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception
// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_GRID_ALBERTAGRID_GEOMETRYREFERENCE_HH
#define DUNE_GRID_ALBERTAGRID_GEOMETRYREFERENCE_HH
 
#include <dune/common/typetraits.hh>
 
#include <dune/geometry/type.hh>
 
#include <dune/grid/common/geometry.hh>
 
namespace Dune
{
 
  // GeometryReference
  // -----------------
 
  template< class Implementation >
  class GeometryReference
  {
    typedef GeometryReference< Implementation > This;
 
  public:
    static const int mydimension = Implementation::mydimension;
    static const int coorddimension = Implementation::coorddimension;
 
    typedef typename Implementation::ctype ctype;
 
    typedef typename Implementation::LocalCoordinate LocalCoordinate;
    typedef typename Implementation::GlobalCoordinate GlobalCoordinate;
 
    typedef typename Implementation::JacobianInverseTransposed JacobianInverseTransposed;
    typedef typename Implementation::JacobianTransposed JacobianTransposed;
 
  private:
 
    template<class Implementation_T>
    using JacobianInverseOfImplementation = decltype(typename Implementation_T::JacobianInverse{std::declval<Implementation_T>().jacobianInverse(std::declval<LocalCoordinate>())});
 
    using JacobianInverseDefault = decltype(transpose(std::declval<JacobianInverseTransposed>()));
 
    template<class Implementation_T>
    using JacobianOfImplementation = decltype(typename Implementation_T::Jacobian{std::declval<Implementation_T>().jacobian(std::declval<LocalCoordinate>())});
 
    using JacobianDefault = decltype(transpose(std::declval<JacobianTransposed>()));
 
 
    template <class I = Implementation>
    [[deprecated("Geometry implementatons are required to provide a jacobian(local) method. The default implementation is deprecated and will be removed after release 2.9")]]
    auto deprecatedDefaultJacobian ( const LocalCoordinate& local ) const {
      return transpose(jacobianTransposed(local));
    }
 
    template <class I = Implementation>
    [[deprecated("Geometry implementatons are required to provide a jacobianInverse(local) method. The default implementation is deprecated and will be removed after release 2.9")]]
    auto deprecatedDefaultJacobianInverse ( const LocalCoordinate& local ) const {
      return transpose(jacobianInverseTransposed(local));
    }
 
  public:
 
    using Jacobian = Std::detected_or_t<JacobianDefault, JacobianOfImplementation, Implementation>;
    using JacobianInverse = Std::detected_or_t<JacobianInverseDefault, JacobianInverseOfImplementation, Implementation>;
 
    explicit GeometryReference ( const Implementation &impl )
      : impl_( &impl )
    {}
 
    GeometryType type () const { return impl().type(); }
 
    bool affine() const { return impl().affine(); }
 
    int corners () const { return impl().corners(); }
    GlobalCoordinate corner ( int i ) const { return impl().corner( i ); }
    GlobalCoordinate center () const { return impl().center(); }
 
    GlobalCoordinate global ( const LocalCoordinate &local ) const
    {
      return impl().global( local );
    }
 
    LocalCoordinate local ( const GlobalCoordinate &global ) const
    {
      return impl().local( global );
    }
 
    ctype integrationElement ( const LocalCoordinate &local ) const
    {
      return impl().integrationElement( local );
    }
 
    ctype volume () const { return impl().volume(); }
 
    JacobianTransposed jacobianTransposed ( const LocalCoordinate &local ) const
    {
      return impl().jacobianTransposed( local );
    }
 
    JacobianInverseTransposed jacobianInverseTransposed ( const LocalCoordinate &local ) const
    {
      return impl().jacobianInverseTransposed( local );
    }
 
    Jacobian jacobian ( const LocalCoordinate& local ) const
    {
      if constexpr(Std::is_detected_v<JacobianOfImplementation, Implementation>)
        return impl().jacobian(local);
      else
        return deprecatedDefaultJacobian(local);
    }
 
    JacobianInverse jacobianInverse ( const LocalCoordinate &local ) const
    {
      if constexpr(Std::is_detected_v<JacobianInverseOfImplementation, Implementation>)
        return impl().jacobianInverse(local);
      else
        return deprecatedDefaultJacobianInverse(local);
    }
 
    const Implementation &impl () const { return *impl_; }
 
  private:
    const Implementation *impl_;
  };
 
 
  // LocalGeometryReference
  // -----------------------
 
  template< int mydim, int cdim, class Grid >
  class LocalGeometryReference
    : public GeometryReference< typename std::remove_const< Grid >::type::Traits::template Codim< std::remove_const< Grid >::type::dimension - mydim >::LocalGeometryImpl >
  {
    typedef typename std::remove_const< Grid >::type::Traits::template Codim< std::remove_const< Grid >::type::dimension - mydim >::LocalGeometryImpl Implementation;
 
  public:
    LocalGeometryReference ( const Implementation &impl )
      : GeometryReference< Implementation >( impl )
    {}
  };
 
 
 
  // Definitions of GeometryReference
  // --------------------------------
 
  template< class Implementation >
  const int GeometryReference< Implementation >::mydimension;
 
  template< class Implementation >
  const int GeometryReference< Implementation >::coorddimension;
 
} // namespace Dune
 
#endif // #ifndef DUNE_GRID_ALBERTAGRID_GEOMETRYREFERENCE_HH