def rephase_standard(UuL, UdL, UuR, UdR):
K = UuL.conj().T @ UdL<EOL>f = mixing_phases(K)<EOL>Fdelta = np.diag(np.exp([<NUM_LIT>*f['<STR_LIT>'], <NUM_LIT>*f['<STR_LIT>'], <NUM_LIT>*f['<STR_LIT>']]))<EOL>Fphi = np.diag(np.exp([-<NUM_LIT>*f['<STR_LIT>']/<NUM_LIT>, -<NUM_LIT>*f['<STR_LIT>']/<NUM_LIT>, <NUM_LIT:0>]))<EOL>return UuL @ Fdelta, UdL @ Fphi.conj(), UuR @ Fdelta, UdR @ Fphi.conj()<EOL>
Function to rephase the quark rotation matrices in order to
    obtain the CKM matrix in standard parametrization.

    The input matrices are assumed to diagonalize the up-type and down-type
    quark matrices like

    ```
    UuL.conj().T @ Mu @ UuR = Mu_diag
    UdL.conj().T @ Md @ UdR = Md_diag
    ```

    The CKM matrix is given as `VCKM = UuL.conj().T @ UdL`.

    Returns a tuple with the rephased versions of the input matrices.