Add a general way to compute the diamond norm distance

Currently rotation_sythesis/Channel have specialized methods that compute the diamond norm distance for certain cases like distance between two unitaries, or channels with known analytical formula. it would be nice to have a method that can work for general channels.