[Kwant] What is the algorithm for calculating the mode decomposition

Dear kwant users， In the latest Kwant manual describing the details about kwant.physics.modes, I find a sentence "This function uses the most stable and efficient algorithm for calculating the mode decomposition that the Kwant authors are aware about. Its details are to be published." Do anyone know what the algorithm is? I look forward some references if possible. Regards, Zhan