]> https://gitweb.dealii.org/ - dealii.git/commitdiff
do not use deprecated access to Quadrature::n_quadrature_points anymore
authorGuido Kanschat <dr.guido.kanschat@gmail.com>
Thu, 11 Nov 2010 20:11:47 +0000 (20:11 +0000)
committerGuido Kanschat <dr.guido.kanschat@gmail.com>
Thu, 11 Nov 2010 20:11:47 +0000 (20:11 +0000)
git-svn-id: https://svn.dealii.org/trunk@22693 0785d39b-7218-0410-832d-ea1e28bc413d

121 files changed:
tests/Makefile
tests/aniso/mesh_3d_21.cc
tests/base/function_derivative.cc
tests/base/gauss_lobatto.cc
tests/base/patches.h
tests/base/polynomial1d.cc
tests/base/polynomial_lagrange.cc
tests/base/polynomials_bdm_01.cc
tests/base/polynomials_rt.cc
tests/base/qprojector.cc
tests/base/quadrature_test.cc
tests/bits/fe_tools_cpfqpm_01.cc
tests/bits/fe_tools_cpfqpm_02.cc
tests/bits/fe_tools_cpfqpm_03.cc
tests/bits/fe_tools_cpfqpm_04.cc
tests/bits/joa_1.cc
tests/bits/mesh_3d_14.cc
tests/bits/mesh_3d_16.cc
tests/bits/mesh_3d_20.cc
tests/bits/mesh_3d_7.cc
tests/bits/normals_1.cc
tests/bits/normals_2.cc
tests/bits/normals_3.cc
tests/bits/normals_4.cc
tests/bits/q_point_sum_1.cc
tests/bits/q_point_sum_2.cc
tests/bits/q_point_sum_3.cc
tests/bits/q_point_sum_4.cc
tests/bits/q_points.cc
tests/bits/rt_2.cc
tests/bits/step-13.cc
tests/bits/step-14.cc
tests/bits/step-15.cc
tests/bits/step-16.cc
tests/bits/step-3.cc
tests/bits/step-4.cc
tests/bits/step-5.cc
tests/bits/step-6.cc
tests/bits/step-7.cc
tests/bits/step-8.cc
tests/bits/volume_1.cc
tests/bits/volume_2.cc
tests/bits/volume_3.cc
tests/bits/volume_4.cc
tests/codim_one/bem.cc
tests/codim_one/gradients.cc
tests/codim_one/mapping_q1.cc
tests/codim_one/mapping_q1_eulerian.cc
tests/deal.II/block_matrices.cc
tests/deal.II/derivatives.cc
tests/deal.II/fe_values_view_05.cc
tests/deal.II/fe_values_view_06.cc
tests/deal.II/fe_values_view_07.cc
tests/deal.II/fe_values_view_08.cc
tests/deal.II/fe_values_view_09.cc
tests/deal.II/fe_values_view_10.cc
tests/deal.II/fe_values_view_10_single_01.cc
tests/deal.II/fe_values_view_10_single_02.cc
tests/deal.II/fe_values_view_10_single_03.cc
tests/deal.II/fe_values_view_10_single_04.cc
tests/deal.II/fe_values_view_11.cc
tests/deal.II/fe_values_view_12.cc
tests/deal.II/fe_values_view_13.cc
tests/deal.II/fe_values_view_14.cc
tests/deal.II/fe_values_view_15.cc
tests/deal.II/fe_values_view_16.cc
tests/deal.II/fe_values_view_17.cc
tests/deal.II/fe_values_view_18.cc
tests/deal.II/fe_values_view_19.cc
tests/deal.II/fe_values_view_20.cc
tests/deal.II/fe_values_view_21.cc
tests/deal.II/filtered_matrix.cc
tests/fail/abf_approximation_01.cc
tests/fail/circular_01.cc
tests/fail/hp-step-14.cc
tests/fail/hp-step-15.cc
tests/fail/rt_distorted_01.cc
tests/fail/rt_distorted_02.cc
tests/fe/abf_01.cc
tests/fe/check_derivatives.cc
tests/fe/deformed_projection.h
tests/fe/derivatives_face.cc
tests/fe/function.cc
tests/fe/interpolate_common.h
tests/fe/mapping.cc
tests/fe/mapping_q1_eulerian.cc
tests/fe/mapping_q_eulerian.cc
tests/fe/nedelec_2.cc
tests/fe/nedelec_3.cc
tests/fe/non_primitive_1.cc
tests/fe/non_primitive_2.cc
tests/fe/rt_10.cc
tests/fe/rt_11.cc
tests/fe/rt_12.cc
tests/fe/rt_13.cc
tests/fe/rt_14.cc
tests/fe/rt_15.cc
tests/fe/rt_7.cc
tests/fe/rt_8.cc
tests/fe/rt_approximation_01.cc
tests/fe/rt_normal_02.cc
tests/fe/rtdiff.cc
tests/fe/rtn_2.cc
tests/fe/shapes.cc
tests/hp/crash_17.cc
tests/hp/crash_17_compressed_set_sparsity.cc
tests/hp/crash_18.cc
tests/hp/crash_18_compressed_set_sparsity.cc
tests/hp/step-13.cc
tests/hp/step-3.cc
tests/hp/step-3a.cc
tests/hp/step-3b.cc
tests/hp/step-3c.cc
tests/hp/step-4.cc
tests/hp/step-5.cc
tests/hp/step-6.cc
tests/hp/step-7.cc
tests/hp/step-8.cc
tests/integrators/cochain_01.cc
tests/integrators/cochain_01/cmp/generic
tests/umfpack/umfpack_04.cc

index dcea3d384c9ee4007c9e2d02c3f2c173c363b042..b1ad499c51e9a1c459f84fc98ec35f212f8539a1 100644 (file)
@@ -11,12 +11,13 @@ nofail-dirs = a-framework       \
               lac               \
               fe                \
               deal.II           \
+              integrators       \
               multigrid         \
               bits              \
               hp                \
-              aniso     \
-              codim_one \
-              serialization \
+              aniso             \
+              codim_one         \
+              serialization     \
              $(if $(DEAL_II_USE_MPI:no=),      $(if $(USE_CONTRIB_P4EST:no=),mpi,),) \
              $(if $(USE_CONTRIB_P4EST:no=),    distributed_grids,) \
              $(if $(USE_CONTRIB_PETSC:no=),    petsc,) \
index 601a257c17655277aadb7e8aeffcea13acdc313f..22c732c78da6e134a70191ce3914f036a9889495 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -71,7 +71,7 @@ void check_this (Triangulation<3> &tria)
             fe_face_values1.reinit (neighbor_child, neighbor_neighbor);
             fe_face_values2.reinit (cell, face_no, subface_no);
 
-            for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<quadrature.size(); ++q)
               {
                 Assert ((fe_face_values1.quadrature_point(q)-
                          fe_face_values2.quadrature_point(q)).square()
index 28171943cde01298aa4a4cc1f74b738be5e7a0ff..89e610dd8ad27c84a0a8dcf06f81368b17a7d6aa 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2007 by the deal.II authors
+//    Copyright (C) 2007, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -35,8 +35,8 @@ check_derivative_order(const std::vector<Tensor<1,dim> >& gradients,
                       const unsigned int direction,
                       const double order)
 {
-  std::vector<double> derivatives(quadrature.n_quadrature_points);
-  std::vector<double> differences(quadrature.n_quadrature_points);
+  std::vector<double> derivatives(quadrature.size());
+  std::vector<double> differences(quadrature.size());
 
                                   // Compute derivatives with one
                                   // step size and store errors
@@ -73,8 +73,8 @@ check_hessian_order(const std::vector<double>& values,
                    const unsigned int direction,
                    const double order)
 {
-  std::vector<Tensor<1,dim> > derivatives(quadrature.n_quadrature_points);
-  std::vector<Tensor<1,dim> > differences(quadrature.n_quadrature_points);
+  std::vector<Tensor<1,dim> > derivatives(quadrature.size());
+  std::vector<Tensor<1,dim> > differences(quadrature.size());
 
   const double h = (order < 3) ? 1.e-2 : 1.e-1 ;
                                   // Compute derivatives with one
@@ -122,8 +122,8 @@ check_sine(unsigned int nquad)
   
   Functions::FourierSineFunction<dim> f(wave_vector);
   
-  std::vector<double> values(quadrature.n_quadrature_points);
-  std::vector<Tensor<1,dim> > gradients(quadrature.n_quadrature_points);
+  std::vector<double> values(quadrature.size());
+  std::vector<Tensor<1,dim> > gradients(quadrature.size());
   
   f.value_list(quadrature.get_points(), values);
   f.gradient_list(quadrature.get_points(), gradients);
index 91c3d24d0b00069afc207cb0b13aba08c7e90fd2..ff64026acab2e5050b26d1cf99fe09a8a970b221 100644 (file)
@@ -72,7 +72,7 @@ check_cells (std::vector<Quadrature<dim>*>& quadratures)
          quadrature_int=0;
                                           // Check the polynomial x^i*y^i
 
-         for (unsigned int x=0; x<quadrature.n_quadrature_points; ++x)
+         for (unsigned int x=0; x<quadrature.size(); ++x)
            {
              double f=1.;
              switch (dim)
@@ -101,13 +101,13 @@ check_cells (std::vector<Quadrature<dim>*>& quadratures)
                                           // check the ordering of
                                           // the quadrature points
          bool in_order=true;
-         for (unsigned int x=1; x<quadrature.n_quadrature_points; ++x)
+         for (unsigned int x=1; x<quadrature.size(); ++x)
            {
              if (points[x](0)<=points[x-1](0))
                in_order=false;
            }
          if (!in_order)
-           for (unsigned int x=0; x<quadrature.n_quadrature_points; ++x)
+           for (unsigned int x=0; x<quadrature.size(); ++x)
              deallog << points[x] << std::endl;
        }
     }
@@ -146,7 +146,7 @@ check_faces (const std::vector<Quadrature<dim-1>*>& quadratures, const bool sub)
                                           // Check the polynomial
                                           // x^i*y^i*z^i
 
-         for (unsigned int x=0; x<quadrature.n_quadrature_points; ++x)
+         for (unsigned int x=0; x<quadrature.size(); ++x)
            {
              long double f=1.;
              switch (dim)
@@ -200,21 +200,21 @@ int main()
       deallog << "QGaussLobatto(" << n << ")" << std::endl;
 
       QGaussLobatto<1> q(n);
-      for (unsigned int i=0; i<q.n_quadrature_points; ++i)
+      for (unsigned int i=0; i<q.size(); ++i)
        deallog << q.point(i) << ' ' << q.weight(i) << std::endl;
 
                                       // the points must be
                                       // symmetrically located around
                                       // 0.5
       double p = 0;
-      for (unsigned int i=0; i<q.n_quadrature_points; ++i)
+      for (unsigned int i=0; i<q.size(); ++i)
        p += (q.point(i)[0] - 0.5);
       Assert (std::fabs(p) < 1e-12, ExcInternalError());
 
                                       // the sum of weights must be
                                       // one
       double w = 0;
-      for (unsigned int i=0; i<q.n_quadrature_points; ++i)
+      for (unsigned int i=0; i<q.size(); ++i)
        w += q.weight(i);
       Assert (std::fabs(w-1) < 1e-12, ExcInternalError());
     }
index eb3e838e899ddec303fbe3dcb7b2ddc72a25e044..779f89c287331981997a2e89dee6d06d009bcc39 100644 (file)
@@ -99,7 +99,7 @@ create_continuous_patches(
        {
          DataOutBase::Patch<dim, dim> patch;
          patch.n_subdivisions = n_sub;
-         for (unsigned int k=0;k<trapez.n_quadrature_points;++k)
+         for (unsigned int k=0;k<trapez.size();++k)
            {
              Point<dim> p = trapez.point(k);
              if (dim>=1) p(0) += i1;
index 9dc57908ef6633caa7c902100ef391d469ee9180..3b31efe5d622f6ea85638608f9a6d93f1537a892 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2008 by the deal.II authors
+//    Copyright (C) 2000, 2001, 2002, 2003, 2004, 2005, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -36,7 +36,7 @@ double scalar_product (const Polynomial<double>& p1,
   QGauss<1> gauss(degree);
 
   double sum = 0.;
-  for (unsigned int i=0;i<gauss.n_quadrature_points;++i)
+  for (unsigned int i=0;i<gauss.size();++i)
     {
       double x = gauss.point(i)(0);
       double P1 = p1.value(x);
index 718bce6bd8948d23b05d7fa9d2e978eaae902a7f..0a0c7b0daf15a3bfbb55fb24c933eb6c4c2bb9fa 100644 (file)
@@ -55,7 +55,7 @@ void check_interpolation (const std::vector<Polynomial<double> >& p,
 void
 check_poly (const Quadrature<1>& q)
 {
-  deallog << "Points: " << q.n_quadrature_points << std::endl;
+  deallog << "Points: " << q.size() << std::endl;
   std::vector<Polynomial<double> > p = generate_complete_Lagrange_basis(q.get_points());
   check_interpolation(p, q.get_points());
 }
index 8e2473f24a1c0cba36a21493311ba1e9320cb245..87768ff35c5ad0b30b1ad73dde5d510ceea67b61 100644 (file)
@@ -36,7 +36,7 @@ void plot(const PolynomialsBDM<dim>& poly)
   std::vector<Tensor<3,dim> > grads2;
 
   
-  for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
+  for (unsigned int k=0;k<quadrature.size();++k)
     {
       if (k%(poly.degree()+4) == 0)
        deallog << "BDM" << poly.degree() << '<' << dim << '>' << std::endl;
index b0923d10f48ab0de460d23941379cdb802d8b19a..cdde066afed6f9d91e12fd2dd5aeb4bced2eaf53 100644 (file)
@@ -35,7 +35,7 @@ void plot(const PolynomialsRaviartThomas<dim>& poly)
   std::vector<Tensor<3,dim> > grads2;
 
   
-  for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
+  for (unsigned int k=0;k<quadrature.size();++k)
     {
       if (k%(poly.degree()+4) == 0)
        deallog << "RT" << poly.degree() << '<' << dim << '>' << std::endl;
index 47e24c8a2506b3f4ec12a991fa2d06bc0eec22d3..46dbf00be1948cbc3f45fb92fafaf216246d9123 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2008, 2009 by the deal.II authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2008, 2009, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -43,7 +43,7 @@ void check_line(Quadrature<1>& quadrature)
   Quadrature<dim> q = QProjector<dim>::project_to_line(quadrature, p1, p2);
   double s = 0.;
   
-  for (unsigned int k=0;k<q.n_quadrature_points;++k)
+  for (unsigned int k=0;k<q.size();++k)
     {
       deallog << k << '\t' << q.point(k) << std::endl;
       s += q.weight(k);
@@ -55,7 +55,7 @@ template<int dim>
 void check_face(Quadrature<1>& q1)
 {
   deallog << "Checking dim " << dim
-         << " 1d-points " << q1.n_quadrature_points
+         << " 1d-points " << q1.size()
          << std::endl;
   
   Quadrature<dim-1> subquadrature(q1);
@@ -66,7 +66,7 @@ void check_face(Quadrature<1>& q1)
       
       Quadrature<dim> quadrature
        = QProjector<dim>::project_to_face(subquadrature, f);
-      for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
+      for (unsigned int k=0;k<quadrature.size();++k)
        deallog << quadrature.point(k) << std::endl;
     }
 
@@ -78,7 +78,7 @@ void check_face(Quadrature<1>& q1)
        
        Quadrature<dim> quadrature
          = QProjector<dim>::project_to_face(subquadrature, f);
-       for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
+       for (unsigned int k=0;k<quadrature.size();++k)
          deallog << quadrature.point(k) << std::endl;
       }
 }
@@ -86,7 +86,7 @@ void check_face(Quadrature<1>& q1)
 template<int dim>
 void check_faces (Quadrature<1>& q1)
 {
-  const unsigned int nq = q1.n_quadrature_points;
+  const unsigned int nq = q1.size();
   
   deallog << "Checking dim " << dim
          << " 1d-points " << nq
@@ -94,7 +94,7 @@ void check_faces (Quadrature<1>& q1)
 
 
   Quadrature<dim-1> subquadrature(q1);
-  const unsigned int nqs = subquadrature.n_quadrature_points;
+  const unsigned int nqs = subquadrature.size();
   
   Quadrature<dim> faces = QProjector<dim>::project_to_all_faces(subquadrature);
 
index b67be092934046988421b36cb895ee76b73570fd..aa1ea1e777016b0e1ab9d357958bf9b5cfd914f7 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007 by the deal.II authors
+//    Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -72,7 +72,7 @@ check_cells (std::vector<Quadrature<dim>*>& quadratures)
          quadrature_int=0;
                                           // Check the polynomial x^i*y^i
 
-         for (unsigned int x=0; x<quadrature.n_quadrature_points; ++x)
+         for (unsigned int x=0; x<quadrature.size(); ++x)
            {
              double f=1.;
              switch (dim)
@@ -101,13 +101,13 @@ check_cells (std::vector<Quadrature<dim>*>& quadratures)
                                           // check the ordering of
                                           // the quadrature points
          bool in_order=true;
-         for (unsigned int x=1; x<quadrature.n_quadrature_points; ++x)
+         for (unsigned int x=1; x<quadrature.size(); ++x)
            {
              if (points[x](0)<=points[x-1](0))
                in_order=false;
            }
          if (!in_order)
-           for (unsigned int x=0; x<quadrature.n_quadrature_points; ++x)
+           for (unsigned int x=0; x<quadrature.size(); ++x)
              deallog << points[x] << std::endl;
        }
     }  
@@ -146,7 +146,7 @@ check_faces (const std::vector<Quadrature<dim-1>*>& quadratures, const bool sub)
                                           // Check the polynomial
                                           // x^i*y^i*z^i
 
-         for (unsigned int x=0; x<quadrature.n_quadrature_points; ++x)
+         for (unsigned int x=0; x<quadrature.size(); ++x)
            {
              long double f=1.;
              switch (dim)
index 2cf444b9c06423db6ca5db84b1f5e2fbf54e085f..29d99bf934b6ec294c3d19c2ae6782c5a307a6c8 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005 by the deal.II authors
+//    Copyright (C) 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -54,7 +54,7 @@ check_this (const FiniteElement<dim> &fe,
   QGauss<dim> q_rhs(fe.degree+1>2 ? fe.degree+1-2 : 1);
       
   FullMatrix<double> X (fe.dofs_per_cell,
-                        q_rhs.n_quadrature_points);
+                        q_rhs.size());
 
   FETools::compute_projection_from_quadrature_points_matrix (fe,
                                                              q_lhs, q_rhs,
index 815bc90e6009b6ca5193a84543189615a6ff2e08..45e98b47a8d67694129ef1cdef81eb41d00b1325 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005 by the deal.II authors
+//    Copyright (C) 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -60,15 +60,15 @@ check_this (const FiniteElement<dim> &fe,
                                    // at least as many degrees of freedom in
                                    // the finite element as there are
                                    // quadrature points
-  if (fe.dofs_per_cell < q_rhs.n_quadrature_points)
+  if (fe.dofs_per_cell < q_rhs.size())
     return;
 
   deallog << "dofs_per_cell=" << fe.dofs_per_cell
-          << ", n_q_points=" << q_rhs.n_quadrature_points
+          << ", n_q_points=" << q_rhs.size()
           << std::endl;
   
   FullMatrix<double> X (fe.dofs_per_cell,
-                        q_rhs.n_quadrature_points);
+                        q_rhs.size());
 
   FETools::compute_projection_from_quadrature_points_matrix (fe,
                                                              q_lhs, q_rhs,
@@ -77,12 +77,12 @@ check_this (const FiniteElement<dim> &fe,
                                    // then compute the matrix that
                                    // interpolates back to the quadrature
                                    // points
-  FullMatrix<double> I_q (q_rhs.n_quadrature_points, fe.dofs_per_cell);
+  FullMatrix<double> I_q (q_rhs.size(), fe.dofs_per_cell);
   FETools::compute_interpolation_to_quadrature_points_matrix (fe, q_rhs,
                                                               I_q);
 
-  FullMatrix<double> product (q_rhs.n_quadrature_points,
-                              q_rhs.n_quadrature_points);
+  FullMatrix<double> product (q_rhs.size(),
+                              q_rhs.size());
   I_q.mmult (product, X);
 
                                    // the product should be the identity
index c8178f98156254a21738f4f5cd9dbe76a067ad7b..9e075a709c32997ee6513a0a3d7a047e4c5aa41d 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005 by the deal.II authors
+//    Copyright (C) 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -56,7 +56,7 @@ check_this (const FiniteElement<dim> &fe,
   QGauss<dim> q_rhs(fe.degree+3);
       
   FullMatrix<double> X (fe.dofs_per_cell,
-                        q_rhs.n_quadrature_points);
+                        q_rhs.size());
 
   FETools::compute_projection_from_quadrature_points_matrix (fe,
                                                              q_lhs, q_rhs,
index 8cbc834e9b1d76478fa5619dcdec239d1f23362c..93bebb68d31791d39a909850ce3b8fcaaed66028 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005 by the deal.II authors
+//    Copyright (C) 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -60,7 +60,7 @@ check_this (const FiniteElement<dim> &fe,
                                             1./fe.dofs_per_cell));
       
   FullMatrix<double> X (fe.dofs_per_cell,
-                        q_rhs.n_quadrature_points);
+                        q_rhs.size());
 
   Assert (X.m() == X.n(), ExcInternalError());
   
index 23348a12b1bc29049105a481ef14ff54914db753..5db3d9bd07c2f997a13dca058ec37fdb037c7fbf 100644 (file)
@@ -2,7 +2,7 @@
 //    joa_1.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2003, 2004, 2005, 2007 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2007, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -555,7 +555,7 @@ void LaplaceProblem<dim>::assemble_system ()
                           update_q_points  |  update_JxW_values);
 
   const unsigned int   dofs_per_cell = fe.dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula.n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula.size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
index 37d5ccc26567dde0cbd96e2688b48f63e1693e46..9b357e2dfdf9006a041d9c54f4cb2bedeb944d5f 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -72,7 +72,7 @@ void check_this (Triangulation<3> &tria)
             fe_face_values1.reinit (neighbor_child, neighbor_neighbor);
             fe_face_values2.reinit (cell, face_no, subface_no);
 
-            for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<quadrature.size(); ++q)
               {
                 Assert ((fe_face_values1.quadrature_point(q)-
                          fe_face_values2.quadrature_point(q)).square()
index e3a42a50a2f2fa3a93c6fa827c90891e117141d9..a5addd44170122cbc4726c370df90ddabf865048 100644 (file)
@@ -61,7 +61,7 @@ void check (Triangulation<3> &tria)
       for (unsigned int f=0; f<GeometryInfo<3>::faces_per_cell; ++f)
         {
           fe_face_values.reinit (cell, f);
-          for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+          for (unsigned int q=0; q<q_face.size(); ++q)
             n1 += fe_face_values.normal_vector(q) *
                   fe_face_values.JxW(q);
         }
@@ -76,7 +76,7 @@ void check (Triangulation<3> &tria)
         for (unsigned int sf=0; sf<GeometryInfo<3>::max_children_per_face; ++sf)
         {
           fe_subface_values.reinit (cell, f, sf);
-          for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+          for (unsigned int q=0; q<q_face.size(); ++q)
             n2 += fe_subface_values.normal_vector(q) *
                   fe_subface_values.JxW(q);
         }
index 4b0ed3f6890b236d81d4e94fe25bc0100ccffebe..62dcfd8bf3ef82e701e7a42aa6f14eb48a11aceb 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -60,7 +60,7 @@ void check_this (Triangulation<3> &tria)
           fe_face_values1.reinit (cell, f);
           fe_face_values2.reinit (cell->neighbor(f), nn);
 
-         for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+         for (unsigned int q=0; q<quadrature.size(); ++q)
             {
                                                // in order to reduce
                                                // output file size,
index d4df19c9244f26ed10780c44dac309284526ba7b..6b6bbe32f128988339fe69136d97d1942a0292b0 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -61,7 +61,7 @@ void check_this (Triangulation<3> &tria)
           deallog << "Cell " << cell << ", face " << f
                   << std::endl;
 
-          for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+          for (unsigned int q=0; q<quadrature.size(); ++q)
             {
                                                // in order to reduce
                                                // output file size,
index 84306165315b1b0e19fc8e6463dc82099ae58919..2d0e3eaa4c59806b8d4b89a5437c843906c9ecdf 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -59,7 +59,7 @@ void check (const Triangulation<dim> &tria)
       for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
         {
           fe_face_values.reinit (cell, f);
-          for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+          for (unsigned int q=0; q<q_face.size(); ++q)
             n1 += fe_face_values.normal_vector(q) *
                   fe_face_values.JxW(q);
         }
@@ -74,7 +74,7 @@ void check (const Triangulation<dim> &tria)
         for (unsigned int sf=0; sf<GeometryInfo<dim>::max_children_per_face; ++sf)
           {
             fe_subface_values.reinit (cell, f, sf);
-            for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<q_face.size(); ++q)
               n2 += fe_subface_values.normal_vector(q) *
                     fe_subface_values.JxW(q);
           }
index 882d6bce899cba67448ac199a87b0ff8833c44ad..9700002baaab88dabebbca18fa954e969657dac6 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -69,7 +69,7 @@ void check (const Triangulation<dim> &tria,
       for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
         {
           fe_face_values.reinit (cell, f);
-          for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+          for (unsigned int q=0; q<q_face.size(); ++q)
             n1 += fe_face_values.normal_vector(q) *
                   fe_face_values.JxW(q);
         }
@@ -84,7 +84,7 @@ void check (const Triangulation<dim> &tria,
         for (unsigned int sf=0; sf<GeometryInfo<dim>::max_children_per_face; ++sf)
           {
             fe_subface_values.reinit (cell, f, sf);
-            for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<q_face.size(); ++q)
               n2 += fe_subface_values.normal_vector(q) *
                     fe_subface_values.JxW(q);
           }
index 36eb7096f4a9f6c6829ed81fb227412728e2e746..e71e47e061a5af3ba4c0ecfb440fa7bf9b2d0509 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -64,7 +64,7 @@ void check (const Triangulation<dim> &tria)
       for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
         {
           fe_face_values.reinit (cell, f);
-          for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+          for (unsigned int q=0; q<q_face.size(); ++q)
             n1 += fe_face_values.normal_vector(q) *
                   fe_face_values.JxW(q);
         }
@@ -79,7 +79,7 @@ void check (const Triangulation<dim> &tria)
         for (unsigned int sf=0; sf<GeometryInfo<dim>::max_children_per_face; ++sf)
           {
             fe_subface_values.reinit (cell, f, sf);
-            for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<q_face.size(); ++q)
               n2 += fe_subface_values.normal_vector(q) *
                     fe_subface_values.JxW(q);
           }
index 321d72df80d3ee0a3e11e7ebf8719b0bb43dd353..f2410dfdb1196f17db0e029ac2c0251085f31757 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -64,7 +64,7 @@ void check (const Triangulation<dim> &tria)
       for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
         {
           fe_face_values.reinit (cell, f);
-          for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+          for (unsigned int q=0; q<q_face.size(); ++q)
             n1 += fe_face_values.normal_vector(q) *
                   fe_face_values.JxW(q);
         }
@@ -79,7 +79,7 @@ void check (const Triangulation<dim> &tria)
         for (unsigned int sf=0; sf<GeometryInfo<dim>::max_children_per_face; ++sf)
           {
             fe_subface_values.reinit (cell, f, sf);
-            for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<q_face.size(); ++q)
               n2 += fe_subface_values.normal_vector(q) *
                     fe_subface_values.JxW(q);
           }
index 3a218af35cb0dab517ecf2abd94a61991aec07f0..d63e37feb5abab694934ddfb8233387066adeba3 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -57,7 +57,7 @@ void check (const Triangulation<dim> &tria)
         if (cell->at_boundary(f))
           {
             fe_face_values.reinit (cell, f);
-            for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<q_face.size(); ++q)
               n1 += fe_face_values.quadrature_point(q) *
                     fe_face_values.JxW(q);
           }
@@ -69,7 +69,7 @@ void check (const Triangulation<dim> &tria)
           for (unsigned int sf=0; sf<GeometryInfo<dim>::max_children_per_face; ++sf)
             {
               fe_subface_values.reinit (cell, f, sf);
-              for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+              for (unsigned int q=0; q<q_face.size(); ++q)
                 n2 += fe_subface_values.quadrature_point(q) *
                       fe_subface_values.JxW(q);
             }
index 749556e59ceeeed3956a801bac34cc11ef1a2817..98085da53f47127e664e5cc2f101564ca33eed78 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -65,7 +65,7 @@ void check (const Triangulation<dim> &tria,
         if (cell->at_boundary(f))
           {
             fe_face_values.reinit (cell, f);
-            for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<q_face.size(); ++q)
               n1 += fe_face_values.quadrature_point(q) *
                     fe_face_values.JxW(q);
           }
@@ -77,7 +77,7 @@ void check (const Triangulation<dim> &tria,
           for (unsigned int sf=0; sf<GeometryInfo<dim>::max_children_per_face; ++sf)
             {
               fe_subface_values.reinit (cell, f, sf);
-              for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+              for (unsigned int q=0; q<q_face.size(); ++q)
                 n2 += fe_subface_values.quadrature_point(q) *
                       fe_subface_values.JxW(q);
             }
index 9382376adc194f31e99eca9eee6e81e681649eec..eb4843e9b5819b000f4ddb9eb31845cfa76ccd38 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -64,7 +64,7 @@ void check (const Triangulation<dim> &tria)
         if (cell->at_boundary(f))
           {
             fe_face_values.reinit (cell, f);
-            for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<q_face.size(); ++q)
               n1 += fe_face_values.quadrature_point(q) *
                     fe_face_values.JxW(q);
           }
@@ -76,7 +76,7 @@ void check (const Triangulation<dim> &tria)
           for (unsigned int sf=0; sf<GeometryInfo<dim>::max_children_per_face; ++sf)
             {
               fe_subface_values.reinit (cell, f, sf);
-              for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+              for (unsigned int q=0; q<q_face.size(); ++q)
                 n2 += fe_subface_values.quadrature_point(q) *
                       fe_subface_values.JxW(q);
             }
index 472765107033761d61eafdbba7b82e3f569798d3..2eea4a76f085f7a3e17e094812a1d080edc07a51 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -64,7 +64,7 @@ void check (const Triangulation<dim> &tria)
         if (cell->at_boundary(f))
           {
             fe_face_values.reinit (cell, f);
-            for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<q_face.size(); ++q)
               n1 += fe_face_values.quadrature_point(q) *
                     fe_face_values.JxW(q);
           }
@@ -76,7 +76,7 @@ void check (const Triangulation<dim> &tria)
           for (unsigned int sf=0; sf<GeometryInfo<dim>::max_children_per_face; ++sf)
             {
               fe_subface_values.reinit (cell, f, sf);
-              for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+              for (unsigned int q=0; q<q_face.size(); ++q)
                 n2 += fe_subface_values.quadrature_point(q) *
                       fe_subface_values.JxW(q);
             }
index 6d387c19a27db2ecdd48d78878689e6453cd8833..9b8010728a008187a5435a77fe1758f522db10e3 100644 (file)
@@ -105,7 +105,7 @@ void check (Triangulation<3> &tria)
           deallog << "Cell " << cell << ", face " << f
                   << std::endl;
 
-          for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+          for (unsigned int q=0; q<quadrature.size(); ++q)
             {
               deallog << "  " << fe_face_values1.quadrature_point(q)
                       << std::endl;
index 555db2c0305eb8f52fb2875762955c6094b27fb2..0e787f824ce1f1d7a81eebd8e35f08faaa53c093 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006 by the deal.II authors and Oliver Kayser-Herold
+//    Copyright (C) 2005, 2006, 2010 by the deal.II authors and Oliver Kayser-Herold
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -67,7 +67,7 @@ void evaluate_normal (DoFHandler<2>  &dof_handler,
                                           update_normal_vectors |
                                           update_JxW_values));
 
-  const unsigned int   n_q_face    = quad.n_quadrature_points;
+  const unsigned int   n_q_face    = quad.size();
   const unsigned int   n_components   = dof_handler.get_fe().n_components();
 
                                   // Cell iterators
index 9bc680901026318b3720d7dccff27d6f41b1d704..ed34157c4a6959ea171cc3ccbe5100d2607fd58d 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -402,7 +402,7 @@ namespace LaplaceSolver
                             update_gradients | update_JxW_values);
 
     const unsigned int   dofs_per_cell = fe->dofs_per_cell;
-    const unsigned int   n_q_points    = quadrature->n_quadrature_points;
+    const unsigned int   n_q_points    = quadrature->size();
 
     FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
 
@@ -525,7 +525,7 @@ namespace LaplaceSolver
                              update_JxW_values);
 
     const unsigned int   dofs_per_cell = this->fe->dofs_per_cell;
-    const unsigned int   n_q_points    = this->quadrature->n_quadrature_points;
+    const unsigned int   n_q_points    = this->quadrature->size();
 
     Vector<double>       cell_rhs (dofs_per_cell);
     std::vector<double>  rhs_values (n_q_points);
index 1c5170ba6af29758c1c324ffdd04f53363321d34..58665b850c9effffd264b7bc111381c7db32929b 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 2005, 2006, 2007, 2008, 2009 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2008, 2009, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -189,7 +189,7 @@ namespace Evaluation
                             vertex_quadrature,
                             update_gradients | update_q_points);
     std::vector<Tensor<1,dim> >
-      solution_gradients (vertex_quadrature.n_quadrature_points);
+      solution_gradients (vertex_quadrature.size());
 
     typename DoFHandler<dim>::active_cell_iterator
       cell = dof_handler.begin_active(),
@@ -489,7 +489,7 @@ namespace LaplaceSolver
                             update_gradients | update_JxW_values);
 
     const unsigned int   dofs_per_cell = fe->dofs_per_cell;
-    const unsigned int   n_q_points    = quadrature->n_quadrature_points;
+    const unsigned int   n_q_points    = quadrature->size();
 
     FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
 
@@ -672,7 +672,7 @@ namespace LaplaceSolver
                              update_JxW_values);
 
     const unsigned int   dofs_per_cell = this->fe->dofs_per_cell;
-    const unsigned int   n_q_points    = this->quadrature->n_quadrature_points;
+    const unsigned int   n_q_points    = this->quadrature->size();
 
     Vector<double>       cell_rhs (dofs_per_cell);
     std::vector<double>  rhs_values (n_q_points);
@@ -1476,7 +1476,7 @@ namespace LaplaceSolver
                  right_hand_side (&right_hand_side)
   {
     const unsigned int n_q_points
-      = quadrature.n_quadrature_points;
+      = quadrature.size();
 
     cell_residual.resize(n_q_points);
     rhs_values.resize(n_q_points);
@@ -1505,7 +1505,7 @@ namespace LaplaceSolver
                                          update_gradients)
   {
     const unsigned int n_face_q_points
-      = face_quadrature.n_quadrature_points;
+      = face_quadrature.size();
 
     jump_residual.resize(n_face_q_points);
     dual_weights.resize(n_face_q_points);
index d9761a95d6dbdd3a0b7c3150a8911abd1dc2a905..c5f1246802916ee8bf413fab42ca59607dadb0a6 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -191,7 +191,7 @@ void MinimizationProblem<dim>::assemble_step ()
                            update_q_points | update_JxW_values);
 
   const unsigned int dofs_per_cell = fe.dofs_per_cell;
-  const unsigned int n_q_points    = quadrature_formula.n_quadrature_points;
+  const unsigned int n_q_points    = quadrature_formula.size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
@@ -400,9 +400,9 @@ void MinimizationProblem<1>::refine_grid ()
   FEValues<dim> neighbor_fe_values (fe, quadrature,
                                     update_gradients);
 
-  std::vector<double> local_values (quadrature.n_quadrature_points);
-  std::vector<Tensor<1,dim> > local_gradients (quadrature.n_quadrature_points);
-  std::vector<Tensor<2,dim> > local_2nd_derivs (quadrature.n_quadrature_points);
+  std::vector<double> local_values (quadrature.size());
+  std::vector<Tensor<1,dim> > local_gradients (quadrature.size());
+  std::vector<Tensor<2,dim> > local_2nd_derivs (quadrature.size());
 
   DoFHandler<dim>::active_cell_iterator
     cell = dof_handler.begin_active (),
@@ -415,7 +415,7 @@ void MinimizationProblem<1>::refine_grid ()
       fe_values.get_function_2nd_derivatives (present_solution, local_2nd_derivs);
 
       double cell_residual_norm = 0;
-      for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+      for (unsigned int q=0; q<quadrature.size(); ++q)
         {
           const double x             = fe_values.quadrature_point(q)[0];
           const double u             = local_values[q];
@@ -518,7 +518,7 @@ MinimizationProblem<dim>::energy (const DoFHandler<dim> &dof_handler,
                           update_values   | update_gradients |
                            update_q_points | update_JxW_values);
 
-  const unsigned int   n_q_points    = quadrature_formula.n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula.size();
 
   std::vector<double>         local_solution_values (n_q_points);
   std::vector<Tensor<1,dim> > local_solution_grads (n_q_points);
index 9e51efaca463dfc8f01117151e8c5a2b88ca51ae..106eca2829a8c050ee44289b8e47a17c39aab377 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -139,7 +139,7 @@ void LaplaceProblem<dim>::assemble_system ()
                            update_q_points | update_JxW_values);
 
   const unsigned int   dofs_per_cell = fe.dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula.n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula.size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
@@ -196,7 +196,7 @@ void LaplaceProblem<dim>::assemble_multigrid ()
                            update_q_points | update_JxW_values);
 
   const unsigned int   dofs_per_cell = fe.dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula.n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula.size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
 
index fe72c0b3d96422ea1ac8cb4370693f9a7cc2c76c..ca323f91f1686b85743feede1f1ea36e0645d3cb 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -120,7 +120,7 @@ void LaplaceProblem::assemble_system ()
                         update_values | update_gradients | update_JxW_values);
   
   const unsigned int   dofs_per_cell = fe.dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula.n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula.size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
index bbb0fe6f641a2e145458cdb7d863b002da0e2f71..af30712ed2e7dac51cfdeb689df24ed475caed25 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -178,7 +178,7 @@ void LaplaceProblem<dim>::assemble_system ()
                            update_q_points | update_JxW_values);
 
   const unsigned int   dofs_per_cell = fe.dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula.n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula.size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
index 29ba8f1114641baee6aba169668936477ba0620e..d012acd0fd4ad60f2137baa523e0065bd306c1e4 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -174,7 +174,7 @@ void LaplaceProblem<dim>::assemble_system ()
                           update_q_points  |  update_JxW_values);
 
   const unsigned int   dofs_per_cell = fe.dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula.n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula.size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
index 2815a27c2972a9c024ceffa84aa06497a214c299..24228a0da3aa4e4df89154beaa112bf428c146fd 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -199,7 +199,7 @@ void LaplaceProblem<dim>::assemble_system ()
                           update_q_points  |  update_JxW_values);
 
   const unsigned int   dofs_per_cell = fe.dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula.n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula.size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
index 3336c0b02d4cc82d988c0142fc51433a9f6b1c51..d784cb6298c4969a102e100fcd0ef9224d574c7f 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -266,8 +266,8 @@ void HelmholtzProblem<dim>::assemble_system ()
   QGauss<dim>   quadrature_formula (3);
   QGauss<dim-1> face_quadrature_formula (3);
 
-  const unsigned int n_q_points    = quadrature_formula.n_quadrature_points;
-  const unsigned int n_face_q_points = face_quadrature_formula.n_quadrature_points;
+  const unsigned int n_q_points    = quadrature_formula.size();
+  const unsigned int n_face_q_points = face_quadrature_formula.size();
 
   const unsigned int dofs_per_cell = (*fe).dofs_per_cell;
 
index 03758a5782ab4de523483746a75e3d9e634a7471..1474af915901b80bab32013bac21cd5dd68a672d 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -205,7 +205,7 @@ void ElasticProblem<dim>::assemble_system ()
                            update_q_points | update_JxW_values);
 
   const unsigned int   dofs_per_cell = fe.dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula.n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula.size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
index b836bac1153a7f374c7646c54ded5fdd69a2a272..30435613cc14a0a90404c263ccd7cb4e1d7e46ca 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -63,7 +63,7 @@ void check (const Triangulation<dim> &tria)
         if (cell->at_boundary(f))
           {
             fe_face_values.reinit (cell, f);
-            for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<q_face.size(); ++q)
               v1 += (fe_face_values.normal_vector(q) *
                      fe_face_values.quadrature_point(q)) *
                     fe_face_values.JxW(q);
@@ -76,7 +76,7 @@ void check (const Triangulation<dim> &tria)
           for (unsigned int sf=0; sf<GeometryInfo<dim>::max_children_per_face; ++sf)
             {
               fe_subface_values.reinit (cell, f, sf);
-              for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+              for (unsigned int q=0; q<q_face.size(); ++q)
                 v2 += (fe_subface_values.normal_vector(q) *
                        fe_subface_values.quadrature_point(q)) *
                       fe_subface_values.JxW(q);
index 437f3413b9efdf9079b2151a6d524f6caf705b11..3bde6a01be79cbb055d21641089ecc4d2c1e045c 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -73,7 +73,7 @@ void check (const Triangulation<dim> &tria,
         if (cell->at_boundary(f))
           {
             fe_face_values.reinit (cell, f);
-            for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<q_face.size(); ++q)
               v1 += (fe_face_values.normal_vector(q) *
                      fe_face_values.quadrature_point(q)) *
                     fe_face_values.JxW(q);
@@ -86,7 +86,7 @@ void check (const Triangulation<dim> &tria,
           for (unsigned int sf=0; sf<GeometryInfo<dim>::max_children_per_face; ++sf)
             {
               fe_subface_values.reinit (cell, f, sf);
-              for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+              for (unsigned int q=0; q<q_face.size(); ++q)
                 v2 += (fe_subface_values.normal_vector(q) *
                        fe_subface_values.quadrature_point(q)) *
                       fe_subface_values.JxW(q);
index a74f771970459a5c2cc17289e8766b38986e02ed..a256679dedd60d5bcb63facde19b669b341b5eb3 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -69,7 +69,7 @@ void check (const Triangulation<dim> &tria)
         if (cell->at_boundary(f))
           {
             fe_face_values.reinit (cell, f);
-            for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<q_face.size(); ++q)
               v1 += (fe_face_values.normal_vector(q) *
                      fe_face_values.quadrature_point(q)) *
                     fe_face_values.JxW(q);
@@ -82,7 +82,7 @@ void check (const Triangulation<dim> &tria)
           for (unsigned int sf=0; sf<GeometryInfo<dim>::max_children_per_face; ++sf)
             {
               fe_subface_values.reinit (cell, f, sf);
-              for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+              for (unsigned int q=0; q<q_face.size(); ++q)
                 v2 += (fe_subface_values.normal_vector(q) *
                        fe_subface_values.quadrature_point(q)) *
                       fe_subface_values.JxW(q);
index e6438a31aa27ee516fdf6908e2be12bb1caf62f7..2d499188006140b788c357fcfd011f2ce6ae935d 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2003, 2004, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2004, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -69,7 +69,7 @@ void check (const Triangulation<dim> &tria)
         if (cell->at_boundary(f))
           {
             fe_face_values.reinit (cell, f);
-            for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+            for (unsigned int q=0; q<q_face.size(); ++q)
               v1 += (fe_face_values.normal_vector(q) *
                      fe_face_values.quadrature_point(q)) *
                     fe_face_values.JxW(q);
@@ -82,7 +82,7 @@ void check (const Triangulation<dim> &tria)
           for (unsigned int sf=0; sf<GeometryInfo<dim>::max_children_per_face; ++sf)
             {
               fe_subface_values.reinit (cell, f, sf);
-              for (unsigned int q=0; q<q_face.n_quadrature_points; ++q)
+              for (unsigned int q=0; q<q_face.size(); ++q)
                 v2 += (fe_subface_values.normal_vector(q) *
                        fe_subface_values.quadrature_point(q)) *
                       fe_subface_values.JxW(q);
index 61ea6ac4394b5fe30cd45b2590842f651f2ddb82..25ccbfc74a627fd7fd7613117b9d7c65b60b3e0b 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2008 by the deal.II authors 
+//    Copyright (C) 2005, 2008, 2010 by the deal.II authors 
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -566,8 +566,8 @@ BEM<spacedim>::solve()
                                                  update_gradients |
                                                  update_cell_normal_vectors );
    
-    std::vector< Point<spacedim> > cell_normals(q_iterated.n_quadrature_points);
-    std::vector< Point<spacedim> > cell_tangentials(q_iterated.n_quadrature_points);
+    std::vector< Point<spacedim> > cell_normals(q_iterated.size());
+    std::vector< Point<spacedim> > cell_tangentials(q_iterated.size());
     std::vector<unsigned int> local_dof_indices (fe_q.dofs_per_cell);
 
     typename DoFHandler<spacedim-1, spacedim>::active_cell_iterator
@@ -593,7 +593,7 @@ BEM<spacedim>::solve()
 
        
        cell_normals = fe_values_q.get_cell_normal_vectors();
-       for(unsigned int i=0; i<q_iterated.n_quadrature_points; ++i)
+       for(unsigned int i=0; i<q_iterated.size(); ++i)
        {
            cell_tangentials[i][0] = cell_normals[i][1]; 
            cell_tangentials[i][1] = -cell_normals[i][0];
@@ -605,8 +605,8 @@ BEM<spacedim>::solve()
                                // factor (smooth_solution..) is taken
                                // so that it is the coefficient of
                                // the fun_th shape function of the cell.
-       std::vector< Tensor<1,spacedim> > gradient(q_iterated.n_quadrature_points);
-       for(unsigned int pnt=0; pnt<q_iterated.n_quadrature_points; ++pnt)
+       std::vector< Tensor<1,spacedim> > gradient(q_iterated.size());
+       for(unsigned int pnt=0; pnt<q_iterated.size(); ++pnt)
            for(unsigned int fun=0; fun<fe_q.dofs_per_cell; ++fun)
            {
            gradient[pnt]+=
@@ -635,7 +635,7 @@ BEM<spacedim>::solve()
 
 //     std::cout
 //         << "gradienti "<<std::endl;
-//     for(unsigned int pnt=0; pnt<q_iterated.n_quadrature_points; ++pnt)
+//     for(unsigned int pnt=0; pnt<q_iterated.size(); ++pnt)
 //         std::cout
 //             << pnt << " - "
 //             << gradient[pnt]
@@ -647,7 +647,7 @@ BEM<spacedim>::solve()
 //         << std::endl;
 
 
-       for(unsigned int pnt=0; pnt<q_iterated.n_quadrature_points; ++pnt)
+       for(unsigned int pnt=0; pnt<q_iterated.size(); ++pnt)
        {
            tangential_derivative(cell->index())=
                contract(
@@ -670,7 +670,7 @@ BEM<spacedim>::solve()
                
        }
 
-// //  for(unsigned int pnt=0; pnt<q_iterated.n_quadrature_points; ++pnt)
+// //  for(unsigned int pnt=0; pnt<q_iterated.size(); ++pnt)
 //         std::cout
 // //          << pnt <<" - "
 //             << tangential_derivative(cell->index())
index 81a0c0406c908c3e8e04feeb9b7457b9272fd66c..9af95b8dd925c5f3d1ebb785bb9f26819bc3ed81 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 2005, 2008, 2009 by the deal.II authors
+//    Copyright (C) 2005, 2008, 2009, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -121,8 +121,8 @@ void test(std::string filename)
 
   FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
 
-  std::vector< Point<spacedim> > cell_normals(q_midpoint.n_quadrature_points);
-  std::vector< Point<spacedim> > cell_tangentials(q_midpoint.n_quadrature_points);
+  std::vector< Point<spacedim> > cell_normals(q_midpoint.size());
+  std::vector< Point<spacedim> > cell_tangentials(q_midpoint.size());
   std::vector<double> shape_directional_derivative(dofs_per_cell);
   Vector<double> projected_directional_derivative(triangulation.n_cells());
 
index dd43a3935ad94ae9a2d1041a0bcf1122af1bc37b..82bc5a943fa45c9c8ca7266edcaa39aa065a2806 100644 (file)
@@ -3,7 +3,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005 by the deal.II authors 
+//    Copyright (C) 2005, 2010 by the deal.II authors 
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -54,7 +54,7 @@ void test(std::string filename) {
   for(;cell!=endc;++cell)
     {
       deallog<<cell<< std::endl;       
-      for(unsigned int q=0; q<quad.n_quadrature_points; ++q)
+      for(unsigned int q=0; q<quad.size(); ++q)
        {
          real = mapping.transform_unit_to_real_cell(cell, quad.point(q));
          // unit = mapping.transform_real_to_unit_cell(cell, real);
index 9e57327c9da5966c36827e377b52d0c2924697c3..e0b37c0b9b45dd54a939124cafd2a71e3b9a950c 100644 (file)
@@ -3,7 +3,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2009 by the deal.II authors 
+//    Copyright (C) 2005, 2009, 2010 by the deal.II authors 
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -69,7 +69,7 @@ void test(std::string filename) {
   for(;cell!=endc;++cell)
     {
       deallog<<cell<< std::endl;       
-      for(unsigned int q=0; q<quad.n_quadrature_points; ++q)
+      for(unsigned int q=0; q<quad.size(); ++q)
        {
          real = mapping.transform_unit_to_real_cell(cell, quad.point(q));
          // unit = mapping.transform_real_to_unit_cell(cell, real);
index d3a43b99b974bf532ca4b365c7a4bfda25e79c83..cc1220af7d6a5d20b9864bb29883b8bf1a8ecd47 100644 (file)
@@ -263,7 +263,7 @@ void LaplaceProblem<Vector,Matrix,Sparsity>::assemble_system ()
                                     update_JxW_values));
 
   const unsigned int   dofs_per_cell = fe.dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula.n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula.size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   ::Vector<double>     cell_rhs (dofs_per_cell);
index 73dbcedcded7c0e724df41d97438c96b95aa31df..2752eb389685e476315fa506926d33df9eb834a9 100644 (file)
@@ -75,7 +75,7 @@ check (const unsigned int level,
                    | update_q_points | update_JxW_values);
 
   std::vector <unsigned int> global_dofs (element.dofs_per_cell);
-  std::vector <double> function (quadrature.n_quadrature_points);
+  std::vector <double> function (quadrature.size());
 
   Vector<double> u (dof.n_dofs ());
   Vector<double> f (dof.n_dofs ());
@@ -96,7 +96,7 @@ check (const unsigned int level,
       cell->get_dof_indices (global_dofs);
       cosine.value_list (fe.get_quadrature_points(), function);
 
-      for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
+      for (unsigned int k=0;k<quadrature.size();++k)
        {
          double dx = fe.JxW (k);
 
@@ -132,11 +132,11 @@ check (const unsigned int level,
   double h1 = 0.;
   double h2 = 0.;
 
-  std::vector<double> u_local (quadrature.n_quadrature_points);
-  std::vector<Tensor<1,dim> > Du (quadrature.n_quadrature_points);
-  std::vector<Tensor<1,dim> > Df (quadrature.n_quadrature_points);
-  std::vector<Tensor<2,dim> > DDu (quadrature.n_quadrature_points);
-  std::vector<Tensor<2,dim> > DDf (quadrature.n_quadrature_points);
+  std::vector<double> u_local (quadrature.size());
+  std::vector<Tensor<1,dim> > Du (quadrature.size());
+  std::vector<Tensor<1,dim> > Df (quadrature.size());
+  std::vector<Tensor<2,dim> > DDu (quadrature.size());
+  std::vector<Tensor<2,dim> > DDf (quadrature.size());
   
   for (cell = dof.begin_active(); cell != end;++cell)
     {
@@ -149,7 +149,7 @@ check (const unsigned int level,
       fe2.get_function_grads (u, Du);
       fe2.get_function_2nd_derivatives (u,DDu);
       
-      for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
+      for (unsigned int k=0;k<quadrature.size();++k)
        {
          const double dx = fe.JxW (k);
          double e =  u_local[k];
index 6a0f06cb1c03ce8ccdfd55ae3888d7456eab0618..14820b78731a305fe495099dde97586afd7f1fd7 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2007, 2008 by the deal.II authors
+//    Copyright (C) 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -52,8 +52,8 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<double> scalar_values (quadrature.n_quadrature_points);
-  std::vector<Vector<double> > vector_values (quadrature.n_quadrature_points,
+  std::vector<double> scalar_values (quadrature.size());
+  std::vector<Vector<double> > vector_values (quadrature.size(),
                                              Vector<double>(fe.n_components()));
 
   fe_values.get_function_values (fe_function, vector_values);
index f3def277196d6bdef6b6c7cf0b4c43a254e3fea8..2af650b40c1b31a90e55be30a4b0176e1beabc8a 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2007, 2008 by the deal.II authors
+//    Copyright (C) 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -52,9 +52,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<Tensor<1,dim> > scalar_values (quadrature.n_quadrature_points);
+  std::vector<Tensor<1,dim> > scalar_values (quadrature.size());
   std::vector<std::vector<Tensor<1,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<1,dim> >(fe.n_components()));
 
   fe_values.get_function_gradients (fe_function, vector_values);
index 0b86d9a67741d65885cf05e9e9ec85377a01b6f5..7ef4e5a8ab634405e3d42a2c021d8029700cbe2a 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2007, 2008 by the deal.II authors
+//    Copyright (C) 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -52,9 +52,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<Tensor<2,dim> > scalar_values (quadrature.n_quadrature_points);
+  std::vector<Tensor<2,dim> > scalar_values (quadrature.size());
   std::vector<std::vector<Tensor<2,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<2,dim> >(fe.n_components()));
 
   fe_values.get_function_hessians (fe_function, vector_values);
index 6f08f0848388a7b4fef0edc802e5a6d9d5f97d96..c1a6af2830f52f629bec5a26e8d36c821e4016b7 100644 (file)
@@ -54,8 +54,8 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<double> scalar_values (quadrature.n_quadrature_points);
-  std::vector<Vector<double> > vector_values (quadrature.n_quadrature_points,
+  std::vector<double> scalar_values (quadrature.size());
+  std::vector<Vector<double> > vector_values (quadrature.size(),
                                              Vector<double>(fe.n_components()));
 
   fe_values.get_function_values (fe_function, vector_values);
index 892ae7959813bcc583fbaa0ad66a274deec1e681..f4ec9b4ed93a54f0ea4dc8e2448ccfac424728d7 100644 (file)
@@ -54,9 +54,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<Tensor<1,dim> > scalar_values (quadrature.n_quadrature_points);
+  std::vector<Tensor<1,dim> > scalar_values (quadrature.size());
   std::vector<std::vector<Tensor<1,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<1,dim> >(fe.n_components()));
 
   fe_values.get_function_gradients (fe_function, vector_values);
index dbbdc7f7a21dc09e10241152595c4a88536eef7d..18456a53dc297ee64710f5f61dbd75d3d9435a75 100644 (file)
@@ -54,9 +54,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<Tensor<2,dim> > scalar_values (quadrature.n_quadrature_points);
+  std::vector<Tensor<2,dim> > scalar_values (quadrature.size());
   std::vector<std::vector<Tensor<2,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<2,dim> >(fe.n_components()));
 
   fe_values.get_function_hessians (fe_function, vector_values);
index 5413dfae1c27e3b3ef9e164fc10ad0a47599b1c2..a6aee7bf700055e9060de4e2287d10afa98c18f1 100644 (file)
@@ -57,9 +57,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<Tensor<2,dim> > scalar_values (quadrature.n_quadrature_points);
+  std::vector<Tensor<2,dim> > scalar_values (quadrature.size());
   std::vector<std::vector<Tensor<2,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<2,dim> >(fe.n_components()));
 
   fe_values.get_function_hessians (fe_function, vector_values);
index 99cf71533594179f9a64759cde31da6d8d94cd40..5fa6702355e1dc25f4da7c42a013458a9ee9e687 100644 (file)
@@ -57,9 +57,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<Tensor<2,dim> > scalar_values (quadrature.n_quadrature_points);
+  std::vector<Tensor<2,dim> > scalar_values (quadrature.size());
   std::vector<std::vector<Tensor<2,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<2,dim> >(fe.n_components()));
 
   fe_values.get_function_hessians (fe_function, vector_values);
index 88c359a8a3f7244d771553482ced90c8bbc7716c..5a5f5b5a2453206fa63b1f055aec89eff2ed68b6 100644 (file)
@@ -57,9 +57,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<Tensor<2,dim> > scalar_values (quadrature.n_quadrature_points);
+  std::vector<Tensor<2,dim> > scalar_values (quadrature.size());
   std::vector<std::vector<Tensor<2,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<2,dim> >(fe.n_components()));
 
   fe_values.get_function_hessians (fe_function, vector_values);
index 254390275095bc8797579f02144eb658d77ba862..071c99b343f71389d7954090143dc8526b3114ac 100644 (file)
@@ -59,9 +59,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<double> scalar_values (quadrature.n_quadrature_points);
+  std::vector<double> scalar_values (quadrature.size());
   std::vector<Vector<double> >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   Vector<double>(fe.n_components()));
 
   fe_values.get_function_values (fe_function, vector_values);
index 8a1aba746cad99e381ff9a8dde9e15261d4cc5eb..ab091d0a10be266f87b031eb6280bb3f6ed92083 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2007, 2008 by the deal.II authors
+//    Copyright (C) 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -52,8 +52,8 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<Tensor<1,dim> > selected_vector_values (quadrature.n_quadrature_points);
-  std::vector<Vector<double> > vector_values (quadrature.n_quadrature_points,
+  std::vector<Tensor<1,dim> > selected_vector_values (quadrature.size());
+  std::vector<Vector<double> > vector_values (quadrature.size(),
                                              Vector<double>(fe.n_components()));
 
   fe_values.get_function_values (fe_function, vector_values);
index 5bfe2c8da6385184198aec7011df283af2f7e018..efdad29d6db61d0b2c4eea600176738012536e72 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2007, 2008 by the deal.II authors
+//    Copyright (C) 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -52,9 +52,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<Tensor<2,dim> > selected_vector_values (quadrature.n_quadrature_points);
+  std::vector<Tensor<2,dim> > selected_vector_values (quadrature.size());
   std::vector<std::vector<Tensor<1,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<1,dim> >(fe.n_components()));
 
   fe_values.get_function_gradients (fe_function, vector_values);
index 8267200fb6998ed7167eb0f3e861548578cf5510..26dde699f6e08b441180158b3e80c86a179b27e3 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2007, 2008 by the deal.II authors
+//    Copyright (C) 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -52,9 +52,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<Tensor<3,dim> > selected_vector_values (quadrature.n_quadrature_points);
+  std::vector<Tensor<3,dim> > selected_vector_values (quadrature.size());
   std::vector<std::vector<Tensor<2,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<2,dim> >(fe.n_components()));
 
   fe_values.get_function_hessians (fe_function, vector_values);
index 53f85f61afbe3649c993b3c03670256b96beb990..e10223c2d010ef38e7a1038915157545a3ae3233 100644 (file)
@@ -54,8 +54,8 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<Tensor<1,dim> > selected_vector_values (quadrature.n_quadrature_points);
-  std::vector<Vector<double> > vector_values (quadrature.n_quadrature_points,
+  std::vector<Tensor<1,dim> > selected_vector_values (quadrature.size());
+  std::vector<Vector<double> > vector_values (quadrature.size(),
                                              Vector<double>(fe.n_components()));
 
   fe_values.get_function_values (fe_function, vector_values);
index db5774b1c3a3ebdb0a2a73c89652d17a68f7b4c1..de50ae94add1c5e1b46b4e361078c973eaff2fa2 100644 (file)
@@ -54,9 +54,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<Tensor<2,dim> > selected_vector_values (quadrature.n_quadrature_points);
+  std::vector<Tensor<2,dim> > selected_vector_values (quadrature.size());
   std::vector<std::vector<Tensor<1,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<1,dim> >(fe.n_components()));
 
   fe_values.get_function_gradients (fe_function, vector_values);
index 2a63ba31ebbd1164e50f99e31d3e4f83f9165e23..fea3184986df802ab4feb0e84e801a61e96fd5cc 100644 (file)
@@ -54,9 +54,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<Tensor<3,dim> > selected_vector_values (quadrature.n_quadrature_points);
+  std::vector<Tensor<3,dim> > selected_vector_values (quadrature.size());
   std::vector<std::vector<Tensor<2,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<2,dim> >(fe.n_components()));
 
   fe_values.get_function_hessians (fe_function, vector_values);
index ca8d025f94cbec953785d27ead26d8d7129bb680..823f27b511cdc77824d0cb47cd066f94d2a60289 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2007, 2008 by the deal.II authors
+//    Copyright (C) 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -52,9 +52,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<SymmetricTensor<2,dim> > selected_vector_values (quadrature.n_quadrature_points);
+  std::vector<SymmetricTensor<2,dim> > selected_vector_values (quadrature.size());
   std::vector<std::vector<Tensor<1,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<1,dim> >(fe.n_components()));
 
   fe_values.get_function_gradients (fe_function, vector_values);
index 7786f9bf2e9d21c747194f51079dd5ee35758c5d..c3dfa2e70a8e7ae9c580ce761485e4a788cdb29c 100644 (file)
@@ -54,9 +54,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<SymmetricTensor<2,dim> > selected_vector_values (quadrature.n_quadrature_points);
+  std::vector<SymmetricTensor<2,dim> > selected_vector_values (quadrature.size());
   std::vector<std::vector<Tensor<1,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<1,dim> >(fe.n_components()));
 
   fe_values.get_function_gradients (fe_function, vector_values);
index bafdab36867fa40d194ce909a46b3965e8c345a0..a9b1b08235aa5992c57b57bd77dddd4cbe9c7a83 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2007, 2008 by the deal.II authors
+//    Copyright (C) 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -52,9 +52,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<double> selected_vector_values (quadrature.n_quadrature_points);
+  std::vector<double> selected_vector_values (quadrature.size());
   std::vector<std::vector<Tensor<1,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<1,dim> >(fe.n_components()));
 
   fe_values.get_function_gradients (fe_function, vector_values);
index fdf288d078d215249c7e8d76c55ebadc2c972cfd..70dc0963df6949147b35e7f56ffd2e8a9799f266 100644 (file)
@@ -54,9 +54,9 @@ void test (const Triangulation<dim>& tr,
                           update_values | update_gradients | update_hessians);
   fe_values.reinit (dof.begin_active());
 
-  std::vector<double> selected_vector_values (quadrature.n_quadrature_points);
+  std::vector<double> selected_vector_values (quadrature.size());
   std::vector<std::vector<Tensor<1,dim> > >
-    vector_values (quadrature.n_quadrature_points,
+    vector_values (quadrature.size(),
                   std::vector<Tensor<1,dim> >(fe.n_components()));
 
   fe_values.get_function_gradients (fe_function, vector_values);
index 738e4832417e581369fc316fbea815b32a1b64ee..2769247fd039c3fb3c344e2f3553eec8efdb67bc 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$
 //
-//    Copyright (C) 2009 by Andrew McBride and the deal.II authors
+//    Copyright (C) 2009, 2010 by Andrew McBride and the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -165,8 +165,7 @@ void MixedElastoPlasticity<dim>::assemble_system()
 
   const unsigned int dofs_per_cell = fe.dofs_per_cell;
   deallog << "dofs_per_cell: " << fe.dofs_per_cell << std::endl;
-  const unsigned int n_q_points = quadrature_formula.size();
-
+  
   FullMatrix<double> local_matrix(dofs_per_cell, dofs_per_cell);
   Vector<double> local_rhs(dofs_per_cell);
 
index f2f3ed4a0d524756092cef00c580953187e8e674..0bc3d37d7fda32952ecd17ea129cbcaeae81d26b 100644 (file)
@@ -122,7 +122,7 @@ check ()
                    | update_q_points | update_JxW_values);
 
   std::vector <unsigned int> global_dofs (element.dofs_per_cell);
-  std::vector <double> function (quadrature.n_quadrature_points);
+  std::vector <double> function (quadrature.size());
 
   Vector<double> f (dof.n_dofs ());
 
@@ -142,7 +142,7 @@ check ()
       cell->get_dof_indices (global_dofs);
       cosine.value_list (fe.get_quadrature_points(), function);
 
-      for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
+      for (unsigned int k=0;k<quadrature.size();++k)
        {
          double dx = fe.JxW (k);
          
index 7f04e514db7d7a2b063a4c979801adbcd9319bd0..544c8364a968507da77d858b9895763d17bb80b5 100644 (file)
@@ -372,7 +372,7 @@ double TestProjection (Mapping<2> &mapping,
                                         update_JxW_values|
                                         update_contravariant_transformation));
 
-      const unsigned int   n_q_points    = quad.n_quadrature_points;
+      const unsigned int   n_q_points    = quad.size();
       const unsigned int   n_components   = dof_handler->get_fe().n_components();
       
                                       // Cell iterators
@@ -452,7 +452,7 @@ int main ()
   QGauss<2> quad_temp(6);
   sprintf (buf, "DoFs per Quad: %i per Line %i per Vert %i\n", fe.dofs_per_quad, fe.dofs_per_line, fe.dofs_per_vertex);
   deallog << buf;
-  sprintf (buf, "QPoints %i\n", quad_temp.n_quadrature_points);
+  sprintf (buf, "QPoints %i\n", quad_temp.size());
   deallog << buf;
 
                                   // Create an deformation object for the Eulerian mapping
index 71ce22a3c6308b38667afbd056fde5f6ee408ab5..ea3783f6f7c9e6bfd924ca88080c91b7a6386678 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -371,7 +371,7 @@ estimate_smoothness (Vector<float> &smoothness_indicators) const
        for (unsigned int i=0; i<fe_collection[fe].dofs_per_cell; ++i)
          {
            std::complex<double> sum = 0;
-           for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+           for (unsigned int q=0; q<quadrature.size(); ++q)
              {
                const Point<dim> x_q = quadrature.point(q);
                sum += std::exp(std::complex<double>(0,1) *
index d3105522906379045e2749d0d544c1feb98bd158..515f683ec441d5120916a5480dc66b056a1c2ebe 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2007 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -189,7 +189,7 @@ namespace Evaluation
                             vertex_quadrature,
                             update_gradients | update_q_points);
     std::vector<Tensor<1,dim> >
-      solution_gradients (vertex_quadrature.n_quadrature_points);
+      solution_gradients (vertex_quadrature.size());
     
     typename hp::DoFHandler<dim>::active_cell_iterator
       cell = dof_handler.begin_active(),
@@ -486,7 +486,7 @@ namespace LaplaceSolver
                             update_gradients | update_JxW_values);
 
     const unsigned int   dofs_per_cell = (*fe)[0].dofs_per_cell;
-    const unsigned int   n_q_points    = (*quadrature)[0].n_quadrature_points;
+    const unsigned int   n_q_points    = (*quadrature)[0].size();
 
     FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
 
@@ -669,7 +669,7 @@ namespace LaplaceSolver
                              update_JxW_values);
 
     const unsigned int   dofs_per_cell = (*this->fe)[0].dofs_per_cell;
-    const unsigned int   n_q_points    = (*this->quadrature)[0].n_quadrature_points;
+    const unsigned int   n_q_points    = (*this->quadrature)[0].size();
 
     Vector<double>       cell_rhs (dofs_per_cell);
     std::vector<double>  rhs_values (n_q_points);
@@ -1473,7 +1473,7 @@ namespace LaplaceSolver
                  right_hand_side (&right_hand_side)
   {  
     const unsigned int n_q_points
-      = quadrature[0].n_quadrature_points;
+      = quadrature[0].size();
   
     cell_residual.resize(n_q_points);
     rhs_values.resize(n_q_points);    
@@ -1502,7 +1502,7 @@ namespace LaplaceSolver
                                          update_gradients)
   {  
     const unsigned int n_face_q_points
-      = face_quadrature[0].n_quadrature_points;
+      = face_quadrature[0].size();
   
     jump_residual.resize(n_face_q_points);
     dual_weights.resize(n_face_q_points);    
index d379f3f1b259622541b568401e503fd33a42d465..511c1ae2336660e73530afc701765a775db7b882 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -192,7 +192,7 @@ void MinimizationProblem<dim>::assemble_step ()
                            update_q_points | update_JxW_values);
 
   const unsigned int dofs_per_cell = fe[0].dofs_per_cell;
-  const unsigned int n_q_points    = quadrature_formula[0].n_quadrature_points;
+  const unsigned int n_q_points    = quadrature_formula[0].size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
@@ -402,9 +402,9 @@ void MinimizationProblem<1>::refine_grid ()
   hp::FEValues<dim> neighbor_fe_values (fe, quadrature,
                                     update_gradients);
 
-  std::vector<double> local_values (quadrature[0].n_quadrature_points);
-  std::vector<Tensor<1,dim> > local_gradients (quadrature[0].n_quadrature_points);
-  std::vector<Tensor<2,dim> > local_2nd_derivs (quadrature[0].n_quadrature_points);
+  std::vector<double> local_values (quadrature[0].size());
+  std::vector<Tensor<1,dim> > local_gradients (quadrature[0].size());
+  std::vector<Tensor<2,dim> > local_2nd_derivs (quadrature[0].size());
 
   hp::DoFHandler<dim>::active_cell_iterator
     cell = dof_handler.begin_active (),
@@ -417,7 +417,7 @@ void MinimizationProblem<1>::refine_grid ()
       fe_values.get_present_fe_values().get_function_2nd_derivatives (present_solution, local_2nd_derivs);
 
       double cell_residual_norm = 0;
-      for (unsigned int q=0; q<quadrature[0].n_quadrature_points; ++q)
+      for (unsigned int q=0; q<quadrature[0].size(); ++q)
         {
           const double x             = fe_values.get_present_fe_values().quadrature_point(q)[0];
           const double u             = local_values[q];
@@ -520,7 +520,7 @@ MinimizationProblem<dim>::energy (const hp::DoFHandler<dim> &dof_handler,
                           update_values   | update_gradients |
                            update_q_points | update_JxW_values);
 
-  const unsigned int   n_q_points    = quadrature_formula[0].n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula[0].size();
 
   std::vector<double>         local_solution_values (n_q_points);
   std::vector<Tensor<1,dim> > local_solution_grads (n_q_points);
index e1b51ed79b9ae51fa9ad2a3c88f5a60e06079ce3..18c76377258d397156593edd53ec75ecb5b8cb6b 100644 (file)
@@ -2,7 +2,7 @@
 //    rt_distorted_01.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2003, 2005, 2006, 2007 by the deal.II authors
+//    Copyright (C) 2003, 2005, 2006, 2007, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -273,7 +273,7 @@ double EvaluateArea (Mapping<2> &mapping,
                                     update_q_points  |
                                     update_JxW_values));
 
-  const unsigned int   n_q_points    = quad.n_quadrature_points;
+  const unsigned int   n_q_points    = quad.size();
   const unsigned int   n_components   = dof_handler->get_fe().n_components();
 
                                   // Cell iterators
index 47c1caa2e5bf346f4b28c5bf4346cc510dd0b283..00a256d0e89bb3893c795d5430a84dd7318c1e3c 100644 (file)
@@ -287,7 +287,7 @@ double EvaluateDiver (Mapping<2> &mapping,
                                     update_JxW_values|
                                     update_contravariant_transformation));
 
-  const unsigned int   n_q_points    = quad.n_quadrature_points;
+  const unsigned int   n_q_points    = quad.size();
   const unsigned int   n_components   = dof_handler.get_fe().n_components();
 
                                   // Cell iterators
@@ -347,7 +347,7 @@ int main ()
   sprintf (buf,"DoFs per Quad: %i per Line %i per Vert %i\n", fe.dofs_per_quad, fe.dofs_per_line, fe.dofs_per_vertex);
   deallog << buf;
   
-  sprintf (buf,"QPoints %i\n", quad_temp.n_quadrature_points);
+  sprintf (buf,"QPoints %i\n", quad_temp.size());
   deallog << buf;
 
                                   // Create an deformation object for the Eulerian mapping
index 4f183ff433ff76c69e5136ff43551b2a34e40546..279c3e62e750e381034cc2fa2a6154e383305c7f 100644 (file)
@@ -2,7 +2,7 @@
 //    abf_01.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2003, 2005, 2006, 2008 by the deal.II authors
+//    Copyright (C) 2003, 2005, 2006, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -73,7 +73,7 @@ void EvaluateDerivative (DoFHandler<2> *dof_handler,
                                       update_gradients |
                                       update_JxW_values));
 
-    const unsigned int   n_q_points    = quad.n_quadrature_points;
+    const unsigned int   n_q_points    = quad.size();
     const unsigned int   n_components   = dof_handler->get_fe().n_components();
     const unsigned int dofs_per_cell = dof_handler->get_fe().dofs_per_cell;
 
index e062974ef8c21158025f80605be5551601734282..8c5938b924aaa6ad0431c4d467156a83a546d7eb 100644 (file)
@@ -41,7 +41,7 @@ void test (const FiniteElement<dim> &fe,
   deallog << fe.get_name() << ' ' << fe.dofs_per_cell << ' ';
   
   for (unsigned int i=0; i<fe.dofs_per_cell; ++i)
-    for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+    for (unsigned int q=0; q<quadrature.size(); ++q)
       for (unsigned int c=0; c<fe.n_components(); ++c)
         {
           const Point<dim> point = quadrature.point(q);
index 149947ffb0a73617a0a128d5df7828dfae11a142..929ac99058cb63ebdeaaaa66f3dfa6a9a5fbbecd 100644 (file)
@@ -117,7 +117,7 @@ void EvaluateDerivative (DoFHandler<2> *dof_handler,
                                     update_gradients |
                                     update_JxW_values));
 
-  const unsigned int   n_q_points    = quad.n_quadrature_points;
+  const unsigned int   n_q_points    = quad.size();
   const unsigned int   n_components   = dof_handler->get_fe().n_components();
   const unsigned int dofs_per_cell = dof_handler->get_fe().dofs_per_cell;
 
index 5237ca2594b005efc492256bdebf8425121fe6e6..32851350cbb735b34bbecb7948cbb623126d129d 100644 (file)
@@ -43,7 +43,7 @@ plot_derivatives(Mapping<dim>& mapping,
     {
       fe.reinit(c, face);
       
-      for (unsigned int k=0; k<q.n_quadrature_points;++k)
+      for (unsigned int k=0; k<q.size();++k)
        {
          deallog << "Face " << face
                  << " Point " << q.point(k) << std::endl;
index 96db890febc8f8f425eea1c45f6a8d83c2d8c4df..ece08bd5da484ded9949a3b10a33481acb2f12e6 100644 (file)
@@ -48,7 +48,7 @@ void vector_values(const FiniteElement<dim>& fe)
     v(i) = i;
   
   FEValues<dim> feval(fe, quadrature, update_values);
-  std::vector<Vector<double> > local(quadrature.n_quadrature_points,
+  std::vector<Vector<double> > local(quadrature.size(),
                                     Vector<double>(fe.n_components()));
   
   typename DoFHandler<dim>::active_cell_iterator cell = dof.begin_active();
@@ -63,7 +63,7 @@ void vector_values(const FiniteElement<dim>& fe)
       for (unsigned int c=0;c<fe.n_components();++c)
        {
          deallog << "Component " << c;
-         for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
+         for (unsigned int k=0;k<quadrature.size();++k)
            deallog << '\t' << (int) local[k](c);
          deallog << std::endl;
        }
index 2acd5293c57713e49ff2376b3eaed378953a8706..ecc0fd43fe875b112930944c5fcbc9d2599928bf 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -32,10 +32,10 @@ double difference(
   double result = 0.;
   QGauss<dim> quadrature(fe.degree+1);
   
-  std::vector<double> f(quadrature.n_quadrature_points);
+  std::vector<double> f(quadrature.size());
   function.value_list(quadrature.get_points(), f);
   
-  for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
+  for (unsigned int k=0;k<quadrature.size();++k)
     {
       double diff = f[k];
       for (unsigned int i=0;i<dofs.size();++i)
@@ -56,12 +56,12 @@ double vector_difference(
   double result = 0.;
   QGauss<dim> quadrature(fe.degree+1);
   
-  std::vector<Vector<double> > f(quadrature.n_quadrature_points,
+  std::vector<Vector<double> > f(quadrature.size(),
                                 Vector<double>(function.n_components));
 
   function.vector_value_list(quadrature.get_points(), f);
   
-  for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
+  for (unsigned int k=0;k<quadrature.size();++k)
     for (unsigned int comp=0;comp<fe.n_components();++comp)
       {
        double diff = f[k](comp+offset);
index 6e8e30ecae9105fadc6c0bc7962601ea188b7b49..c0d6b17d557811da5e5957b896f3a5819bf0cffb 100644 (file)
@@ -79,7 +79,7 @@ plot_faces(Mapping<dim> &mapping,
   deallog.push(name);
 
   QGauss<dim-1> q(4);
-  const unsigned int nq = (unsigned int) (.01 + std::pow(q.n_quadrature_points, 1./(dim-1)));
+  const unsigned int nq = (unsigned int) (.01 + std::pow(q.size(), 1./(dim-1)));
 
   FEFaceValues<dim> fe_values(mapping, fe, q,
                              UpdateFlags(update_q_points
@@ -123,7 +123,7 @@ plot_subfaces(Mapping<dim> &mapping,
   deallog.push(name);
 
   QGauss<dim-1> q(4);
-  const unsigned int nq = (unsigned int) (.01 + std::pow(q.n_quadrature_points, 1./(dim-1)));
+  const unsigned int nq = (unsigned int) (.01 + std::pow(q.size(), 1./(dim-1)));
 
   FESubfaceValues<dim> fe_values(mapping, fe, q,
                                 UpdateFlags(update_q_points
index a78f14cd40282ebfc3032679b733bbcade185308..2755cbf7baab9983b90f266bf15bd328b411372b 100644 (file)
@@ -57,7 +57,7 @@ show_values(FiniteElement<dim>& fe,
   typename DoFHandler<dim>::cell_iterator c = dof.begin();
   fe_values.reinit(c);
 
-  for (unsigned int k=0; k<quadrature_formula.n_quadrature_points; ++k)
+  for (unsigned int k=0; k<quadrature_formula.size(); ++k)
     {
       deallog << quadrature_formula.point(k) << std::endl;
       deallog << "JxW: " << fe_values.JxW(k) << std::endl;
index c9b318988823e2a12d2523bd10771e527f0b24fa..1c55813193ebd7c01aa55d0775a5fe9dce51ff0d 100644 (file)
@@ -2,7 +2,7 @@
 //    mapping_q_eulerian.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2008 by the deal.II authors and Joshua White
+//    Copyright (C) 2008, 2010 by the deal.II authors and Joshua White
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -125,7 +125,7 @@ double MappingTest<dim>::compute_area ()
   FEValues<dim> fe_values (mapping, fe, quadrature_formula, 
                           update_JxW_values);
 
-  const unsigned int   n_q_points = quadrature_formula.n_quadrature_points;
+  const unsigned int   n_q_points = quadrature_formula.size();
 
   long double area = 0.;
 
index ae5c8ee26bb535a20ace475a44716595fade3599..c8580ce85aac5789ee4ce18d00f800472445d88c 100644 (file)
@@ -91,9 +91,9 @@ plot (const Triangulation<dim> &tr, const unsigned int p)
                                    // function on each cell and on
                                    // each quadrature point
   QTrapez<dim-1> quadrature;
-  std::vector<Vector<double> > shape_values1 (quadrature.n_quadrature_points,
+  std::vector<Vector<double> > shape_values1 (quadrature.size(),
                                              Vector<double>(dim));
-  std::vector<Vector<double> > shape_values2 (quadrature.n_quadrature_points,
+  std::vector<Vector<double> > shape_values2 (quadrature.size(),
                                              Vector<double>(dim));
   
   MappingCartesian<dim> mapping;
index d386fedde02f22cf6a58ac8b515360b2f97d2d75..7da6c1fb0744d06e918bd9b8af0c26ac8d5eb2c1 100644 (file)
@@ -69,7 +69,7 @@ check (const unsigned int p)
                                    // quadrature points of all cells
                                    // of the finer grid
   QTrapez<dim> quadrature;
-  std::vector<Vector<double> > shape_values (quadrature.n_quadrature_points,
+  std::vector<Vector<double> > shape_values (quadrature.size(),
                                              Vector<double>(dim));
   FEValues<dim> fe(fe_ned, quadrature,
                    update_values|update_q_points);
@@ -82,7 +82,7 @@ check (const unsigned int p)
       fe.reinit(c);
       fe.get_function_values (values, shape_values);
 
-      for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+      for (unsigned int q=0; q<quadrature.size(); ++q)
         {
           deallog << ", xq=" << fe.quadrature_point(q)
                   << ", f=[";
index daeba60f485a7938a47bdf391dc1228e92cf1538..8c3147763ebaadc773667dd3b5eff84efdf360cc 100644 (file)
@@ -2,7 +2,7 @@
 //    non_primitive_1.cc,v 1.6 2003/04/21 07:31:15 guido Exp
 //    Version: 
 //
-//    Copyright (C) 2001, 2002, 2003, 2005 by the deal.II authors
+//    Copyright (C) 2001, 2002, 2003, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -58,7 +58,7 @@ create_stokes_matrix_1 (const DoFHandler<dim> &dof_handler,
   FullMatrix<double> local_matrix (dofs_per_cell, dofs_per_cell);
 
   QGauss<dim> quadrature (3);
-  const unsigned int n_q_points = quadrature.n_quadrature_points;
+  const unsigned int n_q_points = quadrature.size();
   
   FEValues<dim> fe_values (fe, quadrature,
                           update_values | update_gradients |
@@ -131,7 +131,7 @@ create_stokes_matrix_2 (const DoFHandler<dim> &dof_handler,
   FullMatrix<double> local_matrix (dofs_per_cell, dofs_per_cell);
 
   QGauss<dim> quadrature (3);
-  const unsigned int n_q_points = quadrature.n_quadrature_points;
+  const unsigned int n_q_points = quadrature.size();
   
   FEValues<dim> fe_values (fe, quadrature,
                           update_values | update_gradients |
@@ -207,7 +207,7 @@ create_stokes_matrix_3 (const DoFHandler<dim> &dof_handler,
   FullMatrix<double> local_matrix (dofs_per_cell, dofs_per_cell);
 
   QGauss<dim> quadrature (3);
-  const unsigned int n_q_points = quadrature.n_quadrature_points;
+  const unsigned int n_q_points = quadrature.size();
   
   FEValues<dim> fe_values (fe, quadrature,
                           update_values | update_gradients |
index b32325b7f299d73d062856a426fb5b1d22d29795..652614d0d47901aa31e7331463385692999e48bd 100644 (file)
@@ -2,7 +2,7 @@
 //    non_primitive_2.cc,v 1.5 2003/04/21 07:31:15 guido Exp
 //    Version: 
 //
-//    Copyright (C) 2001, 2002, 2003, 2005, 2007 by the deal.II authors
+//    Copyright (C) 2001, 2002, 2003, 2005, 2007, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -61,7 +61,7 @@ create_stokes_matrix_2 (const DoFHandler<dim> &dof_handler,
   FullMatrix<double> local_matrix (dofs_per_cell, dofs_per_cell);
 
   QGauss<dim> quadrature (3);
-  const unsigned int n_q_points = quadrature.n_quadrature_points;
+  const unsigned int n_q_points = quadrature.size();
   
   FEValues<dim> fe_values (fe, quadrature,
                           update_values | update_gradients |
@@ -137,7 +137,7 @@ create_stokes_matrix_3 (const DoFHandler<dim> &dof_handler,
   FullMatrix<double> local_matrix (dofs_per_cell, dofs_per_cell);
 
   QGauss<dim> quadrature (3);
-  const unsigned int n_q_points = quadrature.n_quadrature_points;
+  const unsigned int n_q_points = quadrature.size();
   
   FEValues<dim> fe_values (fe, quadrature,
                           update_values | update_gradients |
index 25ae558d523b8f5f9abc541c6701245182e8408e..884428c7166c61b35337de89d97758de363c5be6 100644 (file)
@@ -2,7 +2,7 @@
 //    rt_10.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2003, 2005, 2006 by the deal.II authors
+//    Copyright (C) 2003, 2005, 2006, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -64,7 +64,7 @@ test (const unsigned int degree)
 
       FEValues<dim> fe_values (fe_rt, quadrature, update_values);
       fe_values.reinit (dof.begin_active());
-      for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+      for (unsigned int q=0; q<quadrature.size(); ++q)
        {
          deallog << "    Quadrature point " << q << ": ";
          for (unsigned int i=0; i<fe_rt.dofs_per_cell; ++i)
index 8a9efd8416faf358800ecc2723607f3190586d9d..1debc7ebc0d4f06831cdd6f7a9350b60eebce56c 100644 (file)
@@ -2,7 +2,7 @@
 //    rt_11.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2003, 2005, 2006 by the deal.II authors
+//    Copyright (C) 2003, 2005, 2006, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -70,7 +70,7 @@ test (const unsigned int degree)
 
       FEFaceValues<dim> fe_values (fe_rt, quadrature, update_values);
       fe_values.reinit (dof.begin_active(), 0);
-      for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+      for (unsigned int q=0; q<quadrature.size(); ++q)
        {
          deallog << "    Quadrature point " << q << ": ";
          for (unsigned int i=0; i<fe_rt.dofs_per_cell; ++i)
index c68f33de4de868f161c236b361f28d3a0a8fa229..b27a5df1c286ea4222cbce8754898fd49b8213eb 100644 (file)
@@ -2,7 +2,7 @@
 //    rt_12.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2003, 2005, 2006 by the deal.II authors
+//    Copyright (C) 2003, 2005, 2006, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -63,7 +63,7 @@ test (const unsigned int degree)
 
       FEValues<dim> fe_values (fe_rt, quadrature, update_gradients);
       fe_values.reinit (dof.begin_active());
-      for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+      for (unsigned int q=0; q<quadrature.size(); ++q)
        {
          deallog << "    Quadrature point " << q << ": ";
          for (unsigned int i=0; i<fe_rt.dofs_per_cell; ++i)
index 736cf6cd26406923db5a5db9af9e22461c147f07..e8e2417038a314591d3cb06e683c8956187bf19e 100644 (file)
@@ -2,7 +2,7 @@
 //    rt_13.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2003, 2005, 2006 by the deal.II authors
+//    Copyright (C) 2003, 2005, 2006, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -68,7 +68,7 @@ test (const unsigned int degree)
 
       FEFaceValues<dim> fe_values (fe_rt, quadrature, update_gradients);
       fe_values.reinit (dof.begin_active(), 0);
-      for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+      for (unsigned int q=0; q<quadrature.size(); ++q)
        {
          deallog << "    Quadrature point " << q << ": ";
          for (unsigned int i=0; i<fe_rt.dofs_per_cell; ++i)
index e22cf651e27a112c81d9ba270dadb6093a0c15e8..83dcb11cee79407ebc9d77c400f746f52d95314c 100644 (file)
@@ -2,7 +2,7 @@
 //    rt_14.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2003, 2005, 2006 by the deal.II authors
+//    Copyright (C) 2003, 2005, 2006, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -68,7 +68,7 @@ test (const unsigned int degree)
 
       FESubfaceValues<dim> fe_values (fe_rt, quadrature, update_values);
       fe_values.reinit (dof.begin_active(), 0, 0);
-      for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+      for (unsigned int q=0; q<quadrature.size(); ++q)
        {
          deallog << "    Quadrature point " << q << ": ";
          for (unsigned int i=0; i<fe_rt.dofs_per_cell; ++i)
index 38723fd3f446a56aa065ac3f32d8a97758e2a099..c2e42ed6f3e76ee86a47f1e4c1ee554509080dba 100644 (file)
@@ -2,7 +2,7 @@
 //    rt_15.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2003, 2005, 2006 by the deal.II authors
+//    Copyright (C) 2003, 2005, 2006, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -68,7 +68,7 @@ test (const unsigned int degree)
 
       FESubfaceValues<dim> fe_values (fe_rt, quadrature, update_gradients);
       fe_values.reinit (dof.begin_active(), 0, 0);
-      for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+      for (unsigned int q=0; q<quadrature.size(); ++q)
        {
          deallog << "    Quadrature point " << q << ": ";
          for (unsigned int i=0; i<fe_rt.dofs_per_cell; ++i)
index 9af0a16d874dbd1abb5c80a3a53b8b4e7f3bafca..ba987cdae3df5df7efede2a92c778812f5ed9e09 100644 (file)
@@ -2,7 +2,7 @@
 //    rt_7.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2003, 2005 by the deal.II authors
+//    Copyright (C) 2003, 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -53,9 +53,9 @@ plot_shape_functions(const unsigned int degree)
 
   Assert (fe.get_fe().n_components() == dim, ExcInternalError());
   
-  for (unsigned int q_point=0; q_point<q.n_quadrature_points; ++q_point)
+  for (unsigned int q_point=0; q_point<q.size(); ++q_point)
     {
-      if (q_point % QIterated<1>(q_trapez,div).n_quadrature_points == 0)
+      if (q_point % QIterated<1>(q_trapez,div).size() == 0)
         deallog << std::endl;
       
       deallog << fe.quadrature_point(q_point) << " ";
index 1984b7a97d61d512de83dd2ffb0ed60546248eaa..2d1de513e7447494073344a45d2f3786948f36c3 100644 (file)
@@ -62,7 +62,7 @@ test (const unsigned int degree)
   Assert (fe.get_fe().n_components() == dim, ExcInternalError());
 
 
-  for (unsigned int q_point=0; q_point<q.n_quadrature_points; ++q_point)
+  for (unsigned int q_point=0; q_point<q.size(); ++q_point)
     for (unsigned int i=0; i<dofs_per_cell; ++i)
       for (unsigned int j=0; j<dofs_per_cell; ++j)
         for (unsigned int d=0; d<dim; ++d)
index 8e5de4982c46433d6d28b5f951919b4eae24a6bf..5dc83cd12fb4f69ac5a7345827fafbcc138cf13c 100644 (file)
@@ -372,7 +372,7 @@ double TestProjection (Mapping<2> &mapping,
                                         update_JxW_values|
                                         update_contravariant_transformation));
 
-      const unsigned int   n_q_points    = quad.n_quadrature_points;
+      const unsigned int   n_q_points    = quad.size();
       const unsigned int   n_components   = dof_handler->get_fe().n_components();
       
                                       // Cell iterators
@@ -459,7 +459,7 @@ int main ()
          << ", dofs per line: " << fe.dofs_per_line
          << ", dofs per vertex: " <<  fe.dofs_per_vertex
          << std::endl;
-  deallog << "n_q_points=" << quad_temp.n_quadrature_points
+  deallog << "n_q_points=" << quad_temp.size()
          << std::endl;
 
                                   // Create an deformation object for
index c7a575c3ca0b2dffc0810d6eb8c8669ab6cf27df..35f1e6284cad6b4bba2fbc7ea03cefb9748eacd1 100644 (file)
@@ -2,7 +2,7 @@
 //    rt_normal_02.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2003, 2005, 2006 by the deal.II authors
+//    Copyright (C) 2003, 2005, 2006, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -80,8 +80,8 @@ void EvaluateNormal2 (DoFHandler<2> *dof_handler,
                                    update_gradients |
                                    update_JxW_values));
 
-    const unsigned int   n_q_face    = quad.n_quadrature_points;
-    const unsigned int   n_q_proj    = qproject.n_quadrature_points;
+    const unsigned int   n_q_face    = quad.size();
+    const unsigned int   n_q_proj    = qproject.size();
     const unsigned int   n_components   = dof_handler->get_fe().n_components();
     const unsigned int   dofs_per_cell = dof_handler->get_fe().dofs_per_cell;
 
@@ -110,7 +110,7 @@ void EvaluateNormal2 (DoFHandler<2> *dof_handler,
                             cell->face_orientation(f),
                             cell->face_flip(f),
                             cell->face_rotation(f),
-                            quad.n_quadrature_points));
+                            quad.size()));
                fe_v_face.reinit (cell, f);
                
                DoFHandler<2>::active_cell_iterator cell_n = cell->neighbor (f);
@@ -123,7 +123,7 @@ void EvaluateNormal2 (DoFHandler<2> *dof_handler,
                             cell_n->face_orientation(neighbor),
                             cell_n->face_flip(neighbor),
                             cell_n->face_rotation(neighbor),
-                            quad.n_quadrature_points));
+                            quad.size()));
 
                // Get values from solution vector (For Trap.Rule)
                std::vector<Vector<double> > this_value
@@ -191,7 +191,7 @@ void EvaluateNormal (DoFHandler<2> *dof_handler,
                                             update_normal_vectors |
                                             update_JxW_values));
 
-    const unsigned int   n_q_face    = quad.n_quadrature_points;
+    const unsigned int   n_q_face    = quad.size();
     const unsigned int   n_components   = dof_handler->get_fe().n_components();
     const unsigned int   dofs_per_cell = dof_handler->get_fe().dofs_per_cell;
 
index 6fe163256769d8a862da95169cc3e4577844c064..d9b8f38f02e155b3a63511b0d70d649ae6ec56db 100644 (file)
@@ -2,7 +2,7 @@
 //    rt_1.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2005 by the deal.II authors
+//    Copyright (C) 2005, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -86,11 +86,11 @@ compare_shapes (const FiniteElement<dim>& other,
                FullMatrix<double>& M)
 {
   QGauss<dim> quadrature(other.degree+1);
-  Table<3,double> other_values(quadrature.n_quadrature_points, other.dofs_per_cell, dim);
-  Table<3,double> nodes_values(quadrature.n_quadrature_points, other.dofs_per_cell, dim);
-  Table<3,Tensor<1,dim> > other_grads(quadrature.n_quadrature_points, other.dofs_per_cell, dim);
-  Table<3,Tensor<1,dim> > nodes_grads(quadrature.n_quadrature_points, other.dofs_per_cell, dim);
-  for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
+  Table<3,double> other_values(quadrature.size(), other.dofs_per_cell, dim);
+  Table<3,double> nodes_values(quadrature.size(), other.dofs_per_cell, dim);
+  Table<3,Tensor<1,dim> > other_grads(quadrature.size(), other.dofs_per_cell, dim);
+  Table<3,Tensor<1,dim> > nodes_grads(quadrature.size(), other.dofs_per_cell, dim);
+  for (unsigned int k=0;k<quadrature.size();++k)
     for (unsigned int i=0;i<other.dofs_per_cell;++i)
       for (unsigned int d=0;d<dim;++d)
        {
@@ -100,7 +100,7 @@ compare_shapes (const FiniteElement<dim>& other,
          nodes_grads[k][i][d] = nodes.shape_grad_component(i,quadrature.point(k),d);
        }
 
-  for (unsigned int k=0;k<quadrature.n_quadrature_points;++k)
+  for (unsigned int k=0;k<quadrature.size();++k)
     {
       for (unsigned int i=0;i<other.dofs_per_cell;++i)
        for (unsigned int d=0;d<dim;++d)
index a5a5e6a60bae6fa5f05c268fab40af91aa0b8614..7dbcb698e6af3f8a14fd4cf107d1da2f5b22defa 100644 (file)
@@ -2,7 +2,7 @@
 //    rtn_2.cc,v 1.3 2003/06/09 16:00:38 wolf Exp
 //    Version: 
 //
-//    Copyright (C) 2003, 2005, 2006, 2008 by the deal.II authors
+//    Copyright (C) 2003, 2005, 2006, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -52,9 +52,9 @@ plot_shape_functions(const unsigned int degree)
   FEValues<dim> fe(fe_rt, q, update_values|update_gradients|update_q_points);
   fe.reinit(c);
 
-  for (unsigned int q_point=0; q_point<q.n_quadrature_points; ++q_point)
+  for (unsigned int q_point=0; q_point<q.size(); ++q_point)
     {
-      if (q_point % QIterated<1>(q_trapez,div).n_quadrature_points == 0)
+      if (q_point % QIterated<1>(q_trapez,div).size() == 0)
         deallog << std::endl;
       
       deallog << fe.quadrature_point(q_point) << " ";
index f4390738479e9f2951c211374f309aef4d6f3892..a787a495833fe5b237aff1906e325613a3698044 100644 (file)
@@ -259,7 +259,7 @@ check_values_and_derivatives (const FiniteElement<dim> &fe,
                               const Quadrature<dim>    &q)
 {
                                    // check values
-  for (unsigned int x=0; x<q.n_quadrature_points; ++x)
+  for (unsigned int x=0; x<q.size(); ++x)
     for (unsigned int i=0; i<fe.dofs_per_cell; ++i)
       {
         for (unsigned int c=0; c<fe.n_components(); ++c)
@@ -288,7 +288,7 @@ check_values_and_derivatives (const FiniteElement<dim> &fe,
       };
         
                                    // check gradients
-  for (unsigned int x=0; x<q.n_quadrature_points; ++x)
+  for (unsigned int x=0; x<q.size(); ++x)
     for (unsigned int i=0; i<fe.dofs_per_cell; ++i)
       {
         for (unsigned int c=0; c<fe.n_components(); ++c)
@@ -310,7 +310,7 @@ check_values_and_derivatives (const FiniteElement<dim> &fe,
 
                                    // check second derivatives
   double max_diff=0.;
-  for (unsigned int x=0; x<q.n_quadrature_points; ++x)
+  for (unsigned int x=0; x<q.size(); ++x)
     for (unsigned int i=0; i<fe.dofs_per_cell; ++i)
       {
         for (unsigned int c=0; c<fe.n_components(); ++c)
index 2725693b6ab6484ffa3add8893931337404830c8..06485b6e19384cbd47cf89e9b2a77728f5af4828 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2006, 2007, 2008, 2009 by the deal.II authors
+//    Copyright (C) 2006, 2007, 2008, 2009, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -371,7 +371,7 @@ estimate_smoothness (Vector<float> &smoothness_indicators) const
        for (unsigned int i=0; i<fe_collection[fe].dofs_per_cell; ++i)
          {
            std::complex<double> sum = 0;
-           for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+           for (unsigned int q=0; q<quadrature.size(); ++q)
              {
                const Point<dim> x_q = quadrature.point(q);
                sum += std::exp(std::complex<double>(0,1) *
index 0d98839d82a67f9fd1a0d1c525ea8262fa388d45..5adb2ce3e476f19ef597629536981db8fd95517c 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2006, 2007, 2008, 2009 by the deal.II authors
+//    Copyright (C) 2006, 2007, 2008, 2009, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -373,7 +373,7 @@ estimate_smoothness (Vector<float> &smoothness_indicators) const
        for (unsigned int i=0; i<fe_collection[fe].dofs_per_cell; ++i)
          {
            std::complex<double> sum = 0;
-           for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+           for (unsigned int q=0; q<quadrature.size(); ++q)
              {
                const Point<dim> x_q = quadrature.point(q);
                sum += std::exp(std::complex<double>(0,1) *
index 86040c81503186c32264aeb556a00ebcca5344af..5d72ecb3db92350ec321ecd944e4cf2ea0057902 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -376,7 +376,7 @@ estimate_smoothness (Vector<float> &smoothness_indicators) const
        for (unsigned int i=0; i<fe_collection[fe].dofs_per_cell; ++i)
          {
            std::complex<double> sum = 0;
-           for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+           for (unsigned int q=0; q<quadrature.size(); ++q)
              {
                const Point<dim> x_q = quadrature.point(q);
                sum += std::exp(std::complex<double>(0,1) *
index db902cd0e7f11a69c9b6e21d92e97ef15b9a67ae..66c53a2ed1792990a1ed775ec4da26c13c03e662 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -379,7 +379,7 @@ estimate_smoothness (Vector<float> &smoothness_indicators) const
        for (unsigned int i=0; i<fe_collection[fe].dofs_per_cell; ++i)
          {
            std::complex<double> sum = 0;
-           for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+           for (unsigned int q=0; q<quadrature.size(); ++q)
              {
                const Point<dim> x_q = quadrature.point(q);
                sum += std::exp(std::complex<double>(0,1) *
index a5c5915eb633743de3fc8f024c232ed4a3f466d4..8b3522a83909c42aac139cb1b0ab34ccb3a2a79d 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -402,7 +402,7 @@ namespace LaplaceSolver
                             update_gradients | update_JxW_values);
 
     const unsigned int   dofs_per_cell = (*fe)[0].dofs_per_cell;
-    const unsigned int   n_q_points    = (*quadrature)[0].n_quadrature_points;
+    const unsigned int   n_q_points    = (*quadrature)[0].size();
 
     FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
 
@@ -525,7 +525,7 @@ namespace LaplaceSolver
                              update_JxW_values);
 
     const unsigned int   dofs_per_cell = (*this->fe)[0].dofs_per_cell;
-    const unsigned int   n_q_points    = (*this->quadrature)[0].n_quadrature_points;
+    const unsigned int   n_q_points    = (*this->quadrature)[0].size();
 
     Vector<double>       cell_rhs (dofs_per_cell);
     std::vector<double>  rhs_values (n_q_points);
index b254ac30b7400d94256e7c1492f9b7014525e533..b99eb1293c9e05fda1fd437fc58899009dbf56b0 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -120,7 +120,7 @@ void LaplaceProblem::assemble_system ()
                         update_values | update_gradients | update_JxW_values);
   
   const unsigned int   dofs_per_cell = fe[0].dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula[0].n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula[0].size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
index e742e9614b3392b16fa7b9237612acb5226d1625..ad6b53519e332349ce0cb1ac74bea2ebbb7e8107 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -154,7 +154,7 @@ void LaplaceProblem::assemble_system ()
                               update_values | update_gradients | update_JxW_values);
   
   const unsigned int   max_dofs_per_cell = fe.max_dofs_per_cell ();
-  const unsigned int   n_q_points    = quadrature_formula[0].n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula[0].size();
 
   FullMatrix<double>   cell_matrix (max_dofs_per_cell, max_dofs_per_cell);
   Vector<double>       cell_rhs (max_dofs_per_cell);
index c28e052f1aac9d9fbef8f3f1b2c2f12f4c4cf16f..d1e9ff3005ae3e330a0d5b3922a1981de51eaa85 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -151,7 +151,7 @@ void LaplaceProblem::assemble_system ()
                               update_values | update_gradients | update_JxW_values);
   
   const unsigned int   max_dofs_per_cell = fe.max_dofs_per_cell ();
-  const unsigned int   n_q_points    = quadrature_formula[0].n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula[0].size();
 
   FullMatrix<double>   cell_matrix (max_dofs_per_cell, max_dofs_per_cell);
   Vector<double>       cell_rhs (max_dofs_per_cell);
index 2783cf0222638f76e329f5fef0556a497fe7d6a2..196f8c749ad836db799023e33d823fdacd9e4703 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -173,7 +173,7 @@ void LaplaceProblem::assemble_system ()
       cell_rhs = 0;
 
       const unsigned int dofs_per_cell = cell->get_fe ().dofs_per_cell;
-      const unsigned int   n_q_points    = quadrature_formula[cell->active_fe_index ()].n_quadrature_points;
+      const unsigned int   n_q_points    = quadrature_formula[cell->active_fe_index ()].size();
       for (unsigned int i=0; i<dofs_per_cell; ++i)
        for (unsigned int j=0; j<dofs_per_cell; ++j)
          for (unsigned int q_point=0; q_point<n_q_points; ++q_point)
index e29d71d4274e6e4d2612e7749a807e833c35f5db..6a3e7399e30b278d85f595a7482046fcbf961753 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -177,7 +177,7 @@ void LaplaceProblem<dim>::assemble_system ()
                            update_q_points | update_JxW_values);
 
   const unsigned int   dofs_per_cell = fe[0].dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula[0].n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula[0].size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
index 9046856fbbdbdd6b6652b3d423c94f2bd9b2c984..a2213ac9d8dbed4f4ca0c728548607d571cf947b 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -173,7 +173,7 @@ void LaplaceProblem<dim>::assemble_system ()
                           update_q_points  |  update_JxW_values);
 
   const unsigned int   dofs_per_cell = fe[0].dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula[0].n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula[0].size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
index ead4eed476e251a7f815c521b0f3cdf6422b714f..bbc4fb09c6b9a10ab0dc27179eba87efbf42ec05 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -198,7 +198,7 @@ void LaplaceProblem<dim>::assemble_system ()
                           update_q_points  |  update_JxW_values);
 
   const unsigned int   dofs_per_cell = fe[0].dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula[0].n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula[0].size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
index 9085ab590edfdfdbc87294ba4503e07f4f6a08e7..f814069218058d768b08ccf33f3dc080cecb3170 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2007, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -265,8 +265,8 @@ void HelmholtzProblem<dim>::assemble_system ()
   hp::QCollection<dim>   quadrature_formula(QGauss<dim>(3));
   hp::QCollection<dim-1> face_quadrature_formula(QGauss<dim-1>(3));
 
-  const unsigned int n_q_points    = quadrature_formula[0].n_quadrature_points;
-  const unsigned int n_face_q_points = face_quadrature_formula[0].n_quadrature_points;
+  const unsigned int n_q_points    = quadrature_formula[0].size();
+  const unsigned int n_face_q_points = face_quadrature_formula[0].size();
 
   const unsigned int dofs_per_cell = (*fe)[0].dofs_per_cell;
 
index fc28e387e1aa41dedba9e80eaa2805c09a62d5cc..d615dc5643d4360b1d13bde80ab5f30f5cf2e33a 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2005, 2006, 2008 by the deal.II authors
+//    Copyright (C) 2005, 2006, 2008, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -204,7 +204,7 @@ void ElasticProblem<dim>::assemble_system ()
                            update_q_points | update_JxW_values);
 
   const unsigned int   dofs_per_cell = fe[0].dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula[0].n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula[0].size();
 
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);
index c217acda5297246a6608adaf22a9c536c37cbeac..fc5ce290a168bc4e2d420faf6991a1d925ac9bcf 100644 (file)
@@ -148,8 +148,10 @@ test_cochain(const Triangulation<dim>& tr, const FiniteElement<dim>& fe)
   
   MeshWorker::IntegrationInfoBox<dim> info_box;
   UpdateFlags update_flags = update_values | update_gradients;
-  info_box.add_update_flags(update_flags, true, true, true, true);
+  info_box.add_update_flags_cell(update_flags);
   info_box.initialize(fe, mapping, &dof.block_info());
+  if (debugging)
+    deallog << "Infobox ready" << std::endl;
   
   MeshWorker::DoFInfo<dim> dof_info(dof.block_info());
   
@@ -166,7 +168,8 @@ test_cochain(const Triangulation<dim>& tr, const FiniteElement<dim>& fe)
       for (unsigned int i=0;i<matrices.block(b).matrix.m();++i)
        if (matrices.block(b).matrix.diag_element(i) == 0.)
          matrices.block(b).matrix.diag_element(i) = 1.;
-
+  if (debugging)
+    deallog << "Matrices ready" << std::endl;
                                   // Set up vectors
   BlockVector<double> source(dof.block_info().global());
   BlockVector<double> result1(dof.block_info().global());
@@ -195,7 +198,10 @@ test_cochain(const Triangulation<dim>& tr, const FiniteElement<dim>& fe)
       
       matrices.matrix(m+2).vmult(result1.block(d+1), source.block(d));      
 
-      solver.solve(matrices.matrix(m+3), aux.block(d+1), result1.block(d+1), PreconditionIdentity());
+      PreconditionSSOR<SparseMatrix<double> > precondition;
+      precondition.initialize(matrices.matrix(m+3), 1.2);
+      
+      solver.solve(matrices.matrix(m+3), aux.block(d+1), result1.block(d+1), precondition);
       
       matrices.matrix(m+2).Tvmult(result1.block(d), aux.block(d+1));
       matrices.matrix(m+1).vmult(result2.block(d), source.block(d));
@@ -289,6 +295,7 @@ int main()
   const std::string logname = JobIdentifier::base_name(__FILE__) + std::string("/output");
   std::ofstream logfile(logname.c_str());
   deallog.attach(logfile);
+  deallog.log_execution_time(true);
   if (!debugging)
     {
       deallog.depth_console(0);
index 58aefa2e74e28de01bb62a5f7b7e9e1d4af72594..9e3b07d61508683542837273004fa129da5f76b4 100644 (file)
 
-DEAL:d2-p0::Setup
-DEAL:d2-p0::Triangulation hypercube 2D refinement 1 steps 4 active cells 5 total cells 
-DEAL:d2-p0::Form 0
-DEAL:d2-p0::Difference d*d 0
-DEAL:d2-p0::Form 1
-DEAL:d2-p0::d^2            0
-DEAL:d2-p0::d*d^2          0
-DEAL:d2-p0::Difference d*d 0
-DEAL:d2-p0::Triangulation hypercube 2D refinement 2 local  steps 7 active cells 9 total cells 
-DEAL:d2-p0::Form 0
-DEAL:d2-p0::Difference d*d 0
-DEAL:d2-p0::Form 1
-DEAL:d2-p0::d^2            0
-DEAL:d2-p0::d*d^2          0
-DEAL:d2-p0::Difference d*d 0
-DEAL:d2-p0::Triangulation hypercube 2D refinement 3 local  steps 25 active cells 33 total cells 
-DEAL:d2-p0::Form 0
-DEAL:d2-p0::Difference d*d 0
-DEAL:d2-p0::Form 1
-DEAL:d2-p0::d^2            0
-DEAL:d2-p0::d*d^2          0
-DEAL:d2-p0::Difference d*d 0
-DEAL:d2-p1::Setup
-DEAL:d2-p1::Triangulation hypercube 2D refinement 1 steps 4 active cells 5 total cells 
-DEAL:d2-p1::Form 0
-DEAL:d2-p1::Difference d*d 0
-DEAL:d2-p1::Form 1
-DEAL:d2-p1::d^2            0
-DEAL:d2-p1::d*d^2          0
-DEAL:d2-p1::Difference d*d 0
-DEAL:d2-p1::Triangulation hypercube 2D refinement 2 local  steps 7 active cells 9 total cells 
-DEAL:d2-p1::Form 0
-DEAL:d2-p1::Difference d*d 0
-DEAL:d2-p1::Form 1
-DEAL:d2-p1::d^2            0
-DEAL:d2-p1::d*d^2          0
-DEAL:d2-p1::Difference d*d 0
-DEAL:d2-p1::Triangulation hypercube 2D refinement 3 local  steps 25 active cells 33 total cells 
-DEAL:d2-p1::Form 0
-DEAL:d2-p1::Difference d*d 0
-DEAL:d2-p1::Form 1
-DEAL:d2-p1::d^2            0
-DEAL:d2-p1::d*d^2          0
-DEAL:d2-p1::Difference d*d 0
-DEAL:d2-p2::Setup
-DEAL:d2-p2::Triangulation hypercube 2D refinement 1 steps 4 active cells 5 total cells 
-DEAL:d2-p2::Form 0
-DEAL:d2-p2::Difference d*d 0
-DEAL:d2-p2::Form 1
-DEAL:d2-p2::d^2            0
-DEAL:d2-p2::d*d^2          0
-DEAL:d2-p2::Difference d*d 0
-DEAL:d2-p2::Triangulation hypercube 2D refinement 2 local  steps 7 active cells 9 total cells 
-DEAL:d2-p2::Form 0
-DEAL:d2-p2::Difference d*d 0
-DEAL:d2-p2::Form 1
-DEAL:d2-p2::d^2            0
-DEAL:d2-p2::d*d^2          0
-DEAL:d2-p2::Difference d*d 0
-DEAL:d2-p2::Triangulation hypercube 2D refinement 3 local  steps 25 active cells 33 total cells 
-DEAL:d2-p2::Form 0
-DEAL:d2-p2::Difference d*d 0
-DEAL:d2-p2::Form 1
-DEAL:d2-p2::d^2            0
-DEAL:d2-p2::d*d^2          0
-DEAL:d2-p2::Difference d*d 0
-DEAL:d3-p0::Setup
-DEAL:d3-p0::Triangulation hypercube 3D refinement 1 steps 8 active cells 9 total cells 
-DEAL:d3-p0::Form 0
-DEAL:d3-p0::Difference d*d 0
-DEAL:d3-p0::Form 1
-DEAL:d3-p0::d^2            0
-DEAL:d3-p0::d*d^2          0
-DEAL:d3-p0::Difference d*d 0
-DEAL:d3-p0::Form 2
-DEAL:d3-p0::d^2            0
-DEAL:d3-p0::d*d^2          0
-DEAL:d3-p0::Difference d*d 0
-DEAL:d3-p1::Setup
-DEAL:d3-p1::Triangulation hypercube 3D refinement 1 steps 8 active cells 9 total cells 
-DEAL:d3-p1::Form 0
-DEAL:d3-p1::Difference d*d 0
-DEAL:d3-p1::Form 1
-DEAL:d3-p1::d^2            0
-DEAL:d3-p1::d*d^2          0
-DEAL:d3-p1::Difference d*d 0
-DEAL:d3-p1::Form 2
-DEAL:d3-p1::d^2            0
-DEAL:d3-p1::d*d^2          0
-DEAL:d3-p1::Difference d*d 0
+0.0200000:DEAL:d2-p0::Setup
+0.0600000:DEAL:d2-p0::Triangulation hypercube 2D refinement 1 steps 4 active cells 5 total cells 
+0.460000:DEAL:d2-p0::Form 0
+0.480000:DEAL:d2-p0::Difference d*d 0
+0.480000:DEAL:d2-p0::Form 1
+0.480000:DEAL:d2-p0::d^2            0
+0.480000:DEAL:d2-p0::d*d^2          0
+0.480000:DEAL:d2-p0::Difference d*d 0
+0.480000:DEAL:d2-p0::Triangulation hypercube 2D refinement 2 local  steps 7 active cells 9 total cells 
+0.940000:DEAL:d2-p0::Form 0
+0.940000:DEAL:d2-p0::Difference d*d 0
+0.940000:DEAL:d2-p0::Form 1
+0.940000:DEAL:d2-p0::d^2            0
+0.940000:DEAL:d2-p0::d*d^2          0
+0.950000:DEAL:d2-p0::Difference d*d 0
+0.950000:DEAL:d2-p0::Triangulation hypercube 2D refinement 3 local  steps 25 active cells 33 total cells 
+1.37000:DEAL:d2-p0::Form 0
+1.41000:DEAL:d2-p0::Difference d*d 0
+1.41000:DEAL:d2-p0::Form 1
+1.41000:DEAL:d2-p0::d^2            0
+1.41000:DEAL:d2-p0::d*d^2          0
+1.41000:DEAL:d2-p0::Difference d*d 0
+1.41000:DEAL:d2-p1::Setup
+1.49000:DEAL:d2-p1::Triangulation hypercube 2D refinement 1 steps 4 active cells 5 total cells 
+1.95000:DEAL:d2-p1::Form 0
+1.95000:DEAL:d2-p1::Difference d*d 0
+1.95000:DEAL:d2-p1::Form 1
+1.95000:DEAL:d2-p1::d^2            0
+1.95000:DEAL:d2-p1::d*d^2          0
+1.95000:DEAL:d2-p1::Difference d*d 0
+1.95000:DEAL:d2-p1::Triangulation hypercube 2D refinement 2 local  steps 7 active cells 9 total cells 
+2.46000:DEAL:d2-p1::Form 0
+2.48000:DEAL:d2-p1::Difference d*d 0
+2.49000:DEAL:d2-p1::Form 1
+2.49000:DEAL:d2-p1::d^2            0
+2.49000:DEAL:d2-p1::d*d^2          0
+2.49000:DEAL:d2-p1::Difference d*d 0
+2.50000:DEAL:d2-p1::Triangulation hypercube 2D refinement 3 local  steps 25 active cells 33 total cells 
+3.07000:DEAL:d2-p1::Form 0
+3.12000:DEAL:d2-p1::Difference d*d 0
+3.12000:DEAL:d2-p1::Form 1
+3.12000:DEAL:d2-p1::d^2            0
+3.12000:DEAL:d2-p1::d*d^2          0
+3.12000:DEAL:d2-p1::Difference d*d 0
+3.12000:DEAL:d2-p2::Setup
+3.36000:DEAL:d2-p2::Triangulation hypercube 2D refinement 1 steps 4 active cells 5 total cells 
+4.30000:DEAL:d2-p2::Form 0
+4.30000:DEAL:d2-p2::Difference d*d 0
+4.31000:DEAL:d2-p2::Form 1
+4.31000:DEAL:d2-p2::d^2            0
+4.31000:DEAL:d2-p2::d*d^2          0
+4.31000:DEAL:d2-p2::Difference d*d 0
+4.31000:DEAL:d2-p2::Triangulation hypercube 2D refinement 2 local  steps 7 active cells 9 total cells 
+5.27000:DEAL:d2-p2::Form 0
+5.29000:DEAL:d2-p2::Difference d*d 0
+5.29000:DEAL:d2-p2::Form 1
+5.29000:DEAL:d2-p2::d^2            0
+5.29000:DEAL:d2-p2::d*d^2          0
+5.32000:DEAL:d2-p2::Difference d*d 0
+5.32000:DEAL:d2-p2::Triangulation hypercube 2D refinement 3 local  steps 25 active cells 33 total cells 
+6.77000:DEAL:d2-p2::Form 0
+6.90000:DEAL:d2-p2::Difference d*d 0
+6.90000:DEAL:d2-p2::Form 1
+6.90000:DEAL:d2-p2::d^2            0
+6.91000:DEAL:d2-p2::d*d^2          0
+6.95000:DEAL:d2-p2::Difference d*d 0
+6.95000:DEAL:d3-p0::Setup
+7.06000:DEAL:d3-p0::Triangulation hypercube 3D refinement 1 steps 8 active cells 9 total cells 
+8.44000:DEAL:d3-p0::Form 0
+8.44000:DEAL:d3-p0::Difference d*d 0
+8.44000:DEAL:d3-p0::Form 1
+8.45000:DEAL:d3-p0::d^2            0
+8.45000:DEAL:d3-p0::d*d^2          0
+8.46000:DEAL:d3-p0::Difference d*d 0
+8.46000:DEAL:d3-p0::Form 2
+8.46000:DEAL:d3-p0::d^2            0
+8.46000:DEAL:d3-p0::d*d^2          0
+8.46000:DEAL:d3-p0::Difference d*d 0
+8.47000:DEAL:d3-p0::Triangulation hypercube 3D refinement 2 local  steps 15 active cells 17 total cells 
+9.94000:DEAL:d3-p0::Form 0
+9.98000:DEAL:d3-p0::Difference d*d 0
+9.98000:DEAL:d3-p0::Form 1
+9.98000:DEAL:d3-p0::d^2            0
+9.98000:DEAL:d3-p0::d*d^2          0
+9.99000:DEAL:d3-p0::Difference d*d 0
+9.99000:DEAL:d3-p0::Form 2
+9.99000:DEAL:d3-p0::d^2            0
+9.99000:DEAL:d3-p0::d*d^2          0
+9.99000:DEAL:d3-p0::Difference d*d 0
+9.99000:DEAL:d3-p1::Setup
+16.8400:DEAL:d3-p1::Triangulation hypercube 3D refinement 1 steps 8 active cells 9 total cells 
+29.4100:DEAL:d3-p1::Form 0
+29.5100:DEAL:d3-p1::Difference d*d 0
+29.5100:DEAL:d3-p1::Form 1
+29.5100:DEAL:d3-p1::d^2            0
+29.5100:DEAL:d3-p1::d*d^2          0
+29.5500:DEAL:d3-p1::Difference d*d 0
+29.5500:DEAL:d3-p1::Form 2
+29.5500:DEAL:d3-p1::d^2            0
+29.5500:DEAL:d3-p1::d*d^2          0
+29.5500:DEAL:d3-p1::Difference d*d 0
+29.5700:DEAL:d3-p1::Triangulation hypercube 3D refinement 2 local  steps 15 active cells 17 total cells 
+47.7600:DEAL:d3-p1::Form 0
+48.3300:DEAL:d3-p1::Difference d*d 0
+48.3300:DEAL:d3-p1::Form 1
+48.3300:DEAL:d3-p1::d^2            0
+48.3300:DEAL:d3-p1::d*d^2          0
+48.4100:DEAL:d3-p1::Difference d*d 0
+48.4100:DEAL:d3-p1::Form 2
+48.4100:DEAL:d3-p1::d^2            0
+48.4200:DEAL:d3-p1::d*d^2          0
+48.4200:DEAL:d3-p1::Difference d*d 0
index 098c8f3c291135140f02f4c0f87fd8a6a9e1f5dd..2f986def67717c98b63b4d2d3671d09980519803 100644 (file)
@@ -2,7 +2,7 @@
 //    $Id$
 //    Version: $Name$ 
 //
-//    Copyright (C) 2002, 2003, 2004, 2005, 2006, 2009 by the deal.II authors
+//    Copyright (C) 2002, 2003, 2004, 2005, 2006, 2009, 2010 by the deal.II authors
 //
 //    This file is subject to QPL and may not be  distributed
 //    without copyright and license information. Please refer
@@ -49,7 +49,7 @@ void assemble_laplace (MATRIX &B, VECTOR &bb, DoFHandler<dim> &dof_handler, Fini
                           update_values | update_gradients | update_JxW_values);
   
   const unsigned int   dofs_per_cell = fe.dofs_per_cell;
-  const unsigned int   n_q_points    = quadrature_formula.n_quadrature_points;
+  const unsigned int   n_q_points    = quadrature_formula.size();
   
   FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
   Vector<double>       cell_rhs (dofs_per_cell);

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.