pub fn cg<T, A>(
operator: &A,
b: &Array1<T>,
config: &CgConfig<T::Real>,
) -> CgSolution<T>where
T: ComplexField,
A: LinearOperator<T>,Expand description
Solve Ax = b using the Conjugate Gradient method
Note: This method is only correct for symmetric positive definite matrices. For non-symmetric systems, use GMRES or BiCGSTAB instead.