Effective Condition Number and Its Applications

Zi-Cai Li,Hung-Tsai Huang,Jeng-Tzong Chen,Yimin Wei

DOI: https://doi.org/10.1007/s00607-010-0098-8

2010-01-01

Computing

Abstract:Consider the over-determined system Fx = b where \({{\bf F}\in\mathcal{R}^{m \times n}, m \geq n}\) and rank (F) = r ≤ n, the effective condition number is defined by \({{\rm Cond_{-}eff }= \frac {\|{\bf b}\|}{\sigma_r\|{\bf x}\|}}\), where the singular values of F are given as σ max = σ 1 ≥ σ 2 ≥ . . . ≥ σ r > 0 and σ r+1 = . . . = σ n = 0. For the general perturbed system (A+ΔA)(x+Δx) = b+Δb involving both ΔA and Δb, the new error bounds pertinent to Cond_eff are derived. Next, we apply the effective condition number to the solutions of Motz’s problem by the collocation Trefftz methods (CTM). Motz’s problem is the benchmark of singularity problems. We choose the general particular solutions \({v_L = \sum\nolimits_{k=0}^L d_k (\frac {r}{R_p})^{k+\frac 12}}\) \({{\rm cos}(k +\frac 12)\theta}\) with a radius parameter R p . The CTM is used to seek the coefficients d i by satisfying the boundary conditions only. Based on the new effective condition number, the optimal parameter R p = 1 is found. which is completely in accordance with the numerical results. However, if based on the traditional condition number Cond, the optimal choice of R p is misleading. Under the optimal choice R p = 1, the Cond grows exponentially as L increases, but Cond_eff is only linear. The smaller effective condition number explains well the very accurate solutions obtained. The error analysis in [14,15] and the stability analysis in this paper grant the CTM to become the most efficient and competent boundary method.