In inverse_spd 0.5*(a + a') when computing LDLT to ensure that gradients are spread out and update tests.