Statement

Let n(t)∈Rd∖{0}n(t)\in\mathbb R^d\setminus\{0\} be C1C^1, and let A(t),g(t)A(t),g(t) be continuous. Write

Pn=I−n⊗n∣n∣2.P_n=I-\frac{n\otimes n}{|n|^2}.

The constrained equation a′=Aa+g+μna'=Aa+g+\mu n, n⋅a=0n\cdot a=0, is equivalent, for initially transverse data, to

a′=Pn(Aa+g)−n(n′⋅a)∣n∣2,μ=−n′⋅a+n⋅(Aa+g)∣n∣2.a'=P_n(Aa+g)-\frac{n(n'\cdot a)}{|n|^2}, \qquad \mu=-\frac{n'\cdot a+n\cdot(Aa+g)}{|n|^2}.

The alone is insufficient when the normal moves.

Verification

Differentiate n⋅a=0n\cdot a=0 to obtain n′⋅a+n⋅a′=0n'\cdot a+n\cdot a'=0, then solve for μ\mu. Conversely, the displayed evolution makes the derivative of n⋅an\cdot a zero for every solution, so an initially zero constraint remains zero. The same formulas hold for complex a,ga,g by complex linear extension with the real normal nn.