Domain decomposition: 3D full-anisotropy coupling solver