changes in the assembler_mass_matrix_diagonal function to (\varphi_i,1)

git-svn-id: https://svn.dealii.org/trunk@28554 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/examples/step-42/step-42.cc b/deal.II/examples/step-42/step-42.cc
index 21296f7d5ad8fa1b6f286af6e4b5a54c75693b82..454abed00124bbbed49738e8167f7d79b011b8fb 100644 (file)
--- a/deal.II/examples/step-42/step-42.cc
+++ b/deal.II/examples/step-42/step-42.cc
@@ -981,7 +981,7 @@ namespace Step42
   template <int dim>
   void PlasticityContactProblem<dim>::assemble_mass_matrix_diagonal (TrilinosWrappers::SparseMatrix &mass_matrix)
   {
-    QTrapez<dim-1>  face_quadrature_formula;
+       QTrapez<dim-1>  face_quadrature_formula;
 
     FEFaceValues<dim> fe_values_face (fe, face_quadrature_formula,
                                       update_values   |
@@ -992,6 +992,9 @@ namespace Step42
     const unsigned int   n_face_q_points    = face_quadrature_formula.size();
 
     FullMatrix<double>   cell_matrix (dofs_per_cell, dofs_per_cell);
+    Tensor<1,dim,double> ones (dim);
+    for (unsigned i=0; i<dim; i++)
+       ones[i] = 1.0;
 
     std::vector<unsigned int> local_dof_indices (dofs_per_cell);
 
@@ -1013,7 +1016,7 @@ namespace Step42
               for (unsigned int q_point=0; q_point<n_face_q_points; ++q_point)
                 for (unsigned int i=0; i<dofs_per_cell; ++i)
                   cell_matrix(i,i) += (fe_values_face[displacement].value (i, q_point) *
-                                       fe_values_face[displacement].value (i, q_point) *
+                                                          ones *
                                        fe_values_face.JxW (q_point));
 
               cell->get_dof_indices (local_dof_indices);