<p>We propose three floating-point validated algorithms to compute respectively fast inversion, Euclidean division and Hensel lifting over <InlineEquation ID="IEq1"> <EquationSource Format="TEX">\(\mathbb {C}[[x]][y]\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi mathvariant="double-struck">C</mi> <mo stretchy="false">[</mo> <mo stretchy="false">[</mo> <mi>x</mi> <mo stretchy="false">]</mo> <mo stretchy="false">]</mo> <mo stretchy="false">[</mo> <mi>y</mi> <mo stretchy="false">]</mo> </mrow> </math></EquationSource> </InlineEquation>. This is the second step (after&#xa0;Bréhard, Poteaux and Soudant in ISSAC 2023) towards a validated numerical Newton–Puiseux algorithm, and will also be useful towards a validated OM-algorithm over <InlineEquation ID="IEq2"> <EquationSource Format="TEX">\(\mathbb {C}[[x]][y]\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi mathvariant="double-struck">C</mi> <mo stretchy="false">[</mo> <mo stretchy="false">[</mo> <mi>x</mi> <mo stretchy="false">]</mo> <mo stretchy="false">]</mo> <mo stretchy="false">[</mo> <mi>y</mi> <mo stretchy="false">]</mo> </mrow> </math></EquationSource> </InlineEquation>. Our strategy is simply to first compute a floating-point approximation using the classical algorithm, then to a posteriori validate the result using a Newton-like fixed-point operator. We also provide a prototype Julia implementation of these algorithms and several examples.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

Validated Numerical Hensel Lifting

  • Florent Bréhard,
  • Jasmin Krüger,
  • Adrien Poteaux,
  • Arthur Vinciguerra

摘要

We propose three floating-point validated algorithms to compute respectively fast inversion, Euclidean division and Hensel lifting over \(\mathbb {C}[[x]][y]\) C [ [ x ] ] [ y ] . This is the second step (after Bréhard, Poteaux and Soudant in ISSAC 2023) towards a validated numerical Newton–Puiseux algorithm, and will also be useful towards a validated OM-algorithm over \(\mathbb {C}[[x]][y]\) C [ [ x ] ] [ y ] . Our strategy is simply to first compute a floating-point approximation using the classical algorithm, then to a posteriori validate the result using a Newton-like fixed-point operator. We also provide a prototype Julia implementation of these algorithms and several examples.