[dune-pdelab] Header stokesdg_vecfem.hh

[Christian Engwer]

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

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