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

[Markus Blatt]

Hi,

On Mon, Sep 24, 2012 at 06:17:06PM +0200, Christian Engwer wrote:
> Hi Markus,
> 
> I guess, I'm just too dumb, but I don't understand the problem?!
> 

As pointed out, this is a hack. I am under the impression that the
problem is in the Interface. Probably it is not yet finished and still
in the transition from Navier-Stokes with fixed rho to Navier-Stokes
with variable rho (patch r1954). In the NavierStokesParameterTraits it reads: 

typedef Dune::FieldVector<DomainField,dimDomain> Domain;

Unfortunately, I have only a PDELab:LocalVector in jacobian_volume of
NavierStokesMass and  I do not sure which position to take in the
entity for varying rho. I would leave to the original other of patch
r1954 to decide this. (In the meanwhile for pdelab-howto at least the
old functionality with constant rho works.)

Markus

> > ***** 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
> 
> _______________________________________________
> dune-pdelab mailing list
> dune-pdelab at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune-pdelab
> 

-- 
Do you need more support with DUNE or HPC in general? 

Dr. Markus Blatt - HPC-Simulation-Software & Services http://www.dr-blatt.de
Hans-Bunte-Str. 8-10, 69123 Heidelberg, Germany
Tel.: +49 (0) 160 97590858  Fax: +49 (0)322 1108991658