search

aesara-devs / aehmc

An HMC/NUTS implementation in Aesara
MIT License
33 stars 6 forks source link

Change `at.inv` to `at.reciprocal` #12

Closed zoj613 closed 3 years ago

zoj613 commented 3 years ago

Is there a reason why at.inv is used instead of 1 / inverse_mass_matrix when inverse_mass_matrix is 1-dimensional? https://github.com/aesara-devs/aehmc/blob/87b0c4d20c276d191522132c1f3f00fec70239b1/aehmc/metrics.py#L15

Does it have an advantage over the latter? Apologies if this is obvious to some but I was confused by the line when I read it.

brandonwillard commented 3 years ago

at.inv is the old, deprecated name for at.reciprocal.