num::assume_diagonally_dominant

Defined in header <operator/properties.hpp>

template <class T>
with_law<T, law::diagonally_dominant> assume_diagonally_dominant(T A,
                                                                 call_site site = {});

Attach strict diagonal dominance, which Jacobi, Gauss--Seidel and LU without pivoting require. It is checked exactly, since the test reads the matrix once.