From 832e4aec0b675d3b05df529e28fe6fc87370f515 Mon Sep 17 00:00:00 2001 From: frohne Date: Sat, 23 Feb 2013 02:51:58 +0000 Subject: [PATCH] setting the newton iteration index like in the documentation git-svn-id: https://svn.dealii.org/trunk@28536 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-42/step-42.cc | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/deal.II/examples/step-42/step-42.cc b/deal.II/examples/step-42/step-42.cc index a8814b1a5b..21296f7d5a 100644 --- a/deal.II/examples/step-42/step-42.cc +++ b/deal.II/examples/step-42/step-42.cc @@ -572,8 +572,8 @@ namespace Step42 return_value = p(1); if (component == 2) { - return_value = -std::sqrt (0.36 - (p(0)-0.5)*(p(0)-0.5) - (p(1)-0.5)*(p(1)-0.5)) + 1.59; -// return_value = 1.999 - input_obstacle_copy->obstacle_function (p(0), p(1)); + return_value = -std::sqrt (0.36 - (p(0)-0.5)*(p(0)-0.5) - (p(1)-0.5)*(p(1)-0.5)) + 1.59; + // return_value = 1.999 - input_obstacle_copy->obstacle_function (p(0), p(1)); } return return_value; } @@ -1263,14 +1263,14 @@ namespace Step42 additional_data.aggregation_threshold = 1e-2; IndexSet active_set_old (active_set); - unsigned int j = 0; + unsigned int j = 1; unsigned int number_assemble_system = 0; for (; j<=100; j++) { // Solve an elastic problem to obtain a better start solution - if (j == 0) + if (j == 1) plast_lin_hard->set_sigma_0 (1e+10); - else if (j == 1) + else if (j == 2) plast_lin_hard->set_sigma_0 (sigma_hlp); pcout<< " " <