weno_matrices.py 文件源码

def oscillation_indicator():
    """ Generate the oscillation indicator matrix
    """
    ? = zeros([N+1, N+1])
    for a in range(1,N+1):
        ?Dera = ?Der[a]
        for p, m in product(range(N+1), range(N+1)):
            antiderivative = polyint(?Dera[p] * ?Dera[m])
            ?[p,m] += antiderivative(1) - antiderivative(0)
    return ?