### linear-algebra

#### Linear iterative solver vs direct solver stability

Is iterative solver more stable than direct solver based on LU factorization. For LU based solver, we always have cond(A) < cond(L) * cond(U), so factorization amplifies numerical inaccuracy. So in the event of an ill conditioned matrix A, whose condition number is large than 1e10, will it be better off using iterative solver for stability and numerical accuracy?

There are two factors involved into answering your question. 1) The physical system you are analyzing is ill-conditioned by itself (in mechanical terms, the system is pretty "loose", so its equilibrium state may vary greatly depending on just a small variation in the boundary conditions) 2) The physical system is OK, but the matrix has not been scaled properly before the solution process begins. In the first case, there isn't much you can do: the physical system is inherently unstable. Consider applying different boundary conditions, for example. In the second case, a preconditioner should be helpful; for example, the Jacobi preconditioner makes the matrix having all diagonal values equal to 1. In this case, the iterations are more likely to converge.The condition ratio of 1e10 shouldn't represent too much trouble, provided a preconditioning is used.

### Related Links

Using solve and/or linsolve with the symbolic toolbox in R2010b

Discarding null space of matrix A

Solving linear equations with numeric.js

accessing only left eigenvectors in armadillo c

conjugate gradient with a positive semidefinite matrix

Most efficient way to solve SEVERAL linear systems Ax=b with SMALL A (minimum 3x3 maximum 8x8)

Assigning a sparse matrix in Eigen

The significance of the rotation matrix after computing structure from motion

Changing the LAPACK implementation used by IDL linear algebra routines?

approximate polynomial by linear inequalities

Maximization in quadratic programming using CGAL

Efficient recalculation of weighted least squares regression when weights change

How does ATLAS tuning work?

Vector direction on 3d line

hlsl dot function

How do units flow through matrix operations?