Fix bug when using single_particle_density with composite addresses…

… and integer vectors

src/strategies_and_params/poststepstrategy.jl

@@ -268,7 +268,7 @@ julia> single_particle_density(v; component=1)
 """
 function single_particle_density(dvec; component=0)
     K = keytype(dvec)
-    V = valtype(dvec)
+    V = float(valtype(dvec))
     M = num_modes(K)
     N = num_particles(K)
 
@@ -283,12 +283,12 @@ function single_particle_density(dvec; component=0)
 end
 
 function single_particle_density(add::SingleComponentFockAddress; component=0)
-    return Tuple(onr(add))
+    return float.(Tuple(onr(add)))
 end
 function single_particle_density(add::Union{CompositeFS,BoseFS2C}; component=0)
     if component == 0
-        return Float64.(Tuple(sum(onr(add))))
+        return float.(Tuple(sum(onr(add))))
     else
-        return Float64.(Tuple(onr(add)[component]))
+        return float.(Tuple(onr(add)[component]))
     end
 end

test/lomc.jl

@@ -518,6 +518,7 @@ using DataFrames
                 @test single_particle_density(add) == (1, 3, 3)
                 @test single_particle_density(add; component=1) == (1, 2, 3)
                 @test single_particle_density(add; component=2) == (0, 1, 0)
+                @test single_particle_density(DVec(add => 1); component=2) == (0, 1, 0)
             end
         end
     end