Dear all,
I am using the 3D IGA contact (celmt06.f) and cannot understand what exactly is happening in the subroutine cont_projection06.f.
My question is - why do we not add geometrical part of the matrix while updating the trial point in the Newton loop unless the contribution is very small.
183 ! Add geometric part if small enough
184 if(jgmax.lt.jjmax*1.d-03) then
185 jj = jj + jg
186 endif
Thank you in advance for any hints.
Best regards
