This appears to contradict the theory presented in the paper? Moreover, numerous models fail to converge when the parameter ifhardedge=true is set, yet converge when ifhardedge=false. Could this be related to the solution method employed? The straightforward Newton's method may not support solving complex scenarios.
1005266_sf_rem.zip