Function compute_divergence

scripts/geodesic.py:106–119

(mesh, X)

    return np.array(grad);

def compute_divergence(mesh, X):
    cotans = mesh.get_attribute("cotan").reshape((-1, 3));
    edges = mesh.get_attribute("edges").reshape((-1, 3, 3));
    div_X = np.zeros(mesh.num_vertices);

    for fi,f in enumerate(mesh.faces):
        x_fi = X[fi];
        for i in range(3):
            j = (i+1) % 3;
            k = (i+2) % 3;
            div_X[f[i]] += cotans[fi, k] * np.dot(x_fi, edges[fi, i]) / 2.0;
            div_X[f[i]] += cotans[fi, j] * np.dot(x_fi,-edges[fi, k]) / 2.0;

    return div_X;

def main():
    args = parse_args();