fv_fluxes.py 文件源码

def Aint(qL, qR, d):
    """ Returns the Osher-Solomon jump matrix for A, in the dth direction
        NB: eig function should be replaced with analytical function, if known
    """
    ret = zeros(n, dtype=complex128)
    ?q = qR - qL
    for i in range(N+1):
        q = qL + nodes[i] * ?q
        J = jacobian(q, d)
        ?, R = eig(J, overwrite_a=1, check_finite=0)
        b = solve(R, ?q, check_finite=0)
        ret += weights[i] * dot(R, abs(?)*b)
    return ret.real