geometrytraits.hh

// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
#ifndef DUNE_GEOMETRY_GENERICGEOMETRY_GEOMETRYTRAITS_HH
#define DUNE_GEOMETRY_GENERICGEOMETRY_GEOMETRYTRAITS_HH
 
#include "../type.hh"
#include "matrixhelper.hh"
#include "cornermapping.hh"
 
namespace Dune
{
  namespace GenericGeometry
  {
 
    // DuneCoordTraits
    // ---------------
 
    template< class ct >
    struct DuneCoordTraits
    {
      typedef ct ctype;
 
      template< int dim >
      struct Vector
      {
        typedef FieldVector< ctype, dim > type;
      };
 
      template< int rows, int cols >
      struct Matrix
      {
        typedef FieldMatrix< ctype, rows, cols > type;
      };
 
      // This limit is, e.g., used in the termination criterion of the Newton
      // scheme within the generic implementation of the method local
      static const ctype epsilon ()
      {
        return 1e-6;
      }
    };
 
 
 
    // MappingTraits
    // -------------
    template< class CT, unsigned int dim, unsigned int dimW >
    struct MappingTraits
    {
      typedef CT CoordTraits;
 
      static const unsigned int dimension = dim;
      static const unsigned int dimWorld = dimW;
 
      typedef typename CoordTraits :: ctype FieldType;
      typedef typename CoordTraits :: template Vector< dimension > :: type LocalCoordinate;
      typedef typename CoordTraits :: template Vector< dimWorld > :: type GlobalCoordinate;
 
      typedef typename CoordTraits :: template Matrix< dimWorld, dimension > :: type
      JacobianType;
      typedef typename CoordTraits :: template Matrix< dimension, dimWorld > :: type
      JacobianTransposedType;
 
      typedef GenericGeometry :: MatrixHelper< CoordTraits > MatrixHelper;
    };
 
 
 
    enum EvaluationType
    {
      ComputeOnDemand,
      PreCompute
    };
 
 
 
    // DefaultGeometryTraits
    // ---------------------
 
    template< class ctype, int dimG, int dimW, bool alwaysAffine = false >
    struct DefaultGeometryTraits
    {
      typedef DuneCoordTraits< ctype > CoordTraits;
 
      static const int dimWorld = dimW;
 
      template< int dim >
      struct hasSingleGeometryType
      {
        static const bool v = false;
        static const unsigned int topologyId = ~0u;
      };
 
      template< class Topology >
      struct Mapping
      {
        typedef CoordStorage< CoordTraits, Topology, dimWorld > CornerStorage;
        typedef CornerMapping< CoordTraits, Topology, dimWorld, CornerStorage, alwaysAffine > type;
      };
 
      struct Caching
      {
        static const EvaluationType evaluateJacobianTransposed = ComputeOnDemand;
        static const EvaluationType evaluateJacobianInverseTransposed = ComputeOnDemand;
        static const EvaluationType evaluateIntegrationElement = ComputeOnDemand;
      };
 
      struct UserData {};
    };
 
 
 
    template< class Grid >
    struct GlobalGeometryTraits;
 
    template< class Grid >
    struct GlobalGeometryTraits< const Grid >
      : public GlobalGeometryTraits< Grid >
    {};
 
 
 
    template< class Grid >
    struct LocalGeometryTraits;
 
    template< class Grid >
    struct LocalGeometryTraits< const Grid >
      : public LocalGeometryTraits< Grid >
    {};
  }
 
}
 
#endif // #ifndef DUNE_GEOMETRY_GENERICGEOMETRY_GEOMETRYTRAITS_HH