const SYTEN_SRDEF_TYPE SYTEN_BROADENING_EPSILON = std::pow(std::numeric_limits<SYTEN_SRDEF_TYPE>::epsilon(), 12./16.) |
Broadening threshold epsilon
to be used in inexact inversion \(
\frac{1}{x} \approx \frac{x}{x^2 + \epsilon} \).
Potentially multiplied by an overall scale factor.
Referenced by syten::STensorImpl::hadamard_inv_diffdiff(), and syten::invertDiagonal().