Python code for generating random matrices relevant in quantum optics and computing Bloch-Messiah and Williamson decompositions.