Matched pairs approach to set theoretic solutions of the Yang–Baxter equation