Example #1
0
    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
Example #2
0
    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
Example #6
0
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