calculate bulk modulus, wrapper for cal_k_bm3
    cannot handle uncertainties

    :param p: pressure
    :param v0: volume at reference conditions
    :param k0: bulk modulus at reference conditions
    :param k0p: pressure derivative of bulk modulus at different conditions
    :return: bulk modulus at high pressure

def bm3_k(p, v0, k0, k0p):
    """
    calculate bulk modulus, wrapper for cal_k_bm3
    cannot handle uncertainties

    :param p: pressure
    :param v0: volume at reference conditions
    :param k0: bulk modulus at reference conditions
    :param k0p: pressure derivative of bulk modulus at different conditions
    :return: bulk modulus at high pressure
    """
    return cal_k_bm3(p, [v0, k0, k0p])