carstenbauer / StableDQMC.jl

Numerical stabilization routines for determinant quantum Monte Carlo
https://carstenbauer.github.io/StableDQMC.jl/dev/
MIT License
29 stars 5 forks source link

Order fo operations in equation 13 #5

Closed ffreyer closed 5 years ago

ffreyer commented 5 years ago

Shouldn't it be X⁻¹ Dₚ⁻¹ rather than Dₚ⁻¹ X⁻¹ in equation 13? Because Dₚ⁻¹ X⁻¹ Dₚ X ≠ I

ffreyer commented 5 years ago

As far as I can tell the UDT implementation uses X⁻¹ Dₚ⁻¹

carstenbauer commented 5 years ago

Well spotted! Maybe you want to check Eq. 22 as well 😄