]> https://gitweb.dealii.org/ - dealii.git/commitdiff
fix problem with zero quedrature weights in FEPolyTensor
authorGuido Kanschat <dr.guido.kanschat@gmail.com>
Thu, 20 Mar 2008 21:27:04 +0000 (21:27 +0000)
committerGuido Kanschat <dr.guido.kanschat@gmail.com>
Thu, 20 Mar 2008 21:27:04 +0000 (21:27 +0000)
git-svn-id: https://svn.dealii.org/trunk@15922 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/source/fe/fe.cc

index 2b7e92b6db97a3c209dce97c1a7c4672507b59c2..b16f24bddf76665346da8bc2176076d03028707c 100644 (file)
@@ -73,14 +73,14 @@ InternalDataBase::initialize_2nd (const FiniteElement<dim> *element,
                                       // gradients, not more
       for (unsigned int i=0; i<diff_points.size(); ++i)
        diff_points[i] = quadrature.point(i) + shift;
-      const Quadrature<dim> plus_quad (diff_points);
+      const Quadrature<dim> plus_quad (diff_points, quadrature.get_weights());
       differences[d] = new FEValues<dim> (mapping, *element,
                                          plus_quad, update_gradients);
 
                                       // now same in minus-direction
       for (unsigned int i=0; i<diff_points.size(); ++i)
        diff_points[i] = quadrature.point(i) - shift;
-      const Quadrature<dim> minus_quad (diff_points);
+      const Quadrature<dim> minus_quad (diff_points, quadrature.get_weights());
       differences[d+dim] = new FEValues<dim> (mapping, *element,
                                              minus_quad, update_gradients); 
     }

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.