Workaround for single-orbital cases

JuliaAtoms · Jun 13, 2024 · 11c8972 · 11c8972
1 parent 2aa72f7
commit 11c8972
Showing 1 changed file with 3 additions and 0 deletions.
diff --git a/src/multi_configurational_equations.jl b/src/multi_configurational_equations.jl
@@ -164,6 +164,9 @@ function Base.diff(E::EM, orbitals::VO; verbosity=0) where {EM<:EnergyExpression
         isnothing(p) || ProgressMeter.next!(p)
         eq
     end
+    if !(equations isa AbstractVector)
+        equations = [equations]
+    end
 
     vd = collect(values(integrals.data))
     verbosity > 3 && @info "Integrals" integrals.data filter(∉(1:length(integrals.data)), vd)