]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Clarify some things with the update flags for FEValues objects.
authorWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 20 May 1998 08:41:03 +0000 (08:41 +0000)
committerWolfgang Bangerth <bangerth@math.tamu.edu>
Wed, 20 May 1998 08:41:03 +0000 (08:41 +0000)
git-svn-id: https://svn.dealii.org/trunk@320 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/deal.II/Attic/examples/convergence/convergence.cc
deal.II/deal.II/Attic/examples/error-estimation/error-estimation.cc
deal.II/deal.II/Attic/examples/poisson/problem.cc
deal.II/deal.II/source/numerics/error_estimator.cc
deal.II/deal.II/source/numerics/matrices.cc
deal.II/deal.II/source/numerics/vectors.cc
tests/big-tests/convergence/convergence.cc
tests/big-tests/error-estimation/error-estimation.cc
tests/big-tests/poisson/problem.cc

index 7c09fcc3f59490bb45b5e1d38d496e98c89fe9e8..b2424130dc329343c4381a0e99bbf3d77fca5f59 100644 (file)
@@ -252,7 +252,7 @@ void PoissonProblem<dim>::run (const unsigned int level) {
 
   cout << "    Assembling matrices..." << endl;
   UpdateFlags update_flags = UpdateFlags(update_q_points  | update_gradients |
-                                        update_jacobians | update_JxW_values);
+                                        update_JxW_values);
   
   ProblemBase<dim>::FunctionMap dirichlet_bc;
   dirichlet_bc[0] = boundary_values;
index fc2e8078388928986a147650d6780fcc5ffc837d..063f9902c5e6ab59981b6dd09114ef245902f484 100644 (file)
@@ -438,7 +438,7 @@ void PoissonProblem<dim>::run (ParameterHandler &prm) {
 
       cout << "    Assembling matrices..." << endl;
       UpdateFlags update_flags = UpdateFlags(update_q_points  | update_gradients |
-                                            update_jacobians | update_JxW_values);
+                                            update_JxW_values);
   
       ProblemBase<dim>::FunctionMap dirichlet_bc;
       dirichlet_bc[0] = solution_function;
index cd95664c8703c99bacdc63850ae37a4ca167d63d..d8fb3c517fa590d01f760ad08924e61662dfa955 100644 (file)
@@ -459,7 +459,7 @@ void PoissonProblem<dim>::run (ParameterHandler &prm) {
   dirichlet_bc[0] = boundary_values;
   assemble (equation, quadrature, fe,
            UpdateFlags(update_gradients | update_JxW_values |
-                       update_jacobians | update_q_points),
+                       update_q_points),
            dirichlet_bc, *boundary);
 
   cout << "    Solving..." << endl;
index c5be6740d4bc8a7acf9bb7f2f5e74e4bb6026f11..a4cc86f98924d241432a46c6bb2f7190f3be5d79 100644 (file)
@@ -76,14 +76,9 @@ void KellyErrorEstimator<dim>::estimate_error (const DoFHandler<dim>    &dof,
   FEFaceValues<dim> fe_face_values_cell (fe, quadrature,
                                         UpdateFlags(update_gradients  |
                                                     update_JxW_values |
-                                                    update_jacobians  |
                                                     update_normal_vectors)); 
-  FEFaceValues<dim> fe_face_values_neighbor (fe, quadrature,
-                                            UpdateFlags(update_gradients |
-                                                        update_jacobians));
-  FESubfaceValues<dim> fe_subface_values (fe, quadrature,
-                                         UpdateFlags(update_gradients |
-                                                     update_jacobians));
+  FEFaceValues<dim> fe_face_values_neighbor (fe, quadrature, update_gradients);
+  FESubfaceValues<dim> fe_subface_values (fe, quadrature, update_gradients);
   
                                   // loop over all cells
   const active_cell_iterator endc = dof.end();
index 2ba113cd36b12d9641a327a43daf4e4c68b64e55..c2bd7354ece5036bca0a791b0a64b34a47c452be 100644 (file)
@@ -21,8 +21,7 @@ void MatrixCreator<dim>::create_mass_matrix (const DoFHandler<dim>    &dof,
                                             dSMatrix                 &matrix,
                                             const Function<dim> * const a) {
   dVector dummy;    // no entries, should give an error if accessed
-  UpdateFlags update_flags = UpdateFlags(update_jacobians |
-                                        update_JxW_values);
+  UpdateFlags update_flags = update_JxW_values;
   if (a != 0)
     update_flags = UpdateFlags (update_flags | update_q_points);
   const AssemblerData<dim> data (dof,
@@ -55,7 +54,6 @@ void MatrixCreator<dim>::create_mass_matrix (const DoFHandler<dim>    &dof,
                                             dVector                  &rhs_vector,
                                             const Function<dim> * const a) {
   UpdateFlags update_flags = UpdateFlags(update_q_points |
-                                        update_jacobians |
                                         update_JxW_values);
   const AssemblerData<dim> data (dof,
                                 true, true,
@@ -86,7 +84,6 @@ void MatrixCreator<dim>::create_laplace_matrix (const DoFHandler<dim>    &dof,
                                                const Function<dim> * const a) {
   dVector dummy;   // no entries, should give an error if accessed
   UpdateFlags update_flags = UpdateFlags(update_gradients |
-                                        update_jacobians |
                                         update_JxW_values);
   if (a != 0)
     update_flags = UpdateFlags(update_flags | update_q_points);
@@ -120,7 +117,6 @@ void MatrixCreator<dim>::create_laplace_matrix (const DoFHandler<dim>    &dof,
                                                const Function<dim> * const a) {
   UpdateFlags update_flags = UpdateFlags(update_q_points  |
                                         update_gradients |
-                                        update_jacobians |
                                         update_JxW_values);
   const AssemblerData<dim> data (dof,
                                 true, true,
index a9136e875acd53c3a366b6c887db32bbcede0a47..bfffb5a645628d8074da8f0ac5af6baa6f3d570e 100644 (file)
@@ -124,7 +124,6 @@ void VectorTools<dim>::integrate_difference (const DoFHandler<dim>    &dof,
   difference.reinit (dof.get_tria().n_active_cells());
   
   UpdateFlags update_flags = UpdateFlags (update_q_points  |
-                                         update_jacobians |
                                          update_JxW_values);
   if ((norm==H1_seminorm) || (norm==H1_norm))
     update_flags = UpdateFlags (update_flags | update_gradients);
index 7c09fcc3f59490bb45b5e1d38d496e98c89fe9e8..b2424130dc329343c4381a0e99bbf3d77fca5f59 100644 (file)
@@ -252,7 +252,7 @@ void PoissonProblem<dim>::run (const unsigned int level) {
 
   cout << "    Assembling matrices..." << endl;
   UpdateFlags update_flags = UpdateFlags(update_q_points  | update_gradients |
-                                        update_jacobians | update_JxW_values);
+                                        update_JxW_values);
   
   ProblemBase<dim>::FunctionMap dirichlet_bc;
   dirichlet_bc[0] = boundary_values;
index fc2e8078388928986a147650d6780fcc5ffc837d..063f9902c5e6ab59981b6dd09114ef245902f484 100644 (file)
@@ -438,7 +438,7 @@ void PoissonProblem<dim>::run (ParameterHandler &prm) {
 
       cout << "    Assembling matrices..." << endl;
       UpdateFlags update_flags = UpdateFlags(update_q_points  | update_gradients |
-                                            update_jacobians | update_JxW_values);
+                                            update_JxW_values);
   
       ProblemBase<dim>::FunctionMap dirichlet_bc;
       dirichlet_bc[0] = solution_function;
index cd95664c8703c99bacdc63850ae37a4ca167d63d..d8fb3c517fa590d01f760ad08924e61662dfa955 100644 (file)
@@ -459,7 +459,7 @@ void PoissonProblem<dim>::run (ParameterHandler &prm) {
   dirichlet_bc[0] = boundary_values;
   assemble (equation, quadrature, fe,
            UpdateFlags(update_gradients | update_JxW_values |
-                       update_jacobians | update_q_points),
+                       update_q_points),
            dirichlet_bc, *boundary);
 
   cout << "    Solving..." << endl;

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.