1.3.1 Minor Bugfix:

In the Hamiltonian creation routine, the nonlocal (SF and PH) terms were added to the matrix with swapped indices, as a consequence of a legacy bugfix that has been superseded. The indices were only swapped in the NORMAL case. Fixed by swapping the indices back.
QcmPlab · Feb 7, 2024 · c965c20 · c965c20
1 parent fc0c742
commit c965c20
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/src/ED_NORMAL/stored/H_non_local.f90 b/src/ED_NORMAL/stored/H_non_local.f90
@@ -39,9 +39,9 @@
                  !
                  select case(MpiStatus)
                  case (.true.)
-                    call sp_insert_element(MpiComm,spH0nd,htmp,j,i)
+                    call sp_insert_element(MpiComm,spH0nd,htmp,i,j)
                  case (.false.)
-                    call sp_insert_element(spH0nd,htmp,j,i)
+                    call sp_insert_element(spH0nd,htmp,i,j)
                  end select
                  !
               endif
@@ -72,9 +72,9 @@
                  !
                  select case(MpiStatus)
                  case (.true.)
-                    call sp_insert_element(MpiComm,spH0nd,htmp,j,i)
+                    call sp_insert_element(MpiComm,spH0nd,htmp,i,j)
                  case (.false.)
-                    call sp_insert_element(spH0nd,htmp,j,i)
+                    call sp_insert_element(spH0nd,htmp,i,j)
                  end select
                  !
               endif