LCPQ/qp2
9
1
Fork 0
mirror of https://github.com/QuantumPackage/qp2.git synced 2024-07-27 21:07:23 +02:00
Code Releases Activity

fix segfault restore_symmetry if all elem > thresh

Browse Source
This commit is contained in:
Yann Damour 2023-02-02 20:43:16 +01:00
parent 6d5d8701c8
commit ce236f4d24
1 changed files with 5 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
src/utils/linear_algebra.irp.f
View File

@ -1736,7 +1736,7 @@ subroutine restore_symmetry(m,n,A,LDA,thresh)
! Symmetrize ! Symmetrize
i = 1 i = 1
do while( (i < sze).and.(-copy(i) > thresh) ) do while (i < sze)
pi = i pi = i
pf = i pf = i
val = 1d0/copy(i) val = 1d0/copy(i)
@ -1760,6 +1760,10 @@ subroutine restore_symmetry(m,n,A,LDA,thresh)
! Update i ! Update i
i = pf i = pf
endif endif
! Exit if the remaining elements are below thresh
if (-copy(i) <= thresh) exit
! Update i ! Update i
i = i + 1 i = i + 1
enddo enddo