Hi,
Is it possible to iterate contact in the CG solver ??
My contact definition :
n = specialcf.normal(3)
penalty = 1e10
a += penalty * InnerProduct((u * n), (v * n)) * dx(element_boundary=True,
definedon=mesh.Boundaries(contact_faces))
My CGSolver call :
from ngsolve.krylovspace import CGSolver
inv = CGSolver(a.mat, pre, printrates='\r', tol=1e-6, maxiter=500)
gfu.vec.data += inv * (f.vec - a.Apply(gfu.vec))
Thanks