identitygridgeometry.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_IDENTITYGRIDGEOMETRY_HH
#define DUNE_IDENTITYGRIDGEOMETRY_HH

#include <dune/common/fmatrix.hh>
#include <dune/common/typetraits.hh>
#include <dune/grid/common/geometry.hh>

namespace Dune {

  template<int mydim, int coorddim, class GridImp>
  class IdentityGridGeometry :
    public GeometryDefaultImplementation <mydim, coorddim, GridImp, IdentityGridGeometry>
  {
  private:

    typedef typename GridImp::ctype ctype;


  public:

    // The codimension of this entitypointer wrt the host grid
    constexpr static int CodimInHostGrid = GridImp::HostGridType::dimension - mydim;
    constexpr static int DimensionWorld = GridImp::HostGridType::dimensionworld;

    // select appropriate hostgrid geometry via typeswitch
    typedef typename GridImp::HostGridType::Traits::template Codim<CodimInHostGrid>::Geometry HostGridGeometryType;
    typedef typename GridImp::HostGridType::Traits::template Codim<CodimInHostGrid>::Geometry HostGridLocalGeometryType;

    typedef typename std::conditional<coorddim==DimensionWorld, HostGridGeometryType, HostGridLocalGeometryType>::type HostGridGeometry;

    typedef typename HostGridGeometryType::JacobianInverseTransposed JacobianInverseTransposed;
    typedef typename HostGridGeometryType::JacobianTransposed JacobianTransposed;


    IdentityGridGeometry(const HostGridGeometry& hostGeometry)
      : hostGeometry_(hostGeometry)
    {}


    GeometryType type () const {
      return hostGeometry_.type();
    }

    // return whether we have an affine mapping
    bool affine() const {
      return hostGeometry_.affine();
    }

    int corners () const {
      return hostGeometry_.corners();
    }


    const FieldVector<ctype, coorddim> corner (int i) const {
      return hostGeometry_.corner(i);
    }


    FieldVector<ctype, coorddim> global (const FieldVector<ctype, mydim>& local) const {
      return hostGeometry_.global(local);
    }

    JacobianTransposed
    jacobianTransposed ( const FieldVector<ctype, mydim>& local ) const {
      return hostGeometry_.jacobianTransposed(local);
    }

    FieldVector<ctype, mydim> local (const FieldVector<ctype, coorddim>& global) const {
      return hostGeometry_.local(global);
    }


    bool checkInside(const FieldVector<ctype, mydim> &local) const {
      return hostGeometry_.checkInside(local);
    }


    ctype integrationElement (const FieldVector<ctype, mydim>& local) const {
      return hostGeometry_.integrationElement(local);
    }


    JacobianInverseTransposed jacobianInverseTransposed (const FieldVector<ctype, mydim>& local) const {
      return hostGeometry_.jacobianInverseTransposed(local);
    }


    HostGridGeometry hostGeometry_;

  };

}  // namespace Dune

#endif