From 162fdf4e3ccd17a79812797257e6d047175cdd93 Mon Sep 17 00:00:00 2001 From: frohne Date: Mon, 25 Feb 2013 15:45:14 +0000 Subject: [PATCH] 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 --- deal.II/examples/step-42/step-42.cc | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) 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); -- 2.39.5