master/dune-pdelab/a00602_source.html

#ifndef DUNE_PDELAB_LOCALOPERATOR_PERMEABILITY_ADAPTER_HH

#define DUNE_PDELAB_LOCALOPERATOR_PERMEABILITY_ADAPTER_HH


#include <dune/pdelab/common/function.hh>


template<typename T>


class PermeabilityAdapter

  : public Dune::PDELab::GridFunctionBase<Dune::PDELab::GridFunctionTraits<typename T::Traits::GridViewType,

                                       typename T::Traits::RangeFieldType,

                                       1,Dune::FieldVector<typename T::Traits::RangeFieldType,1> >

                      ,PermeabilityAdapter<T> >

{

public:

  typedef Dune::PDELab::GridFunctionTraits<typename T::Traits::GridViewType,

                       typename T::Traits::RangeFieldType,

                       1,Dune::FieldVector<typename T::Traits::RangeFieldType,1> > Traits;


  PermeabilityAdapter (const typename Traits::GridViewType& g_, T& t_)

    : g(g_), t(t_)

  {}


  inline void evaluate (const typename Traits::ElementType& e,

            const typename Traits::DomainType& x,

            typename Traits::RangeType& y) const

  {

    using namespace std;

    y = log(t.A(e,x)[0][0]);

  }


  inline const typename Traits::GridViewType& getGridView () const

  {

    return g;

  }


  inline void setTime(double time_)

  {

    t.setTime(time_);

  }


private:

  const typename Traits::GridViewType& g;

  T& t;

};


template<typename T>


class DiagonalPermeabilityAdapter

  : public Dune::PDELab::GridFunctionBase<Dune::PDELab::GridFunctionTraits<typename T::Traits::GridViewType,

                                       typename T::Traits::RangeFieldType,

                                       T::Traits::dimDomain,Dune::FieldVector<typename T::Traits::RangeFieldType,T::Traits::dimDomain> >

                      ,DiagonalPermeabilityAdapter<T> >

{

public:

  typedef Dune::PDELab::GridFunctionTraits<typename T::Traits::GridViewType,

                       typename T::Traits::RangeFieldType,

                       T::Traits::dimDomain,Dune::FieldVector<typename T::Traits::RangeFieldType,T::Traits::dimDomain> > Traits;


  DiagonalPermeabilityAdapter (const typename Traits::GridViewType& g_, T& t_)

    : g(g_), t(t_)

  {}


  inline void evaluate (const typename Traits::ElementType& e,

            const typename Traits::DomainType& x,

            typename Traits::RangeType& y) const

  {

    using namespace std;

    for (int i=0; i<T::Traits::dimDomain; i++)

      y[i] = log(t.A(e,x)[i][i]);

  }


  inline const typename Traits::GridViewType& getGridView () const

  {

    return g;

  }


  inline void setTime(double time_)

  {

    t.setTime(time_);

  }


private:

  const typename Traits::GridViewType& g;

  T& t;

};


#endif // DUNE_PDELAB_LOCALOPERATOR_PERMEABILITY_ADAPTER_HH

function.hh

std
STL namespace.

Dune::FieldVector

Dune::PDELab::FunctionTraits< GV::Grid::ctype, GV::dimension, Dune::FieldVector< GV::Grid::ctype, GV::dimension >, RF, m, R >::RangeType
R RangeType
range type
Definition function.hh:62

Dune::PDELab::PowerCompositeGridFunctionTraits::ElementType
GV::Traits::template Codim< 0 >::Entity ElementType
codim 0 entity
Definition function.hh:119

Dune::PDELab::PowerCompositeGridFunctionTraits::GridViewType
GV GridViewType
The type of the grid view the function lives on.
Definition function.hh:116

Dune::PDELab::GridFunctionTraits
traits class holding the function signature, same as in local function
Definition function.hh:183

Dune::PDELab::GridFunctionBase
leaf of a function tree
Definition function.hh:302

PermeabilityAdapter
Definition permeability_adapter.hh:16

PermeabilityAdapter::setTime
void setTime(double time_)
Definition permeability_adapter.hh:41

PermeabilityAdapter::evaluate
void evaluate(const typename Traits::ElementType &e, const typename Traits::DomainType &x, typename Traits::RangeType &y) const
Definition permeability_adapter.hh:28

PermeabilityAdapter::PermeabilityAdapter
PermeabilityAdapter(const typename Traits::GridViewType &g_, T &t_)
constructor
Definition permeability_adapter.hh:23

PermeabilityAdapter::getGridView
const Traits::GridViewType & getGridView() const
Definition permeability_adapter.hh:36

PermeabilityAdapter::Traits
Dune::PDELab::GridFunctionTraits< typename T::Traits::GridViewType, typename T::Traits::RangeFieldType, 1, Dune::FieldVector< typename T::Traits::RangeFieldType, 1 > > Traits
Definition permeability_adapter.hh:20

DiagonalPermeabilityAdapter
Definition permeability_adapter.hh:61

DiagonalPermeabilityAdapter::setTime
void setTime(double time_)
Definition permeability_adapter.hh:87

DiagonalPermeabilityAdapter::DiagonalPermeabilityAdapter
DiagonalPermeabilityAdapter(const typename Traits::GridViewType &g_, T &t_)
constructor
Definition permeability_adapter.hh:68

DiagonalPermeabilityAdapter::Traits
Dune::PDELab::GridFunctionTraits< typename T::Traits::GridViewType, typename T::Traits::RangeFieldType, T::Traits::dimDomain, Dune::FieldVector< typename T::Traits::RangeFieldType, T::Traits::dimDomain > > Traits
Definition permeability_adapter.hh:65

DiagonalPermeabilityAdapter::getGridView
const Traits::GridViewType & getGridView() const
Definition permeability_adapter.hh:82

DiagonalPermeabilityAdapter::evaluate
void evaluate(const typename Traits::ElementType &e, const typename Traits::DomainType &x, typename Traits::RangeType &y) const
Definition permeability_adapter.hh:73