[Box] variable switch at boundary
It seems that there is a bug in the variable switch for box at the boundary.
initPriVarSwitch_ method, the states of Dirichlet nodes are corrected.
This method is called in
newtonBegin, which however is called after the initialization
Therefore, when calling the
newtonUpdate method, the corrected states of
We could fix this by either calling
uLastIter is initialized
or by setting this vector after
So we could simply change
if (numSteps_ > 0) uLastIter = uCurrentIter;
if (numSteps_ >= 0) uLastIter = uCurrentIter;