Usually, the first-principles packages like VASP and QE involve many aspects to totally figure out what their subroutines do. For example, to output Berry connection $\mathcal{A}^\alpha_{nn}$ used by these codes to calculate electronic polarization is hard and time-consuming, especially for those...