(git:cb5d5fc)
Loading...
Searching...
No Matches
input_restart_force_eval.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
8MODULE input_restart_force_eval
9
10 USE atomic_kind_list_types, ONLY: atomic_kind_list_type
11 USE atomic_kind_types, ONLY: get_atomic_kind
12 USE cell_types, ONLY: cell_type,&
13 real_to_scaled
14 USE cp_control_types, ONLY: dft_control_type
15 USE cp_files, ONLY: close_file,&
16 open_file
17 USE cp_linked_list_input, ONLY: cp_sll_val_create,&
18 cp_sll_val_get_length,&
19 cp_sll_val_type
20 USE cp_subsys_types, ONLY: cp_subsys_get,&
21 cp_subsys_type
22 USE cp_units, ONLY: cp_unit_from_cp2k
23 USE distribution_1d_types, ONLY: distribution_1d_type
24 USE eip_environment_types, ONLY: eip_env_get
25 USE force_env_types, ONLY: force_env_get,&
26 force_env_type,&
27 multiple_fe_list,&
28 use_eip_force,&
29 use_nnp_force,&
30 use_qmmm,&
31 use_qs_force
32 USE input_constants, ONLY: do_coord_cp2k,&
33 ehrenfest,&
34 mol_dyn_run,&
35 mon_car_run,&
36 pint_run,&
37 use_rt_restart
38 USE input_cp2k_restarts_util, ONLY: section_velocity_val_set
39 USE input_restart_rng, ONLY: section_rng_val_set
40 USE input_section_types, ONLY: &
41 section_get_keyword_index, section_type, section_vals_add_values, section_vals_get, &
42 section_vals_get_subs_vals, section_vals_get_subs_vals3, section_vals_remove_values, &
43 section_vals_type, section_vals_val_get, section_vals_val_set, section_vals_val_unset
44 USE input_val_types, ONLY: val_create,&
45 val_release,&
46 val_type
47 USE kinds, ONLY: default_string_length,&
48 dp
49 USE message_passing, ONLY: mp_para_env_type
50 USE molecule_kind_types, ONLY: get_molecule_kind
51 USE molecule_list_types, ONLY: molecule_list_type
52 USE molecule_types, ONLY: get_molecule,&
53 molecule_type
54 USE multipole_types, ONLY: multipole_type
55 USE nnp_environment_types, ONLY: nnp_env_get
56 USE parallel_rng_types, ONLY: rng_record_length
57 USE particle_list_types, ONLY: particle_list_type
58 USE qmmm_ff_fist, ONLY: qmmm_ff_precond_only_qm
59 USE qs_environment_types, ONLY: get_qs_env
60 USE string_utilities, ONLY: string_to_ascii
61 USE virial_types, ONLY: virial_type
62#include "./base/base_uses.f90"
63
64 IMPLICIT NONE
65
66 PRIVATE
67
68 CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'input_restart_force_eval'
69
70 PUBLIC :: update_force_eval, update_subsys
71
72CONTAINS
73
74! **************************************************************************************************
75!> \brief Updates the force_eval section of the input file
76!> \param force_env ...
77!> \param root_section ...
78!> \param write_binary_restart_file ...
79!> \param respa ...
80!> \par History
81!> 01.2006 created [teo]
82!> \author Teodoro Laino
83! **************************************************************************************************
84 SUBROUTINE update_force_eval(force_env, root_section, &
85 write_binary_restart_file, respa)
86
87 TYPE(force_env_type), POINTER :: force_env
88 TYPE(section_vals_type), POINTER :: root_section
89 LOGICAL, INTENT(IN) :: write_binary_restart_file
90 LOGICAL, OPTIONAL :: respa
91
92 INTEGER :: i, i_subsys, iforce_eval, myid, &
93 nforce_eval
94 INTEGER, DIMENSION(:), POINTER :: i_force_eval
95 LOGICAL :: is_present, multiple_subsys, &
96 skip_vel_section
97 REAL(kind=dp), DIMENSION(:), POINTER :: work
98 TYPE(cell_type), POINTER :: cell
99 TYPE(cp_subsys_type), POINTER :: subsys
100 TYPE(dft_control_type), POINTER :: dft_control
101 TYPE(section_vals_type), POINTER :: cell_section, dft_section, &
102 force_env_sections, qmmm_section, &
103 rng_section, subsys_section, &
104 tmp_section
105 TYPE(virial_type), POINTER :: virial
106
107 NULLIFY (rng_section, subsys_section, cell_section, virial, subsys, cell, dft_control, work, tmp_section)
108 ! If it's not a dynamical run don't update the velocity section
109 CALL section_vals_val_get(root_section, "GLOBAL%RUN_TYPE", i_val=myid)
110 skip_vel_section = ((myid /= mol_dyn_run) .AND. &
111 (myid /= mon_car_run) .AND. &
112 (myid /= pint_run) .AND. &
113 (myid /= ehrenfest))
114
115 ! Go on updatig the force_env_section
116 force_env_sections => section_vals_get_subs_vals(root_section, "FORCE_EVAL")
117 CALL multiple_fe_list(force_env_sections, root_section, i_force_eval, nforce_eval)
118 ! The update of the input MUST be realized only on the main force_eval
119 ! All the others will be left not updated because there is no real need to update them...
120 iforce_eval = 1
121 subsys_section => section_vals_get_subs_vals3(force_env_sections, "SUBSYS", &
122 i_rep_section=i_force_eval(iforce_eval))
123 CALL update_subsys(subsys_section, force_env, skip_vel_section, &
124 write_binary_restart_file)
125
126 ! For RESPA we need to update all subsystems
127 IF (PRESENT(respa)) THEN
128 IF (respa) THEN
129 DO i_subsys = 1, 2
130 iforce_eval = i_subsys
131 subsys_section => section_vals_get_subs_vals3(force_env_sections, "SUBSYS", &
132 i_rep_section=i_force_eval(iforce_eval))
133 CALL update_subsys(subsys_section, force_env, skip_vel_section, &
134 write_binary_restart_file)
135 END DO
136 END IF
137 END IF
138
139 rng_section => section_vals_get_subs_vals(subsys_section, "RNG_INIT")
140 CALL update_rng_particle(rng_section, force_env)
141
142 qmmm_section => section_vals_get_subs_vals3(force_env_sections, "QMMM", &
143 i_rep_section=i_force_eval(iforce_eval))
144 CALL update_qmmm(qmmm_section, force_env)
145
146 ! Non-mixed CDFT calculation: update constraint Lagrangian and counter
147 IF (nforce_eval == 1) THEN
148 IF (ASSOCIATED(force_env%qs_env)) THEN
149 CALL get_qs_env(force_env%qs_env, dft_control=dft_control)
150 IF (dft_control%qs_control%cdft) THEN
151 dft_section => section_vals_get_subs_vals3(force_env_sections, "DFT", &
152 i_rep_section=i_force_eval(iforce_eval))
153 CALL section_vals_val_set(dft_section, "QS%CDFT%COUNTER", &
154 i_val=dft_control%qs_control%cdft_control%ienergy)
155 ALLOCATE (work(SIZE(dft_control%qs_control%cdft_control%strength)))
156 work = dft_control%qs_control%cdft_control%strength
157 CALL section_vals_val_set(dft_section, "QS%CDFT%STRENGTH", &
158 r_vals_ptr=work)
159 END IF
160 END IF
161 END IF
162 ! And now update only the cells of all other force_eval
163 ! This is to make consistent for cell variable runs..
164 IF (nforce_eval > 1) THEN
165 CALL force_env_get(force_env, subsys=subsys, cell=cell)
166 CALL cp_subsys_get(subsys, virial=virial)
167 CALL section_vals_val_get(root_section, "MULTIPLE_FORCE_EVALS%MULTIPLE_SUBSYS", l_val=multiple_subsys)
168 IF (virial%pv_availability .AND. multiple_subsys) THEN
169 DO iforce_eval = 2, nforce_eval
170 subsys_section => section_vals_get_subs_vals3(force_env_sections, "SUBSYS", &
171 i_rep_section=i_force_eval(iforce_eval))
172 cell_section => section_vals_get_subs_vals(subsys_section, "CELL")
173 CALL update_cell_section(cell, cell_section)
174 END DO
175 END IF
176 ! With mixed CDFT, update value of constraint Lagrangian
177 IF (ASSOCIATED(force_env%mixed_env)) THEN
178 IF (force_env%mixed_env%do_mixed_cdft) THEN
179 DO iforce_eval = 2, nforce_eval
180 dft_section => section_vals_get_subs_vals3(force_env_sections, "DFT", &
181 i_rep_section=i_force_eval(iforce_eval))
182 ALLOCATE (work(SIZE(force_env%mixed_env%strength, 2)))
183 work = force_env%mixed_env%strength(iforce_eval - 1, :)
184 CALL section_vals_val_set(dft_section, "QS%CDFT%STRENGTH", &
185 r_vals_ptr=work)
186 CALL section_vals_val_set(dft_section, "QS%CDFT%COUNTER", &
187 i_val=force_env%mixed_env%cdft_control%sim_step)
188 END DO
189 END IF
190 END IF
191 END IF
192
193 IF (ASSOCIATED(force_env%qs_env)) THEN
194 CALL get_qs_env(force_env%qs_env, dft_control=dft_control)
195 dft_section => section_vals_get_subs_vals(force_env_sections, "DFT", &
196 i_rep_section=i_force_eval(iforce_eval))
197 tmp_section => section_vals_get_subs_vals(dft_section, "REAL_TIME_PROPAGATION")
198 CALL section_vals_get(tmp_section, explicit=is_present)
199 IF (is_present) THEN
200 CALL section_vals_val_set(root_section, "FORCE_EVAL%DFT%REAL_TIME_PROPAGATION%INITIAL_WFN", &
201 i_val=use_rt_restart)
202 IF (ASSOCIATED(dft_control%efield_fields)) THEN
203 tmp_section => section_vals_get_subs_vals(dft_section, "EFIELD")
204 ALLOCATE (work(3))
205 DO i = 1, SIZE(dft_control%efield_fields)
206 work = dft_control%efield_fields(1)%efield%vec_pot_initial
207 CALL section_vals_val_set(tmp_section, "VEC_POT_INITIAL", i_rep_section=i, &
208 r_vals_ptr=work)
209 END DO
210 END IF
211 END IF
212 END IF
213 DEALLOCATE (i_force_eval)
214
215 END SUBROUTINE update_force_eval
216
217! **************************************************************************************************
218!> \brief Updates the qmmm section if needed
219!> \param qmmm_section ...
220!> \param force_env ...
221!> \par History
222!> 08.2007 created [teo]
223!> \author Teodoro Laino
224! **************************************************************************************************
225 SUBROUTINE update_qmmm(qmmm_section, force_env)
226
227 TYPE(section_vals_type), POINTER :: qmmm_section
228 TYPE(force_env_type), POINTER :: force_env
229
230 LOGICAL :: explicit
231 REAL(kind=dp), DIMENSION(:), POINTER :: work
232
233 SELECT CASE (force_env%in_use)
234 CASE (use_qmmm)
235 CALL section_vals_get(qmmm_section, explicit=explicit)
236 cpassert(explicit)
237
238 ALLOCATE (work(3))
239 work = force_env%qmmm_env%qm%transl_v
240 CALL section_vals_val_set(qmmm_section, "INITIAL_TRANSLATION_VECTOR", r_vals_ptr=work)
241 END SELECT
242
243 END SUBROUTINE update_qmmm
244
245! **************************************************************************************************
246!> \brief Updates the rng section of the input file
247!> Write current status of the parallel random number generator (RNG)
248!> \param rng_section ...
249!> \param force_env ...
250!> \par History
251!> 01.2006 created [teo]
252!> \author Teodoro Laino
253! **************************************************************************************************
254 SUBROUTINE update_rng_particle(rng_section, force_env)
255
256 TYPE(section_vals_type), POINTER :: rng_section
257 TYPE(force_env_type), POINTER :: force_env
258
259 CHARACTER(LEN=rng_record_length) :: rng_record
260 INTEGER :: iparticle, iparticle_kind, &
261 iparticle_local, nparticle, &
262 nparticle_kind, nparticle_local
263 INTEGER, ALLOCATABLE, DIMENSION(:, :) :: ascii
264 TYPE(atomic_kind_list_type), POINTER :: atomic_kinds
265 TYPE(cp_subsys_type), POINTER :: subsys
266 TYPE(distribution_1d_type), POINTER :: local_particles
267 TYPE(mp_para_env_type), POINTER :: para_env
268 TYPE(particle_list_type), POINTER :: particles
269
270 CALL force_env_get(force_env, subsys=subsys, para_env=para_env)
271 CALL cp_subsys_get(subsys, atomic_kinds=atomic_kinds, local_particles=local_particles, &
272 particles=particles)
273
274 IF (ASSOCIATED(local_particles%local_particle_set)) THEN
275 nparticle_kind = atomic_kinds%n_els
276 nparticle = particles%n_els
277
278 ALLOCATE (ascii(rng_record_length, nparticle))
279 ascii = 0
280
281 DO iparticle = 1, nparticle
282 DO iparticle_kind = 1, nparticle_kind
283 nparticle_local = local_particles%n_el(iparticle_kind)
284 DO iparticle_local = 1, nparticle_local
285 IF (iparticle == local_particles%list(iparticle_kind)%array(iparticle_local)) THEN
286 CALL local_particles%local_particle_set(iparticle_kind)% &
287 rng(iparticle_local)%stream%dump(rng_record=rng_record)
288 CALL string_to_ascii(rng_record, ascii(:, iparticle))
289 END IF
290 END DO
291 END DO
292 END DO
293
294 CALL para_env%sum(ascii)
295
296 CALL section_rng_val_set(rng_section=rng_section, nsize=nparticle, ascii=ascii)
297
298 DEALLOCATE (ascii)
299 END IF
300
301 END SUBROUTINE update_rng_particle
302
303! **************************************************************************************************
304!> \brief Updates the subsys section of the input file
305!> \param subsys_section ...
306!> \param force_env ...
307!> \param skip_vel_section ...
308!> \param write_binary_restart_file ...
309!> \par History
310!> 01.2006 created [teo]
311!> \author Teodoro Laino
312! **************************************************************************************************
313 SUBROUTINE update_subsys(subsys_section, force_env, skip_vel_section, &
314 write_binary_restart_file)
315
316 TYPE(section_vals_type), POINTER :: subsys_section
317 TYPE(force_env_type), POINTER :: force_env
318 LOGICAL, INTENT(IN) :: skip_vel_section, &
319 write_binary_restart_file
320
321 CHARACTER(LEN=*), PARAMETER :: routinen = 'update_subsys'
322
323 CHARACTER(LEN=default_string_length) :: coord_file_name, unit_str
324 INTEGER :: coord_file_format, handle, output_unit
325 INTEGER, DIMENSION(:), POINTER :: multiple_unit_cell
326 LOGICAL :: scale, use_ref_cell
327 REAL(kind=dp) :: conv_factor
328 TYPE(cell_type), POINTER :: cell
329 TYPE(cp_subsys_type), POINTER :: subsys
330 TYPE(molecule_list_type), POINTER :: molecules
331 TYPE(mp_para_env_type), POINTER :: para_env
332 TYPE(multipole_type), POINTER :: multipoles
333 TYPE(particle_list_type), POINTER :: core_particles, particles, &
334 shell_particles
335 TYPE(section_vals_type), POINTER :: work_section
336
337 NULLIFY (work_section, core_particles, particles, shell_particles, &
338 subsys, cell, para_env, multipoles)
339 CALL timeset(routinen, handle)
340 CALL force_env_get(force_env, subsys=subsys, cell=cell, para_env=para_env)
341
342 CALL cp_subsys_get(subsys, particles=particles, molecules=molecules, &
343 shell_particles=shell_particles, core_particles=core_particles, &
344 multipoles=multipoles)
345
346 ! Remove the multiple_unit_cell from the input structure.. at this point we have
347 ! already all the information we need..
348 ALLOCATE (multiple_unit_cell(3))
349 multiple_unit_cell = 1
350 CALL section_vals_val_set(subsys_section, "TOPOLOGY%MULTIPLE_UNIT_CELL", &
351 i_vals_ptr=multiple_unit_cell)
352
353 ! Coordinates and Velocities
354 work_section => section_vals_get_subs_vals(subsys_section, "COORD")
355 CALL section_vals_val_get(work_section, "UNIT", c_val=unit_str)
356 CALL section_vals_val_get(work_section, "SCALED", l_val=scale)
357 conv_factor = cp_unit_from_cp2k(1.0_dp, trim(unit_str))
358 IF (.NOT. write_binary_restart_file) THEN
359 CALL section_vals_val_get(subsys_section, "TOPOLOGY%COORD_FILE_FORMAT", &
360 i_val=coord_file_format)
361 IF (coord_file_format == do_coord_cp2k) THEN
362 CALL section_vals_val_get(subsys_section, "TOPOLOGY%COORD_FILE_NAME", &
363 c_val=coord_file_name)
364 output_unit = 0
365 IF (para_env%is_source()) THEN
366 CALL open_file(file_name=trim(adjustl(coord_file_name)), &
367 file_action="READWRITE", &
368 file_form="FORMATTED", &
369 file_position="REWIND", &
370 file_status="UNKNOWN", &
371 unit_number=output_unit)
372 CALL dump_coordinates_cp2k(particles, molecules, cell, conv_factor, &
373 output_unit=output_unit, &
374 core_or_shell=.false., &
375 scaled_coordinates=scale)
376 CALL close_file(unit_number=output_unit)
377 END IF
378 ELSE
379 ! Remove COORD_FILE_NAME and COORD_FILE_FORMAT keywords from restart
380 CALL section_vals_val_unset(subsys_section, "TOPOLOGY%COORD_FILE_NAME")
381 CALL section_vals_val_unset(subsys_section, "TOPOLOGY%COORD_FILE_FORMAT")
382 CALL section_coord_val_set(work_section, particles, molecules, conv_factor, scale, &
383 cell)
384 END IF
385 END IF
386 CALL section_vals_val_set(subsys_section, "TOPOLOGY%NUMBER_OF_ATOMS", i_val=particles%n_els)
387 work_section => section_vals_get_subs_vals(subsys_section, "VELOCITY")
388 IF (.NOT. skip_vel_section) THEN
389 IF (.NOT. write_binary_restart_file) THEN
390 CALL section_velocity_val_set(work_section, particles, conv_factor=1.0_dp)
391 END IF
392 ELSE
393 CALL section_vals_remove_values(work_section)
394 END IF
395
396 ! Write restart input for core-shell model
397 IF (.NOT. write_binary_restart_file) THEN
398 IF (ASSOCIATED(shell_particles)) THEN
399 work_section => section_vals_get_subs_vals(subsys_section, "SHELL_COORD")
400 CALL section_vals_val_get(work_section, "UNIT", c_val=unit_str)
401 CALL section_vals_val_get(work_section, "SCALED", l_val=scale)
402 conv_factor = cp_unit_from_cp2k(1.0_dp, trim(unit_str))
403 CALL section_coord_val_set(work_section, shell_particles, molecules, &
404 conv_factor, scale, cell, shell=.true.)
405 IF (.NOT. skip_vel_section) THEN
406 work_section => section_vals_get_subs_vals(subsys_section, "SHELL_VELOCITY")
407 CALL section_velocity_val_set(work_section, shell_particles, conv_factor=1.0_dp)
408 END IF
409 END IF
410 IF (ASSOCIATED(core_particles)) THEN
411 work_section => section_vals_get_subs_vals(subsys_section, "CORE_COORD")
412 CALL section_vals_val_get(work_section, "UNIT", c_val=unit_str)
413 CALL section_vals_val_get(work_section, "SCALED", l_val=scale)
414 conv_factor = cp_unit_from_cp2k(1.0_dp, trim(unit_str))
415 CALL section_coord_val_set(work_section, core_particles, molecules, &
416 conv_factor, scale, cell, shell=.true.)
417 IF (.NOT. skip_vel_section) THEN
418 work_section => section_vals_get_subs_vals(subsys_section, "CORE_VELOCITY")
419 CALL section_velocity_val_set(work_section, core_particles, conv_factor=1.0_dp)
420 END IF
421 END IF
422 END IF
423
424 ! Updating cell info
425 CALL force_env_get(force_env, cell=cell)
426 work_section => section_vals_get_subs_vals(subsys_section, "CELL")
427 CALL update_cell_section(cell, cell_section=work_section)
428
429 ! Updating cell_ref info
430 use_ref_cell = .false.
431 SELECT CASE (force_env%in_use)
432 CASE (use_qs_force)
433 CALL get_qs_env(force_env%qs_env, cell_ref=cell, use_ref_cell=use_ref_cell)
434 CASE (use_eip_force)
435 CALL eip_env_get(force_env%eip_env, cell_ref=cell, use_ref_cell=use_ref_cell)
436 CASE (use_nnp_force)
437 CALL nnp_env_get(force_env%nnp_env, cell_ref=cell, use_ref_cell=use_ref_cell)
438 END SELECT
439 IF (use_ref_cell) THEN
440 work_section => section_vals_get_subs_vals(subsys_section, "CELL%CELL_REF")
441 CALL update_cell_section(cell, cell_section=work_section)
442 END IF
443
444 ! Updating multipoles
445 IF (ASSOCIATED(multipoles)) THEN
446 work_section => section_vals_get_subs_vals(subsys_section, "MULTIPOLES")
447 DO
448 IF (SIZE(work_section%values, 2) == 1) EXIT
449 CALL section_vals_add_values(work_section)
450 END DO
451 IF (ASSOCIATED(multipoles%dipoles)) THEN
452 work_section => section_vals_get_subs_vals(subsys_section, "MULTIPOLES%DIPOLES")
453 CALL update_dipoles_section(multipoles%dipoles, work_section)
454 END IF
455 IF (ASSOCIATED(multipoles%quadrupoles)) THEN
456 work_section => section_vals_get_subs_vals(subsys_section, "MULTIPOLES%QUADRUPOLES")
457 CALL update_quadrupoles_section(multipoles%quadrupoles, work_section)
458 END IF
459 END IF
460
461 CALL timestop(handle)
462
463 END SUBROUTINE update_subsys
464
465! **************************************************************************************************
466!> \brief Routine to update a cell section
467!> \param cell ...
468!> \param cell_section ...
469!> \author Ole Schuett
470! **************************************************************************************************
471 SUBROUTINE update_cell_section(cell, cell_section)
472
473 TYPE(cell_type), POINTER :: cell
474 TYPE(section_vals_type), POINTER :: cell_section
475
476 CHARACTER(LEN=*), PARAMETER :: routinen = 'update_cell_section'
477
478 INTEGER :: handle
479 INTEGER, DIMENSION(:), POINTER :: iwork
480 REAL(kind=dp), DIMENSION(:), POINTER :: work
481
482 NULLIFY (work, iwork)
483 CALL timeset(routinen, handle)
484
485 ! CELL VECTORS - A
486 ALLOCATE (work(3))
487 work(1:3) = cell%hmat(1:3, 1)
488 CALL section_vals_val_set(cell_section, "A", r_vals_ptr=work)
489
490 ! CELL VECTORS - B
491 ALLOCATE (work(3))
492 work(1:3) = cell%hmat(1:3, 2)
493 CALL section_vals_val_set(cell_section, "B", r_vals_ptr=work)
494
495 ! CELL VECTORS - C
496 ALLOCATE (work(3))
497 work(1:3) = cell%hmat(1:3, 3)
498 CALL section_vals_val_set(cell_section, "C", r_vals_ptr=work)
499
500 ! MULTIPLE_UNIT_CELL
501 ALLOCATE (iwork(3))
502 iwork = 1
503 CALL section_vals_val_set(cell_section, "MULTIPLE_UNIT_CELL", i_vals_ptr=iwork)
504
505 ! Unset unused or misleading information
506 CALL section_vals_val_unset(cell_section, "ABC")
507 CALL section_vals_val_unset(cell_section, "ALPHA_BETA_GAMMA")
508
509 CALL timestop(handle)
510
511 END SUBROUTINE update_cell_section
512
513! **************************************************************************************************
514!> \brief routine to dump coordinates.. fast implementation
515!> \param coord_section ...
516!> \param particles ...
517!> \param molecules ...
518!> \param conv_factor ...
519!> \param scale ...
520!> \param cell ...
521!> \param shell ...
522!> \par History
523!> 02.2006 created [teo]
524!> \author Teodoro Laino
525! **************************************************************************************************
526 SUBROUTINE section_coord_val_set(coord_section, particles, molecules, conv_factor, &
527 scale, cell, shell)
528
529 TYPE(section_vals_type), POINTER :: coord_section
530 TYPE(particle_list_type), POINTER :: particles
531 TYPE(molecule_list_type), POINTER :: molecules
532 REAL(kind=dp) :: conv_factor
533 LOGICAL, INTENT(IN) :: scale
534 TYPE(cell_type), POINTER :: cell
535 LOGICAL, INTENT(IN), OPTIONAL :: shell
536
537 CHARACTER(LEN=*), PARAMETER :: routinen = 'section_coord_val_set'
538
539 CHARACTER(LEN=2*default_string_length) :: line
540 CHARACTER(LEN=default_string_length) :: my_tag, name
541 INTEGER :: handle, ik, imol, irk, last_atom, nlist
542 LOGICAL :: ldum, molname_generated, my_shell
543 REAL(kind=dp), DIMENSION(3) :: r, s
544 TYPE(cp_sll_val_type), POINTER :: new_pos, vals
545 TYPE(molecule_type), POINTER :: molecule_now
546 TYPE(section_type), POINTER :: section
547 TYPE(val_type), POINTER :: my_val, old_val
548
549 CALL timeset(routinen, handle)
550
551 NULLIFY (my_val, old_val, section, vals)
552 my_shell = .false.
553 IF (PRESENT(shell)) my_shell = shell
554 cpassert(ASSOCIATED(coord_section))
555 cpassert(coord_section%ref_count > 0)
556 section => coord_section%section
557 ik = section_get_keyword_index(section, "_DEFAULT_KEYWORD_")
558 IF (ik == -2) &
559 CALL cp_abort(__location__, &
560 "section "//trim(section%name)//" does not contain keyword "// &
561 "_DEFAULT_KEYWORD_")
562
563 DO
564 IF (SIZE(coord_section%values, 2) == 1) EXIT
565 CALL section_vals_add_values(coord_section)
566 END DO
567 vals => coord_section%values(ik, 1)%list
568 nlist = 0
569 IF (ASSOCIATED(vals)) THEN
570 nlist = cp_sll_val_get_length(vals)
571 END IF
572 imol = 0
573 last_atom = 0
574 DO irk = 1, particles%n_els
575 CALL get_atomic_kind(particles%els(irk)%atomic_kind, name=name)
576 IF (my_shell) THEN
577 s = particles%els(irk)%r
578 IF (scale) THEN
579 r = s
580 CALL real_to_scaled(s, r, cell)
581 ELSE
582 s = s*conv_factor
583 END IF
584 WRITE (unit=line, fmt="(T7,A,3ES25.16E3,1X,I0)") &
585 trim(adjustl(name)), s(1:3), particles%els(irk)%atom_index
586 CALL val_create(my_val, lc_val=line)
587 IF (nlist /= 0) THEN
588 IF (irk == 1) THEN
589 new_pos => vals
590 ELSE
591 new_pos => new_pos%rest
592 END IF
593 old_val => new_pos%first_el
594 CALL val_release(old_val)
595 new_pos%first_el => my_val
596 ELSE
597 IF (irk == 1) THEN
598 NULLIFY (new_pos)
599 CALL cp_sll_val_create(new_pos, first_el=my_val)
600 vals => new_pos
601 ELSE
602 NULLIFY (new_pos%rest)
603 CALL cp_sll_val_create(new_pos%rest, first_el=my_val)
604 new_pos => new_pos%rest
605 END IF
606 END IF
607 NULLIFY (my_val)
608 ELSE
609 IF (last_atom < irk) THEN
610 imol = imol + 1
611 molecule_now => molecules%els(imol)
612 CALL get_molecule(molecule_now, last_atom=last_atom)
613 CALL get_molecule_kind(molecule_now%molecule_kind, molname_generated=molname_generated, &
614 name=my_tag)
615 IF (molname_generated) my_tag = ""
616 END IF
617 ldum = qmmm_ff_precond_only_qm(my_tag)
618 ldum = qmmm_ff_precond_only_qm(name)
619 s = particles%els(irk)%r
620 IF (scale) THEN
621 r = s
622 CALL real_to_scaled(s, r, cell)
623 ELSE
624 s = s*conv_factor
625 END IF
626 IF (len_trim(my_tag) > 0) THEN
627 WRITE (unit=line, fmt="(T7,A,3ES25.16E3,1X,A,1X,I0)") &
628 trim(adjustl(name)), s(1:3), trim(my_tag), imol
629 ELSE
630 WRITE (unit=line, fmt="(T7,A,3ES25.16E3)") &
631 trim(adjustl(name)), s(1:3)
632 END IF
633 CALL val_create(my_val, lc_val=line)
634
635 IF (nlist /= 0) THEN
636 IF (irk == 1) THEN
637 new_pos => vals
638 ELSE
639 new_pos => new_pos%rest
640 END IF
641 old_val => new_pos%first_el
642 CALL val_release(old_val)
643 new_pos%first_el => my_val
644 ELSE
645 IF (irk == 1) THEN
646 NULLIFY (new_pos)
647 CALL cp_sll_val_create(new_pos, first_el=my_val)
648 vals => new_pos
649 ELSE
650 NULLIFY (new_pos%rest)
651 CALL cp_sll_val_create(new_pos%rest, first_el=my_val)
652 new_pos => new_pos%rest
653 END IF
654 END IF
655 NULLIFY (my_val)
656 END IF
657 END DO
658
659 coord_section%values(ik, 1)%list => vals
660
661 CALL timestop(handle)
662
663 END SUBROUTINE section_coord_val_set
664
665! **************************************************************************************************
666!> \brief routine to dump dipoles.. fast implementation
667!> \param dipoles ...
668!> \param dipoles_section ...
669!> \par History
670!> 12.2007 created [teo]
671!> \author Teodoro Laino
672! **************************************************************************************************
673 SUBROUTINE update_dipoles_section(dipoles, dipoles_section)
674
675 REAL(kind=dp), DIMENSION(:, :), POINTER :: dipoles
676 TYPE(section_vals_type), POINTER :: dipoles_section
677
678 CHARACTER(LEN=*), PARAMETER :: routinen = 'update_dipoles_section'
679
680 INTEGER :: handle, ik, irk, nlist, nloop
681 REAL(kind=dp), DIMENSION(:), POINTER :: work
682 TYPE(cp_sll_val_type), POINTER :: new_pos, vals
683 TYPE(section_type), POINTER :: section
684 TYPE(val_type), POINTER :: my_val, old_val
685
686 CALL timeset(routinen, handle)
687 NULLIFY (my_val, old_val, section, vals)
688 cpassert(ASSOCIATED(dipoles_section))
689 cpassert(dipoles_section%ref_count > 0)
690 section => dipoles_section%section
691 ik = section_get_keyword_index(section, "_DEFAULT_KEYWORD_")
692 IF (ik == -2) &
693 CALL cp_abort(__location__, &
694 "section "//trim(section%name)//" does not contain keyword "// &
695 "_DEFAULT_KEYWORD_")
696
697 ! At least one of the two arguments must be present..
698 nloop = SIZE(dipoles, 2)
699 DO
700 IF (SIZE(dipoles_section%values, 2) == 1) EXIT
701 CALL section_vals_add_values(dipoles_section)
702 END DO
703 vals => dipoles_section%values(ik, 1)%list
704 nlist = 0
705 IF (ASSOCIATED(vals)) THEN
706 nlist = cp_sll_val_get_length(vals)
707 END IF
708 DO irk = 1, nloop
709 ALLOCATE (work(3))
710 ! Always stored in A.U.
711 work = dipoles(1:3, irk)
712 CALL val_create(my_val, r_vals_ptr=work)
713
714 IF (nlist /= 0) THEN
715 IF (irk == 1) THEN
716 new_pos => vals
717 ELSE
718 new_pos => new_pos%rest
719 END IF
720 old_val => new_pos%first_el
721 CALL val_release(old_val)
722 new_pos%first_el => my_val
723 ELSE
724 IF (irk == 1) THEN
725 NULLIFY (new_pos)
726 CALL cp_sll_val_create(new_pos, first_el=my_val)
727 vals => new_pos
728 ELSE
729 NULLIFY (new_pos%rest)
730 CALL cp_sll_val_create(new_pos%rest, first_el=my_val)
731 new_pos => new_pos%rest
732 END IF
733 END IF
734 NULLIFY (my_val)
735 END DO
736 dipoles_section%values(ik, 1)%list => vals
737
738 CALL timestop(handle)
739
740 END SUBROUTINE update_dipoles_section
741
742! **************************************************************************************************
743!> \brief routine to dump quadrupoles.. fast implementation
744!> \param quadrupoles ...
745!> \param quadrupoles_section ...
746!> \par History
747!> 12.2007 created [teo]
748!> \author Teodoro Laino
749! **************************************************************************************************
750 SUBROUTINE update_quadrupoles_section(quadrupoles, quadrupoles_section)
751
752 REAL(kind=dp), DIMENSION(:, :, :), POINTER :: quadrupoles
753 TYPE(section_vals_type), POINTER :: quadrupoles_section
754
755 CHARACTER(LEN=*), PARAMETER :: routinen = 'update_quadrupoles_section'
756
757 INTEGER :: handle, i, ik, ind, irk, j, nlist, nloop
758 REAL(kind=dp), DIMENSION(:), POINTER :: work
759 TYPE(cp_sll_val_type), POINTER :: new_pos, vals
760 TYPE(section_type), POINTER :: section
761 TYPE(val_type), POINTER :: my_val, old_val
762
763 CALL timeset(routinen, handle)
764 NULLIFY (my_val, old_val, section, vals)
765 cpassert(ASSOCIATED(quadrupoles_section))
766 cpassert(quadrupoles_section%ref_count > 0)
767 section => quadrupoles_section%section
768 ik = section_get_keyword_index(section, "_DEFAULT_KEYWORD_")
769 IF (ik == -2) &
770 CALL cp_abort(__location__, &
771 "section "//trim(section%name)//" does not contain keyword "// &
772 "_DEFAULT_KEYWORD_")
773
774 ! At least one of the two arguments must be present..
775 nloop = SIZE(quadrupoles, 2)
776 DO
777 IF (SIZE(quadrupoles_section%values, 2) == 1) EXIT
778 CALL section_vals_add_values(quadrupoles_section)
779 END DO
780 vals => quadrupoles_section%values(ik, 1)%list
781 nlist = 0
782 IF (ASSOCIATED(vals)) THEN
783 nlist = cp_sll_val_get_length(vals)
784 END IF
785 DO irk = 1, nloop
786 ALLOCATE (work(6))
787 ! Always stored in A.U.
788 ind = 0
789 DO i = 1, 3
790 DO j = i, 3
791 ind = ind + 1
792 work(ind) = quadrupoles(j, i, irk)
793 END DO
794 END DO
795 CALL val_create(my_val, r_vals_ptr=work)
796
797 IF (nlist /= 0) THEN
798 IF (irk == 1) THEN
799 new_pos => vals
800 ELSE
801 new_pos => new_pos%rest
802 END IF
803 old_val => new_pos%first_el
804 CALL val_release(old_val)
805 new_pos%first_el => my_val
806 ELSE
807 IF (irk == 1) THEN
808 NULLIFY (new_pos)
809 CALL cp_sll_val_create(new_pos, first_el=my_val)
810 vals => new_pos
811 ELSE
812 NULLIFY (new_pos%rest)
813 CALL cp_sll_val_create(new_pos%rest, first_el=my_val)
814 new_pos => new_pos%rest
815 END IF
816 END IF
817 NULLIFY (my_val)
818 END DO
819
820 quadrupoles_section%values(ik, 1)%list => vals
821
822 CALL timestop(handle)
823
824 END SUBROUTINE update_quadrupoles_section
825
826! **************************************************************************************************
827!> \brief Dump atomic, core, or shell coordinates to a file in CP2K &COORD
828!> section format
829!> \param particles ...
830!> \param molecules ...
831!> \param cell ...
832!> \param conv_factor ...
833!> \param output_unit ...
834!> \param core_or_shell ...
835!> \param scaled_coordinates ...
836!> \par History
837!> 07.02.2011 (Creation, MK)
838!> \author Matthias Krack (MK)
839!> \version 1.0
840! **************************************************************************************************
841 SUBROUTINE dump_coordinates_cp2k(particles, molecules, cell, conv_factor, &
842 output_unit, core_or_shell, &
843 scaled_coordinates)
844
845 TYPE(particle_list_type), POINTER :: particles
846 TYPE(molecule_list_type), POINTER :: molecules
847 TYPE(cell_type), POINTER :: cell
848 REAL(kind=dp), INTENT(IN) :: conv_factor
849 INTEGER, INTENT(IN) :: output_unit
850 LOGICAL, INTENT(IN) :: core_or_shell, scaled_coordinates
851
852 CHARACTER(LEN=*), PARAMETER :: routinen = 'dump_coordinates_cp2k'
853
854 CHARACTER(LEN=default_string_length) :: kind_name, tag
855 INTEGER :: handle, imolecule, iparticle, last_atom
856 LOGICAL :: ldum, molname_generated
857 REAL(kind=dp), DIMENSION(3) :: r, s
858 TYPE(molecule_type), POINTER :: molecule
859
860 CALL timeset(routinen, handle)
861
862 cpassert(ASSOCIATED(particles))
863 IF (.NOT. core_or_shell) THEN
864 cpassert(ASSOCIATED(molecules))
865 END IF
866 IF (scaled_coordinates) THEN
867 cpassert(ASSOCIATED(cell))
868 END IF
869
870 kind_name = ""
871 tag = ""
872 imolecule = 0
873 last_atom = 0
874 DO iparticle = 1, particles%n_els
875 CALL get_atomic_kind(particles%els(iparticle)%atomic_kind, name=kind_name)
876 IF (.NOT. core_or_shell) THEN
877 IF (iparticle > last_atom) THEN
878 imolecule = imolecule + 1
879 molecule => molecules%els(imolecule)
880 CALL get_molecule(molecule, last_atom=last_atom)
881 CALL get_molecule_kind(molecule%molecule_kind, &
882 molname_generated=molname_generated, &
883 name=tag)
884 IF (molname_generated) tag = ""
885 END IF
886 ldum = qmmm_ff_precond_only_qm(tag)
887 ldum = qmmm_ff_precond_only_qm(kind_name)
888 END IF
889 IF (scaled_coordinates) THEN
890 CALL real_to_scaled(s, particles%els(iparticle)%r, cell)
891 r(1:3) = s(1:3)
892 ELSE
893 r(1:3) = particles%els(iparticle)%r(1:3)*conv_factor
894 END IF
895 IF (core_or_shell) THEN
896 WRITE (unit=output_unit, fmt="(A,3ES25.16E3,1X,I0)") &
897 trim(adjustl(kind_name)), r(1:3), particles%els(iparticle)%atom_index
898 ELSE
899 IF (len_trim(tag) > 0) THEN
900 WRITE (unit=output_unit, fmt="(A,3ES25.16E3,1X,A,1X,I0)") &
901 trim(adjustl(kind_name)), r(1:3), trim(tag), imolecule
902 ELSE
903 WRITE (unit=output_unit, fmt="(A,3ES25.16E3)") &
904 trim(adjustl(kind_name)), r(1:3)
905 END IF
906 END IF
907 END DO
908
909 CALL timestop(handle)
910
911 END SUBROUTINE dump_coordinates_cp2k
912
913END MODULE input_restart_force_eval
atomic_kind_list_types
represent a simple array based list of the given type
Definition atomic_kind_list_types.F:15
atomic_kind_types
Define the atomic kind types and their sub types.
Definition atomic_kind_types.F:23
atomic_kind_types::get_atomic_kind
subroutine, public get_atomic_kind(atomic_kind, fist_potential, element_symbol, name, mass, kind_number, natom, atom_list, rcov, rvdw, z, qeff, apol, cpol, mm_radius, shell, shell_active, damping)
Get attributes of an atomic kind.
Definition atomic_kind_types.F:138
cell_types
Handles all functions related to the CELL.
Definition cell_types.F:15
cell_types::real_to_scaled
subroutine, public real_to_scaled(s, r, cell)
Transform real to scaled cell coordinates. s=h_inv*r.
Definition cell_types.F:491
cp_control_types
Defines control structures, which contain the parameters and the settings for the DFT-based calculati...
Definition cp_control_types.F:12
cp_files
Utility routines to open and close files. Tracking of preconnections.
Definition cp_files.F:16
cp_files::open_file
subroutine, public open_file(file_name, file_status, file_form, file_action, file_position, file_pad, unit_number, debug, skip_get_unit_number, file_access)
Opens the requested file using a free unit number.
Definition cp_files.F:311
cp_files::close_file
subroutine, public close_file(unit_number, file_status, keep_preconnection)
Close an open file given by its logical unit number. Optionally, keep the file and unit preconnected.
Definition cp_files.F:122
cp_linked_list_input
Definition cp_linked_list_input.F:9
cp_linked_list_input::cp_sll_val_get_length
integer function, public cp_sll_val_get_length(sll)
returns the length of the list
Definition cp_linked_list_input.F:1076
cp_linked_list_input::cp_sll_val_create
subroutine, public cp_sll_val_create(sll, first_el, rest)
allocates and initializes a single linked list
Definition cp_linked_list_input.F:509
cp_subsys_types
types that represent a subsys, i.e. a part of the system
Definition cp_subsys_types.F:16
cp_subsys_types::cp_subsys_get
subroutine, public cp_subsys_get(subsys, ref_count, atomic_kinds, atomic_kind_set, particles, particle_set, local_particles, molecules, molecule_set, molecule_kinds, molecule_kind_set, local_molecules, para_env, colvar_p, shell_particles, core_particles, gci, multipoles, natom, nparticle, ncore, nshell, nkind, atprop, virial, results, cell, cell_ref, use_ref_cell)
returns information about various attributes of the given subsys
Definition cp_subsys_types.F:364
cp_units
unit conversion facility
Definition cp_units.F:30
cp_units::cp_unit_from_cp2k
real(kind=dp) function, public cp_unit_from_cp2k(value, unit_str, defaults, power)
converts from the internal cp2k units to the given unit
Definition cp_units.F:1178
distribution_1d_types
stores a lists of integer that are local to a processor. The idea is that these integers represent ob...
Definition distribution_1d_types.F:24
eip_environment_types
The environment for the empirical interatomic potential methods.
Definition eip_environment_types.F:14
eip_environment_types::eip_env_get
subroutine, public eip_env_get(eip_env, eip_model, eip_energy, eip_energy_var, eip_forces, coord_avg, coord_var, count, subsys, atomic_kind_set, particle_set, local_particles, molecule_kind_set, molecule_set, local_molecules, eip_input, force_env_input, cell, cell_ref, use_ref_cell, eip_kinetic_energy, eip_potential_energy, virial)
Returns various attributes of the eip environment.
Definition eip_environment_types.F:174
force_env_types
Interface for the force calculations.
Definition force_env_types.F:18
force_env_types::use_qmmm
integer, parameter, public use_qmmm
Definition force_env_types.F:83
force_env_types::multiple_fe_list
subroutine, public multiple_fe_list(force_env_sections, root_section, i_force_eval, nforce_eval)
returns the order of the multiple force_env
Definition force_env_types.F:711
force_env_types::use_eip_force
integer, parameter, public use_eip_force
Definition force_env_types.F:83
force_env_types::force_env_get
recursive subroutine, public force_env_get(force_env, in_use, fist_env, qs_env, meta_env, fp_env, subsys, para_env, potential_energy, additional_potential, kinetic_energy, harmonic_shell, kinetic_shell, cell, sub_force_env, qmmm_env, qmmmx_env, eip_env, pwdft_env, globenv, input, force_env_section, method_name_id, root_section, mixed_env, nnp_env, embed_env, ipi_env)
returns various attributes about the force environment
Definition force_env_types.F:342
force_env_types::use_qs_force
integer, parameter, public use_qs_force
Definition force_env_types.F:83
force_env_types::use_nnp_force
integer, parameter, public use_nnp_force
Definition force_env_types.F:83
input_constants
collects all constants needed in input so that they can be used without circular dependencies
Definition input_constants.F:17
input_constants::ehrenfest
integer, parameter, public ehrenfest
Definition input_constants.F:124
input_constants::mon_car_run
integer, parameter, public mon_car_run
Definition input_constants.F:124
input_constants::mol_dyn_run
integer, parameter, public mol_dyn_run
Definition input_constants.F:124
input_constants::use_rt_restart
integer, parameter, public use_rt_restart
Definition input_constants.F:948
input_constants::pint_run
integer, parameter, public pint_run
Definition input_constants.F:124
input_constants::do_coord_cp2k
integer, parameter, public do_coord_cp2k
Definition input_constants.F:419
input_cp2k_restarts_util
Definition input_cp2k_restarts_util.F:8
input_cp2k_restarts_util::section_velocity_val_set
subroutine, public section_velocity_val_set(velocity_section, particles, velocity, conv_factor)
routine to dump velocities.. fast implementation
Definition input_cp2k_restarts_util.F:44
input_restart_force_eval
Definition input_restart_force_eval.F:8
input_restart_force_eval::update_subsys
subroutine, public update_subsys(subsys_section, force_env, skip_vel_section, write_binary_restart_file)
Updates the subsys section of the input file.
Definition input_restart_force_eval.F:315
input_restart_force_eval::update_force_eval
subroutine, public update_force_eval(force_env, root_section, write_binary_restart_file, respa)
Updates the force_eval section of the input file.
Definition input_restart_force_eval.F:86
input_restart_rng
Definition input_restart_rng.F:8
input_restart_rng::section_rng_val_set
subroutine, public section_rng_val_set(rng_section, nsize, ascii)
routine to dump rngs.. fast implementation
Definition input_restart_rng.F:44
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_val_unset
subroutine, public section_vals_val_unset(section_vals, keyword_name, i_rep_section, i_rep_val)
unsets (removes) the requested value (if it is a keyword repetitions removes the repetition,...
Definition input_section_types.F:1354
input_section_types::section_vals_val_set
subroutine, public section_vals_val_set(section_vals, keyword_name, i_rep_section, i_rep_val, val, l_val, i_val, r_val, c_val, l_vals_ptr, i_vals_ptr, r_vals_ptr, c_vals_ptr)
sets the requested value
Definition input_section_types.F:1223
input_section_types::section_get_keyword_index
integer function, public section_get_keyword_index(section, keyword_name)
returns the index of the requested keyword (or -2 if not found)
Definition input_section_types.F:436
input_section_types::section_vals_remove_values
subroutine, public section_vals_remove_values(section_vals)
removes the values of a repetition of the section
Definition input_section_types.F:896
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:731
input_section_types::section_vals_get
subroutine, public section_vals_get(section_vals, ref_count, n_repetition, n_subs_vals_rep, section, explicit)
returns various attributes about the section_vals
Definition input_section_types.F:704
input_section_types::section_vals_get_subs_vals3
type(section_vals_type) function, pointer, public section_vals_get_subs_vals3(section_vals, subsection_name, i_rep_section)
returns the values of the n-th non default subsection (null if no such section exists (not so many no...
Definition input_section_types.F:824
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:1047
input_section_types::section_vals_add_values
subroutine, public section_vals_add_values(section_vals)
adds the place to store the values of a repetition of the section
Definition input_section_types.F:848
input_val_types
a wrapper for basic fortran types.
Definition input_val_types.F:14
input_val_types::val_create
subroutine, public val_create(val, l_val, l_vals, l_vals_ptr, i_val, i_vals, i_vals_ptr, r_val, r_vals, r_vals_ptr, c_val, c_vals, c_vals_ptr, lc_val, lc_vals, lc_vals_ptr, enum)
creates a keyword value
Definition input_val_types.F:103
input_val_types::val_release
subroutine, public val_release(val)
releases the given val
Definition input_val_types.F:262
kinds
Defines the basic variable types.
Definition kinds.F:23
kinds::dp
integer, parameter, public dp
Definition kinds.F:34
kinds::default_string_length
integer, parameter, public default_string_length
Definition kinds.F:57
message_passing
Interface to the message passing library MPI.
Definition message_passing.F:23
molecule_kind_types
Define the molecule kind structure types and the corresponding functionality.
Definition molecule_kind_types.F:16
molecule_kind_types::get_molecule_kind
subroutine, public get_molecule_kind(molecule_kind, atom_list, bond_list, bend_list, ub_list, impr_list, opbend_list, colv_list, fixd_list, g3x3_list, g4x6_list, vsite_list, torsion_list, shell_list, name, mass, charge, kind_number, natom, nbend, nbond, nub, nimpr, nopbend, nconstraint, nconstraint_fixd, nfixd, ncolv, ng3x3, ng4x6, nvsite, nfixd_restraint, ng3x3_restraint, ng4x6_restraint, nvsite_restraint, nrestraints, nmolecule, nsgf, nshell, ntorsion, molecule_list, nelectron, nelectron_alpha, nelectron_beta, bond_kind_set, bend_kind_set, ub_kind_set, impr_kind_set, opbend_kind_set, torsion_kind_set, molname_generated)
Get informations about a molecule kind.
Definition molecule_kind_types.F:512
molecule_list_types
represent a simple array based list of the given type
Definition molecule_list_types.F:15
molecule_types
Define the data structure for the molecule information.
Definition molecule_types.F:16
molecule_types::get_molecule
subroutine, public get_molecule(molecule, molecule_kind, lmi, lci, lg3x3, lg4x6, lcolv, first_atom, last_atom, first_shell, last_shell)
Get components from a molecule data set.
Definition molecule_types.F:326
multipole_types
Multipole structure: for multipole (fixed and induced) in FF based MD.
Definition multipole_types.F:12
nnp_environment_types
Data types for neural network potentials.
Definition nnp_environment_types.F:13
nnp_environment_types::nnp_env_get
subroutine, public nnp_env_get(nnp_env, nnp_forces, subsys, atomic_kind_set, particle_set, local_particles, molecule_kind_set, molecule_set, local_molecules, nnp_input, force_env_input, cell, cell_ref, use_ref_cell, nnp_potential_energy, virial)
Returns various attributes of the nnp environment.
Definition nnp_environment_types.F:418
parallel_rng_types
Parallel (pseudo)random number generator (RNG) for multiple streams and substreams of random numbers.
Definition parallel_rng_types.F:51
parallel_rng_types::rng_record_length
integer, parameter, public rng_record_length
Definition parallel_rng_types.F:65
particle_list_types
represent a simple array based list of the given type
Definition particle_list_types.F:15
qmmm_ff_fist
Definition qmmm_ff_fist.F:11
qmmm_ff_fist::qmmm_ff_precond_only_qm
logical function, public qmmm_ff_precond_only_qm(id1, id2, id3, id4, is_link)
This function handles the atom names and modifies the "_QM_" prefix, in order to find the parameters ...
Definition qmmm_ff_fist.F:39
qs_environment_types
Definition qs_environment_types.F:14
qs_environment_types::get_qs_env
subroutine, public get_qs_env(qs_env, atomic_kind_set, qs_kind_set, cell, super_cell, cell_ref, use_ref_cell, kpoints, dft_control, mos, sab_orb, sab_all, qmmm, qmmm_periodic, mimic, sac_ae, sac_ppl, sac_lri, sap_ppnl, sab_vdw, sab_scp, sap_oce, sab_lrc, sab_se, sab_xtbe, sab_tbe, sab_core, sab_xb, sab_xtb_pp, sab_xtb_nonbond, sab_almo, sab_kp, sab_kp_nosym, sab_cneo, particle_set, energy, force, matrix_h, matrix_h_im, matrix_ks, matrix_ks_im, matrix_vxc, run_rtp, rtp, matrix_h_kp, matrix_h_im_kp, matrix_ks_kp, matrix_ks_im_kp, matrix_vxc_kp, kinetic_kp, matrix_s_kp, matrix_w_kp, matrix_s_ri_aux_kp, matrix_s, matrix_s_ri_aux, matrix_w, matrix_p_mp2, matrix_p_mp2_admm, rho, rho_xc, pw_env, ewald_env, ewald_pw, active_space, mpools, input, para_env, blacs_env, scf_control, rel_control, kinetic, qs_charges, vppl, xcint_weights, rho_core, rho_nlcc, rho_nlcc_g, ks_env, ks_qmmm_env, wf_history, scf_env, local_particles, local_molecules, distribution_2d, dbcsr_dist, molecule_kind_set, molecule_set, subsys, cp_subsys, oce, local_rho_set, rho_atom_set, task_list, task_list_soft, rho0_atom_set, rho0_mpole, rhoz_set, rhoz_cneo_set, ecoul_1c, rho0_s_rs, rho0_s_gs, rhoz_cneo_s_rs, rhoz_cneo_s_gs, do_kpoints, has_unit_metric, requires_mo_derivs, mo_derivs, mo_loc_history, nkind, natom, nelectron_total, nelectron_spin, efield, neighbor_list_id, linres_control, xas_env, virial, cp_ddapc_env, cp_ddapc_ewald, outer_scf_history, outer_scf_ihistory, x_data, et_coupling, dftb_potential, results, se_taper, se_store_int_env, se_nddo_mpole, se_nonbond_env, admm_env, lri_env, lri_density, exstate_env, ec_env, harris_env, dispersion_env, gcp_env, vee, rho_external, external_vxc, mask, mp2_env, bs_env, kg_env, wanniercentres, atprop, ls_scf_env, do_transport, transport_env, v_hartree_rspace, s_mstruct_changed, rho_changed, potential_changed, forces_up_to_date, mscfg_env, almo_scf_env, gradient_history, variable_history, embed_pot, spin_embed_pot, polar_env, mos_last_converged, eeq, rhs, do_rixs, tb_tblite)
Get the QUICKSTEP environment.
Definition qs_environment_types.F:530
string_utilities
Utilities for string manipulations.
Definition string_utilities.F:16
string_utilities::string_to_ascii
subroutine, public string_to_ascii(string, nascii)
Convert a string to sequence of integer numbers.
Definition string_utilities.F:3255
virial_types
Definition virial_types.F:13
atomic_kind_list_types::atomic_kind_list_type
represent a list of objects
Definition atomic_kind_list_types.F:45
cell_types::cell_type
Type defining parameters related to the simulation cell.
Definition cell_types.F:60
cp_control_types::dft_control_type
Definition cp_control_types.F:633
cp_linked_list_input::cp_sll_val_type
represent a single linked list that stores pointers to the elements
Definition cp_linked_list_input.F:341
cp_subsys_types::cp_subsys_type
represents a system: atoms, molecules, their pos,vel,...
Definition cp_subsys_types.F:89
distribution_1d_types::distribution_1d_type
structure to store local (to a processor) ordered lists of integers.
Definition distribution_1d_types.F:62
force_env_types::force_env_type
wrapper to abstract the force evaluation of the various methods
Definition force_env_types.F:145
input_section_types::section_type
represent a section of the input file
Definition input_section_types.F:102
input_section_types::section_vals_type
stores the values of a section
Definition input_section_types.F:127
input_val_types::val_type
a type to have a wrapper that stores any basic fortran type
Definition input_val_types.F:63
message_passing::mp_para_env_type
stores all the informations relevant to an mpi environment
Definition message_passing.F:743
molecule_list_types::molecule_list_type
represent a list of objects
Definition molecule_list_types.F:46
molecule_types::molecule_type
Definition molecule_types.F:159
multipole_types::multipole_type
Define multipole type.
Definition multipole_types.F:44
particle_list_types::particle_list_type
represent a list of objects
Definition particle_list_types.F:46
virial_types::virial_type
Definition virial_types.F:26