1 #ifndef DUNE_ALBERTA_TRANSFORMATION_HH
2 #define DUNE_ALBERTA_TRANSFORMATION_HH
4 #include <dune/common/fvector.hh>
26 : matrix_( (trafo != NULL ? trafo->M :
GlobalSpace::identityMatrix()) ),
27 shift_( (trafo != NULL ? trafo->t :
GlobalSpace::nullVector()) )
44 y[ i ] += row[ j ] * x[ j ];
56 const ctype v = y[ i ] - shift_[ i ];
58 x[ j ] += row[ j ] * v;
70 #endif // #if HAVE_ALBERTA
72 #endif // #ifndef DUNE_ALBERTA_TRANSFORMATION_HH