adjoint.py 文件源码

python
阅读 26 收藏 0 点赞 0 评论 0

项目:SeisFlows_tunnel 作者: DmBorisov 项目源码 文件源码
def Envelope(syn, obs, nt, dt, eps=0.05):
    # envelope difference
    # (Yuan et al 2015, eq 16)
    esyn = abs(_analytic(syn))
    eobs = abs(_analytic(obs))
    if esyn.all() == 0.0:
        etmp = 0.0
    else:
        etmp = (esyn - eobs)/(esyn + eps*esyn.max())

    wadj = etmp*syn - _np.imag(_analytic(etmp*_np.imag(_analytic(syn))))
    return wadj
评论列表
文章目录


问题


面经


文章

微信
公众号

扫码关注公众号