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

Merge pull request #240 from Ydrnan/qp2_add
Some checks failed
continuous-integration/drone/push Build is failing
Details

Browse Source 
fix segfault restore_symmetry if all elem > thresh
This commit is contained in:
Anthony Scemama 2023-02-03 02:01:38 +01:00 committed by GitHub
commit 1273c19897
No known key found for this signature in database
GPG Key ID: 4AEE18F83AFDEB23
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