LCPQ/qp2
9
1
Fork 0
mirror of https://github.com/QuantumPackage/qp2.git synced 2024-12-21 11:03:29 +01:00
Code Releases Activity

Fixed bug in deterministic PT2 for buffer size

Browse Source
This commit is contained in:
Anthony Scemama 2020-06-13 00:05:11 +02:00
parent b4bbd01574
commit ebf49ce789
3 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
src/cipsi/cipsi.irp.f
View File

@ -88,7 +88,7 @@ subroutine run_cipsi
threshold_generators = threshold_generators_save threshold_generators = threshold_generators_save
SOFT_TOUCH threshold_generators SOFT_TOUCH threshold_generators
else else
call ZMQ_selection(to_select, pt2, variance, norm) call ZMQ_selection(to_select, pt2, variance, norm)
endif endif
do k=1,N_states do k=1,N_states

3
src/cipsi/run_selection_slave.irp.f
View File

@ -59,8 +59,7 @@ subroutine run_selection_slave(thread,iproc,energy)
read(task,*) subset, i_generator, N read(task,*) subset, i_generator, N
if(buf%N == 0) then if(buf%N == 0) then
! Only first time ! Only first time
bsize = min(N, (elec_alpha_num * (mo_num-elec_alpha_num))**2) call create_selection_buffer(N, N*2, buf)
call create_selection_buffer(bsize, bsize*2, buf)
buffer_ready = .True. buffer_ready = .True.
else else
if (N /= buf%N) then if (N /= buf%N) then

2
src/cipsi/zmq_selection.irp.f
View File

@ -16,6 +16,7 @@ subroutine ZMQ_selection(N_in, pt2, variance, norm)
! PROVIDE psi_det psi_coef N_det qp_max_mem N_states pt2_F s2_eig N_det_generators ! PROVIDE psi_det psi_coef N_det qp_max_mem N_states pt2_F s2_eig N_det_generators
N = max(N_in,1) N = max(N_in,1)
N = min(N, (elec_alpha_num * (mo_num-elec_alpha_num))**2)
if (.True.) then if (.True.) then
PROVIDE pt2_e0_denominator nproc PROVIDE pt2_e0_denominator nproc
PROVIDE psi_bilinear_matrix_columns_loc psi_det_alpha_unique psi_det_beta_unique PROVIDE psi_bilinear_matrix_columns_loc psi_det_alpha_unique psi_det_beta_unique
@ -78,6 +79,7 @@ subroutine ZMQ_selection(N_in, pt2, variance, norm)
stop 'Unable to add task to task server' stop 'Unable to add task to task server'
endif endif
endif endif
N = max(N_in,1)
ASSERT (associated(b%det)) ASSERT (associated(b%det))