2015-12-09 23:20:34 +01:00
|
|
|
program qp_ao_ints
|
2016-02-19 00:20:28 +01:00
|
|
|
use omp_lib
|
2015-12-09 23:20:34 +01:00
|
|
|
implicit none
|
|
|
|
BEGIN_DOC
|
|
|
|
! Increments a running calculation to compute AO integrals
|
|
|
|
END_DOC
|
2016-02-19 00:20:28 +01:00
|
|
|
integer :: i
|
2015-12-09 23:20:34 +01:00
|
|
|
|
2016-02-19 00:20:28 +01:00
|
|
|
call switch_qp_run_to_master
|
|
|
|
|
2016-04-19 20:34:43 +02:00
|
|
|
zmq_context = f77_zmq_ctx_new ()
|
|
|
|
|
2015-12-09 23:20:34 +01:00
|
|
|
! Set the state of the ZMQ
|
|
|
|
zmq_state = 'ao_integrals'
|
|
|
|
|
|
|
|
! Provide everything needed
|
|
|
|
double precision :: integral, ao_bielec_integral
|
|
|
|
integral = ao_bielec_integral(1,1,1,1)
|
|
|
|
|
2016-07-26 18:29:52 +02:00
|
|
|
character*(64) :: state
|
|
|
|
call wait_for_state(zmq_state,state)
|
|
|
|
do while (state /= 'Stopped')
|
|
|
|
!$OMP PARALLEL DEFAULT(PRIVATE) PRIVATE(i)
|
|
|
|
i = omp_get_thread_num()
|
|
|
|
call ao_bielec_integrals_in_map_slave_tcp(i)
|
|
|
|
!$OMP END PARALLEL
|
|
|
|
call wait_for_state(zmq_state,state)
|
|
|
|
enddo
|
2015-12-09 23:20:34 +01:00
|
|
|
|
|
|
|
print *, 'Done'
|
|
|
|
end
|
|
|
|
|