randomized trace estimation is a popular and well studied technique that
approximates the trace of a large-scale matrix $B$ by computing the average of
$x^T Bx$ for many samples of a random vector $X$. Often, $B$ is symmetric
positive definite (SPD) but a number of applications give ri