[dune-pdelab] Header stokesdg_vecfem.hh

[Christian Engwer]

On Fri, Aug 30, 2013 at 02:27:04PM +0200, Christian Engwer wrote:
> On Fri, Aug 30, 2013 at 11:16:09AM +0200, Christoph Grüninger wrote:
> > Hi PDELab,
> > currently only one header fails headercheck in PDELab for me.
> > Attached patch fixes the remaining one but I am not sure whether
> > deleting the PressureDirichlet part is a good idea.
> 
> It is not... 
> 
> I'll look into it

Thanks for pointing it out! It should work now :-)

> 
> Christian
> 
> > Bye
> > Christoph
> > 
> > -- 
> > Digital information lasts forever – or five years, whichever comes first.
> > -- Jeff Rothenberg, 1997
> > *********************************************
> > CMWR 2014: 10th - 13th June 2014 in Stuttgart
> >          Please visit www.cmwr14.de
> > *********************************************
> 
> > From e8805cbf16c743e8a86943c532864b7d4ce1a347 Mon Sep 17 00:00:00 2001
> > From: =?UTF-8?q?Christoph=20Gr=C3=BCninger?= <gruenich at dune-project.org>
> > Date: Fri, 30 Aug 2013 11:11:36 +0200
> > Subject: [PATCH] [stokesdg_vecfem.hh] Fix headercheck.
> > 
> > ---
> >  dune/pdelab/localoperator/stokesdg_vecfem.hh | 18 ++++--------------
> >  1 file changed, 4 insertions(+), 14 deletions(-)
> > 
> > diff --git a/dune/pdelab/localoperator/stokesdg_vecfem.hh b/dune/pdelab/localoperator/stokesdg_vecfem.hh
> > index ffd1a8a..537c000 100644
> > --- a/dune/pdelab/localoperator/stokesdg_vecfem.hh
> > +++ b/dune/pdelab/localoperator/stokesdg_vecfem.hh
> > @@ -546,16 +546,6 @@ namespace Dune {
> >  
> >            } // DirichletVelocity
> >  
> > -          if (bctype == BC::PressureDirichlet){
> > -            typename P::Traits::RangeType p0;
> > -            p.evaluateGlobal(global,p0);
> > -
> > -            for (size_t i=0; i<lfsv_v_s.size(); i++){
> > -              const RF val = p0 * (normal*phi_v_s[i]) * weight;
> > -              r_s.accumulate(lfsv_v_s, i, val);
> > -            }
> > -          } // PressureDirichlet
> > -
> >          } // it - quadrature
> >  
> >        }
> > @@ -576,10 +566,10 @@ namespace Dune {
> >        template<class M, class R>
> >        static void add_compute_flux(const M & du, const R & n, R & result)
> >        {
> > -        const int N = du.N();
> > -        const int M = du.M();
> > -        for(int r=0; r<N; ++r)
> > -          for(int c=0; c<M; ++c)
> > +        const int duN = du.N();
> > +        const int duM = du.M();
> > +        for(int r=0; r<duN; ++r)
> > +          for(int c=0; c<duM; ++c)
> >              result[r] += du[r][c] * n[c];
> >        }
> >  
> > -- 
> > 1.8.1.4
> > 
> 
> > _______________________________________________
> > dune-pdelab mailing list
> > dune-pdelab at dune-project.org
> > http://lists.dune-project.org/mailman/listinfo/dune-pdelab
> 
> 
> -- 
> Prof. Dr. Christian Engwer 
> Institut für Numerische und Angewandte Mathematik
> Fachbereich Mathematik und Informatik der Universität Münster
> Einsteinstrasse 62
> 48149 Münster
> 
> E-Mail	christian.engwer at uni-muenster.de
> Telefon	+49 251 83-35067
> FAX		+49 251 83-32729
> 
> _______________________________________________
> dune-pdelab mailing list
> dune-pdelab at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune-pdelab

-- 
Prof. Dr. Christian Engwer 
Institut für Numerische und Angewandte Mathematik
Fachbereich Mathematik und Informatik der Universität Münster
Einsteinstrasse 62
48149 Münster

E-Mail	christian.engwer at uni-muenster.de
Telefon	+49 251 83-35067
FAX		+49 251 83-32729