3 #ifndef DUNE_ALBERTA_GEOMETRYCACHE_HH
4 #define DUNE_ALBERTA_GEOMETRYCACHE_HH
23 static const unsigned int flagIntegrationElement = (1 << 0);
24 static const unsigned int flagJacobianTransposed = (1 << 1);
25 static const unsigned int flagJacobianInverseTransposed = (1 << 2);
37 if( (flags_ & flagIntegrationElement) == 0 )
40 assert( integrationElement_ > 1e-14 );
41 flags_ |= flagIntegrationElement;
43 return integrationElement_;
48 if( (flags_ & flagJacobianTransposed) == 0 )
52 for(
int i = 0; i < dim; ++i )
56 jacobianTransposed_[ i ][ j ] = y[ j ] - x[ j ];
58 flags_ |= flagJacobianTransposed;
60 return jacobianTransposed_;
66 if( (flags_ & flagJacobianInverseTransposed) == 0 )
69 assert( integrationElement_ > 1e-14 );
70 flags_ |= flagIntegrationElement | flagJacobianInverseTransposed;
72 return jacobianInverseTransposed_;
77 Real integrationElement_;
78 FieldMatrix< Real, dim, dimWorld > jacobianTransposed_;
79 FieldMatrix< Real, dimWorld, dim > jacobianInverseTransposed_;
94 : geometryCache_( geometryCache ),
100 return geometryCache_.integrationElement( elInfo_ );
105 return geometryCache_.jacobianTransposed( elInfo_ );
110 return geometryCache_.jacobianInverseTransposed( elInfo_ );
115 const ALBERTA EL_INFO &elInfo_;
122 #endif // #if HAVE_ALBERTA
124 #endif // #ifndef DUNE_ALBERTA_GEOMETRYCACHE_HH