Very often, the primary computational bottleneck in the solution of a problem in applied geometry and other areas of applications is the numerical solution of one or several linear systems that arise throughout the computation:

Ax = b,

where A ∈ Rn×n and x,b ∈ Rn. Even a modest understanding of the properties of the underlying matrix may help in making a good choice of a solution method.