Skip to content
Snippets Groups Projects
Commit c057615b authored by Philipp Rüssmann's avatar Philipp Rüssmann
Browse files

Fix a possible divide by zero in LDAU writeout

parent 1934962e
No related branches found
No related tags found
No related merge requests found
......@@ -207,6 +207,13 @@ contains
write (1337, '(12X,A,2F8.4,A,/)') 'Ueff and Jeff = ', ueff(i1), jeff(i1), ' Ry'
write (1337, '(12X,"Scaling factor for F^n :",F10.6,/)') scl
write (1337, '(12X," n F^n calculated F^n scaled ")')
if (scl==0) then
! to prevent division by zero
! we can do this here because scl is not use anymore in the remainder
scl = 1.0
end if
do lf = 2, lfmax, 2
write (1337, '(12X,I3,2(2X,F12.8," Ry"))') lf, fclmb(lf)/scl, fclmb(lf)
end do
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment