Use the BiConjugateGradientStabilizedSolver class to solve large, non-symmetric sparse
systems of equations using the Stabilized Bi-Conjugate Gradient method (BiCGStab).
Convergence is more stable than the Conjugate Gradient Squared (CGS) method.

To solve a system of equations using this class, construct a new instance, passing the matrix of
the system to the constructor.
You can optionally set a preconditioner by setting the
Preconditioner property. By default, a JacobiPreconditioner
is used. Finally, the Solve(Vector) method solves the system.

Once the algorithm finished, the Status property indicates
whether the algorithm completed successfully.