Loop subdivision surface

On a regular triangular grid, Loop subdivision retains an old vertex with weight $10/16$ and adds its six neighbors with weight $1/16$ each. A new edge vertex has weights $3/8$ on the endpoints and $1/8$ on the two opposite vertices. Extraordinary vertices require valence-dependent rules.

// Target: analysis.bigb

For the three unit-grid edge vectors $v_1,v_2,v_3$ and a multivariate polynomial $f$ of total degree at most three, sampled regular-grid Loop subdivision generates the graph of $f+\frac{1}{12}{\sum }_{r=1}^3(v_r\cdot \nabla {)}^{2}f$. Thus affine polynomials are reproduced exactly and all polynomials through degree three retain their degree. Fourth-degree data constant along one edge direction reduce to a Cardinal cubic B-spline, whose limit need not be a polynomial.

On an infinite regular triangular grid, height data constant along one edge direction remain constant along that direction under Loop subdivision. Across rows the induced rule is $g_{2j}^{\prime }=(g_{j-1}+6g_j+g_{j+1})/8$, $g_{2j+1}^{\prime }=(g_j+g_{j+1})/2$, the centered Cardinal cubic B-spline rule.

// Target: analysis.bigb