Skip to content

Commit

Permalink
Fix partial charge in the dense eigensolver (#229)
Browse files Browse the repository at this point in the history
  • Loading branch information
william-dawson committed Jan 22, 2024
1 parent 2f07c7c commit b1ca2f5
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion Source/Fortran/FermiOperatorModule.F90
Original file line number Diff line number Diff line change
Expand Up @@ -156,7 +156,7 @@ SUBROUTINE ComputeDenseFOE(H, ISQ, trace, K, inv_temp_in, &
energy_value = energy_value + tlist%DATA(II)%point_value
tlist%DATA(II)%point_value = 1.0_NTREAL
ELSE IF (tlist%DATA(II)%index_column .EQ. CEILING(trace)) THEN
occ_temp = CEILING(trace) - trace
occ_temp = trace - FLOOR(trace)
energy_value = energy_value + &
& occ_temp * tlist%DATA(II)%point_value
tlist%DATA(II)%point_value = occ_temp
Expand Down

0 comments on commit b1ca2f5

Please sign in to comment.