Dear community, I am studying with pleasure step-42 and I got a bit confused in the method PlasticityContactProblem::solve_newton. At the very end of it, we find these instructions:
old_active_set = active_set; previous_residual_norm = residual_norm; if (Utilities::MPI::sum <http://www.dealii.org/8.5.0/doxygen/deal.II/namespaceUtilities_1_1MPI.html#ab544a3bf3301a6dd3e705ee352c5551b>((active_set == old_active_set) ? 0 : 1, mpi_communicator) == 0) { pcout << " Active set did not change!" << std::endl; if (residual_norm < 1e-10) break; } Now, I am a bit confused by the if statement and I am clearly missing something. The first instruction makes old_active_set equal to active_set. Accordingly, wouldn't active_set == old_active_set always be true? What am I missing? Thanks for your help Alberto -- Informativa sulla Privacy: http://www.unibs.it/node/8155 -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. For more options, visit https://groups.google.com/d/optout.