“…At each subdomain interface, the Dirichlet and Neumann boundary conditions should be satisfied, which can be written as where and are the unit normal vectors of subdomain i and j , respectively. The Neumann boundary condition contains an unknown vector Λ , often referred to as the ‘dual’ unknown, which must be solved to recover the field values inside each subdomain . In order to enforce the continuity condition along corner edges, one set of global corner DOFs, A c , often referred to as the ‘primal’ unknowns , is used within the entire computational domain.…”