From ae757db66f1002c42b24a2005a74c65c0c00824a Mon Sep 17 00:00:00 2001 From: Timo Heister Date: Mon, 27 Mar 2017 15:32:56 -0400 Subject: [PATCH] fix step-41 warning --- examples/step-41/step-41.cc | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/examples/step-41/step-41.cc b/examples/step-41/step-41.cc index b14511c069..948caed7af 100644 --- a/examples/step-41/step-41.cc +++ b/examples/step-41/step-41.cc @@ -439,7 +439,10 @@ namespace Step41 TrilinosWrappers::MPI::Vector lambda (complete_index_set(dof_handler.n_dofs())); complete_system_matrix.residual (lambda, solution, complete_system_rhs); - contact_force.ratio (lambda, diagonal_of_mass_matrix); + + // compute contact_force[i] = - lambda[i] * diagonal_of_mass_matrix[i] + contact_force = lambda; + contact_force.scale (diagonal_of_mass_matrix); contact_force *= -1; // The next step is to reset the active set and constraints objects and to -- 2.39.5