PTEROSOR/QuantumPackage
10
0
Fork 0
mirror of https://github.com/QuantumPackage/qp2.git synced 2024-12-23 04:43:45 +01:00
Code Releases Activity

csf fix segfault

Browse Source
This commit is contained in:
Yann Damour 2022-03-22 10:49:39 +01:00
parent b70ae8dc58
commit 22b28fc774
1 changed files with 18 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
src/csf/sigma_vector.irp.f
View File

@ -52,12 +52,17 @@
ncfgpersomo = cfg_seniority_index(i+2) ncfgpersomo = cfg_seniority_index(i+2)
else else
k = 0 k = 0
if ((i+2+k) < size(cfg_seniority_index,1)) then
do while(cfg_seniority_index(i+2+k) < ncfgpersomo) do while(cfg_seniority_index(i+2+k) < ncfgpersomo)
k = k + 2 k = k + 2
if ((i+2+k) >= size(cfg_seniority_index,1)) then
exit
endif
ncfgpersomo = cfg_seniority_index(i+2+k) ncfgpersomo = cfg_seniority_index(i+2+k)
enddo enddo
endif endif
endif endif
endif
ncfg = ncfgpersomo - ncfgprev ncfg = ncfgpersomo - ncfgprev
if(i .EQ. 0 .OR. i .EQ. 1) then if(i .EQ. 0 .OR. i .EQ. 1) then
dimcsfpercfg = 1 dimcsfpercfg = 1
@ -75,11 +80,16 @@
ncfgprev = cfg_seniority_index(i+2) ncfgprev = cfg_seniority_index(i+2)
else else
k = 0 k = 0
if ((i+2+k) < size(cfg_seniority_index,1)) then
do while(cfg_seniority_index(i+2+k) < ncfgprev) do while(cfg_seniority_index(i+2+k) < ncfgprev)
k = k + 2 k = k + 2
if ((i+2+k) >= size(cfg_seniority_index,1)) then
exit
endif
ncfgprev = cfg_seniority_index(i+2+k) ncfgprev = cfg_seniority_index(i+2+k)
enddo enddo
endif endif
endif
enddo enddo
END_PROVIDER END_PROVIDER