def add_scaling(bonds, means, stds):
    """Add the scaling information to the bonds dataframe.

    Args:
        bonds (pd.DataFrame): The dataframe of the bonds, after :func:`enhance_bonds`.
        means (dict): Output of :func:`get_scaling`.
        stds (dict): Output of :func:`get_scaling`.

    Returns:
        pd.DataFrame: Same dataframe, with added columns.

    """
    # Add mean/std scaling factors to bonds dataframe
    bonds["sc_mean"] = bonds["labeled_type"].apply(lambda x: means[x])
    bonds["sc_std"] = bonds["labeled_type"].apply(lambda x: stds[x])
    if "scalar_coupling_constant" in bonds.columns:
        bonds["sc_scaled"] = (bonds["scalar_coupling_constant"] - bonds["sc_mean"]) / bonds["sc_std"]
    return bonds