Gradient should have zero curl; harmonic has zero local curl but nonzero periods; the positive control has nonzero curl.
One smooth field · several discrete addresses
Curl is circulation.
But around which cell?
Run the same manufactured tangent field through a primal triangle loop and a barycentric dual loop. Then separate that measurement question from the choice of tangent connection.
Triangle and dual-cell circulation return scalars at different mesh locations. “Truth error” compares each with the analytic curl.
One resolution is a picture. The bottom chart is the actual consistency experiment.
(Σ αᵢⱼ) / Af- truth error
- —
- curl RMS
- —
- samples
- —
(∮ X♭) / A★v- truth error
- —
- curl RMS
- —
- samples
- —
- u period
- —
- v period
- —
- analytic curl RMS
- —
The analytic torus Levi–Civita angle is the manufactured reference. These are connection errors, not curl values.
- extrinsic RMS
- —
- intrinsic RMS
- —
- edge family
- —
| grid | triangle curl error | dual curl error | extrinsic angle error | intrinsic angle error |
|---|
αᵢⱼ = ½(Xᵢ + Xⱼ) · (pⱼ − pᵢ)Endpoint trapezoids turn tangent samples into oriented line integrals.
curlf = Σ∂f α / AfThree boundary values produce one face-centered scalar.
curlv = ∮∂v★ X♭ / A★vMidpoints and barycenters bound one positive vertex dual cell.
DᵢⱼX = Xⱼ − QᵢⱼXᵢTransport compares vectors across tangent planes; circulation compares scalar integrals.