num::sqrt_lanczos

Defined in header <linear/eigen/lanczos.hpp>

template <class Op>
requires spd_operator<Op, vec<real>>
lanczos_action_result sqrt_lanczos(const Op &A, const vec<real> &right_hand_side,
                                   real tolerance = 1e-8, idx max_steps = 0);

Approximate \(A^{1/2}b\) for a symmetric positive-definite operator.