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