I can do that for you. Do I understand correctly that the only change
that needs to happen is to move
active_set_old = active_set;
to *after* the if-statement?
Exactly.
Great, thanks for confirming, Joerg!
The patch is now here:
https://www.dealii.org/8.5.1/doxygen/deal.II/Tutorial.htm
>>On 01/17/2018 02:45 PM, Frohne, Joerg wrote:
>> obviously this is a mistake. I checked the source files which we have used
>> for results in corresponding paper.
>> There I have found the following coding:
>>
>>resid_old = resid;
>>
>>resid_vector = system_rhs_newton;
>>