From: frohne <frohne@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Mon, 25 Feb 2013 15:45:14 +0000 (+0000)
Subject: changes in the assembler_mass_matrix_diagonal function to (\varphi_i,1)
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=162fdf4e3ccd17a79812797257e6d047175cdd93;p=dealii-svn.git

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 21296f7d5a..454abed001 100644
--- 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);