.. _davidson: .. program:: davidson .. default-role:: option ============== Davidson_Utils ============== Abstract module for Davidson's diagonalization. It contains everything required for the Davidson algorithm, dressed or not. If a dressing is used, the dressing column should be defined and the :ref:`davidsondressed` module should be used. If no dressing is required, the :ref:`davidson` module should be used, and it has a default zero dressing vector. EZFIO parameters ---------------- .. option:: threshold_davidson Thresholds of Davidson's algorithm Default: 1.e-10 .. option:: n_states_diag Number of states to consider during the Davdison diagonalization Default: 4 .. option:: davidson_sze_max Number of micro-iterations before re-contracting Default: 8 .. option:: state_following If |true|, the states are re-ordered to match the input states Default: False .. option:: disk_based_davidson If |true|, disk space is used to store the vectors Default: False .. option:: distributed_davidson If |true|, use the distributed algorithm Default: True Providers --------- .. c:var:: ci_eigenvectors .. code:: text double precision, allocatable :: ci_electronic_energy (N_states_diag) double precision, allocatable :: ci_eigenvectors (N_det,N_states_diag) double precision, allocatable :: ci_eigenvectors_s2 (N_states_diag) File: :file:`diagonalize_CI.irp.f` Eigenvectors/values of the CI matrix .. c:var:: ci_eigenvectors_s2 .. code:: text double precision, allocatable :: ci_electronic_energy (N_states_diag) double precision, allocatable :: ci_eigenvectors (N_det,N_states_diag) double precision, allocatable :: ci_eigenvectors_s2 (N_states_diag) File: :file:`diagonalize_CI.irp.f` Eigenvectors/values of the CI matrix .. c:var:: ci_electronic_energy .. code:: text double precision, allocatable :: ci_electronic_energy (N_states_diag) double precision, allocatable :: ci_eigenvectors (N_det,N_states_diag) double precision, allocatable :: ci_eigenvectors_s2 (N_states_diag) File: :file:`diagonalize_CI.irp.f` Eigenvectors/values of the CI matrix .. c:var:: ci_energy .. code:: text double precision, allocatable :: ci_energy (N_states_diag) File: :file:`diagonalize_CI.irp.f` N_states lowest eigenvalues of the CI matrix .. c:var:: davidson_criterion .. code:: text character(64) :: davidson_criterion File: :file:`parameters.irp.f` Can be : [ energy | residual | both | wall_time | cpu_time | iterations ] .. c:var:: dressed_column_idx .. code:: text integer, allocatable :: dressed_column_idx (N_states) File: :file:`diagonalization_hs2_dressed.irp.f` Index of the dressed columns .. c:var:: n_states_diag .. code:: text integer :: n_states_diag File: :file:`input.irp.f` Number of states to consider during the Davdison diagonalization .. c:var:: nthreads_davidson .. code:: text integer :: nthreads_davidson File: :file:`davidson_parallel.irp.f` Number of threads for Davdison .. c:var:: psi_energy .. code:: text double precision, allocatable :: psi_energy (N_states) File: :file:`u0Hu0.irp.f` Energy of the current wave function .. c:var:: psi_energy_with_nucl_rep .. code:: text double precision, allocatable :: psi_energy_with_nucl_rep (N_states) File: :file:`u0Hu0.irp.f` Energy of the wave function with the nuclear repulsion energy. Subroutines / functions ----------------------- .. c:function:: davidson_collector .. code:: text subroutine davidson_collector(zmq_to_qp_run_socket, zmq_socket_pull, v0, s0, sze, N_st) File: :file:`davidson_parallel.irp.f` .. c:function:: davidson_converged .. code:: text subroutine davidson_converged(energy,residual,wall,iterations,cpu,N_st,converged) File: :file:`parameters.irp.f` True if the Davidson algorithm is converged .. c:function:: davidson_diag_hjj_sjj .. code:: text subroutine davidson_diag_hjj_sjj(dets_in,u_in,H_jj,s2_out,energies,dim_in,sze,N_st,N_st_diag,Nint,dressing_state,converged) File: :file:`diagonalization_hs2_dressed.irp.f` Davidson diagonalization with specific diagonal elements of the H matrix H_jj : specific diagonal H matrix elements to diagonalize de Davidson S2_out : Output : s^2 dets_in : bitmasks corresponding to determinants u_in : guess coefficients on the various states. Overwritten on exit dim_in : leftmost dimension of u_in sze : Number of determinants N_st : Number of eigenstates N_st_diag : Number of states in which H is diagonalized. Assumed > sze Initial guess vectors are not necessarily orthonormal .. c:function:: davidson_diag_hs2 .. code:: text subroutine davidson_diag_hs2(dets_in,u_in,s2_out,dim_in,energies,sze,N_st,N_st_diag,Nint,dressing_state) File: :file:`diagonalization_hs2_dressed.irp.f` Davidson diagonalization. dets_in : bitmasks corresponding to determinants u_in : guess coefficients on the various states. Overwritten on exit dim_in : leftmost dimension of u_in sze : Number of determinants N_st : Number of eigenstates Initial guess vectors are not necessarily orthonormal .. c:function:: davidson_pull_results .. code:: text subroutine davidson_pull_results(zmq_socket_pull, v_t, s_t, imin, imax, task_id) File: :file:`davidson_parallel.irp.f` .. c:function:: davidson_push_results .. code:: text subroutine davidson_push_results(zmq_socket_push, v_t, s_t, imin, imax, task_id) File: :file:`davidson_parallel.irp.f` .. c:function:: davidson_run_slave .. code:: text subroutine davidson_run_slave(thread,iproc) File: :file:`davidson_parallel.irp.f` .. c:function:: davidson_slave_inproc .. code:: text subroutine davidson_slave_inproc(i) File: :file:`davidson_parallel.irp.f` .. c:function:: davidson_slave_tcp .. code:: text subroutine davidson_slave_tcp(i) File: :file:`davidson_parallel.irp.f` .. c:function:: davidson_slave_work .. code:: text subroutine davidson_slave_work(zmq_to_qp_run_socket, zmq_socket_push, N_st, sze, worker_id) File: :file:`davidson_parallel.irp.f` .. c:function:: diagonalize_ci .. code:: text subroutine diagonalize_CI File: :file:`diagonalize_CI.irp.f` Replace the coefficients of the CI states by the coefficients of the eigenstates of the CI matrix .. c:function:: find_reference .. code:: text subroutine find_reference(thresh,n_ref,result) File: :file:`find_reference.irp.f` .. c:function:: h_s2_u_0_nstates_openmp .. code:: text subroutine H_S2_u_0_nstates_openmp(v_0,s_0,u_0,N_st,sze) File: :file:`u0Hu0.irp.f` Computes v_0 = H|u_0> and s_0 = S^2 |u_0> Assumes that the determinants are in psi_det istart, iend, ishift, istep are used in ZMQ parallelization. .. c:function:: h_s2_u_0_nstates_openmp_work .. code:: text subroutine H_S2_u_0_nstates_openmp_work(v_t,s_t,u_t,N_st,sze,istart,iend,ishift,istep) File: :file:`u0Hu0.irp.f` Computes v_t = H|u_t> and s_t = S^2 |u_t> Default should be 1,N_det,0,1 .. c:function:: h_s2_u_0_nstates_openmp_work_1 .. code:: text subroutine H_S2_u_0_nstates_openmp_work_1(v_t,s_t,u_t,N_st,sze,istart,iend,ishift,istep) File: :file:`u0Hu0.irp.f_template_468` Computes v_t = H|u_t> and s_t = S^2 |u_t> Default should be 1,N_det,0,1 .. c:function:: h_s2_u_0_nstates_openmp_work_2 .. code:: text subroutine H_S2_u_0_nstates_openmp_work_2(v_t,s_t,u_t,N_st,sze,istart,iend,ishift,istep) File: :file:`u0Hu0.irp.f_template_468` Computes v_t = H|u_t> and s_t = S^2 |u_t> Default should be 1,N_det,0,1 .. c:function:: h_s2_u_0_nstates_openmp_work_3 .. code:: text subroutine H_S2_u_0_nstates_openmp_work_3(v_t,s_t,u_t,N_st,sze,istart,iend,ishift,istep) File: :file:`u0Hu0.irp.f_template_468` Computes v_t = H|u_t> and s_t = S^2 |u_t> Default should be 1,N_det,0,1 .. c:function:: h_s2_u_0_nstates_openmp_work_4 .. code:: text subroutine H_S2_u_0_nstates_openmp_work_4(v_t,s_t,u_t,N_st,sze,istart,iend,ishift,istep) File: :file:`u0Hu0.irp.f_template_468` Computes v_t = H|u_t> and s_t = S^2 |u_t> Default should be 1,N_det,0,1 .. c:function:: h_s2_u_0_nstates_openmp_work_n_int .. code:: text subroutine H_S2_u_0_nstates_openmp_work_N_int(v_t,s_t,u_t,N_st,sze,istart,iend,ishift,istep) File: :file:`u0Hu0.irp.f_template_468` Computes v_t = H|u_t> and s_t = S^2 |u_t> Default should be 1,N_det,0,1 .. c:function:: h_s2_u_0_nstates_zmq .. code:: text subroutine H_S2_u_0_nstates_zmq(v_0,s_0,u_0,N_st,sze) File: :file:`davidson_parallel.irp.f` Computes v_0 = H|u_0> and s_0 = S^2 |u_0> n : number of determinants H_jj : array of S2_jj : array of .. c:function:: u_0_h_u_0 .. code:: text subroutine u_0_H_u_0(e_0,u_0,n,keys_tmp,Nint,N_st,sze) File: :file:`u0Hu0.irp.f` Computes e_0 = / n : number of determinants .. c:function:: zmq_get_n_states_diag .. code:: text integer function zmq_get_N_states_diag(zmq_to_qp_run_socket, worker_id) File: :file:`davidson_parallel.irp.f` Get N_states_diag from the qp_run scheduler .. c:function:: zmq_put_n_states_diag .. code:: text integer function zmq_put_N_states_diag(zmq_to_qp_run_socket,worker_id) File: :file:`davidson_parallel.irp.f` Put N_states_diag on the qp_run scheduler