solve_cgs

sectionproperties.analysis.solver.solve_cgs(k: csc_matrix, f: npt.NDArray[np.float64], m: LinearOperator, tol: float = 1e-05) → npt.NDArray[np.float64][source]

Solves a linear system using the CGS iterative method.

Parameters:

  • k (csc_matrix) – N x N matrix of the linear system

  • f (npt.NDArray[np.float64]) – N x 1 right hand side of the linear system

  • m (LinearOperator) – Preconditioner for the linear matrix approximating the inverse of k

  • tol (float) – Relative tolerance for the solver to achieve. Defaults to 1e-5.

Returns:

The solution vector to the linear system of equations

Raises:

RuntimeError – If the CGS iterative method does not converge

Return type:

npt.NDArray[np.float64]