Documentation

TauCeti.Analysis.Polynomial.Puiseux.RootDifference.Nonmonic

Laurent unit forms of nonmonic root differences #

Let the leading coefficient and discriminant of a nonmonic polynomial family be centered powers of a distinguished variable times analytic units. If its leading-coefficient-scaled roots extend analytically across the distinguished hyperplane, every difference of distinct root labels is a Laurent power times an analytic unit. The exponent is independent of the other parameters. Negative exponents permit poles; positive exponents permit collisions.

The normalized roots are supplied by a splitting of the integral normalization off the hyperplane. No splitting of the original family on the hyperplane is required, where its degree may drop or it may be nullified. Parameter-independent orders of root differences are used in analytic preparation of nonmonic polynomial families, including those whose branches have poles. The original coefficient and discriminant formulas are required only off the hyperplane: recomputing a discriminant after a degree drop need not preserve the formal discriminant. No constant-coefficient hypothesis is needed for differences.

References #

S. McCallum, A. ParusiΕ„ski, L. Paunescu, Validity proof of Lazard's method for CAD construction, J. Symbolic Comput. 92 (2019), Β§4, Corollary 4.2.

theorem TauCeti.exists_root_sub_eq_zpow_mul_unit {π•œ : Type u_1} {E : Type u_2} [RCLike π•œ] [NormedAddCommGroup E] [NormedSpace π•œ E] {d a c : β„•} {P : E Γ— π•œ β†’ Polynomial π•œ} {r s : Fin d β†’ E Γ— π•œ β†’ π•œ} {xβ‚€ : E} {yβ‚€ : π•œ} {u v : E Γ— π•œ β†’ π•œ} (hs : βˆ€ (i : Fin d), AnalyticAt π•œ (s i) (xβ‚€, yβ‚€)) (hu : AnalyticAt π•œ u (xβ‚€, yβ‚€)) (hu0 : u (xβ‚€, yβ‚€) β‰  0) (hv : AnalyticAt π•œ v (xβ‚€, yβ‚€)) (hv0 : v (xβ‚€, yβ‚€) β‰  0) (hnorm : βˆ€αΆ  (p : E Γ— π•œ) in nhds (xβ‚€, yβ‚€), p.2 β‰  yβ‚€ β†’ (P p).integralNormalization = ∏ i : Fin d, (Polynomial.X - Polynomial.C (s i p))) (hlead : βˆ€αΆ  (p : E Γ— π•œ) in nhds (xβ‚€, yβ‚€), p.2 β‰  yβ‚€ β†’ (P p).coeff d = (p.2 - yβ‚€) ^ c * v p) (hdiscr : βˆ€αΆ  (p : E Γ— π•œ) in nhds (xβ‚€, yβ‚€), p.2 β‰  yβ‚€ β†’ (P p).discr = (p.2 - yβ‚€) ^ a * u p) {i j : Fin d} (hscale : βˆ€αΆ  (p : E Γ— π•œ) in nhds (xβ‚€, yβ‚€), p.2 β‰  yβ‚€ β†’ (P p).coeff d * r i p = s i p ∧ (P p).coeff d * r j p = s j p) (hij : i β‰  j) :
βˆƒ (b : β„€) (w : E Γ— π•œ β†’ π•œ), AnalyticAt π•œ w (xβ‚€, yβ‚€) ∧ w (xβ‚€, yβ‚€) β‰  0 ∧ βˆ€αΆ  (p : E Γ— π•œ) in nhds (xβ‚€, yβ‚€), w p β‰  0 ∧ (p.2 β‰  yβ‚€ β†’ r i p - r j p = (p.2 - yβ‚€) ^ b * w p)

Differences of nonmonic root branches are Laurent powers times analytic units when the leading coefficient and discriminant are centered powers times analytic units. The scaled roots s i must extend analytically and split the integral normalization off the hyperplane. Only the selected pair of original roots must satisfy the scaling identities, and they may have poles there. The exponent is constant in the parameters, and the unit is nonzero on a neighborhood including the hyperplane.