columnobject.hh

#ifndef DUNE_FEM_COLUMNOBJECT_HH
#define DUNE_FEM_COLUMNOBJECT_HH
 
namespace Dune
{
 
  namespace Fem
  {
 
    template< class LinearOperator >
    struct ColumnObject
    {
      typedef typename LinearOperator::ColumnEntityType ColumnEntityType;
      typedef typename LinearOperator::RowEntityType RowEntityType;
 
      typedef typename LinearOperator::DomainSpaceType DomainSpaceType;
      typedef typename LinearOperator::RangeSpaceType RangeSpaceType;
 
      typedef typename LinearOperator::LocalMatrixType LocalMatrixType;
 
      ColumnObject( const LinearOperator &linOp, const ColumnEntityType &colEntity )
      :
        linOp_( linOp ),
        colEntity_( colEntity )
      {}
 
      inline LocalMatrixType localMatrix( const RowEntityType &rowEntity ) const
      {
        return linOp_.localMatrix( rowEntity, colEntity_ );
      }
 
      const DomainSpaceType& domainSpace() const { return linOp_.domainSpace(); }
 
      const RangeSpaceType& rangeSpace() const { return linOp_.rangeSpace(); }
 
    private:
      const LinearOperator &linOp_;
      const ColumnEntityType &colEntity_;
    };
 
  } // namespace Fem
} // namespace Dune
 
#endif //#ifndef DUNE_FEM_COLUMNOBJECT_HH