RimuQMC / Rimu.jl / 12442692251

94%

push

github

Fix normalisation in `single_particle_density` (#300)

The `single_particle_density` was not normalized properly when computing
densities of specific components. This PR fixes it.

Before:
```julia
julia> addr = FermiFS2C((1, 1, 1, 1), (1, 0, 0, 0))
julia> single_particle_density(addr; component=2)
(1, 0, 0, 0)
julia> DVec(addr => 1)
julia> single_particle_density(dvec; component=2)
(5, 0, 0, 0) # <- multiplied by the total number of particles for some reason
```
after:
```julia
julia> addr = FermiFS2C((1, 1, 1, 1), (1, 0, 0, 0))
julia> single_particle_density(addr; component=2)
(1, 0, 0, 0)
julia> DVec(addr => 1)
julia> single_particle_density(dvec; component=2)
(1, 0, 0, 0)
```

---------

Co-authored-by: Joachim Brand <joachim.brand@gmail.com>

3 of 3 new or added lines in 1 file covered. (100.0%)

6682 of 7066 relevant lines covered (94.57%)

8679731.22 hits per line