Line data Source code
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 Defines control structures, which contain the parameters and the
10 : !> settings for the DFT-based calculations.
11 : ! **************************************************************************************************
12 : MODULE cp_control_types
13 : USE cp_fm_types, ONLY: cp_fm_release,&
14 : cp_fm_type
15 : USE eeq_input, ONLY: eeq_solver_type
16 : USE input_constants, ONLY: &
17 : do_full_density, rtp_bse_ham_G0W0, rtp_method_tddft, tblite_cli_born_kernel_auto, &
18 : tblite_cli_solution_state_gsolv, tblite_cli_solvation_none, tblite_guess_ceh, &
19 : tblite_guess_sad, tblite_mixer_damping_default, tblite_mixer_iterations_default, &
20 : tblite_mixer_max_weight_default, tblite_mixer_min_weight_default, &
21 : tblite_mixer_omega0_default, tblite_mixer_weight_factor_default, tblite_scc_mixer_auto, &
22 : tblite_solver_gvd, use_mom_ref_coac
23 : USE kinds, ONLY: default_path_length,&
24 : default_string_length,&
25 : dp
26 : USE pair_potential_types, ONLY: pair_potential_p_release,&
27 : pair_potential_p_type
28 : USE qs_cdft_types, ONLY: cdft_control_create,&
29 : cdft_control_release,&
30 : cdft_control_type
31 : USE smeagol_control_types, ONLY: smeagol_control_create,&
32 : smeagol_control_release,&
33 : smeagol_control_type
34 : USE xas_control, ONLY: xas_control_release,&
35 : xas_control_type
36 : USE xas_tdp_types, ONLY: xas_tdp_control_create,&
37 : xas_tdp_control_release,&
38 : xas_tdp_control_type
39 : #include "./base/base_uses.f90"
40 :
41 : IMPLICIT NONE
42 :
43 : PRIVATE
44 :
45 : !***************************************************************************************************
46 : !\brief Control parameters for hairy-probes DFT
47 : !***************************************************************************************************
48 : TYPE hairy_probes_type
49 : REAL(KIND=dp) :: alpha = 0.0_dp ! solution probes parameter
50 : REAL(KIND=dp) :: mu = 0.0_dp ! chemical potenatial of electrons in reservoir
51 : REAL(KIND=dp) :: T = 0.0_dp ! temperature of electrons in reservoir
52 : REAL(KIND=dp) :: eps_hp = 0.0_dp ! tolerance for accuracy checks on occupation numbers
53 : INTEGER :: natoms = 0, &
54 : last_ao = 0, &
55 : first_ao = 0
56 : INTEGER, DIMENSION(:), POINTER :: atom_ids => NULL() ! atom ids to which the probes are attached
57 : END TYPE hairy_probes_type
58 :
59 : ! **************************************************************************************************
60 : ! \brief Control parameters for pw grids
61 : ! **************************************************************************************************
62 : TYPE pw_grid_option
63 : LOGICAL :: spherical = .FALSE.
64 : LOGICAL :: fullspace = .FALSE.
65 : INTEGER, DIMENSION(2) :: distribution_layout = 0
66 : INTEGER :: blocked = 0
67 : END TYPE pw_grid_option
68 :
69 : ! **************************************************************************************************
70 : ! \brief parameters for EMD/RTP calculations involving MO projections
71 : ! **************************************************************************************************
72 : TYPE proj_mo_type
73 : INTEGER, DIMENSION(:), ALLOCATABLE :: ref_mo_index
74 : INTEGER :: ref_mo_spin = 1
75 : INTEGER :: ref_nlumo = 0
76 : LOGICAL :: sum_on_all_ref = .FALSE.
77 : INTEGER, DIMENSION(:), ALLOCATABLE :: td_mo_index
78 : REAL(dp), DIMENSION(:), ALLOCATABLE :: td_mo_occ
79 : INTEGER :: td_mo_spin = 1
80 : LOGICAL :: sum_on_all_td = .FALSE.
81 : CHARACTER(LEN=default_path_length) :: ref_mo_file_name = ""
82 : LOGICAL :: propagate_ref = .FALSE.
83 : TYPE(cp_fm_type), DIMENSION(:), &
84 : ALLOCATABLE :: mo_ref
85 : END TYPE proj_mo_type
86 :
87 : TYPE proj_mo_p_type
88 : TYPE(proj_mo_type), POINTER :: proj_mo => NULL()
89 : END TYPE proj_mo_p_type
90 :
91 : ! **************************************************************************************************
92 : ! \brief Control parameters for REAL_TIME_PROPAGATION calculations
93 : ! **************************************************************************************************
94 : TYPE rtp_control_type
95 : LOGICAL :: converged = .FALSE.
96 : REAL(KIND=dp) :: eps_ener = 0.0_dp
97 : INTEGER :: max_iter = 0
98 : INTEGER :: mat_exp = 0
99 : INTEGER :: propagator = 0
100 : LOGICAL :: fixed_ions = .FALSE.
101 : INTEGER :: rtp_method = rtp_method_tddft
102 : INTEGER :: rtbse_ham = rtp_bse_ham_G0W0
103 : INTEGER :: initial_wfn = 0
104 : REAL(dp) :: eps_exp = 0.0_dp
105 : LOGICAL :: initial_step = .FALSE.
106 : LOGICAL :: hfx_redistribute = .FALSE.
107 : INTEGER :: aspc_order = 0
108 : INTEGER :: sc_check_start = 0
109 : LOGICAL :: apply_wfn_mix_init_restart = .FALSE.
110 : LOGICAL :: apply_delta_pulse = .FALSE.
111 : LOGICAL :: apply_delta_pulse_mag = .FALSE.
112 : LOGICAL :: periodic = .FALSE.
113 : LOGICAL :: linear_scaling = .FALSE.
114 : LOGICAL :: write_restart = .FALSE.
115 : INTEGER :: mcweeny_max_iter = 0
116 : INTEGER :: acc_ref = 0
117 : REAL(dp) :: mcweeny_eps = 0.0_dp
118 : INTEGER, DIMENSION(3) :: delta_pulse_direction = 0
119 : REAL(KIND=dp) :: delta_pulse_scale = 0.0_dp
120 : LOGICAL :: velocity_gauge = .FALSE.
121 : REAL(KIND=dp), DIMENSION(3) :: field = 0.0_dp
122 : REAL(KIND=dp), DIMENSION(3) :: vec_pot = 0.0_dp
123 : LOGICAL :: nl_gauge_transform = .FALSE.
124 : LOGICAL :: is_proj_mo = .FALSE.
125 : TYPE(proj_mo_p_type), DIMENSION(:), &
126 : POINTER :: proj_mo_list => NULL()
127 : ! Switch to turn on moments trace saving
128 : LOGICAL :: save_local_moments = .FALSE.
129 : INTEGER :: moment_trace_ref_type = use_mom_ref_coac
130 : REAL(dp), DIMENSION(:), POINTER :: moment_trace_user_ref_point => NULL()
131 : REAL(dp) :: ft_damping = -1.0_dp
132 : REAL(dp) :: ft_t0 = 0.0_dp
133 : ! Index 1 : number of the element, Index 2 : coordinates
134 : ! e.g. if xx and xz elements are to be printed
135 : ! print_pol_elements(1,1) = 1
136 : ! print_pol_elements(1,2) = 1
137 : ! print_pol_elements(2,1) = 1
138 : ! print_pol_elements(2,2) = 3
139 : INTEGER, DIMENSION(:, :), POINTER :: print_pol_elements => NULL()
140 : LOGICAL :: pade_requested = .FALSE.
141 : REAL(dp) :: pade_e_min = 0.0_dp
142 : REAL(dp) :: pade_e_step = 0.02_dp
143 : REAL(dp) :: pade_e_max = 100.0_dp
144 : REAL(dp) :: pade_fit_e_min = 0.0_dp
145 : REAL(dp) :: pade_fit_e_max = 300.0_dp
146 : END TYPE rtp_control_type
147 :
148 : ! **************************************************************************************************
149 : ! \brief Control parameters for DFTB calculations
150 : ! **************************************************************************************************
151 : TYPE dftb_control_type
152 : LOGICAL :: self_consistent = .FALSE.
153 : LOGICAL :: orthogonal_basis = .FALSE.
154 : LOGICAL :: dispersion = .FALSE.
155 : INTEGER :: tblite_scc_mixer = tblite_scc_mixer_auto
156 : INTEGER :: dispersion_type = 0
157 : INTEGER :: tblite_mixer_iterations = tblite_mixer_iterations_default
158 : INTEGER :: tblite_mixer_memory = tblite_mixer_iterations_default
159 : INTEGER :: tblite_mixer_solver = tblite_solver_gvd
160 : LOGICAL :: dftb3_diagonal = .FALSE.
161 : LOGICAL :: hb_sr_damp = .FALSE.
162 : REAL(KIND=dp) :: hb_sr_para = 0.0_dp
163 : REAL(KIND=dp) :: tblite_mixer_damping = tblite_mixer_damping_default
164 : REAL(KIND=dp) :: tblite_mixer_max_weight = tblite_mixer_max_weight_default
165 : REAL(KIND=dp) :: tblite_mixer_min_weight = tblite_mixer_min_weight_default
166 : REAL(KIND=dp) :: tblite_mixer_omega0 = tblite_mixer_omega0_default
167 : REAL(KIND=dp) :: tblite_mixer_weight_factor = tblite_mixer_weight_factor_default
168 : REAL(KIND=dp) :: eps_disp = 0.0_dp
169 : REAL(KIND=dp) :: epscn = 0.0_dp
170 : REAL(KIND=dp) :: exp_pre = 0.0_dp
171 : REAL(KIND=dp) :: scaling = 0.0_dp
172 : REAL(KIND=dp) :: rcdisp = 0.0_dp
173 : REAL(KIND=dp), DIMENSION(3) :: sd3 = 0.0_dp
174 : REAL(KIND=dp), DIMENSION(4) :: sd3bj = 0.0_dp
175 : LOGICAL :: do_ewald = .FALSE.
176 : CHARACTER(LEN=default_path_length) :: sk_file_path = ""
177 : CHARACTER(LEN=default_path_length) :: sk_file_list = ""
178 : CHARACTER(LEN=default_string_length), &
179 : DIMENSION(:, :), POINTER :: sk_pair_list => NULL()
180 : CHARACTER(LEN=default_path_length) :: uff_force_field = ""
181 : CHARACTER(LEN=default_path_length) :: dispersion_parameter_file = ""
182 : END TYPE dftb_control_type
183 :
184 : ! **************************************************************************************************
185 : ! \brief Control parameters for xTB calculations
186 : ! **************************************************************************************************
187 : TYPE xtb_reference_cli_guess_type
188 : LOGICAL :: enabled = .FALSE.
189 : LOGICAL :: efield_active = .FALSE.
190 : LOGICAL :: grad = .FALSE.
191 : INTEGER :: method = tblite_guess_ceh
192 : INTEGER :: solver = tblite_solver_gvd
193 : REAL(KIND=dp), DIMENSION(3) :: efield = 0.0_dp
194 : REAL(KIND=dp) :: electronic_temperature_guess = 0.0_dp
195 : CHARACTER(LEN=default_path_length) :: input_file = ""
196 : CHARACTER(LEN=default_path_length) :: json_file = ""
197 : CHARACTER(LEN=default_string_length) :: input_format = "gen"
198 : END TYPE xtb_reference_cli_guess_type
199 :
200 : TYPE xtb_reference_cli_param_type
201 : LOGICAL :: enabled = .FALSE.
202 : LOGICAL :: method_explicit = .FALSE.
203 : INTEGER :: method = 0
204 : CHARACTER(LEN=default_path_length) :: input_file = ""
205 : CHARACTER(LEN=default_path_length) :: output_file = ""
206 : END TYPE xtb_reference_cli_param_type
207 :
208 : TYPE xtb_reference_cli_fit_type
209 : LOGICAL :: enabled = .FALSE.
210 : LOGICAL :: dry_run = .FALSE.
211 : CHARACTER(LEN=default_path_length) :: copy_file = ""
212 : CHARACTER(LEN=default_path_length) :: input_file = ""
213 : CHARACTER(LEN=default_path_length) :: param_file = ""
214 : END TYPE xtb_reference_cli_fit_type
215 :
216 : TYPE xtb_reference_cli_tagdiff_type
217 : LOGICAL :: enabled = .FALSE.
218 : LOGICAL :: fit = .FALSE.
219 : CHARACTER(LEN=default_path_length) :: actual_file = ""
220 : CHARACTER(LEN=default_path_length) :: reference_file = ""
221 : END TYPE xtb_reference_cli_tagdiff_type
222 :
223 : TYPE xtb_reference_cli_type
224 : LOGICAL :: enabled = .FALSE.
225 : LOGICAL :: check_energy = .TRUE.
226 : LOGICAL :: check_forces = .TRUE.
227 : LOGICAL :: check_virial = .FALSE.
228 : LOGICAL :: keep_files = .FALSE.
229 : LOGICAL :: stop_on_error = .FALSE.
230 : LOGICAL :: efield_active = .FALSE.
231 : LOGICAL :: solvation_active = .FALSE.
232 : INTEGER :: guess = tblite_guess_sad
233 : INTEGER :: solvation_born_kernel = tblite_cli_born_kernel_auto
234 : INTEGER :: solvation_model = tblite_cli_solvation_none
235 : INTEGER :: solvation_state = tblite_cli_solution_state_gsolv
236 : REAL(KIND=dp), DIMENSION(3) :: efield = 0.0_dp
237 : REAL(KIND=dp) :: electronic_temperature_guess = 0.0_dp
238 : REAL(KIND=dp) :: error_limit = 1.0E-8_dp
239 : CHARACTER(LEN=default_path_length) :: program_name = "tblite"
240 : CHARACTER(LEN=default_path_length) :: grad_file = ""
241 : CHARACTER(LEN=default_path_length) :: json_file = ""
242 : CHARACTER(LEN=default_path_length) :: post_processing = ""
243 : CHARACTER(LEN=default_path_length) :: post_processing_output_file = ""
244 : CHARACTER(LEN=default_path_length) :: restart_file = ""
245 : CHARACTER(LEN=default_path_length) :: solvation_solvent = ""
246 : CHARACTER(LEN=default_path_length) :: work_directory = "."
247 : CHARACTER(LEN=default_path_length) :: prefix = "tblite-reference"
248 : CHARACTER(LEN=default_string_length) :: input_format = "gen"
249 : TYPE(xtb_reference_cli_guess_type) :: guess_cli
250 : TYPE(xtb_reference_cli_param_type) :: param_cli
251 : TYPE(xtb_reference_cli_fit_type) :: fit_cli
252 : TYPE(xtb_reference_cli_tagdiff_type) :: tagdiff_cli
253 : END TYPE xtb_reference_cli_type
254 :
255 : TYPE xtb_control_type
256 : !
257 : INTEGER :: gfn_type = 1
258 : !
259 : LOGICAL :: do_ewald = .FALSE.
260 : LOGICAL :: do_tblite = .FALSE.
261 : LOGICAL :: tblite_mixer_damping_explicit = .FALSE.
262 : !
263 : INTEGER :: sto_ng = 0
264 : INTEGER :: h_sto_ng = 0
265 : INTEGER :: tblite_method = 0
266 : INTEGER :: tblite_scc_mixer = tblite_scc_mixer_auto
267 : INTEGER :: tblite_mixer_iterations = tblite_mixer_iterations_default
268 : INTEGER :: tblite_mixer_memory = tblite_mixer_iterations_default
269 : INTEGER :: tblite_mixer_solver = tblite_solver_gvd
270 : REAL(KIND=dp) :: tblite_accuracy = 1.0_dp
271 : REAL(KIND=dp) :: tblite_mixer_damping = tblite_mixer_damping_default
272 : REAL(KIND=dp) :: tblite_mixer_max_weight = tblite_mixer_max_weight_default
273 : REAL(KIND=dp) :: tblite_mixer_min_weight = tblite_mixer_min_weight_default
274 : REAL(KIND=dp) :: tblite_mixer_omega0 = tblite_mixer_omega0_default
275 : REAL(KIND=dp) :: tblite_mixer_weight_factor = tblite_mixer_weight_factor_default
276 : CHARACTER(LEN=default_path_length) :: tblite_param_file = ""
277 : !
278 : INTEGER :: vdw_type = -1
279 : CHARACTER(LEN=default_path_length) :: parameter_file_path = ""
280 : CHARACTER(LEN=default_path_length) :: parameter_file_name = ""
281 : CHARACTER(LEN=default_path_length) :: spinpol_param_file_name = ""
282 : !
283 : CHARACTER(LEN=default_path_length) :: dispersion_parameter_file = ""
284 : REAL(KIND=dp) :: epscn = 0.0_dp
285 : REAL(KIND=dp) :: rcdisp = 0.0_dp
286 : REAL(KIND=dp) :: s6 = 0.0_dp, s8 = 0.0_dp
287 : REAL(KIND=dp) :: a1 = 0.0_dp, a2 = 0.0_dp
288 : !
289 : REAL(KIND=dp) :: ks = 0.0_dp, kp = 0.0_dp, kd = 0.0_dp, ksp = 0.0_dp, k2sh = 0.0_dp
290 : REAL(KIND=dp) :: kg = 0.0_dp, kf = 0.0_dp
291 : REAL(KIND=dp) :: kcns = 0.0_dp, kcnp = 0.0_dp, kcnd = 0.0_dp
292 : REAL(KIND=dp) :: ken = 0.0_dp
293 : REAL(KIND=dp) :: ksen = 0.0_dp, kpen = 0.0_dp, kden = 0.0_dp
294 : REAL(KIND=dp) :: ben = 0.0_dp
295 : REAL(KIND=dp) :: kxr = 0.0_dp, kx2 = 0.0_dp
296 : REAL(KIND=dp) :: enscale = 0.0_dp
297 : !
298 : LOGICAL :: xb_interaction = .FALSE.
299 : LOGICAL :: do_nonbonded = .FALSE.
300 : LOGICAL :: do_spinpol = .FALSE.
301 : LOGICAL :: coulomb_interaction = .FALSE.
302 : LOGICAL :: coulomb_lr = .FALSE.
303 : LOGICAL :: tb3_interaction = .FALSE.
304 : LOGICAL :: check_atomic_charges = .FALSE.
305 : LOGICAL :: var_dipole = .FALSE.
306 : !
307 : REAL(KIND=dp) :: xb_radius = 0.0_dp
308 : REAL(KIND=dp) :: coulomb_sr_cut = 0.0_dp
309 : REAL(KIND=dp) :: coulomb_sr_eps = 0.0_dp
310 : !
311 : CHARACTER(LEN=default_string_length), &
312 : DIMENSION(:, :), POINTER :: kab_param => NULL()
313 : INTEGER, DIMENSION(:, :), POINTER :: kab_types => NULL()
314 : INTEGER :: kab_nval = 0
315 : REAL(KIND=dp), DIMENSION(:), POINTER :: kab_vals => NULL()
316 : !
317 : INTEGER, DIMENSION(:), POINTER :: spinpol_type => NULL()
318 : REAL(KIND=dp), DIMENSION(:, :), &
319 : POINTER :: spinpol_vals => NULL()
320 : !
321 : TYPE(pair_potential_p_type), POINTER :: nonbonded => NULL()
322 : REAL(KIND=dp) :: eps_pair = 0.0_dp
323 : REAL(KIND=dp), DIMENSION(:, :), &
324 : POINTER :: rcpair => NULL()
325 : !
326 : ! SRB terms
327 : REAL(KIND=dp) :: ksrb = 0.0_dp, esrb = 0.0_dp, gscal = 0.0_dp
328 : REAL(KIND=dp) :: c1srb = 0.0_dp, c2srb = 0.0_dp, shift = 0.0_dp
329 : !
330 : ! EN shift in EEQ (molecular=1 or crystaline=2)
331 : INTEGER :: enshift_type = 1
332 : TYPE(eeq_solver_type) :: eeq_sparam ! parameters for EEQ solver
333 : TYPE(xtb_reference_cli_type) :: reference_cli
334 : END TYPE xtb_control_type
335 :
336 : ! **************************************************************************************************
337 : ! \brief Control parameters for semi empirical calculations
338 : ! **************************************************************************************************
339 : TYPE semi_empirical_control_type
340 : LOGICAL :: orthogonal_basis = .FALSE.
341 : LOGICAL :: analytical_gradients = .FALSE.
342 : LOGICAL :: force_kdsod_EX = .FALSE.
343 : LOGICAL :: do_ewald = .FALSE., do_ewald_r3 = .FALSE., do_ewald_gks = .FALSE.
344 : INTEGER :: integral_screening = 0, periodic_type = 0
345 : INTEGER :: max_multipole = 0
346 : INTEGER :: ga_ncells = 0
347 : REAL(KIND=dp) :: delta = 0.0_dp
348 : ! Dispersion pair potential
349 : LOGICAL :: dispersion = .FALSE.
350 : REAL(KIND=dp) :: rcdisp = 0.0_dp
351 : REAL(KIND=dp) :: epscn = 0.0_dp
352 : REAL(KIND=dp), DIMENSION(3) :: sd3 = 0.0_dp
353 : CHARACTER(LEN=default_path_length) :: dispersion_parameter_file = ""
354 : ! Parameters controlling the evaluation of the integrals
355 : REAL(KIND=dp) :: cutoff_lrc = 0.0_dp, taper_lrc = 0.0_dp, range_lrc = 0.0_dp
356 : REAL(KIND=dp) :: cutoff_cou = 0.0_dp, taper_cou = 0.0_dp, range_cou = 0.0_dp
357 : REAL(KIND=dp) :: cutoff_exc = 0.0_dp, taper_exc = 0.0_dp, range_exc = 0.0_dp
358 : REAL(KIND=dp) :: taper_scr = 0.0_dp, range_scr = 0.0_dp
359 : END TYPE semi_empirical_control_type
360 :
361 : ! **************************************************************************************************
362 : ! \brief Control parameters for GAPW method within QUICKSTEP ***
363 : ! **************************************************************************************************
364 : TYPE gapw_control_type
365 : INTEGER :: basis_1c = 0
366 : REAL(KIND=dp) :: eps_fit = 0.0_dp, &
367 : eps_iso = 0.0_dp, &
368 : eps_Vrho0 = 0.0_dp, &
369 : eps_svd = 0.0_dp, &
370 : eps_cpc = 0.0_dp
371 : INTEGER :: ladd_rho0 = 0, &
372 : lmax_rho0 = 0, &
373 : lmax_sphere = 0, &
374 : quadrature = 0
375 : LOGICAL :: accurate_xcint = .FALSE.
376 : REAL(KIND=dp) :: aweights = 0.0_dp
377 : REAL(KIND=dp), DIMENSION(:), POINTER :: aw => NULL()
378 : LOGICAL :: alpha0_hard_from_input = .FALSE., &
379 : force_paw = .FALSE., &
380 : non_paw_atoms = .FALSE., &
381 : nopaw_as_gpw = .FALSE.
382 : REAL(KIND=dp) :: alpha0_hard = 0.0_dp
383 : REAL(KIND=dp) :: max_rad_local = 0.0_dp
384 : END TYPE gapw_control_type
385 :
386 : ! **************************************************************************************************
387 : ! \brief parameters for calculations involving a time dependent electric field
388 : ! **************************************************************************************************
389 : TYPE efield_type
390 : REAL(KIND=dp) :: actual_time = 0.0_dp
391 : REAL(KIND=dp), DIMENSION(:), POINTER :: polarisation => NULL()
392 : INTEGER :: envelop_id = 0
393 : REAL(KIND=dp), DIMENSION(:), POINTER :: envelop_r_vars => NULL()
394 : INTEGER, DIMENSION(:), POINTER :: envelop_i_vars => NULL()
395 : REAL(KIND=dp) :: strength = 0.0_dp
396 : REAL(KIND=dp) :: amplitude = 0.0_dp
397 : REAL(KIND=dp) :: phase_offset = 0.0_dp
398 : REAL(KIND=dp) :: wavelength = 0.0_dp
399 : REAL(KIND=dp), DIMENSION(3) :: vec_pot_initial = 0.0_dp
400 : END TYPE efield_type
401 :
402 : TYPE efield_p_type
403 : TYPE(efield_type), POINTER :: efield => NULL()
404 : END TYPE efield_p_type
405 :
406 : ! **************************************************************************************************
407 : ! \brief parameters for calculations involving a time dependent electric field
408 : ! **************************************************************************************************
409 : TYPE period_efield_type
410 : LOGICAL :: displacement_field = .FALSE.
411 : REAL(KIND=dp), DIMENSION(3) :: polarisation = 0.0_dp
412 : REAL(KIND=dp), DIMENSION(3) :: d_filter = 0.0_dp
413 : REAL(KIND=dp) :: strength = 0.0_dp
414 : REAL(KIND=dp), ALLOCATABLE, DIMENSION(:) :: strength_list
415 : INTEGER :: start_frame = 0
416 : INTEGER :: end_frame = -1
417 : END TYPE period_efield_type
418 :
419 : ! **************************************************************************************************
420 : ! \brief some parameters useful for mulliken_restraints
421 : ! **************************************************************************************************
422 : TYPE mulliken_restraint_type
423 : REAL(KIND=dp) :: strength = 0.0_dp
424 : REAL(KIND=dp) :: TARGET = 0.0_dp
425 : INTEGER :: natoms = 0
426 : INTEGER, POINTER, DIMENSION(:) :: atoms => NULL()
427 : END TYPE mulliken_restraint_type
428 :
429 : ! **************************************************************************************************
430 : ! \brief some parameters useful for ddapc_restraints
431 : ! **************************************************************************************************
432 : TYPE ddapc_restraint_type
433 : INTEGER :: ref_count = 0
434 : REAL(KIND=dp) :: strength = 0.0_dp
435 : REAL(KIND=dp) :: TARGET = 0.0_dp
436 : REAL(KIND=dp) :: ddapc_order_p = 0.0_dp
437 : INTEGER :: functional_form = 0
438 : INTEGER :: natoms = 0
439 : INTEGER, POINTER, DIMENSION(:) :: atoms => NULL()
440 : REAL(KIND=dp), POINTER, DIMENSION(:) :: coeff => NULL()
441 : INTEGER :: density_type = 0
442 : END TYPE ddapc_restraint_type
443 :
444 : ! **************************************************************************************************
445 : ! \brief some parameters useful for s2_restraints
446 : ! **************************************************************************************************
447 : TYPE s2_restraint_type
448 : REAL(KIND=dp) :: strength = 0.0_dp
449 : REAL(KIND=dp) :: TARGET = 0.0_dp
450 : REAL(KIND=dp) :: s2_order_p = 0.0_dp
451 : INTEGER :: functional_form = 0
452 : END TYPE s2_restraint_type
453 :
454 : ! **************************************************************************************************
455 : ! \brief some parameters useful for auxiliary density matrix method
456 : ! **************************************************************************************************
457 : TYPE admm_block_type
458 : INTEGER, DIMENSION(:), ALLOCATABLE :: list
459 : END TYPE admm_block_type
460 :
461 : TYPE admm_control_type
462 : REAL(KIND=dp) :: eps_filter = 0.0_dp
463 : INTEGER :: admm_type = 0
464 : INTEGER :: purification_method = 0
465 : INTEGER :: method = 0
466 : LOGICAL :: charge_constrain = .FALSE.
467 : INTEGER :: scaling_model = 0
468 : INTEGER :: aux_exch_func = 0
469 : LOGICAL :: aux_exch_func_param = .FALSE.
470 : REAL(KIND=dp), DIMENSION(3) :: aux_x_param = 0.0_dp
471 : TYPE(admm_block_type), DIMENSION(:), &
472 : ALLOCATABLE :: blocks
473 : END TYPE admm_control_type
474 :
475 : ! **************************************************************************************************
476 : ! \brief Parameters for external potential
477 : ! **************************************************************************************************
478 : TYPE expot_control_type
479 : LOGICAL :: read_from_cube = .FALSE.
480 : LOGICAL :: maxwell_solver = .FALSE.
481 : LOGICAL :: static = .FALSE.
482 : REAL(KIND=dp) :: scaling_factor = 0.0_dp
483 : END TYPE expot_control_type
484 :
485 : ! **************************************************************************************************
486 : ! \brief Parameters useful for Maxwell equation evaluation of external potential
487 : ! **************************************************************************************************
488 : TYPE maxwell_control_type
489 : LOGICAL :: log_test = .FALSE.
490 : INTEGER :: int_test = 0
491 : REAL(KIND=dp) :: real_test = 0.0_dp
492 : END TYPE maxwell_control_type
493 :
494 : ! **************************************************************************************************
495 : ! \brief Control parameters for a QUICKSTEP and KIM-GORDON calculation ***
496 : ! eps_pgf_orb: Cutoff value for the interaction of the primitive
497 : ! Gaussian-type functions (primitive basis functions).
498 : ! **************************************************************************************************
499 : TYPE qs_control_type
500 : INTEGER :: method_id = 0
501 : REAL(KIND=dp) :: eps_core_charge = 0.0_dp, &
502 : eps_kg_orb = 0.0_dp, &
503 : eps_pgf_orb = 0.0_dp, &
504 : eps_ppl = 0.0_dp, &
505 : eps_ppnl = 0.0_dp, &
506 : eps_rho_gspace = 0.0_dp, &
507 : eps_rho_rspace = 0.0_dp, &
508 : eps_filter_matrix = 0.0_dp, &
509 : eps_gvg_rspace = 0.0_dp, &
510 : progression_factor = 0.0_dp, &
511 : relative_cutoff = 0.0_dp
512 : LOGICAL :: do_almo_scf = .FALSE.
513 : LOGICAL :: do_ls_scf = .FALSE.
514 : LOGICAL :: do_kg = .FALSE.
515 : LOGICAL :: commensurate_mgrids = .FALSE.
516 : LOGICAL :: realspace_mgrids = .FALSE.
517 : LOGICAL :: gapw = .FALSE., gapw_xc = .FALSE., gpw = .FALSE., pao = .FALSE.
518 : LOGICAL :: lrigpw = .FALSE., rigpw = .FALSE.
519 : LOGICAL :: lri_optbas = .FALSE.
520 : LOGICAL :: ofgpw = .FALSE.
521 : LOGICAL :: dftb = .FALSE.
522 : LOGICAL :: xtb = .FALSE.
523 : LOGICAL :: semi_empirical = .FALSE.
524 : LOGICAL :: mulliken_restraint = .FALSE.
525 : LOGICAL :: ddapc_restraint = .FALSE.
526 : LOGICAL :: ddapc_restraint_is_spin = .FALSE.
527 : LOGICAL :: ddapc_explicit_potential = .FALSE.
528 : LOGICAL :: cdft = .FALSE.
529 : LOGICAL :: et_coupling_calc = .FALSE.
530 : LOGICAL :: s2_restraint = .FALSE.
531 : INTEGER :: do_ppl_method = 0
532 : INTEGER :: wf_interpolation_method_nr = 0
533 : INTEGER :: wf_extrapolation_order = 0
534 : INTEGER :: periodicity = 0
535 : REAL(KIND=dp) :: pairlist_radius = 0.0_dp
536 : REAL(KIND=dp) :: cutoff = 0.0_dp
537 : REAL(KIND=dp), DIMENSION(:), POINTER :: e_cutoff => NULL()
538 : TYPE(mulliken_restraint_type), &
539 : POINTER :: mulliken_restraint_control => NULL()
540 : TYPE(ddapc_restraint_type), &
541 : DIMENSION(:), POINTER :: ddapc_restraint_control => NULL()
542 : TYPE(cdft_control_type), POINTER :: cdft_control => NULL()
543 : TYPE(s2_restraint_type), POINTER :: s2_restraint_control => NULL()
544 : TYPE(dftb_control_type), POINTER :: dftb_control => NULL()
545 : TYPE(xtb_control_type), POINTER :: xtb_control => NULL()
546 : TYPE(semi_empirical_control_type), &
547 : POINTER :: se_control => NULL()
548 : TYPE(gapw_control_type), POINTER :: gapw_control => NULL()
549 : TYPE(pw_grid_option) :: pw_grid_opt = pw_grid_option()
550 : LOGICAL :: skip_load_balance_distributed = .FALSE.
551 : ! Types of subsystems for embedding
552 : LOGICAL :: ref_embed_subsys = .FALSE.
553 : LOGICAL :: cluster_embed_subsys = .FALSE.
554 : LOGICAL :: high_level_embed_subsys = .FALSE.
555 : LOGICAL :: dfet_embedded = .FALSE.
556 : LOGICAL :: dmfet_embedded = .FALSE.
557 : END TYPE qs_control_type
558 :
559 : ! **************************************************************************************************
560 : ! \brief Control parameters for the SCCS models
561 : ! **************************************************************************************************
562 : TYPE sccs_control_type
563 : LOGICAL :: sccs_activated = .FALSE., &
564 : saa_activated = .FALSE.
565 : INTEGER :: derivative_method = 0, &
566 : max_iter = 0, &
567 : method_id = 0
568 : REAL(KIND=dp) :: alpha_solvent = 0.0_dp, &
569 : beta = 0.0_dp, &
570 : beta_solvent = 0.0_dp, &
571 : delta_rho = 0.0_dp, &
572 : eps_sccs = 0.0_dp, &
573 : eps_scf = 0.0_dp, &
574 : epsilon_solvent = 0.0_dp, &
575 : gamma_solvent = 0.0_dp, &
576 : mixing = 0.0_dp, &
577 : rho_zero = 0.0_dp, &
578 : rho_max = 0.0_dp, &
579 : rho_min = 0.0_dp, &
580 : f0 = 0.0_dp, &
581 : delta_eta = 0.0_dp, &
582 : alpha_zeta = 0.0_dp, &
583 : delta_zeta = 0.0_dp, &
584 : R_solv = 0.0_dp
585 : END TYPE sccs_control_type
586 :
587 : ! **************************************************************************************************
588 : ! \brief Control parameters for simplified Tamm Dancoff approximation (sTDA)
589 : ! \par ATTRIBUTES
590 : ! \par NOTES
591 : ! **************************************************************************************************
592 : TYPE stda_control_type
593 : LOGICAL :: do_ewald = .FALSE.
594 : LOGICAL :: do_exchange = .FALSE.
595 : REAL(KIND=dp) :: hfx_fraction = 0.0_dp
596 : REAL(KIND=dp) :: eps_td_filter = 0.0_dp
597 : REAL(KIND=dp) :: mn_alpha = 0.0_dp
598 : REAL(KIND=dp) :: mn_beta = 0.0_dp
599 : REAL(KIND=dp) :: coulomb_sr_cut = 0.0_dp
600 : REAL(KIND=dp) :: coulomb_sr_eps = 0.0_dp
601 : END TYPE stda_control_type
602 :
603 : ! **************************************************************************************************
604 : ! \brief Control parameters for smeared occupation
605 : ! \par ATTRIBUTES
606 : ! \par NOTES
607 : ! **************************************************************************************************
608 : TYPE smeared_type
609 : REAL(KIND=dp), DIMENSION(:), POINTER :: fermia => NULL()
610 : REAL(KIND=dp), DIMENSION(:, :), POINTER :: fermib => NULL()
611 : END TYPE smeared_type
612 :
613 : ! **************************************************************************************************
614 : ! \brief Control parameters for the planar averaged Hartree potential
615 : ! **************************************************************************************************
616 : TYPE paep_control_type
617 : INTEGER :: surf_normal = 0
618 : END TYPE paep_control_type
619 :
620 : ! **************************************************************************************************
621 : ! \brief Control parameters for the planar counter charge density
622 : ! **************************************************************************************************
623 : TYPE pcc_control_type
624 : INTEGER :: surf_normal = 0
625 : REAL(KIND=dp) :: dist_edge = 0.0_dp, &
626 : gau_a = 0.0_dp, &
627 : gau_c = 0.0_dp, &
628 : charge = 0.0_dp
629 : END TYPE pcc_control_type
630 :
631 : ! **************************************************************************************************
632 : ! \brief Control parameters for a Time-Dependent DFT calculation.
633 : ! **************************************************************************************************
634 : TYPE tddfpt2_control_type
635 : !> compute TDDFPT excitation energies and oscillator strengths
636 : LOGICAL :: enabled = .FALSE.
637 : !> number of excited states to converge
638 : INTEGER :: nstates = 0
639 : !> maximal number of iterations to be performed
640 : INTEGER :: niters = 0
641 : !> maximal number of Krylov space vectors
642 : INTEGER :: nkvs = 0
643 : !> number of unoccupied (virtual) molecular orbitals to consider
644 : INTEGER :: nlumo = 0
645 : !> minimal number of MPI processes to be used per excited state
646 : INTEGER :: nprocs = 0
647 : !> type of kernel function/approximation to use
648 : INTEGER :: kernel = 0
649 : !> type of spin excitations to compute
650 : INTEGER :: spinflip = 0
651 : !> for full kernel, do we have HFX/ADMM
652 : LOGICAL :: do_hfx = .FALSE.
653 : LOGICAL :: do_admm = .FALSE.
654 : !> for full kernel, do we have short-range/long-range HFX and/or Kxc potential
655 : LOGICAL :: do_hfxsr = .FALSE.
656 : LOGICAL :: hfxsr_re_int = .TRUE.
657 : INTEGER :: hfxsr_primbas = 0
658 : LOGICAL :: do_hfxlr = .FALSE.
659 : REAL(KIND=dp) :: hfxlr_rcut = 0.0_dp, hfxlr_scale = 0.0_dp
660 : LOGICAL :: do_exck = .FALSE.
661 : !> options used in sTDA calculation (Kernel)
662 : TYPE(stda_control_type) :: stda_control = stda_control_type()
663 : !> algorithm to correct orbital energies
664 : INTEGER :: oe_corr = 0
665 : !> eigenvalue shifts
666 : REAL(KIND=dp) :: ev_shift = 0.0_dp, eos_shift = 0.0_dp
667 : !> active orbitals
668 : INTEGER, DIMENSION(2) :: nactive = -1
669 : !> target accuracy
670 : REAL(kind=dp) :: conv = 0.0_dp
671 : !> the smallest excitation amplitude to print
672 : REAL(kind=dp) :: min_excitation_amplitude = 0.0_dp
673 : !> threshold which controls when two wave functions considered to be orthogonal:
674 : !> maxabs(Ci^T * S * Cj) <= orthogonal_eps
675 : REAL(kind=dp) :: orthogonal_eps = 0.0_dp
676 : !> read guess wave functions from restart file if exists
677 : LOGICAL :: is_restart = .FALSE.
678 : !> compute triplet excited states using spin-unpolarised molecular orbitals
679 : LOGICAL :: rks_triplets = .FALSE.
680 : !> local resolution of identity for Coulomb contribution
681 : LOGICAL :: do_lrigpw = .FALSE.
682 : !> smeared occupation
683 : LOGICAL :: do_smearing = .FALSE.
684 : !> dynamical correlation
685 : LOGICAL :: do_bse = .FALSE.
686 : LOGICAL :: do_bse_w_only = .FALSE.
687 : LOGICAL :: do_bse_gw_only = .FALSE.
688 : ! automatic generation of auxiliary basis for LRI-TDDFT
689 : INTEGER :: auto_basis_p_lri_aux = 1
690 : !> use symmetric definition of ADMM Kernel correction
691 : LOGICAL :: admm_symm = .FALSE.
692 : !> Use/Ignore possible ADMM Kernel XC correction
693 : LOGICAL :: admm_xc_correction = .FALSE.
694 : ! Compute exciton descriptors
695 : LOGICAL :: do_exciton_descriptors = .FALSE.
696 : LOGICAL :: do_directional_exciton_descriptors = .FALSE.
697 : !
698 : ! DIPOLE_MOMENTS subsection
699 : !
700 : ! form of the dipole operator used to compute oscillator strengths
701 : INTEGER :: dipole_form = 0
702 : !> type of the reference point used for calculation of electrostatic dipole moments
703 : INTEGER :: dipole_reference = 0
704 : !> user-defined reference point
705 : REAL(kind=dp), DIMENSION(:), POINTER :: dipole_ref_point => NULL()
706 : !
707 : ! SOC subsection
708 : LOGICAL :: do_soc = .FALSE.
709 : !
710 : ! MGRID subsection
711 : !
712 : !> number of plain-wave grids
713 : INTEGER :: mgrid_ngrids = 0
714 : !> create commensurate grids (progression factor and cutoff values of sub-grids will be ignored)
715 : LOGICAL :: mgrid_commensurate_mgrids = .FALSE.
716 : !> signals that MGRID section has been explicitly given. Other mgrid_* variables
717 : !> are not initialised when it is equal to .FALSE. as in this case the default
718 : !> set of plain-wave grids will be used
719 : LOGICAL :: mgrid_is_explicit = .FALSE.
720 : !> same as qs_control%realspace_mgrids
721 : LOGICAL :: mgrid_realspace_mgrids = .FALSE.
722 : !> do not perform load balancing
723 : LOGICAL :: mgrid_skip_load_balance = .FALSE.
724 : !> cutoff value at the finest grid level
725 : REAL(kind=dp) :: mgrid_cutoff = 0.0_dp
726 : !> cutoff at the next grid level will be smaller then the cutoff
727 : !> at the current grid by this number of times
728 : REAL(kind=dp) :: mgrid_progression_factor = 0.0_dp
729 : !> cutoff that determines to which grid a particular Gaussian function will be mapped
730 : REAL(kind=dp) :: mgrid_relative_cutoff = 0.0_dp
731 : !> manually provided the list of cutoff values for each grid level
732 : !> (when it is null(), the cutoff values will be assigned automatically)
733 : REAL(kind=dp), DIMENSION(:), POINTER :: mgrid_e_cutoff => NULL()
734 : !> Parameter for smeared occupation TDA
735 : TYPE(smeared_type), DIMENSION(:), POINTER :: smeared_occup => NULL()
736 : END TYPE tddfpt2_control_type
737 :
738 : ! **************************************************************************************************
739 : !> \brief
740 : ! **************************************************************************************************
741 : TYPE rixs_control_type
742 :
743 : LOGICAL :: enabled = .FALSE.
744 : INTEGER :: core_states = 0
745 : INTEGER :: valence_states = 0
746 :
747 : TYPE(tddfpt2_control_type), POINTER :: tddfpt2_control => NULL()
748 : TYPE(xas_tdp_control_type), POINTER :: xas_tdp_control => NULL()
749 :
750 : END TYPE rixs_control_type
751 :
752 : ! **************************************************************************************************
753 : ! \brief Control parameters for a DFT calculation
754 : ! \par History
755 : ! 10.2019 added variables related to surface dipole correction [Soumya Ghosh]
756 : ! **************************************************************************************************
757 : TYPE dft_control_type
758 : TYPE(admm_control_type), POINTER :: admm_control => NULL()
759 : TYPE(period_efield_type), POINTER :: period_efield => NULL()
760 : TYPE(qs_control_type), POINTER :: qs_control => NULL()
761 : TYPE(rtp_control_type), POINTER :: rtp_control => NULL()
762 : TYPE(sccs_control_type), POINTER :: sccs_control => NULL()
763 : TYPE(tddfpt2_control_type), POINTER :: tddfpt2_control => NULL()
764 : TYPE(rixs_control_type), POINTER :: rixs_control => NULL()
765 : TYPE(xas_control_type), POINTER :: xas_control => NULL()
766 : TYPE(expot_control_type), POINTER :: expot_control => NULL()
767 : TYPE(maxwell_control_type), POINTER :: maxwell_control => NULL()
768 : TYPE(smeagol_control_type), POINTER :: smeagol_control => NULL()
769 : TYPE(paep_control_type), POINTER :: paep_control => NULL()
770 : TYPE(pcc_control_type), POINTER :: pcc_control => NULL()
771 : TYPE(efield_p_type), POINTER, &
772 : DIMENSION(:) :: efield_fields => NULL()
773 : TYPE(hairy_probes_type), POINTER, &
774 : DIMENSION(:) :: probe => NULL()
775 : INTEGER :: nspins = 0, &
776 : charge = 0, &
777 : multiplicity = 0, &
778 : sic_method_id = 0, &
779 : plus_u_method_id = 0, &
780 : dir_surf_dip = 0, &
781 : nimages = 1
782 : INTEGER :: sic_list_id = 0
783 : INTEGER :: auto_basis_ri_aux = 1, &
784 : auto_basis_aux_fit = 1, &
785 : auto_basis_lri_aux = 1, &
786 : auto_basis_p_lri_aux = 1, &
787 : auto_basis_ri_hxc = 1, &
788 : auto_basis_ri_xas = 1, &
789 : auto_basis_ri_hfx = 1, &
790 : localize_each = 1, &
791 : mtlr_ikind = 1, &
792 : max_mtlr_iter = 10
793 : REAL(KIND=dp) :: eps_u_j_loop = 0.01_dp, &
794 : relax_multiplicity = 0.0_dp, &
795 : sic_scaling_a = 0.0_dp, &
796 : sic_scaling_b = 0.0_dp, &
797 : pos_dir_surf_dip = 0.0_dp, &
798 : perturbation_strength = 0.0_dp, &
799 : trq(2) = 0.0_dp, &
800 : vhxc(2) = 0.0_dp
801 : LOGICAL :: do_xas_calculation = .FALSE., &
802 : do_xas_tdp_calculation = .FALSE., &
803 : drho_by_collocation = .FALSE., &
804 : use_kinetic_energy_density = .FALSE., &
805 : restricted = .FALSE., &
806 : roks = .FALSE., &
807 : uks = .FALSE., &
808 : lsd = .FALSE., &
809 : dft_plus_u = .FALSE., &
810 : mtlr_u_j = .FALSE., &
811 : apply_efield = .FALSE., &
812 : apply_efield_field = .FALSE., &
813 : apply_vector_potential = .FALSE., &
814 : apply_period_efield = .FALSE., &
815 : apply_external_potential = .FALSE., &
816 : eval_external_potential = .FALSE., &
817 : do_admm = .FALSE., &
818 : do_admm_dm = .FALSE., &
819 : do_admm_mo = .FALSE., &
820 : smear = .FALSE., &
821 : low_spin_roks = .FALSE., &
822 : apply_external_density = .FALSE., &
823 : read_external_density = .FALSE., &
824 : apply_external_vxc = .FALSE., &
825 : read_external_vxc = .FALSE., &
826 : correct_surf_dip = .FALSE., &
827 : surf_dip_correct_switch = .FALSE., &
828 : switch_surf_dip = .FALSE., &
829 : correct_el_density_dip = .FALSE., &
830 : do_sccs = .FALSE., &
831 : apply_embed_pot = .FALSE., &
832 : apply_dmfet_pot = .FALSE., &
833 : hairy_probes = .FALSE., &
834 : do_paep = .FALSE., &
835 : do_pcc = .FALSE., &
836 : ! If true, then no other XC functionals are allowed
837 : use_gauxc = .FALSE., &
838 : mtlr_dft_with_perturbation = .FALSE.
839 : END TYPE dft_control_type
840 :
841 : CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'cp_control_types'
842 :
843 : ! Public data types
844 :
845 : PUBLIC :: dft_control_type, &
846 : qs_control_type, &
847 : gapw_control_type, &
848 : tddfpt2_control_type, &
849 : rixs_control_type, &
850 : proj_mo_type, &
851 : efield_type, &
852 : mulliken_restraint_type, &
853 : ddapc_restraint_type, &
854 : dftb_control_type, &
855 : xtb_control_type, &
856 : xtb_reference_cli_type, &
857 : semi_empirical_control_type, &
858 : s2_restraint_type, &
859 : admm_control_type, &
860 : maxwell_control_type, &
861 : expot_control_type, &
862 : rtp_control_type, &
863 : sccs_control_type, &
864 : stda_control_type, &
865 : smeared_type, &
866 : hairy_probes_type, &
867 : paep_control_type, &
868 : pcc_control_type
869 :
870 : ! Public subroutines
871 :
872 : PUBLIC :: dft_control_release, &
873 : dft_control_create, &
874 : admm_control_create, &
875 : admm_control_release, &
876 : maxwell_control_create, &
877 : expot_control_create, &
878 : ddapc_control_create, &
879 : rixs_control_create, &
880 : rixs_control_release
881 :
882 : CONTAINS
883 :
884 : ! **************************************************************************************************
885 : !> \brief create the mulliken_restraint_type
886 : !> \param mulliken_restraint_control ...
887 : !> \par History
888 : !> 02.2005 created [Joost VandeVondele]
889 : ! **************************************************************************************************
890 8800 : SUBROUTINE mulliken_control_create(mulliken_restraint_control)
891 : TYPE(mulliken_restraint_type), INTENT(OUT) :: mulliken_restraint_control
892 :
893 8800 : mulliken_restraint_control%strength = 0.1_dp
894 8800 : mulliken_restraint_control%target = 1.0_dp
895 : mulliken_restraint_control%natoms = 0
896 : NULLIFY (mulliken_restraint_control%atoms)
897 8800 : END SUBROUTINE mulliken_control_create
898 :
899 : ! **************************************************************************************************
900 : !> \brief release the mulliken_restraint_type
901 : !> \param mulliken_restraint_control ...
902 : !> \par History
903 : !> 02.2005 created [Joost VandeVondele]
904 : ! **************************************************************************************************
905 8800 : SUBROUTINE mulliken_control_release(mulliken_restraint_control)
906 : TYPE(mulliken_restraint_type), INTENT(INOUT) :: mulliken_restraint_control
907 :
908 8800 : IF (ASSOCIATED(mulliken_restraint_control%atoms)) THEN
909 2 : DEALLOCATE (mulliken_restraint_control%atoms)
910 : END IF
911 8800 : mulliken_restraint_control%strength = 0.0_dp
912 8800 : mulliken_restraint_control%target = 0.0_dp
913 8800 : mulliken_restraint_control%natoms = 0
914 8800 : END SUBROUTINE mulliken_control_release
915 :
916 : ! **************************************************************************************************
917 : !> \brief create the ddapc_restraint_type
918 : !> \param ddapc_restraint_control ...
919 : !> \par History
920 : !> 02.2006 created [Joost VandeVondele]
921 : ! **************************************************************************************************
922 18 : SUBROUTINE ddapc_control_create(ddapc_restraint_control)
923 : TYPE(ddapc_restraint_type), INTENT(OUT) :: ddapc_restraint_control
924 :
925 : ddapc_restraint_control%density_type = do_full_density
926 18 : ddapc_restraint_control%strength = 0.1_dp
927 : ddapc_restraint_control%ddapc_order_p = 0.0_dp
928 18 : ddapc_restraint_control%functional_form = -1
929 18 : ddapc_restraint_control%target = 1.0_dp
930 : ddapc_restraint_control%natoms = 0
931 : NULLIFY (ddapc_restraint_control%atoms)
932 : NULLIFY (ddapc_restraint_control%coeff)
933 :
934 18 : END SUBROUTINE ddapc_control_create
935 :
936 : ! **************************************************************************************************
937 : !> \brief release the ddapc_restraint_type
938 : !> \param ddapc_restraint_control ...
939 : !> \par History
940 : !> 02.2006 created [Joost VandeVondele]
941 : ! **************************************************************************************************
942 18 : SUBROUTINE ddapc_control_release(ddapc_restraint_control)
943 : TYPE(ddapc_restraint_type), INTENT(INOUT) :: ddapc_restraint_control
944 :
945 18 : IF (ASSOCIATED(ddapc_restraint_control%atoms)) THEN
946 18 : DEALLOCATE (ddapc_restraint_control%atoms)
947 : END IF
948 18 : IF (ASSOCIATED(ddapc_restraint_control%coeff)) THEN
949 18 : DEALLOCATE (ddapc_restraint_control%coeff)
950 : END IF
951 18 : ddapc_restraint_control%strength = 0.0_dp
952 18 : ddapc_restraint_control%target = 0.0_dp
953 18 : ddapc_restraint_control%natoms = 0
954 18 : END SUBROUTINE ddapc_control_release
955 :
956 : ! **************************************************************************************************
957 : !> \brief create the s2_restraint_type
958 : !> \param s2_restraint_control ...
959 : !> \par History
960 : !> 03.2006 created [Joost VandeVondele]
961 : ! **************************************************************************************************
962 8800 : SUBROUTINE s2_control_create(s2_restraint_control)
963 : TYPE(s2_restraint_type), INTENT(OUT) :: s2_restraint_control
964 :
965 8800 : s2_restraint_control%strength = 0.1_dp
966 : s2_restraint_control%s2_order_p = 0.0_dp
967 8800 : s2_restraint_control%functional_form = -1
968 8800 : s2_restraint_control%target = 1.0_dp
969 8800 : END SUBROUTINE s2_control_create
970 :
971 : ! **************************************************************************************************
972 : !> \brief release the s2_restraint_type
973 : !> \param s2_restraint_control ...
974 : !> \par History
975 : !> 03.2006 created [Joost VandeVondele]
976 : ! **************************************************************************************************
977 8800 : SUBROUTINE s2_control_release(s2_restraint_control)
978 : TYPE(s2_restraint_type), INTENT(INOUT) :: s2_restraint_control
979 :
980 8800 : s2_restraint_control%strength = 0.0_dp
981 8800 : s2_restraint_control%target = 0.0_dp
982 8800 : END SUBROUTINE s2_control_release
983 :
984 : ! **************************************************************************************************
985 : !> \brief allocates and perform a very basic initialization
986 : !> \param dft_control the object to create
987 : !> \par History
988 : !> 02.2003 created [fawzi]
989 : !> \author fawzi
990 : ! **************************************************************************************************
991 52800 : SUBROUTINE dft_control_create(dft_control)
992 : TYPE(dft_control_type), INTENT(OUT) :: dft_control
993 :
994 : NULLIFY (dft_control%xas_control)
995 : NULLIFY (dft_control%qs_control)
996 : NULLIFY (dft_control%tddfpt2_control)
997 : NULLIFY (dft_control%rixs_control)
998 : NULLIFY (dft_control%efield_fields)
999 : NULLIFY (dft_control%period_efield)
1000 : NULLIFY (dft_control%admm_control)
1001 : NULLIFY (dft_control%expot_control)
1002 : NULLIFY (dft_control%maxwell_control)
1003 : NULLIFY (dft_control%smeagol_control)
1004 : NULLIFY (dft_control%rtp_control)
1005 : NULLIFY (dft_control%sccs_control)
1006 : NULLIFY (dft_control%probe)
1007 : NULLIFY (dft_control%paep_control)
1008 : NULLIFY (dft_control%pcc_control)
1009 : dft_control%do_sccs = .FALSE.
1010 : dft_control%apply_embed_pot = .FALSE.
1011 : dft_control%apply_dmfet_pot = .FALSE.
1012 : dft_control%hairy_probes = .FALSE.
1013 : dft_control%do_pcc = .FALSE.
1014 : dft_control%do_paep = .FALSE.
1015 8800 : CALL qs_control_create(dft_control%qs_control)
1016 8800 : CALL tddfpt2_control_create(dft_control%tddfpt2_control)
1017 8800 : CALL rixs_control_create(dft_control%rixs_control)
1018 8800 : CALL smeagol_control_create(dft_control%smeagol_control)
1019 8800 : END SUBROUTINE dft_control_create
1020 :
1021 : ! **************************************************************************************************
1022 : !> \brief ...
1023 : !> \param dft_control ...
1024 : !> \par History
1025 : !> 02.2003 created [fawzi]
1026 : !> \author fawzi
1027 : ! **************************************************************************************************
1028 8800 : SUBROUTINE dft_control_release(dft_control)
1029 : TYPE(dft_control_type), INTENT(INOUT) :: dft_control
1030 :
1031 : INTEGER :: i
1032 :
1033 8800 : CALL qs_control_release(dft_control%qs_control)
1034 8800 : CALL tddfpt2_control_release(dft_control%tddfpt2_control)
1035 8800 : CALL rixs_control_release(dft_control%rixs_control) ! maybe check first if allocated
1036 8800 : IF (ASSOCIATED(dft_control%xas_control)) THEN
1037 42 : CALL xas_control_release(dft_control%xas_control)
1038 42 : DEALLOCATE (dft_control%xas_control)
1039 : END IF
1040 8800 : CALL admm_control_release(dft_control%admm_control)
1041 8800 : CALL expot_control_release(dft_control%expot_control)
1042 8800 : CALL maxwell_control_release(dft_control%maxwell_control)
1043 8800 : CALL smeagol_control_release(dft_control%smeagol_control)
1044 8800 : CALL efield_fields_release(dft_control%efield_fields)
1045 8800 : IF (ASSOCIATED(dft_control%probe)) THEN
1046 12 : DO i = 1, SIZE(dft_control%probe)
1047 12 : DEALLOCATE (dft_control%probe(i)%atom_ids)
1048 : END DO
1049 4 : DEALLOCATE (dft_control%probe)
1050 : END IF
1051 8800 : IF (ASSOCIATED(dft_control%sccs_control)) DEALLOCATE (dft_control%sccs_control)
1052 8800 : IF (ASSOCIATED(dft_control%period_efield)) THEN
1053 76 : DEALLOCATE (dft_control%period_efield)
1054 : END IF
1055 8800 : IF (ASSOCIATED(dft_control%rtp_control)) THEN
1056 308 : CALL proj_mo_list_release(dft_control%rtp_control%proj_mo_list)
1057 308 : DEALLOCATE (dft_control%rtp_control)
1058 : END IF
1059 8800 : IF (ASSOCIATED(dft_control%pcc_control)) DEALLOCATE (dft_control%pcc_control)
1060 8800 : IF (ASSOCIATED(dft_control%paep_control)) DEALLOCATE (dft_control%paep_control)
1061 :
1062 8800 : END SUBROUTINE dft_control_release
1063 :
1064 : ! **************************************************************************************************
1065 : !> \brief ...
1066 : !> \param qs_control ...
1067 : ! **************************************************************************************************
1068 8800 : SUBROUTINE qs_control_create(qs_control)
1069 : TYPE(qs_control_type), POINTER :: qs_control
1070 :
1071 8800 : CPASSERT(.NOT. ASSOCIATED(qs_control))
1072 35200 : ALLOCATE (qs_control)
1073 :
1074 : NULLIFY (qs_control%e_cutoff)
1075 : NULLIFY (qs_control%gapw_control)
1076 : NULLIFY (qs_control%mulliken_restraint_control)
1077 : NULLIFY (qs_control%ddapc_restraint_control)
1078 : NULLIFY (qs_control%s2_restraint_control)
1079 : NULLIFY (qs_control%se_control)
1080 : NULLIFY (qs_control%dftb_control)
1081 : NULLIFY (qs_control%xtb_control)
1082 : NULLIFY (qs_control%cdft_control)
1083 : NULLIFY (qs_control%ddapc_restraint_control)
1084 :
1085 8800 : ALLOCATE (qs_control%mulliken_restraint_control)
1086 8800 : CALL mulliken_control_create(qs_control%mulliken_restraint_control)
1087 8800 : ALLOCATE (qs_control%s2_restraint_control)
1088 8800 : CALL s2_control_create(qs_control%s2_restraint_control)
1089 8800 : ALLOCATE (qs_control%gapw_control)
1090 8800 : CALL se_control_create(qs_control%se_control)
1091 8800 : CALL dftb_control_create(qs_control%dftb_control)
1092 8800 : CALL xtb_control_create(qs_control%xtb_control)
1093 26400 : ALLOCATE (qs_control%cdft_control)
1094 8800 : CALL cdft_control_create(qs_control%cdft_control)
1095 8800 : END SUBROUTINE qs_control_create
1096 :
1097 : ! **************************************************************************************************
1098 : !> \brief ...
1099 : !> \param qs_control ...
1100 : ! **************************************************************************************************
1101 8800 : SUBROUTINE qs_control_release(qs_control)
1102 : TYPE(qs_control_type), POINTER :: qs_control
1103 :
1104 : INTEGER :: i
1105 :
1106 8800 : IF (ASSOCIATED(qs_control)) THEN
1107 8800 : CALL mulliken_control_release(qs_control%mulliken_restraint_control)
1108 8800 : DEALLOCATE (qs_control%mulliken_restraint_control)
1109 8800 : CALL s2_control_release(qs_control%s2_restraint_control)
1110 8800 : DEALLOCATE (qs_control%s2_restraint_control)
1111 8800 : CALL se_control_release(qs_control%se_control)
1112 8800 : CALL dftb_control_release(qs_control%dftb_control)
1113 8800 : CALL xtb_control_release(qs_control%xtb_control)
1114 8800 : IF (ASSOCIATED(qs_control%cdft_control)) THEN
1115 8800 : CALL cdft_control_release(qs_control%cdft_control)
1116 8800 : DEALLOCATE (qs_control%cdft_control)
1117 : END IF
1118 :
1119 8800 : IF (ASSOCIATED(qs_control%e_cutoff)) THEN
1120 8800 : DEALLOCATE (qs_control%e_cutoff)
1121 : END IF
1122 8800 : IF (ASSOCIATED(qs_control%gapw_control)) THEN
1123 8800 : IF (ASSOCIATED(qs_control%gapw_control%aw)) THEN
1124 294 : DEALLOCATE (qs_control%gapw_control%aw)
1125 : END IF
1126 8800 : DEALLOCATE (qs_control%gapw_control)
1127 : END IF
1128 8800 : IF (ASSOCIATED(qs_control%ddapc_restraint_control)) THEN
1129 32 : DO i = 1, SIZE(qs_control%ddapc_restraint_control)
1130 32 : CALL ddapc_control_release(qs_control%ddapc_restraint_control(i))
1131 : END DO
1132 14 : DEALLOCATE (qs_control%ddapc_restraint_control)
1133 : END IF
1134 8800 : DEALLOCATE (qs_control)
1135 : END IF
1136 8800 : END SUBROUTINE qs_control_release
1137 :
1138 : ! **************************************************************************************************
1139 : !> \brief allocate control options for Time-Dependent Density Functional Theory calculation
1140 : !> \param tddfpt_control an object to create
1141 : !> \par History
1142 : !> * 05.2016 created [Sergey Chulkov]
1143 : ! **************************************************************************************************
1144 35232 : SUBROUTINE tddfpt2_control_create(tddfpt_control)
1145 : TYPE(tddfpt2_control_type), POINTER :: tddfpt_control
1146 :
1147 : CHARACTER(len=*), PARAMETER :: routineN = 'tddfpt2_control_create'
1148 :
1149 : INTEGER :: handle
1150 :
1151 17616 : CPASSERT(.NOT. ASSOCIATED(tddfpt_control))
1152 17616 : CALL timeset(routineN, handle)
1153 :
1154 70464 : ALLOCATE (tddfpt_control)
1155 : tddfpt_control%do_soc = .FALSE.
1156 :
1157 17616 : CALL timestop(handle)
1158 17616 : END SUBROUTINE tddfpt2_control_create
1159 :
1160 : ! **************************************************************************************************
1161 : !> \brief release memory allocated for TDDFT control options
1162 : !> \param tddfpt_control an object to release
1163 : !> \par History
1164 : !> * 05.2016 created [Sergey Chulkov]
1165 : ! **************************************************************************************************
1166 17616 : SUBROUTINE tddfpt2_control_release(tddfpt_control)
1167 : TYPE(tddfpt2_control_type), POINTER :: tddfpt_control
1168 :
1169 : CHARACTER(len=*), PARAMETER :: routineN = 'tddfpt2_control_release'
1170 :
1171 : INTEGER :: handle
1172 :
1173 17616 : CALL timeset(routineN, handle)
1174 :
1175 17616 : IF (ASSOCIATED(tddfpt_control)) THEN
1176 17616 : DEALLOCATE (tddfpt_control)
1177 : END IF
1178 :
1179 17616 : CALL timestop(handle)
1180 17616 : END SUBROUTINE tddfpt2_control_release
1181 :
1182 : ! **************************************************************************************************
1183 : !> \brief Creates and initializes the rixs_control_type
1184 : !> \param rixs_control the type to initialize
1185 : ! **************************************************************************************************
1186 8816 : SUBROUTINE rixs_control_create(rixs_control)
1187 : TYPE(rixs_control_type), POINTER :: rixs_control
1188 :
1189 8816 : CPASSERT(.NOT. ASSOCIATED(rixs_control))
1190 8816 : ALLOCATE (rixs_control)
1191 :
1192 : NULLIFY (rixs_control%tddfpt2_control)
1193 8816 : CALL tddfpt2_control_create(rixs_control%tddfpt2_control)
1194 8816 : NULLIFY (rixs_control%xas_tdp_control)
1195 8816 : CALL xas_tdp_control_create(rixs_control%xas_tdp_control)
1196 :
1197 8816 : END SUBROUTINE rixs_control_create
1198 :
1199 : ! **************************************************************************************************
1200 : !> \brief Releases the rixs_control_type
1201 : !> \param rixs_control ...
1202 : ! **************************************************************************************************
1203 8816 : SUBROUTINE rixs_control_release(rixs_control)
1204 : TYPE(rixs_control_type), POINTER :: rixs_control
1205 :
1206 8816 : IF (ASSOCIATED(rixs_control)) THEN
1207 8816 : CALL tddfpt2_control_release(rixs_control%tddfpt2_control)
1208 8816 : CALL xas_tdp_control_release(rixs_control%xas_tdp_control)
1209 :
1210 8816 : DEALLOCATE (rixs_control)
1211 : END IF
1212 :
1213 8816 : END SUBROUTINE rixs_control_release
1214 :
1215 : ! **************************************************************************************************
1216 : !> \brief ...
1217 : !> \param proj_mo_list ...
1218 : ! **************************************************************************************************
1219 308 : SUBROUTINE proj_mo_list_release(proj_mo_list)
1220 : TYPE(proj_mo_p_type), DIMENSION(:), POINTER :: proj_mo_list
1221 :
1222 : INTEGER :: i, mo_ref_nbr
1223 :
1224 308 : IF (ASSOCIATED(proj_mo_list)) THEN
1225 50 : DO i = 1, SIZE(proj_mo_list)
1226 50 : IF (ASSOCIATED(proj_mo_list(i)%proj_mo)) THEN
1227 46 : IF (ALLOCATED(proj_mo_list(i)%proj_mo%ref_mo_index)) THEN
1228 46 : DEALLOCATE (proj_mo_list(i)%proj_mo%ref_mo_index)
1229 : END IF
1230 46 : IF (ALLOCATED(proj_mo_list(i)%proj_mo%mo_ref)) THEN
1231 192 : DO mo_ref_nbr = 1, SIZE(proj_mo_list(i)%proj_mo%mo_ref)
1232 192 : CALL cp_fm_release(proj_mo_list(i)%proj_mo%mo_ref(mo_ref_nbr))
1233 : END DO
1234 46 : DEALLOCATE (proj_mo_list(i)%proj_mo%mo_ref)
1235 : END IF
1236 46 : IF (ALLOCATED(proj_mo_list(i)%proj_mo%td_mo_index)) THEN
1237 46 : DEALLOCATE (proj_mo_list(i)%proj_mo%td_mo_index)
1238 : END IF
1239 46 : IF (ALLOCATED(proj_mo_list(i)%proj_mo%td_mo_occ)) THEN
1240 46 : DEALLOCATE (proj_mo_list(i)%proj_mo%td_mo_occ)
1241 : END IF
1242 46 : DEALLOCATE (proj_mo_list(i)%proj_mo)
1243 : END IF
1244 : END DO
1245 4 : DEALLOCATE (proj_mo_list)
1246 : END IF
1247 308 : END SUBROUTINE proj_mo_list_release
1248 :
1249 : ! **************************************************************************************************
1250 : !> \brief ...
1251 : !> \param efield_fields ...
1252 : ! **************************************************************************************************
1253 8800 : SUBROUTINE efield_fields_release(efield_fields)
1254 : TYPE(efield_p_type), DIMENSION(:), POINTER :: efield_fields
1255 :
1256 : INTEGER :: i
1257 :
1258 8800 : IF (ASSOCIATED(efield_fields)) THEN
1259 636 : DO i = 1, SIZE(efield_fields)
1260 636 : IF (ASSOCIATED(efield_fields(i)%efield)) THEN
1261 318 : IF (ASSOCIATED(efield_fields(i)%efield%envelop_r_vars)) THEN
1262 10 : DEALLOCATE (efield_fields(i)%efield%envelop_r_vars)
1263 : END IF
1264 318 : IF (ASSOCIATED(efield_fields(i)%efield%envelop_i_vars)) THEN
1265 308 : DEALLOCATE (efield_fields(i)%efield%envelop_i_vars)
1266 : END IF
1267 318 : IF (ASSOCIATED(efield_fields(i)%efield%polarisation)) THEN
1268 318 : DEALLOCATE (efield_fields(i)%efield%polarisation)
1269 : END IF
1270 318 : DEALLOCATE (efield_fields(i)%efield)
1271 : END IF
1272 : END DO
1273 318 : DEALLOCATE (efield_fields)
1274 : END IF
1275 8800 : END SUBROUTINE efield_fields_release
1276 :
1277 : ! **************************************************************************************************
1278 : !> \brief ...
1279 : !> \param dftb_control ...
1280 : ! **************************************************************************************************
1281 8800 : SUBROUTINE dftb_control_create(dftb_control)
1282 : TYPE(dftb_control_type), POINTER :: dftb_control
1283 :
1284 8800 : CPASSERT(.NOT. ASSOCIATED(dftb_control))
1285 88000 : ALLOCATE (dftb_control)
1286 :
1287 : NULLIFY (dftb_control%sk_pair_list)
1288 8800 : END SUBROUTINE dftb_control_create
1289 :
1290 : ! **************************************************************************************************
1291 : !> \brief ...
1292 : !> \param dftb_control ...
1293 : ! **************************************************************************************************
1294 8800 : SUBROUTINE dftb_control_release(dftb_control)
1295 : TYPE(dftb_control_type), POINTER :: dftb_control
1296 :
1297 8800 : IF (ASSOCIATED(dftb_control)) THEN
1298 8800 : IF (ASSOCIATED(dftb_control%sk_pair_list)) THEN
1299 294 : DEALLOCATE (dftb_control%sk_pair_list)
1300 : END IF
1301 8800 : DEALLOCATE (dftb_control)
1302 : END IF
1303 8800 : END SUBROUTINE dftb_control_release
1304 :
1305 : ! **************************************************************************************************
1306 : !> \brief ...
1307 : !> \param xtb_control ...
1308 : ! **************************************************************************************************
1309 8800 : SUBROUTINE xtb_control_create(xtb_control)
1310 : TYPE(xtb_control_type), POINTER :: xtb_control
1311 :
1312 8800 : CPASSERT(.NOT. ASSOCIATED(xtb_control))
1313 61600 : ALLOCATE (xtb_control)
1314 :
1315 : NULLIFY (xtb_control%kab_param)
1316 : NULLIFY (xtb_control%kab_vals)
1317 : NULLIFY (xtb_control%kab_types)
1318 : NULLIFY (xtb_control%nonbonded)
1319 : NULLIFY (xtb_control%rcpair)
1320 : NULLIFY (xtb_control%spinpol_type)
1321 : NULLIFY (xtb_control%spinpol_vals)
1322 :
1323 8800 : END SUBROUTINE xtb_control_create
1324 :
1325 : ! **************************************************************************************************
1326 : !> \brief ...
1327 : !> \param xtb_control ...
1328 : ! **************************************************************************************************
1329 8800 : SUBROUTINE xtb_control_release(xtb_control)
1330 : TYPE(xtb_control_type), POINTER :: xtb_control
1331 :
1332 8800 : IF (ASSOCIATED(xtb_control)) THEN
1333 8800 : IF (ASSOCIATED(xtb_control%kab_param)) THEN
1334 2 : DEALLOCATE (xtb_control%kab_param)
1335 : END IF
1336 8800 : IF (ASSOCIATED(xtb_control%kab_vals)) THEN
1337 2 : DEALLOCATE (xtb_control%kab_vals)
1338 : END IF
1339 8800 : IF (ASSOCIATED(xtb_control%kab_types)) THEN
1340 2 : DEALLOCATE (xtb_control%kab_types)
1341 : END IF
1342 8800 : IF (ASSOCIATED(xtb_control%rcpair)) THEN
1343 1030 : DEALLOCATE (xtb_control%rcpair)
1344 : END IF
1345 8800 : IF (ASSOCIATED(xtb_control%nonbonded)) THEN
1346 6 : CALL pair_potential_p_release(xtb_control%nonbonded)
1347 : END IF
1348 8800 : IF (ASSOCIATED(xtb_control%spinpol_type)) THEN
1349 2 : DEALLOCATE (xtb_control%spinpol_type)
1350 : END IF
1351 8800 : IF (ASSOCIATED(xtb_control%spinpol_vals)) THEN
1352 2 : DEALLOCATE (xtb_control%spinpol_vals)
1353 : END IF
1354 8800 : DEALLOCATE (xtb_control)
1355 : END IF
1356 8800 : END SUBROUTINE xtb_control_release
1357 :
1358 : ! **************************************************************************************************
1359 : !> \brief ...
1360 : !> \param se_control ...
1361 : ! **************************************************************************************************
1362 8800 : SUBROUTINE se_control_create(se_control)
1363 : TYPE(semi_empirical_control_type), POINTER :: se_control
1364 :
1365 8800 : CPASSERT(.NOT. ASSOCIATED(se_control))
1366 44000 : ALLOCATE (se_control)
1367 8800 : END SUBROUTINE se_control_create
1368 :
1369 : ! **************************************************************************************************
1370 : !> \brief ...
1371 : !> \param se_control ...
1372 : ! **************************************************************************************************
1373 8800 : SUBROUTINE se_control_release(se_control)
1374 : TYPE(semi_empirical_control_type), POINTER :: se_control
1375 :
1376 8800 : IF (ASSOCIATED(se_control)) THEN
1377 8800 : DEALLOCATE (se_control)
1378 : END IF
1379 8800 : END SUBROUTINE se_control_release
1380 :
1381 : ! **************************************************************************************************
1382 : !> \brief ...
1383 : !> \param admm_control ...
1384 : ! **************************************************************************************************
1385 520 : SUBROUTINE admm_control_create(admm_control)
1386 : TYPE(admm_control_type), POINTER :: admm_control
1387 :
1388 520 : CPASSERT(.NOT. ASSOCIATED(admm_control))
1389 2600 : ALLOCATE (admm_control)
1390 :
1391 520 : END SUBROUTINE admm_control_create
1392 :
1393 : ! **************************************************************************************************
1394 : !> \brief ...
1395 : !> \param admm_control ...
1396 : ! **************************************************************************************************
1397 8804 : SUBROUTINE admm_control_release(admm_control)
1398 : TYPE(admm_control_type), POINTER :: admm_control
1399 :
1400 8804 : IF (ASSOCIATED(admm_control)) THEN
1401 560 : DEALLOCATE (admm_control)
1402 : END IF
1403 8804 : END SUBROUTINE admm_control_release
1404 :
1405 : ! **************************************************************************************************
1406 : !> \brief ...
1407 : !> \param expot_control ...
1408 : ! **************************************************************************************************
1409 16 : SUBROUTINE expot_control_create(expot_control)
1410 : TYPE(expot_control_type), POINTER :: expot_control
1411 :
1412 16 : CPASSERT(.NOT. ASSOCIATED(expot_control))
1413 16 : ALLOCATE (expot_control)
1414 : expot_control%read_from_cube = .FALSE.
1415 : expot_control%maxwell_solver = .FALSE.
1416 16 : expot_control%static = .TRUE.
1417 16 : expot_control%scaling_factor = 1.0_dp
1418 :
1419 16 : END SUBROUTINE expot_control_create
1420 :
1421 : ! **************************************************************************************************
1422 : !> \brief ...
1423 : !> \param expot_control ...
1424 : ! **************************************************************************************************
1425 8800 : SUBROUTINE expot_control_release(expot_control)
1426 : TYPE(expot_control_type), POINTER :: expot_control
1427 :
1428 8800 : IF (ASSOCIATED(expot_control)) THEN
1429 16 : DEALLOCATE (expot_control)
1430 : END IF
1431 :
1432 8800 : END SUBROUTINE expot_control_release
1433 :
1434 : ! **************************************************************************************************
1435 : !> \brief ...
1436 : !> \param maxwell_control ...
1437 : ! **************************************************************************************************
1438 0 : SUBROUTINE maxwell_control_create(maxwell_control)
1439 : TYPE(maxwell_control_type), POINTER :: maxwell_control
1440 :
1441 0 : CPASSERT(.NOT. ASSOCIATED(maxwell_control))
1442 0 : ALLOCATE (maxwell_control)
1443 :
1444 0 : END SUBROUTINE maxwell_control_create
1445 :
1446 : ! **************************************************************************************************
1447 : !> \brief ...
1448 : !> \param maxwell_control ...
1449 : ! **************************************************************************************************
1450 8800 : SUBROUTINE maxwell_control_release(maxwell_control)
1451 : TYPE(maxwell_control_type), POINTER :: maxwell_control
1452 :
1453 8800 : IF (ASSOCIATED(maxwell_control)) THEN
1454 0 : DEALLOCATE (maxwell_control)
1455 : END IF
1456 :
1457 8800 : END SUBROUTINE maxwell_control_release
1458 :
1459 0 : END MODULE cp_control_types
|
