carstenbauer / MonteCarlo.jl

Classical and quantum Monte Carlo simulations in Julia
https://carstenbauer.github.io/MonteCarlo.jl/dev/
Other
186 stars 18 forks source link

Fallbacks for complex matrices + optimize T inversion #79

Closed ffreyer closed 3 years ago

ffreyer commented 3 years ago

In relation to #78.

In isolated benchmarks with real 64x64 martices, rdivp! takes 8-9µs. The old method took 61-69µs and should cause GC to trigger more frequently. rdivp! may also be more accurate, as we avoid an extra lu decomposition and explicit inversion.