You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
If we want to follow the original algorithm, option for Rademacher distribution should be available. It's possible to use xnp.sign(xnp.randn(...)), but adding xnp.rand to backends would make it less ambiguous.
馃殌 Feature Request
An implementation of the Hutch++ trace and diagonal estimator
Motivation
The$1/\sqrt{n}$ convergence of the Hutchinson estimator can be very slow and it would be nice to have some alternative algorithms here
Pitch
Implementation going in here:
https://github.com/wilson-labs/cola/tree/main/cola/algorithms
and to be wired up into the
cola.linalg.trace
andcola.linalg.diag
function via the method argumentsThe text was updated successfully, but these errors were encountered: