[dune-pdelab] [dune-pdelab-commit] dune-pdelab r2233 - trunk/dune/pdelab/localoperator

[Christian Engwer]

Hi Markus,

I guess, I'm just too dumb, but I don't understand the problem?!

> ***** HACK ALARM *******
> I had to change the interface of
> TaylorHoodNavierStokesDefaultParameters::rho(eg,x). Namely the type of
> x is not fixed anymore, but can be arbitrary. This was needed as there
> is no way of finding out the global coordinate within
> NavierStokesMass::jacobian_volume (or at least it is unclear which
> value to take).

why would you need the type for a global coordinate? You only pass
local coordinates!

>          /** \brief Density evaluated on a codim 0 geometry. */
> -        template<typename EG>
> -        RF rho(const EG & eg, const typename Traits::Domain & x) const
> +        template<typename EG, typename Domain>
> +        RF rho(const EG & eg, const Domain & x) const
>          { return rho_; }

Christian