]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
fix step-41: n_quadrature_points got removed.
authorheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 16 Nov 2010 19:48:00 +0000 (19:48 +0000)
committerheister <heister@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 16 Nov 2010 19:48:00 +0000 (19:48 +0000)
git-svn-id: https://svn.dealii.org/trunk@22767 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-41/include/base.templates.h
deal.II/examples/step-41/include/camclay.templates.h

index 8664990d80df5cd606838d747a5443447cf549c9..5f6f4c63edbe529c4586f3f5eb381b54dab7661a 100644 (file)
@@ -236,9 +236,9 @@ void Base<dim>::write_stress_strain (const unsigned int &cc,
   QIterated<dim> qf(QTrapez<1>(),2);
   FEValues<dim> fe_v(vspace.get_fe(), qf,
                     update_gradients | update_quadrature_points);
-  std::vector< std::vector< Tensor<1,dim> > > total_grads(qf.n_quadrature_points,
+  std::vector< std::vector< Tensor<1,dim> > > total_grads(qf.size(),
                                                          std::vector<Tensor<1,dim> >(dim));
-  std::vector< Point<dim> > points(qf.n_quadrature_points);
+  std::vector< Point<dim> > points(qf.size());
 
   output << "SCALARS uyy double 1" << std::endl;
   output << "LOOKUP_TABLE default" << std::endl;
@@ -253,7 +253,7 @@ void Base<dim>::write_stress_strain (const unsigned int &cc,
     points = fe_v.get_quadrature_points();
 
     //loop over the quadrature points---------------------------------
-    for (unsigned int qp = 0; qp<qf.n_quadrature_points; ++qp) {
+    for (unsigned int qp = 0; qp<qf.size(); ++qp) {
 
       output << total_grads[qp][1][1] << " ";
 
@@ -300,12 +300,12 @@ void Base<dim>::write_plot_values(const unsigned int &cc,
   QTrapez<dim> qf;
   FEValues<dim> fe_v(vspace.get_fe(), qf,
                     update_values | update_gradients | update_quadrature_points);
-  std::vector< std::vector< Tensor<1,dim> > > elastic_grads(qf.n_quadrature_points,
+  std::vector< std::vector< Tensor<1,dim> > > elastic_grads(qf.size(),
                                                            std::vector<Tensor<1,dim> >(dim));
-  std::vector< std::vector< Tensor<1,dim> > > plastic_grads(qf.n_quadrature_points,
+  std::vector< std::vector< Tensor<1,dim> > > plastic_grads(qf.size(),
                                                            std::vector<Tensor<1,dim> >(dim));
-  std::vector< Vector<double> > hard_values(qf.n_quadrature_points, Vector<double>(dim));
-  std::vector< Point<dim> > points(qf.n_quadrature_points);
+  std::vector< Vector<double> > hard_values(qf.size(), Vector<double>(dim));
+  std::vector< Point<dim> > points(qf.size());
 
   //std::vector< std::vector<double> > strain(dim, std::vector<double>(dim) );
 
index da3c49385fc0e83eef338ae061a5abb382c3d50e..ca9c2bad9e41296dc2b707f9b3d7c9a111b5f9b1 100644 (file)
@@ -78,7 +78,7 @@ void CamClay<dim>::initial_conditions(VectorSpace<dim> &vspace)
                      update_values   | update_gradients |
                      update_quadrature_points | update_JxW_values);
   
-  const unsigned int n_qp_v = qf_v.n_quadrature_points;
+  const unsigned int n_qp_v = qf_v.size();
   
 
   typename MGDoFHandler<dim>::active_cell_iterator cell = vspace.get_dh().begin_active(),
@@ -114,7 +114,7 @@ void CamClay<dim>::update_internal_variables(VectorSpace<dim> &vspace,
                      update_values   | update_gradients |
                      update_quadrature_points | update_JxW_values);
   
-  const unsigned int n_qp_v = qf_v.n_quadrature_points;
+  const unsigned int n_qp_v = qf_v.size();
 
   std::vector<Point<dim> > points(n_qp_v);
   

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.