From: Timo Heister Date: Mon, 27 Mar 2017 19:32:56 +0000 (-0400) Subject: fix step-41 warning X-Git-Tag: v9.0.0-rc1~1776^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=ae757db66f1002c42b24a2005a74c65c0c00824a;p=dealii.git fix step-41 warning --- 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