When dealing with Boson system, we often meet Bogoliubov transformation, which is different from the unitary transformation. The matrix may be large that needs the numerical calculation. My question is how to perform the Bogoliubov transformation numerically?