StabilizedModes.sqrt_hop scaling and off-diagonal entries change when adding lead_conservation laws
I multiply StabilizedModes.vecs by StabilizedModes.sqrt_hop (H.C.) to get a matrix related to the incoming lead modes. But sqrt_hop behaves very differently when I add a lead_conservation law like spin or particle-hole. The scaling and off-diagonal elements are very different to what I expect. See this minimal example: sqrt_hop_snip.py
What is should I expect from sqrt_hop upon adding a lead_conservation law?