LCPQ/quantum_package
10
0
Fork 0
mirror of https://github.com/LCPQ/quantum_package synced 2025-01-03 18:16:12 +01:00
Code Issues Releases Wiki Activity

fragmentation bug in selection.irp.f

Browse Source
This commit is contained in:
Yann Garniron 2017-02-10 19:52:44 +01:00
parent 3b935080b4
commit 3fd719ed7e
1 changed files with 23 additions and 24 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

47
plugins/Full_CI_ZMQ/selection.irp.f
View File

@ -448,32 +448,31 @@ subroutine select_doubles(i_generator,hole_mask,particle_mask,fock_diag_tmp,E0,p
logical :: banned(mo_tot_num, mo_tot_num,2) logical :: banned(mo_tot_num, mo_tot_num,2)
logical :: bannedOrb(mo_tot_num, 2) logical :: bannedOrb(mo_tot_num, 2)
maskInd += 1 h2 = hole_list(i2,s2)
if(subset == 0 .or. mod(maskInd, fragment_count) == (subset-1)) then call apply_hole(pmask, s2,h2, mask, ok, N_int)
h2 = hole_list(i2,s2) banned = .false.
call apply_hole(pmask, s2,h2, mask, ok, N_int) do j=1,mo_tot_num
banned = .false. bannedOrb(j, 1) = .true.
do j=1,mo_tot_num bannedOrb(j, 2) = .true.
bannedOrb(j, 1) = .true. enddo
bannedOrb(j, 2) = .true. do s3=1,2
do i=1,N_particles(s3)
bannedOrb(particle_list(i,s3), s3) = .false.
enddo enddo
do s3=1,2 enddo
do i=1,N_particles(s3) if(s1 /= s2) then
bannedOrb(particle_list(i,s3), s3) = .false. if(monoBdo) then
enddo bannedOrb(h1,s1) = .false.
enddo
if(s1 /= s2) then
if(monoBdo) then
bannedOrb(h1,s1) = .false.
end if
if(monoAdo) then
bannedOrb(h2,s2) = .false.
monoAdo = .false.
end if
end if end if
if(monoAdo) then
bannedOrb(h2,s2) = .false.
monoAdo = .false.
end if
end if
maskInd += 1
if(subset == 0 .or. mod(maskInd, fragment_count) == (subset-1)) then
call spot_isinwf(mask, fullminilist, i_generator, fullinteresting(0), banned, fullMatch, fullinteresting) call spot_isinwf(mask, fullminilist, i_generator, fullinteresting(0), banned, fullMatch, fullinteresting)
if(fullMatch) cycle if(fullMatch) cycle