<p>Preconditioning is essential in iterative methods for solving linear systems. It is also the implicit objective in updating approximations of Jacobians in optimization methods, e.g.,&#xa0;in quasi-Newton methods. We study a nonclassic matrix condition number, the <InlineEquation ID="IEq4"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10589_2025_669_Article_IEq4.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="14" /> </InlineMediaObject> <EquationSource Format="TEX">\(\omega \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>ω</mi> </math></EquationSource> </InlineEquation>-condition number, <InlineEquation ID="IEq5"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10589_2025_669_Article_IEq5.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="14" /> </InlineMediaObject> <EquationSource Format="TEX">\(\omega \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>ω</mi> </math></EquationSource> </InlineEquation> for short. <InlineEquation ID="IEq6"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10589_2025_669_Article_IEq6.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="14" /> </InlineMediaObject> <EquationSource Format="TEX">\(\omega \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>ω</mi> </math></EquationSource> </InlineEquation> is the ratio of: the arithmetic and geometric means of the singular values, rather than the largest and smallest for the classical <InlineEquation ID="IEq7"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10589_2025_669_Article_IEq7.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="13" /> </InlineMediaObject> <EquationSource Format="TEX">\(\kappa \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>κ</mi> </math></EquationSource> </InlineEquation>-condition number. The simple functions in <InlineEquation ID="IEq8"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10589_2025_669_Article_IEq8.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="14" /> </InlineMediaObject> <EquationSource Format="TEX">\(\omega \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>ω</mi> </math></EquationSource> </InlineEquation> allow one to exploit first order optimality conditions. We use this fact to derive <i>explicit</i> formulae for (i) <InlineEquation ID="IEq9"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10589_2025_669_Article_IEq9.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="14" /> </InlineMediaObject> <EquationSource Format="TEX">\(\omega \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>ω</mi> </math></EquationSource> </InlineEquation>-optimal low rank updating of generalized Jacobians arising in the context of nonsmooth Newton methods; and (ii) <InlineEquation ID="IEq10"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10589_2025_669_Article_IEq10.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="14" /> </InlineMediaObject> <EquationSource Format="TEX">\(\omega \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>ω</mi> </math></EquationSource> </InlineEquation>-optimal preconditioners of special structure for iterative methods for linear systems. In the latter context, we analyze the benefits of <InlineEquation ID="IEq11"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10589_2025_669_Article_IEq11.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="14" /> </InlineMediaObject> <EquationSource Format="TEX">\(\omega \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>ω</mi> </math></EquationSource> </InlineEquation> for (a) improving the clustering of eigenvalues; (b) reducing the number of iterations; and (c) estimating the actual condition of a linear system. Moreover we show strong theoretical connections between the <InlineEquation ID="IEq12"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10589_2025_669_Article_IEq12.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="14" /> </InlineMediaObject> <EquationSource Format="TEX">\(\omega \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>ω</mi> </math></EquationSource> </InlineEquation>-optimal preconditioners and incomplete Cholesky factorizations, and highlight the misleading effects arising from the inverse invariance of <InlineEquation ID="IEq13"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10589_2025_669_Article_IEq13.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="13" /> </InlineMediaObject> <EquationSource Format="TEX">\(\kappa \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>κ</mi> </math></EquationSource> </InlineEquation>. Our results confirm the efficacy of using the <InlineEquation ID="IEq14"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10589_2025_669_Article_IEq14.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="14" /> </InlineMediaObject> <EquationSource Format="TEX">\(\omega \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>ω</mi> </math></EquationSource> </InlineEquation>-condition number compared to the <InlineEquation ID="IEq15"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10589_2025_669_Article_IEq15.gif" Format="GIF" Height="10" Rendition="HTML" Resolution="72" Type="Linedraw" Width="13" /> </InlineMediaObject> <EquationSource Format="TEX">\(\kappa \)</EquationSource> <EquationSource Format="MATHML"><math> <mi>κ</mi> </math></EquationSource> </InlineEquation>-condition number.</p>

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

The \(\omega \)-condition number: applications to preconditioning and low rank generalized Jacobian updating

  • Woosuk L. Jung,
  • David Torregrosa-Belén,
  • Henry Wolkowicz

摘要

Preconditioning is essential in iterative methods for solving linear systems. It is also the implicit objective in updating approximations of Jacobians in optimization methods, e.g., in quasi-Newton methods. We study a nonclassic matrix condition number, the \(\omega \) ω -condition number, \(\omega \) ω for short. \(\omega \) ω is the ratio of: the arithmetic and geometric means of the singular values, rather than the largest and smallest for the classical \(\kappa \) κ -condition number. The simple functions in \(\omega \) ω allow one to exploit first order optimality conditions. We use this fact to derive explicit formulae for (i) \(\omega \) ω -optimal low rank updating of generalized Jacobians arising in the context of nonsmooth Newton methods; and (ii) \(\omega \) ω -optimal preconditioners of special structure for iterative methods for linear systems. In the latter context, we analyze the benefits of \(\omega \) ω for (a) improving the clustering of eigenvalues; (b) reducing the number of iterations; and (c) estimating the actual condition of a linear system. Moreover we show strong theoretical connections between the \(\omega \) ω -optimal preconditioners and incomplete Cholesky factorizations, and highlight the misleading effects arising from the inverse invariance of \(\kappa \) κ . Our results confirm the efficacy of using the \(\omega \) ω -condition number compared to the \(\kappa \) κ -condition number.