(git:f2099e5)
Loading...
Searching...
No Matches
mp2_setup.F
Go to the documentation of this file.
1!--------------------------------------------------------------------------------------------------!
2! CP2K: A general program to perform molecular dynamics simulations !
3! Copyright 2000-2026 CP2K developers group <https://cp2k.org> !
4! !
5! SPDX-License-Identifier: GPL-2.0-or-later !
6!--------------------------------------------------------------------------------------------------!
7
8! **************************************************************************************************
9!> \brief Types needed for MP2 calculations
10!> \par History
11!> 2011.05 created [Mauro Del Ben]
12!> \author MDB
13! **************************************************************************************************
14MODULE mp2_setup
15 USE bse_input, ONLY: read_bse_section
16 USE cp_eri_mme_interface, ONLY: cp_eri_mme_init_read_input
17 USE cp_log_handling, ONLY: cp_get_default_logger,&
18 cp_logger_type
19 USE cp_output_handling, ONLY: cp_p_file,&
20 cp_print_key_finished_output,&
21 cp_print_key_should_output,&
22 cp_print_key_unit_nr
23 USE cp_parser_methods, ONLY: read_float_object
24 USE input_constants, ONLY: &
25 do_eri_mme, do_potential_short, evgw, evgw0, mp2_method_direct, mp2_method_gpw, &
26 mp2_method_none, mp2_ri_optimize_basis, ri_mp2_laplace, ri_mp2_method_gpw, &
27 ri_rpa_method_gpw
28 USE input_section_types, ONLY: section_vals_get_subs_vals,&
29 section_vals_type,&
30 section_vals_val_get
31 USE kinds, ONLY: default_string_length,&
32 dp,&
33 max_line_length
34 USE machine, ONLY: m_flush
35 USE mathlib, ONLY: erfc_cutoff
36 USE mp2_types, ONLY: mp2_method_direct,&
37 mp2_method_gpw,&
38 mp2_method_none,&
39 mp2_ri_optimize_basis,&
40 mp2_type,&
41 ri_mp2_laplace,&
42 ri_mp2_method_gpw,&
43 ri_rpa_method_gpw
44#include "./base/base_uses.f90"
45
46 IMPLICIT NONE
47
48 PRIVATE
49
50 CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'mp2_setup'
51
52 PUBLIC :: read_mp2_section
53
54CONTAINS
55
56! **************************************************************************************************
57!> \brief ...
58!> \param input ...
59!> \param mp2_env ...
60! **************************************************************************************************
61 SUBROUTINE read_mp2_section(input, mp2_env)
62 TYPE(section_vals_type), POINTER :: input
63 TYPE(mp2_type), INTENT(INOUT) :: mp2_env
64
65 CHARACTER(len=*), PARAMETER :: routinen = 'read_mp2_section'
66
67 CHARACTER(LEN=default_string_length), &
68 DIMENSION(:), POINTER :: string_pointer
69 CHARACTER(LEN=max_line_length) :: error_message
70 INTEGER :: gw_sc, handle, i, i_special_kp, ival, &
71 unit_nr
72 INTEGER, DIMENSION(:), POINTER :: tmplist
73 LOGICAL :: do_mp2, do_opt_ri_basis, do_ri_mp2, &
74 do_ri_sos_mp2, do_rpa
75 REAL(kind=dp), DIMENSION(:), POINTER :: r_vals
76 TYPE(cp_logger_type), POINTER :: logger
77 TYPE(section_vals_type), POINTER :: bse_section, cphf_section, &
78 eri_mme_section, gw_section, &
79 low_scaling_section, mp2_section
80
81 CALL timeset(routinen, handle)
82 logger => cp_get_default_logger()
83
84 mp2_section => section_vals_get_subs_vals(input, "DFT%XC%WF_CORRELATION")
85
86 mp2_env%method = mp2_method_none
87
88 ! should come from input
89 CALL section_vals_val_get(mp2_section, "MEMORY", r_val=mp2_env%mp2_memory)
90 CALL section_vals_val_get(mp2_section, "SCALE_S", r_val=mp2_env%scale_S)
91 CALL section_vals_val_get(mp2_section, "SCALE_T", r_val=mp2_env%scale_T)
92 CALL section_vals_val_get(mp2_section, "GROUP_SIZE", i_val=mp2_env%mp2_num_proc)
93 CALL section_vals_val_get(mp2_section, "E_GAP", r_val=mp2_env%e_gap)
94 CALL section_vals_val_get(mp2_section, "E_RANGE", r_val=mp2_env%e_range)
95
96 CALL section_vals_val_get(mp2_section, "MP2%_SECTION_PARAMETERS_", l_val=do_mp2)
97 CALL section_vals_val_get(mp2_section, "MP2%BIG_SEND", l_val=mp2_env%direct_canonical%big_send)
98 IF (do_mp2) THEN
99 CALL check_method(mp2_env%method)
100 CALL section_vals_val_get(mp2_section, "MP2%METHOD", i_val=mp2_env%method)
101 END IF
102
103 CALL section_vals_val_get(mp2_section, "RI_RPA%_SECTION_PARAMETERS_", l_val=do_rpa)
104
105 !CALL section_vals_val_get(mp2_section, "RI_RPA%SIGMA_FUNCTIONAL",l_val=mp2_env%ri_rpa%do_sigma)
106 CALL section_vals_val_get(mp2_section, "RI_RPA%SIGMA_FUNCTIONAL", i_val=mp2_env%ri_rpa%sigma_param)
107
108 IF (do_rpa) THEN
109 CALL check_method(mp2_env%method)
110 mp2_env%method = ri_rpa_method_gpw
111 END IF
112 CALL section_vals_val_get(mp2_section, "RI_RPA%RPA_NUM_QUAD_POINTS", i_val=mp2_env%ri_rpa%rpa_num_quad_points)
113 CALL section_vals_val_get(mp2_section, "RI_RPA%NUM_INTEG_GROUPS", i_val=mp2_env%ri_rpa%rpa_num_integ_groups)
114 CALL section_vals_val_get(mp2_section, "RI_RPA%MM_STYLE", i_val=mp2_env%ri_rpa%mm_style)
115 CALL section_vals_val_get(mp2_section, "RI_RPA%MINIMAX_QUADRATURE", l_val=mp2_env%ri_rpa%minimax_quad)
116 CALL section_vals_val_get(mp2_section, "RI_RPA%ADMM", l_val=mp2_env%ri_rpa%do_admm)
117 CALL section_vals_val_get(mp2_section, "RI_RPA%SCALE_RPA", r_val=mp2_env%ri_rpa%scale_rpa)
118 mp2_env%ri_rpa%reuse_hfx = .false.
119
120 CALL section_vals_val_get(mp2_section, "RI_RPA%EXCHANGE_CORRECTION%_SECTION_PARAMETERS_", &
121 i_val=mp2_env%ri_rpa%exchange_correction)
122 CALL section_vals_val_get(mp2_section, "RI_RPA%EXCHANGE_CORRECTION%BLOCK_SIZE", &
123 i_val=mp2_env%ri_rpa%exchange_block_size)
124 CALL section_vals_val_get(mp2_section, "RI_RPA%EXCHANGE_CORRECTION%USE_HFX_IMPLEMENTATION", &
125 l_val=mp2_env%ri_rpa%use_hfx_implementation)
126
127 CALL section_vals_val_get(mp2_section, "RI_RPA%RSE", l_val=mp2_env%ri_rpa%do_rse)
128
129 CALL section_vals_val_get(mp2_section, "RI_RPA%PRINT_DGEMM_INFO", l_val=mp2_env%ri_rpa%print_dgemm_info)
130
131 NULLIFY (gw_section)
132 gw_section => section_vals_get_subs_vals(mp2_section, "RI_RPA%GW")
133 CALL section_vals_val_get(gw_section, "_SECTION_PARAMETERS_", &
134 l_val=mp2_env%ri_rpa%do_ri_g0w0)
135 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%SELF_CONSISTENCY", &
136 i_val=gw_sc)
137 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%CORR_MOS_OCC", &
138 i_val=mp2_env%ri_g0w0%corr_mos_occ)
139 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%CORR_MOS_VIRT", &
140 i_val=mp2_env%ri_g0w0%corr_mos_virt)
141 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%NUMB_POLES", &
142 i_val=mp2_env%ri_g0w0%num_poles)
143 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%NPARAM_PADE", &
144 i_val=mp2_env%ri_g0w0%nparam_pade)
145 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%ANALYTIC_CONTINUATION", &
146 i_val=mp2_env%ri_g0w0%analytic_continuation)
147 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%OMEGA_MAX_FIT", &
148 r_val=mp2_env%ri_g0w0%omega_max_fit)
149 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%CROSSING_SEARCH", &
150 i_val=mp2_env%ri_g0w0%crossing_search)
151 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%FERMI_LEVEL_OFFSET", &
152 r_val=mp2_env%ri_g0w0%fermi_level_offset)
153 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%HEDIN_SHIFT", &
154 l_val=mp2_env%ri_g0w0%do_hedin_shift)
155 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%EV_GW_ITER", &
156 i_val=mp2_env%ri_g0w0%iter_evGW)
157 IF (gw_sc == evgw) mp2_env%ri_g0w0%iter_evGW = 10
158 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%SC_GW0_ITER", &
159 i_val=mp2_env%ri_g0w0%iter_sc_GW0)
160 IF (gw_sc == evgw0) mp2_env%ri_g0w0%iter_sc_GW0 = 10
161 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%EPS_ITER", &
162 r_val=mp2_env%ri_g0w0%eps_iter)
163 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT_EXX", &
164 i_val=mp2_env%ri_g0w0%print_exx)
165 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT_SELF_ENERGY", &
166 l_val=mp2_env%ri_g0w0%print_self_energy)
167 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%RI_SIGMA_X", &
168 l_val=mp2_env%ri_g0w0%do_ri_Sigma_x)
169
170 NULLIFY (r_vals)
171 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%IC_CORR_LIST", &
172 r_vals=r_vals)
173 ALLOCATE (mp2_env%ri_g0w0%ic_corr_list(1)%array(SIZE(r_vals)))
174 mp2_env%ri_g0w0%ic_corr_list(1)%array(:) = r_vals(:)
175 NULLIFY (r_vals)
176 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%IC_CORR_LIST_BETA", &
177 r_vals=r_vals)
178 ALLOCATE (mp2_env%ri_g0w0%ic_corr_list(2)%array(SIZE(r_vals)))
179 mp2_env%ri_g0w0%ic_corr_list(2)%array(:) = r_vals(:)
180 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%GAMMA_ONLY_SIGMA", &
181 l_val=mp2_env%ri_g0w0%do_gamma_only_sigma)
182 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%UPDATE_XC_ENERGY", &
183 l_val=mp2_env%ri_g0w0%update_xc_energy)
184 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%REGULARIZATION_MINIMAX", &
185 r_val=mp2_env%ri_g0w0%regularization_minimax)
186 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%SOC", &
187 i_val=mp2_env%ri_g0w0%soc_type)
188 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%SOC_ENERGY_WINDOW", &
189 r_val=mp2_env%ri_g0w0%soc_energy_window)
190
191 NULLIFY (bse_section)
192 bse_section => section_vals_get_subs_vals(mp2_section, "RI_RPA%GW%BSE")
193 CALL read_bse_section(bse_section, mp2_env%bse)
194
195 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%IMAGE_CHARGE_MODEL", &
196 l_val=mp2_env%ri_g0w0%do_ic_model)
197 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%IC%PRINT_IC_LIST", &
198 l_val=mp2_env%ri_g0w0%print_ic_values)
199 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%IC%EPS_DIST", &
200 r_val=mp2_env%ri_g0w0%eps_dist)
201
202 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PERIODIC_CORRECTION", &
203 l_val=mp2_env%ri_g0w0%do_periodic)
204 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PERIODIC_CORRECTION%KPOINTS", &
205 i_vals=mp2_env%ri_g0w0%kp_grid)
206 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PERIODIC_CORRECTION%NUM_KP_GRIDS", &
207 i_val=mp2_env%ri_g0w0%num_kp_grids)
208 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PERIODIC_CORRECTION%EPS_KPOINT", &
209 r_val=mp2_env%ri_g0w0%eps_kpoint)
210 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PERIODIC_CORRECTION%MO_COEFF_GAMMA", &
211 l_val=mp2_env%ri_g0w0%do_mo_coeff_gamma)
212 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PERIODIC_CORRECTION%AVERAGE_DEGENERATE_LEVELS", &
213 l_val=mp2_env%ri_g0w0%do_average_deg_levels)
214 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PERIODIC_CORRECTION%EPS_EIGENVAL", &
215 r_val=mp2_env%ri_g0w0%eps_eigenval)
216 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PERIODIC_CORRECTION%EXTRAPOLATE_KPOINTS", &
217 l_val=mp2_env%ri_g0w0%do_extra_kpoints)
218 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PERIODIC_CORRECTION%DO_AUX_BAS_GW", &
219 l_val=mp2_env%ri_g0w0%do_aux_bas_gw)
220 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PERIODIC_CORRECTION%FRACTION_AUX_MOS", &
221 r_val=mp2_env%ri_g0w0%frac_aux_mos)
222 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PERIODIC_CORRECTION%NUM_OMEGA_POINTS", &
223 i_val=mp2_env%ri_g0w0%num_omega_points)
224
225 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%KPOINT_SET%NPOINTS", &
226 i_val=mp2_env%ri_g0w0%n_kp_in_kp_line)
227 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%KPOINT_SET%SPECIAL_POINT", &
228 n_rep_val=mp2_env%ri_g0w0%n_special_kp)
229 ALLOCATE (mp2_env%ri_g0w0%xkp_special_kp(3, mp2_env%ri_g0w0%n_special_kp))
230 DO i_special_kp = 1, mp2_env%ri_g0w0%n_special_kp
231 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%KPOINT_SET%SPECIAL_POINT", &
232 i_rep_val=i_special_kp, c_vals=string_pointer)
233 cpassert(SIZE(string_pointer(:), 1) == 3)
234 DO i = 1, 3
235 CALL read_float_object(string_pointer(i), &
236 mp2_env%ri_g0w0%xkp_special_kp(i, i_special_kp), &
237 error_message)
238 IF (len_trim(error_message) > 0) cpabort(trim(error_message))
239 END DO
240 END DO
241
242 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%KPOINTS_SELF_ENERGY", &
243 i_vals=mp2_env%ri_g0w0%kp_grid_Sigma)
244
245 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT%GW_DOS%LOWER_BOUND", &
246 r_val=mp2_env%ri_g0w0%dos_lower)
247 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT%GW_DOS%UPPER_BOUND", &
248 r_val=mp2_env%ri_g0w0%dos_upper)
249 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT%GW_DOS%STEP", &
250 r_val=mp2_env%ri_g0w0%dos_prec)
251 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT%GW_DOS%MIN_LEVEL_SPECTRAL", &
252 i_val=mp2_env%ri_g0w0%dos_min)
253 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT%GW_DOS%MAX_LEVEL_SPECTRAL", &
254 i_val=mp2_env%ri_g0w0%dos_max)
255 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT%GW_DOS%MIN_LEVEL_SELF_ENERGY", &
256 i_val=mp2_env%ri_g0w0%min_level_self_energy)
257 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT%GW_DOS%MAX_LEVEL_SELF_ENERGY", &
258 i_val=mp2_env%ri_g0w0%max_level_self_energy)
259 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT%GW_DOS%BROADENING", &
260 r_val=mp2_env%ri_g0w0%dos_eta)
261
262 mp2_env%ri_g0w0%do_kpoints_Sigma = mp2_env%ri_g0w0%n_special_kp > 0 .OR. &
263 mp2_env%ri_g0w0%kp_grid_Sigma(1)* &
264 mp2_env%ri_g0w0%kp_grid_Sigma(2)* &
265 mp2_env%ri_g0w0%kp_grid_Sigma(3) > 0
266
267 mp2_env%ri_g0w0%print_local_bandgap = btest(cp_print_key_should_output(logger%iter_info, &
268 mp2_section, "RI_RPA%GW%PRINT%LOCAL_BANDGAP"), &
269 cp_p_file)
270 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT%LOCAL_BANDGAP%ENERGY_WINDOW", &
271 r_val=mp2_env%ri_g0w0%energy_window_print_loc_bandgap)
272 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT%LOCAL_BANDGAP%ENERGY_SPACING", &
273 r_val=mp2_env%ri_g0w0%energy_spacing_print_loc_bandgap)
274 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT%LOCAL_BANDGAP%LDOS_THRESHOLD_GAP", &
275 r_val=mp2_env%ri_g0w0%ldos_thresh_print_loc_bandgap)
276 CALL section_vals_val_get(mp2_section, "RI_RPA%GW%PRINT%LOCAL_BANDGAP%STRIDE", &
277 i_vals=mp2_env%ri_g0w0%stride_loc_bandgap)
278
279 NULLIFY (low_scaling_section)
280 low_scaling_section => section_vals_get_subs_vals(mp2_section, "LOW_SCALING")
281 CALL section_vals_val_get(low_scaling_section, "_SECTION_PARAMETERS_", &
282 l_val=mp2_env%do_im_time)
283
284 CALL section_vals_val_get(low_scaling_section, "MEMORY_CUT", i_val=mp2_env%ri_rpa_im_time%cut_memory)
285 CALL section_vals_val_get(low_scaling_section, "MEMORY_INFO", l_val=mp2_env%ri_rpa_im_time%memory_info)
286 CALL section_vals_val_get(low_scaling_section, "EPS_FILTER", r_val=mp2_env%ri_rpa_im_time%eps_filter)
287 CALL section_vals_val_get(low_scaling_section, "EPS_STORAGE_SCALING", r_val=mp2_env%ri_rpa_im_time%eps_compress)
288 mp2_env%ri_rpa_im_time%eps_compress = mp2_env%ri_rpa_im_time%eps_compress*mp2_env%ri_rpa_im_time%eps_filter
289 mp2_env%ri_rpa_im_time%eps_compress = max(mp2_env%ri_rpa_im_time%eps_compress, 1.0e-16_dp)
290 CALL section_vals_val_get(low_scaling_section, "EPS_FILTER_FACTOR", r_val=mp2_env%ri_rpa_im_time%eps_filter_factor)
291
292 CALL section_vals_val_get(low_scaling_section, "DO_KPOINTS", &
293 l_val=mp2_env%ri_rpa_im_time%do_im_time_kpoints)
294 CALL section_vals_val_get(low_scaling_section, "KPOINTS", &
295 i_vals=mp2_env%ri_rpa_im_time%kp_grid)
296 mp2_env%ri_rpa_im_time%do_kpoints_from_Gamma = sum(mp2_env%ri_rpa_im_time%kp_grid) > 0
297 IF (mp2_env%ri_rpa_im_time%do_kpoints_from_Gamma) THEN
298 cpassert(mp2_env%ri_g0w0%do_kpoints_Sigma)
299 END IF
300 CALL section_vals_val_get(low_scaling_section, "KPOINT_WEIGHTS_W", &
301 i_val=mp2_env%ri_rpa_im_time%kpoint_weights_W_method)
302 CALL section_vals_val_get(low_scaling_section, "EXPONENT_TAILORED_WEIGHTS", &
303 r_val=mp2_env%ri_rpa_im_time%exp_tailored_weights)
304 CALL section_vals_val_get(low_scaling_section, "REGULARIZATION_RI", &
305 r_val=mp2_env%ri_rpa_im_time%regularization_RI)
306 CALL section_vals_val_get(low_scaling_section, "EPS_EIGVAL_S", &
307 r_val=mp2_env%ri_rpa_im_time%eps_eigval_S)
308 CALL section_vals_val_get(low_scaling_section, "EPS_EIGVAL_S_GAMMA", &
309 r_val=mp2_env%ri_rpa_im_time%eps_eigval_S_Gamma)
310 CALL section_vals_val_get(low_scaling_section, "MAKE_CHI_POS_DEFINITE", &
311 l_val=mp2_env%ri_rpa_im_time%make_chi_pos_definite)
312 CALL section_vals_val_get(low_scaling_section, "MAKE_OVERLAP_MAT_AO_POS_DEFINITE", &
313 l_val=mp2_env%ri_rpa_im_time%make_overlap_mat_ao_pos_definite)
314 CALL section_vals_val_get(low_scaling_section, "TRUNC_COULOMB_RI_X", &
315 l_val=mp2_env%ri_rpa_im_time%trunc_coulomb_ri_x)
316 CALL section_vals_val_get(low_scaling_section, "DO_EXTRAPOLATE_KPOINTS", &
317 l_val=mp2_env%ri_rpa_im_time%do_extrapolate_kpoints)
318 CALL section_vals_val_get(low_scaling_section, "REL_CUTOFF_TRUNC_COULOMB_RI_X", &
319 r_val=mp2_env%ri_rpa_im_time%rel_cutoff_trunc_coulomb_ri_x)
320 CALL section_vals_val_get(low_scaling_section, "K_MESH_G_FACTOR", &
321 i_val=mp2_env%ri_rpa_im_time%k_mesh_g_factor)
322
323 CALL section_vals_val_get(low_scaling_section, "KEEP_QUADRATURE", &
324 l_val=mp2_env%ri_rpa_im_time%keep_quad)
325
326 CALL section_vals_val_get(low_scaling_section, "MIN_BLOCK_SIZE", &
327 i_val=mp2_env%ri_rpa_im_time%min_bsize)
328
329 CALL section_vals_val_get(low_scaling_section, "MIN_BLOCK_SIZE_MO", &
330 i_val=mp2_env%ri_rpa_im_time%min_bsize_mo)
331
332 CALL section_vals_val_get(mp2_section, "RI_SOS_MP2%_SECTION_PARAMETERS_", l_val=do_ri_sos_mp2)
333 IF (do_ri_sos_mp2) THEN
334 CALL check_method(mp2_env%method)
335 mp2_env%method = ri_mp2_laplace
336 END IF
337 CALL section_vals_val_get(mp2_section, "RI_SOS_MP2%QUADRATURE_POINTS", i_val=mp2_env%ri_laplace%n_quadrature)
338 CALL section_vals_val_get(mp2_section, "RI_SOS_MP2%NUM_INTEG_GROUPS", i_val=mp2_env%ri_laplace%num_integ_groups)
339
340 CALL section_vals_val_get(mp2_section, "RI_MP2%_SECTION_PARAMETERS_", l_val=do_ri_mp2)
341 IF (do_ri_mp2) THEN
342 CALL check_method(mp2_env%method)
343 mp2_env%method = ri_mp2_method_gpw
344 END IF
345 CALL section_vals_val_get(mp2_section, "RI_MP2%BLOCK_SIZE", i_val=mp2_env%ri_mp2%block_size)
346 CALL section_vals_val_get(mp2_section, "RI_MP2%NUMBER_INTEGRATION_GROUPS", i_val=mp2_env%ri_mp2%number_integration_groups)
347 CALL section_vals_val_get(mp2_section, "RI_MP2%PRINT_DGEMM_INFO", l_val=mp2_env%ri_mp2%print_dgemm_info)
348
349 CALL section_vals_val_get(mp2_section, "RI%ROW_BLOCK", i_val=mp2_env%block_size_row)
350 CALL section_vals_val_get(mp2_section, "RI%COL_BLOCK", i_val=mp2_env%block_size_col)
351 CALL section_vals_val_get(mp2_section, "RI%CALC_COND_NUM", l_val=mp2_env%calc_PQ_cond_num)
352 CALL section_vals_val_get(mp2_section, "RI%DO_SVD", l_val=mp2_env%do_svd)
353 CALL section_vals_val_get(mp2_section, "RI%EPS_SVD", r_val=mp2_env%eps_svd)
354 CALL section_vals_val_get(mp2_section, "RI%ERI_BLKSIZE", i_vals=mp2_env%eri_blksize)
355 CALL section_vals_val_get(mp2_section, "RI%RI_METRIC%POTENTIAL_TYPE", i_val=mp2_env%ri_metric%potential_type)
356 CALL section_vals_val_get(mp2_section, "RI%RI_METRIC%OMEGA", r_val=mp2_env%ri_metric%omega)
357 CALL section_vals_val_get(mp2_section, "RI%RI_METRIC%EPS_RANGE", r_val=mp2_env%eps_range)
358 CALL section_vals_val_get(mp2_section, "RI%RI_METRIC%CUTOFF_RADIUS", r_val=mp2_env%ri_metric%cutoff_radius)
359 CALL section_vals_val_get(mp2_section, "RI%RI_METRIC%T_C_G_DATA", c_val=mp2_env%ri_metric%filename)
360 IF (mp2_env%ri_metric%potential_type == do_potential_short) THEN
361 CALL erfc_cutoff(mp2_env%eps_range, mp2_env%ri_metric%omega, mp2_env%ri_metric%cutoff_radius)
362 END IF
363
364 CALL section_vals_val_get(mp2_section, "RI%OPT_RI_BASIS%_SECTION_PARAMETERS_", l_val=do_opt_ri_basis)
365 IF (do_opt_ri_basis) THEN
366 CALL check_method(mp2_env%method)
367 mp2_env%method = mp2_ri_optimize_basis
368 END IF
369 CALL section_vals_val_get(mp2_section, "RI%OPT_RI_BASIS%DELTA_I_REL", &
370 r_val=mp2_env%ri_opt_param%DI_rel)
371 CALL section_vals_val_get(mp2_section, "RI%OPT_RI_BASIS%DELTA_RI", &
372 r_val=mp2_env%ri_opt_param%DRI)
373 CALL section_vals_val_get(mp2_section, "RI%OPT_RI_BASIS%EPS_DERIV", &
374 r_val=mp2_env%ri_opt_param%eps_step)
375 CALL section_vals_val_get(mp2_section, "RI%OPT_RI_BASIS%MAX_ITER", &
376 i_val=mp2_env%ri_opt_param%max_num_iter)
377 CALL section_vals_val_get(mp2_section, "RI%OPT_RI_BASIS%BASIS_SIZE", &
378 i_val=mp2_env%ri_opt_param%basis_quality)
379 NULLIFY (tmplist)
380 CALL section_vals_val_get(mp2_section, "RI%OPT_RI_BASIS%NUM_FUNC", &
381 i_vals=tmplist)
382 IF (tmplist(1) > 0) THEN
383 ALLOCATE (mp2_env%ri_opt_param%RI_nset_per_l(0:SIZE(tmplist) - 1))
384 mp2_env%ri_opt_param%RI_nset_per_l = 0
385 DO ival = 1, SIZE(tmplist)
386 mp2_env%ri_opt_param%RI_nset_per_l(ival - 1) = tmplist(ival)
387 END DO
388 END IF
389
390 CALL section_vals_val_get(mp2_section, "INTEGRALS%ERI_METHOD", i_val=mp2_env%eri_method)
391 CALL section_vals_val_get(mp2_section, "INTEGRALS%SIZE_LATTICE_SUM", i_val=mp2_env%mp2_gpw%size_lattice_sum)
392 CALL section_vals_val_get(mp2_section, "INTEGRALS%WFC_GPW%EPS_FILTER", r_val=mp2_env%mp2_gpw%eps_filter)
393 CALL section_vals_val_get(mp2_section, "INTEGRALS%WFC_GPW%EPS_GRID", r_val=mp2_env%mp2_gpw%eps_grid)
394 CALL section_vals_val_get(mp2_section, "INTEGRALS%WFC_GPW%CUTOFF", r_val=mp2_env%mp2_gpw%cutoff)
395 CALL section_vals_val_get(mp2_section, "INTEGRALS%WFC_GPW%REL_CUTOFF", r_val=mp2_env%mp2_gpw%relative_cutoff)
396 CALL section_vals_val_get(mp2_section, "INTEGRALS%WFC_GPW%PRINT_LEVEL", i_val=mp2_env%mp2_gpw%print_level)
397 CALL section_vals_val_get(mp2_section, "INTEGRALS%WFC_GPW%EPS_PGF_ORB_S", r_val=mp2_env%mp2_gpw%eps_pgf_orb_S)
398 CALL section_vals_val_get(mp2_section, "INTEGRALS%INTERACTION_POTENTIAL%POTENTIAL_TYPE", &
399 i_val=mp2_env%potential_parameter%potential_type)
400 CALL section_vals_val_get(mp2_section, "INTEGRALS%INTERACTION_POTENTIAL%TRUNCATION_RADIUS", &
401 r_val=mp2_env%potential_parameter%cutoff_radius)
402 CALL section_vals_val_get(mp2_section, "INTEGRALS%INTERACTION_POTENTIAL%POTENTIAL_DATA", &
403 c_val=mp2_env%potential_parameter%filename)
404 CALL section_vals_val_get(mp2_section, "INTEGRALS%INTERACTION_POTENTIAL%OMEGA", &
405 r_val=mp2_env%potential_parameter%omega)
406 CALL section_vals_val_get(mp2_section, "INTEGRALS%INTERACTION_POTENTIAL%SCALE_COULOMB", &
407 r_val=mp2_env%potential_parameter%scale_coulomb)
408 CALL section_vals_val_get(mp2_section, "INTEGRALS%INTERACTION_POTENTIAL%SCALE_LONGRANGE", &
409 r_val=mp2_env%potential_parameter%scale_longrange)
410
411 NULLIFY (mp2_env%eri_mme_param)
412 ALLOCATE (mp2_env%eri_mme_param)
413
414 IF (mp2_env%eri_method == do_eri_mme) THEN
415 eri_mme_section => section_vals_get_subs_vals(mp2_section, "INTEGRALS%ERI_MME")
416 CALL cp_eri_mme_init_read_input(eri_mme_section, mp2_env%eri_mme_param)
417 END IF
418
419 ! Set some parameters in case of P screening
420 mp2_env%not_last_hfx = .true.
421 mp2_env%p_screen = .true.
422
423 ! Set the CPHF section
424 CALL section_vals_val_get(mp2_section, "CANONICAL_GRADIENTS%FREE_HFX_BUFFER", l_val=mp2_env%ri_grad%free_hfx_buffer)
425 CALL section_vals_val_get(mp2_section, "CANONICAL_GRADIENTS%EPS_CANONICAL", r_val=mp2_env%ri_grad%eps_canonical)
426 CALL section_vals_val_get(mp2_section, "CANONICAL_GRADIENTS%DOT_PRODUCT_BLKSIZE", i_val=mp2_env%ri_grad%dot_blksize)
427 CALL section_vals_val_get(mp2_section, "CANONICAL_GRADIENTS%MAX_PARALLEL_COMM", i_val=mp2_env%ri_grad%max_parallel_comm)
428 cphf_section => section_vals_get_subs_vals(mp2_section, "CANONICAL_GRADIENTS%CPHF")
429 IF (ASSOCIATED(cphf_section)) THEN
430 CALL section_vals_val_get(cphf_section, "MAX_ITER", i_val=mp2_env%ri_grad%cphf_max_num_iter)
431 CALL section_vals_val_get(cphf_section, "EPS_CONV", r_val=mp2_env%ri_grad%cphf_eps_conv)
432 CALL section_vals_val_get(cphf_section, "SCALE_STEP_SIZE", r_val=mp2_env%ri_grad%scale_step_size)
433 CALL section_vals_val_get(cphf_section, "SOLVER_METHOD", i_val=mp2_env%ri_grad%z_solver_method)
434 CALL section_vals_val_get(cphf_section, "RESTART_EVERY", i_val=mp2_env%ri_grad%cphf_restart)
435 CALL section_vals_val_get(cphf_section, "ENFORCE_DECREASE", l_val=mp2_env%ri_grad%enforce_decrease)
436 CALL section_vals_val_get(cphf_section, "RECALC_RESIDUAL", l_val=mp2_env%ri_grad%recalc_residual)
437 CALL section_vals_val_get(cphf_section, "DO_POLAK_RIBIERE", l_val=mp2_env%ri_grad%polak_ribiere)
438 END IF
439
440 ! print some info about the MP2 parameters
441 unit_nr = cp_print_key_unit_nr(logger, mp2_section, "PRINT", &
442 extension=".mp2Log")
443 IF ((mp2_env%method /= mp2_method_none) .AND. unit_nr > 0) THEN
444 WRITE (unit_nr, '(T2,A)') ""
445 SELECT CASE (mp2_env%method)
446 CASE (mp2_method_direct)
447 WRITE (unit_nr, '(T2,A)') "MP2| using direct canonical MP2"
448 CASE (mp2_method_gpw)
449 WRITE (unit_nr, '(T2,A)') "MP2| using MP2 GPW style"
450 CASE (ri_mp2_method_gpw)
451 WRITE (unit_nr, '(T2,A)') "MP2| using RI-MP2-GPW"
452 CASE (ri_rpa_method_gpw)
453 WRITE (unit_nr, '(T2,A)') "RI-RPA| using GPW style"
454 CASE (ri_mp2_laplace)
455 WRITE (unit_nr, '(T2,A)') "RI-SOS-Laplace-MP2| using GPW style"
456 CASE (mp2_ri_optimize_basis)
457 WRITE (unit_nr, '(T2,A)') "MP2| Optimize RI auxiliary basis"
458 CASE DEFAULT
459 cpabort("Unknown MP2 method")
460 END SELECT
461 WRITE (unit_nr, '(T2,A)') ""
462 CALL m_flush(unit_nr)
463 END IF
464
465 CALL cp_print_key_finished_output(unit_nr, logger, mp2_section, &
466 "PRINT")
467
468 CALL timestop(handle)
469
470 END SUBROUTINE read_mp2_section
471
472! **************************************************************************************************
473!> \brief ...
474!> \param method ...
475! **************************************************************************************************
476 SUBROUTINE check_method(method)
477 INTEGER, INTENT(IN) :: method
478
479 CHARACTER(len=*), PARAMETER :: routinen = 'check_method'
480
481 INTEGER :: handle
482
483 CALL timeset(routinen, handle)
484
485 IF (method /= mp2_method_none) THEN
486 cpabort("Please use not more than one method to compute the correlation energy.")
487 END IF
488
489 CALL timestop(handle)
490
491 END SUBROUTINE check_method
492END MODULE mp2_setup
bse_input
Reads the &BSE input section into the BSE environment; called by read_mp2_section with the section un...
Definition bse_input.F:15
bse_input::read_bse_section
subroutine, public read_bse_section(bse_section, bse_env)
Reads every keyword of the &BSE section into the settings of the BSE environment.
Definition bse_input.F:38
cp_eri_mme_interface
Interface to Minimax-Ewald method for periodic ERI's to be used in CP2K.
Definition cp_eri_mme_interface.F:12
cp_eri_mme_interface::cp_eri_mme_init_read_input
subroutine, public cp_eri_mme_init_read_input(mme_section, param)
Read input and initialize parameter type.
Definition cp_eri_mme_interface.F:249
cp_log_handling
various routines to log and control the output. The idea is that decisions about where to log should ...
Definition cp_log_handling.F:41
cp_log_handling::cp_get_default_logger
type(cp_logger_type) function, pointer, public cp_get_default_logger()
returns the default logger
Definition cp_log_handling.F:234
cp_output_handling
routines to handle the output, The idea is to remove the decision of wheter to output and what to out...
Definition cp_output_handling.F:25
cp_output_handling::cp_print_key_unit_nr
integer function, public cp_print_key_unit_nr(logger, basis_section, print_key_path, extension, middle_name, local, log_filename, ignore_should_output, file_form, file_position, file_action, file_status, do_backup, on_file, is_new_file, mpi_io, fout)
...
Definition cp_output_handling.F:853
cp_output_handling::cp_print_key_finished_output
subroutine, public cp_print_key_finished_output(unit_nr, logger, basis_section, print_key_path, local, ignore_should_output, on_file, mpi_io)
should be called after you finish working with a unit obtained with cp_print_key_unit_nr,...
Definition cp_output_handling.F:1073
cp_output_handling::cp_p_file
integer, parameter, public cp_p_file
Definition cp_output_handling.F:103
cp_output_handling::cp_print_key_should_output
integer function, public cp_print_key_should_output(iteration_info, basis_section, print_key_path, used_print_key, first_time)
returns what should be done with the given property if btest(res,cp_p_store) then the property should...
Definition cp_output_handling.F:345
cp_parser_methods
Utility routines to read data from files. Kept as close as possible to the old parser because.
Definition cp_parser_methods.F:20
cp_parser_methods::read_float_object
elemental subroutine, public read_float_object(string, object, error_message)
Returns a floating point number read from a string including fraction like z1/z2.
Definition cp_parser_methods.F:1259
input_constants
collects all constants needed in input so that they can be used without circular dependencies
Definition input_constants.F:18
input_constants::evgw
integer, parameter, public evgw
Definition input_constants.F:1150
input_constants::mp2_method_direct
integer, parameter, public mp2_method_direct
Definition input_constants.F:1129
input_constants::mp2_ri_optimize_basis
integer, parameter, public mp2_ri_optimize_basis
Definition input_constants.F:1129
input_constants::do_eri_mme
integer, parameter, public do_eri_mme
Definition input_constants.F:1332
input_constants::ri_rpa_method_gpw
integer, parameter, public ri_rpa_method_gpw
Definition input_constants.F:1129
input_constants::ri_mp2_method_gpw
integer, parameter, public ri_mp2_method_gpw
Definition input_constants.F:1129
input_constants::mp2_method_gpw
integer, parameter, public mp2_method_gpw
Definition input_constants.F:1129
input_constants::mp2_method_none
integer, parameter, public mp2_method_none
Definition input_constants.F:1129
input_constants::ri_mp2_laplace
integer, parameter, public ri_mp2_laplace
Definition input_constants.F:1129
input_constants::evgw0
integer, parameter, public evgw0
Definition input_constants.F:1150
input_constants::do_potential_short
integer, parameter, public do_potential_short
Definition input_constants.F:855
input_section_types
objects that represent the structure of input sections and the data contained in an input section
Definition input_section_types.F:15
input_section_types::section_vals_get_subs_vals
recursive type(section_vals_type) function, pointer, public section_vals_get_subs_vals(section_vals, subsection_name, i_rep_section, can_return_null)
returns the values of the requested subsection
Definition input_section_types.F:735
input_section_types::section_vals_val_get
subroutine, public section_vals_val_get(section_vals, keyword_name, i_rep_section, i_rep_val, n_rep_val, val, l_val, i_val, r_val, c_val, l_vals, i_vals, r_vals, c_vals, explicit)
returns the requested value
Definition input_section_types.F:1052
kinds
Defines the basic variable types.
Definition kinds.F:23
kinds::max_line_length
integer, parameter, public max_line_length
Definition kinds.F:59
kinds::dp
integer, parameter, public dp
Definition kinds.F:34
kinds::default_string_length
integer, parameter, public default_string_length
Definition kinds.F:57
machine
Machine interface based on Fortran 2003 and POSIX.
Definition machine.F:17
machine::m_flush
subroutine, public m_flush(lunit)
flushes units if the &GLOBAL flag is set accordingly
Definition machine.F:124
mathlib
Collection of simple mathematical functions and subroutines.
Definition mathlib.F:15
mathlib::erfc_cutoff
subroutine, public erfc_cutoff(eps, omg, r_cutoff)
compute a truncation radius for the shortrange operator
Definition mathlib.F:1827
mp2_setup
Types needed for MP2 calculations.
Definition mp2_setup.F:14
mp2_setup::read_mp2_section
subroutine, public read_mp2_section(input, mp2_env)
...
Definition mp2_setup.F:62
mp2_types
Types needed for MP2 calculations.
Definition mp2_types.F:14
cp_log_handling::cp_logger_type
type of a logger, at the moment it contains just a print level starting at which level it should be l...
Definition cp_log_handling.F:140
input_section_types::section_vals_type
stores the values of a section
Definition input_section_types.F:127
mp2_types::mp2_type
Definition mp2_types.F:303