def METHOD_NAME(ref_rho, K_0, K_prime, G_0, G_prime):
    rock = burnman.Mineral()
    rock.params["V_0"] = 10.0e-6
    rock.params["molar_mass"] = ref_rho * rock.params["V_0"]
    rock.params["K_0"] = K_0
    rock.params["Kprime_0"] = K_prime
    rock.params["G_0"] = G_0
    rock.params["Gprime_0"] = G_prime
    rock.set_method("bm3")
    temperature = np.empty_like(seis_p)
    mat_rho, mat_vphi, mat_vs = rock.evaluate(
        ["density", "v_phi", "v_s"], seis_p, temperature
    )
    return mat_rho, mat_vphi, mat_vs