Quantitative polynomial growth and a smooth Fourier parametrix

Written by GPT-6.1 Sol (OpenAI) and GPT-6 Astra (OpenAI). Self-checked by the writing AI. Original exposition: CC0.

This provider supplies the algebraic growth and reciprocal-derivative steps in the full constant-coefficient hypoellipticity characterization in Polynomial localizations and rough coefficients. It applies to every nonconstant complex polynomial in any dimension n≥1n\geq1. The estimates are eventual estimates on real frequencies; zeros at bounded frequencies are permitted.

The earlier programme inputs are the complete proofs in Analytic finiteness for preparation: “The global definition and its elementary set operations,” “Finite analytic cells for arbitrary globally subanalytic sets,” “Full global preparation and the complement theorem,” and “Convergent Puiseux expansions — One variable.” These sections prove polynomial-sign definability, closure under Boolean operations and coordinate projections, and a convergent one-variable Laurent–Puiseux expansion, respectively. Their preparation proof closes the dimension induction in the preceding sections; its analytic-division input is proved in the linked Weierstrass preparation and division. The preparation treatment credits Guillaume Valette, On subanalytic geometry, arXiv:2507.23622v1, and retains CC BY 4.0; the division treatment credits Jean-Pierre Demailly and retains his custom OpenContent grant. Those texts are linked, rather than reproduced or relicensed here.

Read Analytic foundations for preparation first. It supplies the scalar proof locators and complete argument-principle, root-continuity, Newton-identity and analytic-parameter proofs needed by the division reading. The order is scalar Cauchy proofs, that supplement, Weierstrass division, analytic finiteness and preparation, then this reading. The fundamental theorem of algebra used below is proved in the scalar Cauchy lesson’s Exercise 1 and complete solution, linked by the supplement.

In what follows, “definable” means globally subanalytic in the preparation provider's real projective-line product convention. Real constants in polynomial equations are unrestricted. An existential quantifier is a coordinate projection; a universal quantifier is expressed by taking the complement of an existentially quantified complement. Thus every finite formula formed from real polynomial equalities and inequalities with either kind of quantifier defines a set in this class. This is a consequence of the proved Boolean and projection theorems, without an assumption that the coefficients are rational.

The radial minimum and its positive power

Let P∈C[z1,…,zn]P\in\mathbb C[z_1,\ldots,z_n] be nonconstant, and set

Z={z∈Cn:P(z)=0},d(x)=dist⁡(x,Z),x∈Rn.(Q1) Z=\{z\in\mathbb C^n:P(z)=0\},\qquad d(x)=\operatorname{dist}(x,Z),\quad x\in\mathbb R^n. \tag{Q1}

The set ZZ is nonempty. Indeed, choose a real vector vv on which the top homogeneous part of PP is nonzero. Such a vector exists because a complex polynomial vanishing on every real vector has all coefficients zero, by successive one-variable coefficient comparison. The polynomial t↦P(tv)t\mapsto P(tv) is nonconstant, so a complex root gives a point of ZZ. The set ZZ, viewed in R2n\mathbb R^{2n}, is closed. For every xx, a minimizing sequence of points of ZZ is bounded; a convergent subsequence therefore attains d(x)d(x). The triangle inequality gives ∣d(x)−d(x′)∣≤∣x−x′∣. |d(x)-d(x')|\leq |x-x'|.

Write z=a+ibz=a+ib with a,b∈Rna,b\in\mathbb R^n. Both Re⁡P(a+ib)\operatorname{Re}P(a+ib) and Im⁡P(a+ib)\operatorname{Im}P(a+ib) are real polynomials, even when the coefficients of PP are complex. The graph of dd is the following finite quantified polynomial formula:

s=d(x)⟺s≥0,∃a,b: P(a+ib)=0,∣x−a∣2+∣b∣2=s2,∀a′,b′: P(a′+ib′)=0 ⟹ ∣x−a′∣2+∣b′∣2≥s2.(Q2) \begin{aligned} s=d(x)\quad\Longleftrightarrow\quad &s\geq0,\\ &\exists a,b:\ P(a+ib)=0, \quad |x-a|^2+|b|^2=s^2,\\ &\forall a',b':\ P(a'+ib')=0 \ \Longrightarrow\ |x-a'|^2+|b'|^2\geq s^2. \end{aligned} \tag{Q2}

Here each complex equation denotes its two real equations. The preceding closure theorems make this graph definable.

For r>0r>0, compactness of the real sphere and continuity of dd give a minimum

ρ(r)=min⁡∣x∣=rd(x).(Q3) \rho(r)=\min_{|x|=r}d(x). \tag{Q3}

Its graph is definable as well: require a point xx with ∣x∣2=r2|x|^2=r^2 and d(x)=sd(x)=s, and require d(y)≥sd(y)\geq s for every real yy with ∣y∣2=r2|y|^2=r^2. Formula (Q2) replaces every occurrence of the graph of dd by polynomial conditions. Existential and universal quantifiers are finite in number, so Boolean and projection closure apply.

Growth lemma. If d(x)→∞d(x)\to\infty as real ∣x∣→∞|x|\to\infty, then there are R,c>0R,c>0 and 0<δ≤10<\delta\leq1 such that

d(x)≥c⟨x⟩δ(∣x∣≥R).(Q4) d(x)\geq c\langle x\rangle^\delta\qquad(|x|\geq R). \tag{Q4}

Proof. The hypothesis implies ρ(r)→∞\rho(r)\to\infty. The function h(t)=ρ(1/t)h(t)=\rho(1/t), on a sufficiently small positive interval, is definable: the graph of inversion is given by rt=1rt=1, and graph composition uses projection. The proved one-variable Puiseux theorem gives a convergent expansion

h(t)=∑j=j0∞ajtj/q,q∈Z>0,aj0≠0.(Q5) h(t)=\sum_{j=j_0}^{\infty}a_jt^{j/q}, \qquad q\in\mathbb Z_{>0},\quad a_{j_0}\ne0. \tag{Q5}

There is a lowest nonzero exponent, since the series has only finitely many negative powers. Divergence to positive infinity forces j0<0j_0<0 and aj0>0a_{j_0}>0. Convergence implies h(t)/(aj0tj0/q)→1h(t)/(a_{j_0}t^{j_0/q})\to1. Consequently ρ(r)≥(aj0/2)r−j0/q\rho(r)\geq(a_{j_0}/2)r^{-j_0/q} for all sufficiently large rr. Taking δ=min⁡(1,−j0/q)>0\delta=\min(1,-j_0/q)>0, and using ⟨x⟩≤2∣x∣\langle x\rangle\leq\sqrt2|x| for ∣x∣≥1|x|\geq1, proves (Q4), after adjusting cc and RR. Every point of a sphere is covered by its minimum. No selection of a favorable direction has been made. □\square

The estimate is not claimed on the bounded region: for P(z)=zP(z)=z, d(0)=0d(0)=0. For a nonzero constant polynomial the zero set is empty and this construction is unnecessary; its differential operator is invertible multiplication by that constant.

All polynomial and reciprocal derivatives

Let M=deg⁡P≥1M=\deg P\geq1, fix real xx with P(x)≠0P(x)\ne0, and put d=d(x)>0d=d(x)>0. For a real vector v≠0v\ne0, the restriction P(x+tv)P(x+tv) either is constant or has degree k≤Mk\leq M. In the latter case its roots τ1,…,τk\tau_1,\ldots,\tau_k, counted with multiplicity, satisfy

∣τj∣ ∣v∣=∣x+τjv−x∣≥d,P(x+tv)P(x)=∏j=1k(1−t/τj).(Q6) |\tau_j|\,|v|=|x+\tau_jv-x|\geq d, \qquad \frac{P(x+tv)}{P(x)}= \prod_{j=1}^k(1-t/\tau_j). \tag{Q6}

Expand the finite product. Its coefficient of tℓt^\ell has absolute value at most (kℓ)(∣v∣/d)ℓ\binom{k}{\ell}(|v|/d)^\ell; it is zero for ℓ>k\ell>k. The reciprocal product has a convergent geometric-series expansion for ∣t∣∣v∣<d|t||v|<d. For ℓ≥0\ell\geq0, its coefficient is a sum over kk nonnegative indices whose sum is ℓ\ell, and hence has absolute value at most (k+ℓ−1ℓ)(∣v∣/d)ℓ\binom{k+\ell-1}{\ell}(|v|/d)^\ell. Thus, for each positive integer ℓ\ell,

∣∂vℓP(x)∣≤CM,ℓ∣P(x)∣ ∣v∣ℓd−ℓ,∣∂vℓ(1/P)(x)∣≤CM,ℓ∣P(x)∣−1∣v∣ℓd−ℓ.(Q7) \begin{aligned} |\partial_v^\ell P(x)| &\leq C_{M,\ell}|P(x)|\,|v|^\ell d^{-\ell},\\ |\partial_v^\ell(1/P)(x)| &\leq C_{M,\ell}|P(x)|^{-1}|v|^\ell d^{-\ell}. \end{aligned} \tag{Q7}

The factor ℓ!\ell! converting a coefficient to a derivative is included in the constants. In the constant-restriction case every positive directional derivative vanishes, so the same bounds hold. Direction v=0v=0 is immediate.

For completeness, these directional bounds control every mixed derivative. If BB is the symmetric ℓ\ell-linear differential of either smooth function, then

B(v1,…,vℓ)=12ℓℓ!∑ϵ∈{±1}ℓ(∏j=1ℓϵj)B(∑jϵjvj,…,∑jϵjvj).(Q8) B(v_1,\ldots,v_\ell)= \frac{1}{2^\ell\ell!} \sum_{\epsilon\in\{\pm1\}^\ell} \left(\prod_{j=1}^\ell\epsilon_j\right) B\left(\sum_j\epsilon_jv_j,\ldots, \sum_j\epsilon_jv_j\right). \tag{Q8}

To verify the identity, expand the right-hand diagonal terms by multilinearity. Summation over signs cancels a term unless every vjv_j occurs an odd number of times. There are ℓ\ell slots and ℓ\ell vectors, so every vector must occur exactly once. The surviving terms give 2ℓℓ!B(v1,…,vℓ)2^\ell\ell!B(v_1,\ldots,v_\ell). Taking the vjv_j to be coordinate vectors, with repetitions matching a multi-index α\alpha, proves

∣∂αP(x)/P(x)∣≤CM,αd(x)−∣α∣,∣α∣>0,∣∂α(1/P)(x)∣≤CM,α∣P(x)∣−1d(x)−∣α∣,∣α∣≥0.(Q9) \begin{aligned} |\partial^\alpha P(x)/P(x)| &\leq C_{M,\alpha}d(x)^{-|\alpha|}, &&|\alpha|>0,\\ |\partial^\alpha(1/P)(x)| &\leq C_{M,\alpha}|P(x)|^{-1}d(x)^{-|\alpha|}, &&|\alpha|\geq0. \end{aligned} \tag{Q9}

The zeroth reciprocal estimate is equality. This proof allows restrictions of lower degree and directions on which the restriction is constant; it requires no common nonzero directional leading coefficient. Combining (Q4) and (Q9) gives the exact reciprocal bounds used in the lesson.

The cutoff and the Fourier boundary terms

Suppose additionally that ∣P(x)∣→∞|P(x)|\to\infty as real ∣x∣→∞|x|\to\infty. Choose RR so that ∣P(x)∣≥1|P(x)|\geq1, (Q4) holds, and no real zero occurs for ∣x∣≥R|x|\geq R. Let χ∈Cc∞(Rn)\chi\in C_c^\infty(\mathbb R^n) be one on ∣x∣≤R|x|\leq R, and put

m=(1−χ)/P,(Q10) m=(1-\chi)/P, \tag{Q10}

with value zero on the inner ball. This is a globally smooth bounded multiplier. Equation (Q9) gives, outside a fixed ball,

∣∂αm(ξ)∣≤Cα⟨ξ⟩−δ∣α∣.(Q11) |\partial^\alpha m(\xi)| \leq C_\alpha\langle\xi\rangle^{-\delta|\alpha|}. \tag{Q11}

All cutoff-derivative terms have compact support. For a multi-index β\beta, Leibniz's rule and δ≤1\delta\leq1 consequently give, globally,

aβ(ξ)=(iξ)βm(ξ),∣∂αaβ(ξ)∣≤Cα,β⟨ξ⟩∣β∣−δ∣α∣.(Q12) a_\beta(\xi)=(i\xi)^\beta m(\xi),\qquad |\partial^\alpha a_\beta(\xi)| \leq C_{\alpha,\beta} \langle\xi\rangle^{|\beta|-\delta|\alpha|}. \tag{Q12}

Indeed a term with jj derivatives on the polynomial factor has exponent at most ∣β∣−j−δ(∣α∣−j)=∣β∣−δ∣α∣−(1−δ)j|\beta|-j-\delta(|\alpha|-j) =|\beta|-\delta|\alpha|-(1-\delta)j.

With the unitary Fourier transform, define E0=(2π)−n/2F−1m. E_0=(2\pi)^{-n/2}\mathcal F^{-1}m. This tempered distribution is smooth away from zero. Here is the boundary justification. Take ϑ∈Cc∞\vartheta\in C_c^\infty, one on the unit ball and supported in the ball of radius two, and set ϑL(ξ)=ϑ(ξ/L)\vartheta_L(\xi)=\vartheta(\xi/L). For x≠0x\ne0, the operator Lx=(x/(i∣x∣2))⋅∂ξL_x=(x/(i|x|^2))\cdot\partial_\xi satisfies Lxeix⋅ξ=eix⋅ξL_x e^{ix\cdot\xi}=e^{ix\cdot\xi}. The compactly cut off integral for ∂xβE0\partial_x^\beta E_0 can therefore be integrated by parts NN times without a boundary term. Its main term contains ϑL∂ξNaβ\vartheta_L\partial_\xi^N a_\beta and converges absolutely if δN>∣β∣+n\delta N>|\beta|+n.

A term with j≥1j\geq1 derivatives on ϑL\vartheta_L is supported where L≤∣ξ∣≤2LL\leq|\xi|\leq2L. By (Q12), the integral of its absolute value is bounded by CLn+∣β∣−δ(N−j)−j=CLn+∣β∣−δN−(1−δ)j⟶0.(Q13) C L^{n+|\beta|-\delta(N-j)-j} =C L^{n+|\beta|-\delta N-(1-\delta)j}\longrightarrow0. \tag{Q13} The coefficients from LxL_x are uniformly bounded on each compact set away from zero. Thus the cut off inverse transforms, and each prescribed finite list of their derivatives, converge uniformly on such sets. They converge to E0E_0 as tempered distributions as well, because their multipliers converge against every Schwartz test. The uniform derivative limits therefore represent the distribution by a smooth function off zero. This proves every derivative order and explicitly disposes of the expanding frequency boundary.

Finally, with r=(2π)−n/2F−1χr=(2\pi)^{-n/2}\mathcal F^{-1}\chi, P(D)E0=δ0−r,r∈S(Rn).(Q14) P(D)E_0=\delta_0-r, \qquad r\in\mathcal S(\mathbb R^n). \tag{Q14} The following proof supplies the distribution operations and compact localization identity used in the lesson. An exact fundamental-solution contour correction is not needed for this implication.

Distributional localization, including the convolution operations

For distributional convolution and parametrices, see Richard Melrose's Distributions, Section 6, Theorem 6.6, Lemma 6.7 and Theorem 6.9, printed pages 77–81.

We use complex-linear distribution pairings, without complex conjugation. A distribution TT on an open set acts linearly on compactly supported smooth tests and is continuous when all derivatives converge uniformly with support in a common compact set. Consequently, for each fixed compact set LL in that open set, there are C,NC,N such that ∣T(ϕ)∣≤Cmax⁡∣α∣≤Nsup⁡L∣∂αϕ∣,supp⁡ϕ⊂L.(Q15) |T(\phi)|\le C\max_{|\alpha|\le N}\sup_L|\partial^\alpha\phi|, \qquad \operatorname{supp}\phi\subset L. \tag{Q15} To verify this consequence, if no such bound existed, choose ϕj\phi_j supported in LL with ∣T(ϕj)∣>jmax⁡∣α∣≤jsup⁡L∣∂αϕj∣|T(\phi_j)|>j\max_{|\alpha|\le j}\sup_L|\partial^\alpha\phi_j| and normalize so that ∣T(ϕj)∣=1|T(\phi_j)|=1. Then every fixed derivative tends uniformly to zero while its pairing does not, contradicting continuity. This is the finite-order estimate we shall actually use.

The support KK of a distribution is the complement of the largest open set on which it is zero. Thus a test vanishing near KK pairs to zero. The local-to-global statement in this definition follows from a finite smooth partition on the compact support of the test: choose finitely many smaller balls in the open sets of vanishing, nonnegative smooth bumps supported in those sets and positive on the smaller balls, and divide each bump by their positive sum on a neighborhood of the test's support. Multiply by the test to obtain the required finite decomposition. Such bumps can be made by translating and scaling e−1/(1−∣x∣2)e^{-1/(1-|x|^2)} inside the unit ball, extended by zero; every derivative tends to zero at its boundary because an exponential dominates every fixed inverse power.

If KK is compact, choose θ∈Cc∞\theta\in C_c^\infty equal to one on a neighborhood of KK. Extend the pairing to every smooth function aa by T(a)=T(θa)T(a)=T(\theta a). Two choices differ by a compact test vanishing near KK, so the extension is well-defined. Estimate (Q15), applied to θa\theta a and the finite product rule, bounds this extension by finitely many derivatives of aa on the fixed compact support of θ\theta. In particular a compactly supported distribution is tempered, because those seminorms are bounded by Schwartz seminorms.

Multiplication and differentiation are defined by (aT)(ϕ)=T(aϕ)(aT)(\phi)=T(a\phi) and (∂jT)(ϕ)=−T(∂jϕ)(\partial_jT)(\phi)=-T(\partial_j\phi). The ordinary product rule on tests immediately gives ∂j(aT)=(∂ja)T+a∂jT\partial_j(aT)=(\partial_j a)T+a\partial_jT. Iteration gives the full finite Leibniz formula. It also proves that differentiation does not enlarge support. The same derivative identity holds for the extended pairing of a compactly supported distribution with any smooth function: derivatives of the auxiliary cutoff vanish near its support and hence contribute zero. A distribution multiplied by a compact smooth cutoff inside its domain extends by zero to all of Rn\mathbb R^n, by pairing with that cutoff times each global test. Its support is compact, and (Q15) proves continuity of this extension.

Convolution with one compact factor. Let VV be supported in a compact set KK and let EE be any distribution on Rn\mathbb R^n. For a test ϕ\phi, define Bϕ(x)=Vy(ϕ(x+y)),(E∗V)(ϕ)=Ex(Bϕ(x)).(Q16) B_\phi(x)=V_y\bigl(\phi(x+y)\bigr),\qquad (E*V)(\phi)=E_x(B_\phi(x)). \tag{Q16} The inner pairing is the extension just defined. It is smooth in xx: the difference quotient of θ(y)ϕ(x+y)\theta(y)\phi(x+y) converges, with every derivative in yy uniformly on its fixed compact support, to the corresponding xx derivative, by the integral form of the ordinary Taylor remainder. Estimate (Q15) permits this limit inside VV. Iterating proves ∂xαBϕ=Vy(∂αϕ(x+y))\partial_x^\alpha B_\phi=V_y(\partial^\alpha\phi(x+y)). Its support is contained in the compact set supp⁡ϕ−K\operatorname{supp}\phi-K. Indeed outside that set the function of yy vanishes on a neighborhood of KK, and the same remains true for nearby xx.

For tests supported in one fixed compact set, this support bound is uniform, and the finite-order bound for VV bounds every derivative of BϕB_\phi by finitely many derivatives of ϕ\phi. Applying (Q15) for EE on the resulting fixed compact set proves that (Q16) is a distribution. Its definition involves no interchange of two unbounded distribution pairings.

Since ∂xjϕ(x+y)=∂yjϕ(x+y)\partial_{x_j}\phi(x+y)=\partial_{y_j}\phi(x+y), the definitions and the extended derivative identity give ∂j(E∗V)=(∂jE)∗V=E∗(∂jV),δ0∗V=V.(Q17) \partial_j(E*V)=(\partial_jE)*V=E*(\partial_jV), \qquad \delta_0*V=V. \tag{Q17} For example BϕB_\phi formed with ∂jV\partial_jV equals −B∂jϕ-B_{\partial_j\phi}, which proves the second equality after applying EE. The first follows from distributional differentiation in xx. Thus (Q17) also holds with any constant-coefficient polynomial P(D)P(D) in place of ∂j\partial_j, including complex coefficients.

The smooth terms. If f∈Cc∞f\in C_c^\infty, the function h(z)=Ex(f(z−x))h(z)=E_x(f(z-x)) is smooth. For zz in a compact set, all these tests have support in one compact set, and Taylor remainders converge in every test seminorm; (Q15) again passes each derivative through EE. This function represents E∗fE*f as defined in (Q16). To check the identity, pair with a compact test ϕ(z)\phi(z), pass the integral through EE using uniform Riemann sums in every test seminorm, and change variables y=z−xy=z-x. The resulting inner function is exactly ∫f(y)ϕ(x+y) dy\int f(y)\phi(x+y)\,dy.

More generally, suppose EE agrees with a smooth function ee on a neighborhood of O−KO-K, where OO is an open output set. Then E∗VE*V is smooth on OO and there has the formula (E∗V)(z)=Vy(e(z−y)),z∈O.(Q18) (E*V)(z)=V_y\bigl(e(z-y)\bigr),\qquad z\in O. \tag{Q18} For each compact subset O0⊂OO_0\subset O, choose the auxiliary cutoff in yy supported so close to KK that ee is smooth near O0−supp⁡θO_0-\operatorname{supp}\theta. The existence of such a neighborhood follows by compactness of O0−KO_0-K. The same finite-order and Taylor argument proves smoothness of the right side. To identify it as a distribution, insert a test supported in O0O_0 into (Q16): its inner function is supported where E=eE=e. Replace the outer pairing by integration, pass that integral through VV using compact-support Riemann sums and (Q15), and change variables z=x+yz=x+y. This gives (Q18) paired with the test. Smoothness is local, so the argument covers OO. In particular a globally smooth function convolved with a compactly supported distribution is smooth. If EE is smooth away from zero and OO is disjoint from KK, (Q18) applies after restricting to each compact O0⊂OO_0\subset O.

The Fourier distribution identity. Fourier inversion on Schwartz functions is proved in the earlier Fourier reading. Its unitary version has forward factor (2π)−n/2(2\pi)^{-n/2} and the same inverse factor. It preserves the Schwartz space continuously: differentiating under the integral inserts powers of the input variable; multiplication of the transform by a monomial is integration by parts in the input. Each resulting integral is bounded by finitely many weighted suprema of the input derivatives, since ⟨x⟩−n−1\langle x\rangle^{-n-1} is integrable. Applying the same calculation to the inverse gives the continuous inverse. Transposition therefore defines the Fourier transform and inverse on the continuous dual S′\mathcal S'. Integration by parts on tests proves F(P(D)T)=PFT\mathcal F(P(D)T)=P\mathcal F T and Fδ0=(2π)−n/2\mathcal F\delta_0=(2\pi)^{-n/2}. The bounded function mm defines a tempered distribution by integration, because the integral of a Schwartz test is bounded by a weighted supremum. These facts prove (Q14), including its normalization. Repeated integration by parts in the compact smooth χ\chi shows that r=(2π)−n/2F−1χr=(2\pi)^{-n/2}\mathcal F^{-1}\chi is Schwartz. The same test-space estimates justify the tempered limits of the cutoffs used above.

Local regularity for arbitrary distributions. Suppose P(D)u=fP(D)u=f is smooth near x0x_0, with no integrability assumption on uu. Choose ψ∈Cc∞\psi\in C_c^\infty supported in that neighborhood and equal to one on a smaller neighborhood of x0x_0. The compactly supported distribution v=ψuv=\psi u, extended by zero, satisfies P(D)v=ψf+[P(D),ψ]uP(D)v=\psi f+[P(D),\psi]u. Every term of the commutator has a positive-order derivative on ψ\psi, so its compact support is separated from a still smaller neighborhood of x0x_0. From (Q14), (Q16) and (Q17), v=E0∗(ψf)+E0∗([P(D),ψ]u)+r∗v.(Q19) v=E_0*(\psi f)+E_0*([P(D),\psi]u)+r*v. \tag{Q19} The first term is smooth by the compact smooth convolution proof. The second is smooth near x0x_0 by (Q18), since E0E_0 is smooth away from zero. The third is smooth because rr is smooth and vv has compact support. Thus u=vu=v is smooth near x0x_0. This proves the distributional step of the full characterization, rather than restricting it to square-integrable null solutions.

Further reading

Lars Hörmander, On the theory of general partial differential operators, Acta Mathematica 94 (1955), printed pages 225–228, Lemmas 3.9–3.11: quantitative polynomial growth and reciprocal derivatives.