The condensed interface matrices are assembled to yield the global interface equations. These are solved for the interface variables.
Introduction of the interface variables into the subdomain equations give then the internal variables. Since the problem is non linear this procedure is connected with a Newton-Raphson procedure.