Statement

For Q(a)=H(aa)Q(a)=H(a\odot a) with HH invertible and all aj>0a_j>0, its derivative is DQ(a)[h]=2H(ah)DQ(a)[h]=2H(a\odot h). Therefore the linear map

LΣ=12diag(a)1H1ΣL\Sigma=\frac12\operatorname{diag}(a)^{-1}H^{-1}\Sigma

is an inverse for this :

DQ(a)[LΣ]=Σ.DQ(a)[L\Sigma]=\Sigma.

It accepts arbitrary signed increments Σ\Sigma.

Exact remaining error

The complete update satisfies

Q(a+LΣ)=Q(a)+Σ+H((LΣ)(LΣ)).Q(a+L\Sigma)=Q(a)+\Sigma+H\bigl((L\Sigma)\odot(L\Sigma)\bigr).

The final quadratic term must be retained. The linear identity does not assert that an arbitrary signed finite target is exactly realizable by nonnegative squared amplitudes. In repeated linear corrections the matrix and base amplitudes in LL must be kept at the values used to define that operator, unless a new linearization is explicitly chosen.