the real result is available in closed form: sum(log(diag))
Then the test runs faster