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.

Reply via email to