/usr/include/dune/grid/albertagrid/macroelement.hh is in libdune-grid-dev 2.2.1-2.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 | #ifndef DUNE_ALBERTA_MACROELEMENT_HH
#define DUNE_ALBERTA_MACROELEMENT_HH
#include <dune/grid/albertagrid/misc.hh>
#if HAVE_ALBERTA
namespace Dune
{
namespace Alberta
{
// MacroElement
// ------------
template< int dim >
struct MacroElement
: public ALBERTA MACRO_EL
{
const GlobalVector &coordinate ( const int vertex ) const;
int boundaryId ( const int face ) const;
bool isBoundary ( const int face ) const;
const MacroElement< dim > *neighbor ( const int face ) const;
};
#if DUNE_ALBERTA_VERSION >= 0x300
template< int dim >
inline const GlobalVector &MacroElement< dim >::coordinate ( const int vertex ) const
{
assert( (vertex >= 0) && (vertex < N_VERTICES_MAX) );
return *coord[ vertex ];
}
#endif // #if DUNE_ALBERTA_VERSION >= 0x300
#if DUNE_ALBERTA_VERSION < 0x300
template< int dim >
inline const GlobalVector &MacroElement< dim >::coordinate ( const int vertex ) const
{
assert( (vertex >= 0) && (vertex < N_VERTICES_MAX) );
return *((const GlobalVector *)coord[ vertex ]);
}
#endif // #if DUNE_ALBERTA_VERSION < 0x300
template< int dim >
inline bool MacroElement< dim >::isBoundary ( const int face ) const
{
return (boundaryId( face ) != InteriorBoundary);
}
#if DUNE_ALBERTA_VERSION >= 0x300
template< int dim >
inline int MacroElement< dim >::boundaryId ( const int face ) const
{
return wall_bound[ face ];
}
#endif // #if DUNE_ALBERTA_VERSION >= 0x300
#if DUNE_ALBERTA_VERSION < 0x300
template< int dim >
inline int MacroElement< dim >::boundaryId ( const int face ) const
{
switch( dim )
{
case 1:
assert( (face >= 0) && (face < N_VERTICES_MAX) );
return vertex_bound[ face ];
case 2:
assert( (face >= 0) && (face < N_EDGES_MAX) );
return edge_bound[ face ];
case 3:
assert( (face >= 0) && (face < N_FACES_MAX) );
return face_bound[ face ];
}
}
#endif // #if DUNE_ALBERTA_VERSION < 0x300
template< int dim >
const MacroElement< dim > *MacroElement< dim >::neighbor ( const int face ) const
{
assert( (face >= 0) && (face < N_NEIGH_MAX) );
return static_cast< const MacroElement * >( neigh[ face ] );
}
}
}
#endif // #if HAVE_ALBERTA
#endif // #ifndef DUNE_ALBERTA_MACROELEMENT_HH
|