From: frohne 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 void PlasticityContactProblem::assemble_mass_matrix_diagonal (TrilinosWrappers::SparseMatrix &mass_matrix) { - QTrapez face_quadrature_formula; + QTrapez face_quadrature_formula; FEFaceValues 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 cell_matrix (dofs_per_cell, dofs_per_cell); + Tensor<1,dim,double> ones (dim); + for (unsigned i=0; i local_dof_indices (dofs_per_cell); @@ -1013,7 +1016,7 @@ namespace Step42 for (unsigned int q_point=0; q_pointget_dof_indices (local_dof_indices);