def error(a, b, c): mat_vp, mat_vs, mat_rho = calc_velocities(a, b, c) vs_err = burnman.l2(depths, mat_vs, seis_vs) / 1e9 vp_err = burnman.l2(depths, mat_vp, seis_vp) / 1e9 den_err = burnman.l2(depths, mat_rho, seis_rho) / 1e9 return vs_err + vp_err + den_err
def error(a, b, c): mat_vp, mat_vs, mat_rho = calc_velocities(a, b, c) vs_err = burnman.l2(depths, mat_vs, seis_vs) / 1e9 vp_err = burnman.l2(depths, mat_vp, seis_vp) / 1e9 den_err = burnman.l2(depths, mat_rho, seis_rho) / 1e9 return vs_err + vp_err + den_err
def error(mg_pv_K,mg_pv_K_prime,mg_pv_G,mg_pv_G_prime,fe_pv_K,fe_pv_K_prime,fe_pv_G,fe_pv_G_prime): mat_vp, mat_vs, mat_rho = calc_velocities(mg_pv_K,mg_pv_K_prime,mg_pv_G,mg_pv_G_prime,fe_pv_K,fe_pv_K_prime,fe_pv_G,fe_pv_G_prime) vs_err = burnman.l2(depths, mat_vs, seis_vs) vp_err = burnman.l2(depths, mat_vp, seis_vp) den_err = burnman.l2(depths, mat_rho, seis_rho) return vs_err #+ vp_err #+ den_err
def error(amount_pv, iron_pv, iron_fp): mat_vp, mat_vs, mat_rho = calc_velocities(amount_pv, iron_pv, iron_fp) vs_err = burnman.l2(depths, mat_vs, seis_vs) / 1e9 vp_err = burnman.l2(depths, mat_vp, seis_vp) / 1e9 den_err = burnman.l2(depths, mat_rho, seis_rho) / 1e9 # print vs_err, vp_err, den_err return vs_err + vp_err + den_err
def error(amount_pv, iron_pv, iron_fp): mat_vp, mat_vs, mat_rho = calc_velocities(amount_pv, iron_pv, iron_fp) vs_err = burnman.l2(depths, mat_vs, seis_vs) / 1e9 vp_err = burnman.l2(depths, mat_vp, seis_vp) / 1e9 den_err = burnman.l2(depths, mat_rho, seis_rho) / 1e9 # print vs_err, vp_err, den_err return vs_err + vp_err + den_err
def error(mg_pv_K, mg_pv_K_prime, mg_pv_G, mg_pv_G_prime, fe_pv_K, fe_pv_K_prime, fe_pv_G, fe_pv_G_prime): mat_vp, mat_vs, mat_rho = calc_velocities(mg_pv_K, mg_pv_K_prime, mg_pv_G, mg_pv_G_prime, fe_pv_K, fe_pv_K_prime, fe_pv_G, fe_pv_G_prime) vs_err = burnman.l2(depths, mat_vs, seis_vs) vp_err = burnman.l2(depths, mat_vp, seis_vp) den_err = burnman.l2(depths, mat_rho, seis_rho) return vs_err #+ vp_err #+ den_err