Add the missing case of time_zero, to preserve the correct behavior.

parent 23215288
Pipeline #147138244 failed with stage
in 0 seconds
......@@ -165,7 +165,7 @@ module hamiltonian_elec_oct_m
type(lda_u_t) :: lda_u
integer :: lda_u_level
logical, private :: time_zero
logical, public :: time_zero
type(exchange_operator_t), public :: exxop
type(namespace_t), pointer :: namespace
......
......@@ -1236,7 +1236,7 @@ contains
hm%ep%vdw_forces(1:ks%gr%sb%dim, 1:ks%calc%geo%natoms) = CNST(0.0)
end if
if(ks%calc%time_present) then
if(ks%calc%time_present .or. hm%time_zero) then
call hamiltonian_elec_update(hm, ks%gr%mesh, namespace, time = ks%calc%time)
else
call hamiltonian_elec_update_pot(hm, ks%gr%mesh)
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment