fix for case fock = nothing

pablosanjose · Nov 21, 2024 · 17cf46b · 17cf46b
1 parent 75e66fa
commit 17cf46b
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/meanfield.jl b/src/meanfield.jl
@@ -33,6 +33,7 @@ function meanfield(g::GreenFunction{T,E}, args...;
     Vf = sanitize_potential(fock)
     U = onsite === missing ? T(Vh(zero(SVector{E,T}))) : T(onsite)
     Uf = onsite === missing ? T(Vf(zero(SVector{E,T}))) : T(onsite)
+    Uf = fock === nothing ? zero(Uf) : Uf
     B = blocktype(hamiltonian(g))
     is_nambu_rotated´ = sanitize_nambu_rotated(namburotation, nambu, B)
     Q = sanitize_charge(charge, B, nambu, is_nambu_rotated´)