(git:1155b05)
Loading...
Searching...
No Matches
xc_atom.F
Go to the documentation of this file.
1!--------------------------------------------------------------------------------------------------!
2! CP2K: A general program to perform molecular dynamics simulations !
3! Copyright 2000-2026 CP2K developers group <https://cp2k.org> !
4! !
5! SPDX-License-Identifier: GPL-2.0-or-later !
6!--------------------------------------------------------------------------------------------------!
7
8! **************************************************************************************************
9MODULE xc_atom
10
11 USE cp_linked_list_xc_deriv, ONLY: cp_sll_xc_deriv_next,&
12 cp_sll_xc_deriv_type
13 USE input_section_types, ONLY: section_vals_get_subs_vals,&
14 section_vals_type
15 USE kinds, ONLY: dp
16 USE pw_pool_types, ONLY: pw_pool_type
17 USE pw_types, ONLY: pw_r3d_rs_type
18 USE xc, ONLY: divide_by_norm_drho,&
19 xc_calc_2nd_deriv_analytical
20 USE xc_derivative_desc, ONLY: &
21 deriv_norm_drho, deriv_norm_drhoa, deriv_norm_drhob, deriv_rho, deriv_rhoa, deriv_rhob, &
22 deriv_tau, deriv_tau_a, deriv_tau_b
23 USE xc_derivative_set_types, ONLY: xc_derivative_set_type,&
24 xc_dset_get_derivative
25 USE xc_derivative_types, ONLY: xc_derivative_get,&
26 xc_derivative_type
27 USE xc_derivatives, ONLY: xc_functionals_eval
28 USE xc_rho_cflags_types, ONLY: xc_rho_cflags_type
29 USE xc_rho_set_types, ONLY: xc_rho_set_get,&
30 xc_rho_set_type
31#include "../base/base_uses.f90"
32
33 IMPLICIT NONE
34
35 PRIVATE
36
37 CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'xc_atom'
38
39 PUBLIC :: vxc_of_r_new, vxc_of_r_epr, xc_rho_set_atom_update, xc_2nd_deriv_of_r, fill_rho_set
40
41CONTAINS
42
43! **************************************************************************************************
44!> \brief ...
45!> \param xc_fun_section ...
46!> \param rho_set ...
47!> \param deriv_set ...
48!> \param deriv_order ...
49!> \param needs ...
50!> \param w ...
51!> \param lsd ...
52!> \param na ...
53!> \param nr ...
54!> \param exc ...
55!> \param vxc ...
56!> \param vxg ...
57!> \param vtau ...
58!> \param energy_only ...
59!> \param adiabatic_rescale_factor ...
60! **************************************************************************************************
61 SUBROUTINE vxc_of_r_new(xc_fun_section, rho_set, deriv_set, deriv_order, needs, w, &
62 lsd, na, nr, exc, vxc, vxg, vtau, &
63 energy_only, adiabatic_rescale_factor)
64
65! This routine updates rho_set by giving to it the rho and drho that are needed.
66! Since for the local densities rho1_h and rho1_s local grids are used it is not possible
67! to call xc_rho_set_update.
68! As input of this routine one gets rho and drho on a one dimensional grid.
69! The grid is the angular grid corresponding to a given point ir_pnt on the radial grid.
70! The derivatives are calculated on this one dimensional grid, the results are stored in
71! exc, vxc(1:na,ir_pnt,ispin), vxg(1:na,ir_pnt,ispin), vxg_cross(1:na,ir_pnt,ispin)
72! Afterwords the arrays containing the derivatives are put to zero so that the routine
73! can safely be called for the next radial point ir_pnt
74
75 TYPE(section_vals_type), POINTER :: xc_fun_section
76 TYPE(xc_rho_set_type), INTENT(IN) :: rho_set
77 TYPE(xc_derivative_set_type), INTENT(IN) :: deriv_set
78 INTEGER, INTENT(in) :: deriv_order
79 TYPE(xc_rho_cflags_type), INTENT(IN) :: needs
80 REAL(dp), DIMENSION(:, :), INTENT(IN) :: w
81 LOGICAL, INTENT(IN) :: lsd
82 INTEGER, INTENT(in) :: na, nr
83 REAL(dp) :: exc
84 REAL(dp), DIMENSION(:, :, :), POINTER :: vxc
85 REAL(dp), DIMENSION(:, :, :, :), POINTER :: vxg
86 REAL(dp), DIMENSION(:, :, :), POINTER :: vtau
87 LOGICAL, INTENT(IN), OPTIONAL :: energy_only
88 REAL(dp), INTENT(IN), OPTIONAL :: adiabatic_rescale_factor
89
90 CHARACTER(LEN=*), PARAMETER :: routinen = 'vxc_of_r_new'
91
92 INTEGER :: handle, ia, idir, ir
93 LOGICAL :: gradient_f, my_only_energy
94 REAL(dp) :: my_adiabatic_rescale_factor
95 REAL(dp), DIMENSION(:, :, :), POINTER :: deriv_data
96 REAL(kind=dp) :: drho_cutoff
97 TYPE(xc_derivative_type), POINTER :: deriv_att
98
99 CALL timeset(routinen, handle)
100 my_only_energy = .false.
101 IF (PRESENT(energy_only)) my_only_energy = energy_only
102
103 IF (PRESENT(adiabatic_rescale_factor)) THEN
104 my_adiabatic_rescale_factor = adiabatic_rescale_factor
105 ELSE
106 my_adiabatic_rescale_factor = 1.0_dp
107 END IF
108
109 gradient_f = (needs%drho_spin .OR. needs%norm_drho_spin .OR. &
110 needs%drho .OR. needs%norm_drho)
111
112 ! Calculate the derivatives
113 CALL xc_functionals_eval(xc_fun_section, &
114 lsd=lsd, &
115 rho_set=rho_set, &
116 deriv_set=deriv_set, &
117 deriv_order=deriv_order)
118
119 CALL xc_rho_set_get(rho_set, drho_cutoff=drho_cutoff)
120
121 NULLIFY (deriv_data)
122
123 ! EXC energy
124 deriv_att => xc_dset_get_derivative(deriv_set, [INTEGER::])
125 exc = 0.0_dp
126 IF (ASSOCIATED(deriv_att)) THEN
127 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
128 DO ir = 1, nr
129 DO ia = 1, na
130 exc = exc + deriv_data(ia, ir, 1)*w(ia, ir)
131 END DO
132 END DO
133 NULLIFY (deriv_data)
134 END IF
135 ! Calculate the potential only if needed
136 IF (.NOT. my_only_energy) THEN
137 ! Derivative with respect to the density
138 IF (lsd) THEN
139 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_rhoa])
140 IF (ASSOCIATED(deriv_att)) THEN
141 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
142 vxc(:, :, 1) = deriv_data(:, :, 1)*w(:, :)*my_adiabatic_rescale_factor
143 NULLIFY (deriv_data)
144 END IF
145 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_rhob])
146 IF (ASSOCIATED(deriv_att)) THEN
147 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
148 vxc(:, :, 2) = deriv_data(:, :, 1)*w(:, :)*my_adiabatic_rescale_factor
149 NULLIFY (deriv_data)
150 END IF
151 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_rho])
152 IF (ASSOCIATED(deriv_att)) THEN
153 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
154 vxc(:, :, 1) = vxc(:, :, 1) + deriv_data(:, :, 1)*w(:, :)*my_adiabatic_rescale_factor
155 vxc(:, :, 2) = vxc(:, :, 2) + deriv_data(:, :, 1)*w(:, :)*my_adiabatic_rescale_factor
156 NULLIFY (deriv_data)
157 END IF
158 ELSE
159 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_rho])
160 IF (ASSOCIATED(deriv_att)) THEN
161 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
162 vxc(:, :, 1) = deriv_data(:, :, 1)*w(:, :)*my_adiabatic_rescale_factor
163 NULLIFY (deriv_data)
164 END IF
165 END IF ! lsd
166
167 ! Derivatives with respect to the gradient
168 IF (lsd) THEN
169 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_norm_drhoa])
170 IF (ASSOCIATED(deriv_att)) THEN
171 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
172 DO ir = 1, nr
173 DO ia = 1, na
174 DO idir = 1, 3
175 IF (rho_set%norm_drhoa(ia, ir, 1) > drho_cutoff) THEN
176 vxg(idir, ia, ir, 1) = rho_set%drhoa(idir)%array(ia, ir, 1)* &
177 deriv_data(ia, ir, 1)*w(ia, ir)/ &
178 rho_set%norm_drhoa(ia, ir, 1)*my_adiabatic_rescale_factor
179 ELSE
180 vxg(idir, ia, ir, 1) = 0.0_dp
181 END IF
182 END DO
183 END DO
184 END DO
185 NULLIFY (deriv_data)
186 END IF
187 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_norm_drhob])
188 IF (ASSOCIATED(deriv_att)) THEN
189 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
190 DO ir = 1, nr
191 DO ia = 1, na
192 DO idir = 1, 3
193 IF (rho_set%norm_drhob(ia, ir, 1) > drho_cutoff) THEN
194 vxg(idir, ia, ir, 2) = rho_set%drhob(idir)%array(ia, ir, 1)* &
195 deriv_data(ia, ir, 1)*w(ia, ir)/ &
196 rho_set%norm_drhob(ia, ir, 1)*my_adiabatic_rescale_factor
197 ELSE
198 vxg(idir, ia, ir, 2) = 0.0_dp
199 END IF
200 END DO
201 END DO
202 END DO
203 NULLIFY (deriv_data)
204 END IF
205 ! Cross Terms
206 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_norm_drho])
207 IF (ASSOCIATED(deriv_att)) THEN
208 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
209 DO ir = 1, nr
210 DO ia = 1, na
211 DO idir = 1, 3
212 IF (rho_set%norm_drho(ia, ir, 1) > drho_cutoff) THEN
213 vxg(idir, ia, ir, 1:2) = &
214 vxg(idir, ia, ir, 1:2) + ( &
215 rho_set%drhoa(idir)%array(ia, ir, 1) + &
216 rho_set%drhob(idir)%array(ia, ir, 1))* &
217 deriv_data(ia, ir, 1)*w(ia, ir)/rho_set%norm_drho(ia, ir, 1)* &
218 my_adiabatic_rescale_factor
219 END IF
220 END DO
221 END DO
222 END DO
223 NULLIFY (deriv_data)
224 END IF
225 ELSE
226 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_norm_drho])
227 IF (ASSOCIATED(deriv_att)) THEN
228 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
229 DO ir = 1, nr
230 DO ia = 1, na
231 IF (rho_set%norm_drho(ia, ir, 1) > drho_cutoff) THEN
232 DO idir = 1, 3
233 vxg(idir, ia, ir, 1) = rho_set%drho(idir)%array(ia, ir, 1)* &
234 deriv_data(ia, ir, 1)*w(ia, ir)/ &
235 rho_set%norm_drho(ia, ir, 1)*my_adiabatic_rescale_factor
236 END DO
237 ELSE
238 vxg(1:3, ia, ir, 1) = 0.0_dp
239 END IF
240 END DO
241 END DO
242 NULLIFY (deriv_data)
243 END IF
244 END IF ! lsd
245 ! Derivative with respect to tau
246 IF (lsd) THEN
247 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_tau_a])
248 IF (ASSOCIATED(deriv_att)) THEN
249 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
250 vtau(:, :, 1) = deriv_data(:, :, 1)*w(:, :)*my_adiabatic_rescale_factor
251 NULLIFY (deriv_data)
252 END IF
253 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_tau_b])
254 IF (ASSOCIATED(deriv_att)) THEN
255 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
256 vtau(:, :, 2) = deriv_data(:, :, 1)*w(:, :)*my_adiabatic_rescale_factor
257 NULLIFY (deriv_data)
258 END IF
259 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_tau])
260 IF (ASSOCIATED(deriv_att)) THEN
261 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
262 vtau(:, :, 1) = vtau(:, :, 1) + deriv_data(:, :, 1)*w(:, :)*my_adiabatic_rescale_factor
263 vtau(:, :, 2) = vtau(:, :, 2) + deriv_data(:, :, 1)*w(:, :)*my_adiabatic_rescale_factor
264 NULLIFY (deriv_data)
265 END IF
266 ELSE
267 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_tau])
268 IF (ASSOCIATED(deriv_att)) THEN
269 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
270 vtau(:, :, 1) = deriv_data(:, :, 1)*w(:, :)*my_adiabatic_rescale_factor
271 NULLIFY (deriv_data)
272 END IF
273 END IF ! lsd
274 END IF ! only_energy
275
276 CALL timestop(handle)
277
278 END SUBROUTINE vxc_of_r_new
279
280! **************************************************************************************************
281!> \brief Specific EPR version of vxc_of_r_new
282!> \param xc_fun_section ...
283!> \param rho_set ...
284!> \param deriv_set ...
285!> \param needs ...
286!> \param w ...
287!> \param lsd ...
288!> \param na ...
289!> \param nr ...
290!> \param exc ...
291!> \param vxc ...
292!> \param vxg ...
293!> \param vtau ...
294! **************************************************************************************************
295 SUBROUTINE vxc_of_r_epr(xc_fun_section, rho_set, deriv_set, needs, w, &
296 lsd, na, nr, exc, vxc, vxg, vtau)
297
298 TYPE(section_vals_type), POINTER :: xc_fun_section
299 TYPE(xc_rho_set_type), INTENT(IN) :: rho_set
300 TYPE(xc_derivative_set_type), INTENT(IN) :: deriv_set
301 TYPE(xc_rho_cflags_type), INTENT(IN) :: needs
302 REAL(dp), DIMENSION(:, :), INTENT(IN) :: w
303 LOGICAL, INTENT(IN) :: lsd
304 INTEGER, INTENT(in) :: na, nr
305 REAL(dp) :: exc
306 REAL(dp), DIMENSION(:, :, :), POINTER :: vxc
307 REAL(dp), DIMENSION(:, :, :, :), POINTER :: vxg
308 REAL(dp), DIMENSION(:, :, :), POINTER :: vtau
309
310 CHARACTER(LEN=*), PARAMETER :: routinen = 'vxc_of_r_epr'
311
312 INTEGER :: handle, ia, idir, ir, my_deriv_order
313 LOGICAL :: gradient_f
314 REAL(dp) :: my_adiabatic_rescale_factor
315 REAL(dp), DIMENSION(:, :, :), POINTER :: deriv_data
316 REAL(kind=dp) :: drho_cutoff
317 TYPE(xc_derivative_type), POINTER :: deriv_att
318
319 CALL timeset(routinen, handle)
320
321 mark_used(vxc)
322 mark_used(vtau)
323
324 my_adiabatic_rescale_factor = 1.0_dp
325 my_deriv_order = 2
326
327 gradient_f = (needs%drho_spin .OR. needs%norm_drho_spin .OR. &
328 needs%drho .OR. needs%norm_drho)
329
330 ! Calculate the derivatives
331 CALL xc_functionals_eval(xc_fun_section, &
332 lsd=lsd, &
333 rho_set=rho_set, &
334 deriv_set=deriv_set, &
335 deriv_order=my_deriv_order)
336
337 CALL xc_rho_set_get(rho_set, drho_cutoff=drho_cutoff)
338
339 NULLIFY (deriv_data)
340
341 ! nabla v_xc (using the vxg arrays)
342 ! there's no point doing this when lsd = false
343 IF (lsd) THEN
344 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_rhoa, deriv_rhoa])
345 IF (ASSOCIATED(deriv_att)) THEN
346 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
347 DO ir = 1, nr
348 DO ia = 1, na
349 DO idir = 1, 3
350 vxg(idir, ia, ir, 1) = rho_set%drhoa(idir)%array(ia, ir, 1)* &
351 deriv_data(ia, ir, 1)
352 END DO !idir
353 END DO !ia
354 END DO !ir
355 NULLIFY (deriv_data)
356 END IF
357 deriv_att => xc_dset_get_derivative(deriv_set, [deriv_rhob, deriv_rhob])
358 IF (ASSOCIATED(deriv_att)) THEN
359 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
360 DO ir = 1, nr
361 DO ia = 1, na
362 DO idir = 1, 3
363 vxg(idir, ia, ir, 2) = rho_set%drhob(idir)%array(ia, ir, 1)* &
364 deriv_data(ia, ir, 1)
365 END DO !idir
366 END DO !ia
367 END DO !ir
368 NULLIFY (deriv_data)
369 END IF
370 END IF
371 ! EXC energy ! is that needed for epr?
372 deriv_att => xc_dset_get_derivative(deriv_set, [INTEGER::])
373 exc = 0.0_dp
374 IF (ASSOCIATED(deriv_att)) THEN
375 CALL xc_derivative_get(deriv_att, deriv_data=deriv_data)
376 DO ir = 1, nr
377 DO ia = 1, na
378 exc = exc + deriv_data(ia, ir, 1)*w(ia, ir)
379 END DO
380 END DO
381 NULLIFY (deriv_data)
382 END IF
383
384 CALL timestop(handle)
385
386 END SUBROUTINE vxc_of_r_epr
387
388! **************************************************************************************************
389!> \brief ...
390!> \param rho_set ...
391!> \param rho1_set ...
392!> \param xc_section ...
393!> \param deriv_set ...
394!> \param w ...
395!> \param vxc ...
396!> \param vxg ...
397!> \param do_triplet ...
398!> \param do_sf ...
399! **************************************************************************************************
400 SUBROUTINE xc_2nd_deriv_of_r(rho_set, rho1_set, xc_section, &
401 deriv_set, w, vxc, vxg, do_triplet, do_sf)
402
403! As input of this routine one gets rho and drho on a one dimensional grid.
404! The grid is the angular grid corresponding to a given point ir on the radial grid.
405! The derivatives are calculated on this one dimensional grid, the results are stored in
406! vxc(1:na,ir,ispin), vxg(1:na,ir,ispin), vxg_cross(1:na,ir,ispin)
407! Afterwords the arrays containing the derivatives are put to zero so that the routine
408! can safely be called for the next radial point ir
409
410 TYPE(xc_rho_set_type), INTENT(IN) :: rho_set, rho1_set
411 TYPE(section_vals_type), POINTER :: xc_section
412 TYPE(xc_derivative_set_type), INTENT(INOUT) :: deriv_set
413 REAL(dp), DIMENSION(:, :), INTENT(IN) :: w
414 REAL(dp), CONTIGUOUS, DIMENSION(:, :, :), POINTER :: vxc
415 REAL(dp), DIMENSION(:, :, :, :), POINTER :: vxg
416 LOGICAL, INTENT(IN), OPTIONAL :: do_triplet, do_sf
417
418 CHARACTER(LEN=*), PARAMETER :: routinen = 'xc_2nd_deriv_of_r'
419
420 INTEGER :: handle, ispin, nspins
421 LOGICAL :: lsd, my_do_sf
422 REAL(dp) :: drho_cutoff, my_fac_triplet
423 TYPE(cp_sll_xc_deriv_type), POINTER :: pos
424 TYPE(pw_pool_type), POINTER :: pw_pool
425 TYPE(pw_r3d_rs_type), DIMENSION(:), POINTER :: vxc_pw, vxc_tau_pw
426 TYPE(section_vals_type), POINTER :: xc_fun_section
427 TYPE(xc_derivative_type), POINTER :: deriv_att
428
429 CALL timeset(routinen, handle)
430
431 nspins = SIZE(vxc, 3)
432 lsd = (nspins == 2)
433 IF (ASSOCIATED(rho_set%rhoa)) THEN
434 lsd = .true.
435 END IF
436 my_fac_triplet = 1.0_dp
437 IF (PRESENT(do_triplet)) THEN
438 IF (do_triplet) my_fac_triplet = -1.0_dp
439 END IF
440
441 my_do_sf = .false.
442 IF (PRESENT(do_sf)) my_do_sf = do_sf
443
444 CALL xc_rho_set_get(rho_set, drho_cutoff=drho_cutoff)
445 xc_fun_section => section_vals_get_subs_vals(xc_section, &
446 "XC_FUNCTIONAL")
447
448 ! Calculate the derivatives
449 CALL xc_functionals_eval(xc_fun_section, &
450 lsd=lsd, &
451 rho_set=rho_set, &
452 deriv_set=deriv_set, &
453 deriv_order=2)
454
455 CALL divide_by_norm_drho(deriv_set, rho_set, lsd)
456
457 ! multiply by w
458 pos => deriv_set%derivs
459 DO WHILE (cp_sll_xc_deriv_next(pos, el_att=deriv_att))
460 deriv_att%deriv_data(:, :, 1) = w(:, :)*deriv_att%deriv_data(:, :, 1)
461 END DO
462
463 NULLIFY (pw_pool)
464 ALLOCATE (vxc_pw(nspins))
465 DO ispin = 1, nspins
466 vxc_pw(ispin)%array => vxc(:, :, ispin:ispin)
467 END DO
468
469 NULLIFY (vxc_tau_pw)
470
471 CALL xc_calc_2nd_deriv_analytical(vxc_pw, vxc_tau_pw, deriv_set, rho_set, rho1_set, pw_pool, &
472 xc_section, gapw=.true., vxg=vxg, &
473 tddfpt_fac=my_fac_triplet, spinflip=do_sf)
474
475 DEALLOCATE (vxc_pw)
476
477 ! zero the derivative data for the next call
478 pos => deriv_set%derivs
479 DO WHILE (cp_sll_xc_deriv_next(pos, el_att=deriv_att))
480 deriv_att%deriv_data = 0.0_dp
481 END DO
482
483 CALL timestop(handle)
484
485 END SUBROUTINE xc_2nd_deriv_of_r
486
487! **************************************************************************************************
488!> \brief ...
489!> \param rho_set ...
490!> \param needs ...
491!> \param nspins ...
492!> \param bo ...
493! **************************************************************************************************
494 SUBROUTINE xc_rho_set_atom_update(rho_set, needs, nspins, bo)
495
496! This routine allocates the storage arrays for rho and drho
497! In calculate_vxc_atom this is called once for each atomic_kind,
498! After the loop over all the atoms of the kind and over all the points
499! of the radial grid for each atom, rho_set is deallocated.
500! Within the same kind, at each new point on the radial grid, the rho_set
501! arrays rho and drho are overwritten.
502
503 TYPE(xc_rho_set_type), INTENT(INOUT) :: rho_set
504 TYPE(xc_rho_cflags_type), INTENT(IN) :: needs
505 INTEGER, INTENT(IN) :: nspins
506 INTEGER, DIMENSION(2, 3), INTENT(IN) :: bo
507
508 INTEGER :: idir
509
510 SELECT CASE (nspins)
511 CASE (1)
512! What is this for?
513 IF (needs%rho_1_3) THEN
514 NULLIFY (rho_set%rho_1_3)
515 ALLOCATE (rho_set%rho_1_3(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
516 rho_set%owns%rho_1_3 = .true.
517 rho_set%has%rho_1_3 = .false.
518 END IF
519! Allocate the storage space for the density
520 IF (needs%rho) THEN
521 NULLIFY (rho_set%rho)
522 ALLOCATE (rho_set%rho(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
523 rho_set%owns%rho = .true.
524 rho_set%has%rho = .false.
525 END IF
526! Allocate the storage space for the norm of the gradient of the density
527 IF (needs%norm_drho) THEN
528 NULLIFY (rho_set%norm_drho)
529 ALLOCATE (rho_set%norm_drho(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
530 rho_set%owns%norm_drho = .true.
531 rho_set%has%norm_drho = .false.
532 END IF
533! Allocate the storage space for the three components of the gradient of the density
534 IF (needs%drho) THEN
535 DO idir = 1, 3
536 NULLIFY (rho_set%drho(idir)%array)
537 ALLOCATE (rho_set%drho(idir)%array(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
538 END DO
539 rho_set%owns%drho = .true.
540 rho_set%has%drho = .false.
541 END IF
542 CASE (2)
543! Allocate the storage space for the total density
544 IF (needs%rho) THEN
545 ! this should never be the case unless you use LDA functionals with LSD
546 NULLIFY (rho_set%rho)
547 ALLOCATE (rho_set%rho(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
548 rho_set%owns%rho = .true.
549 rho_set%has%rho = .false.
550 END IF
551! What is this for?
552 IF (needs%rho_1_3) THEN
553 NULLIFY (rho_set%rho_1_3)
554 ALLOCATE (rho_set%rho_1_3(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
555 rho_set%owns%rho_1_3 = .true.
556 rho_set%has%rho_1_3 = .false.
557 END IF
558! What is this for?
559 IF (needs%rho_spin_1_3) THEN
560 NULLIFY (rho_set%rhoa_1_3, rho_set%rhob_1_3)
561 ALLOCATE (rho_set%rhoa_1_3(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
562 ALLOCATE (rho_set%rhob_1_3(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
563 rho_set%owns%rho_spin_1_3 = .true.
564 rho_set%has%rho_spin_1_3 = .false.
565 END IF
566! Allocate the storage space for the spin densities rhoa and rhob
567 IF (needs%rho_spin) THEN
568 NULLIFY (rho_set%rhoa, rho_set%rhob)
569 ALLOCATE (rho_set%rhoa(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
570 ALLOCATE (rho_set%rhob(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
571 rho_set%owns%rho_spin = .true.
572 rho_set%has%rho_spin = .false.
573 END IF
574! Allocate the storage space for the norm of the gradient of the total density
575 IF (needs%norm_drho) THEN
576 NULLIFY (rho_set%norm_drho)
577 ALLOCATE (rho_set%norm_drho(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
578 rho_set%owns%norm_drho = .true.
579 rho_set%has%norm_drho = .false.
580 END IF
581! Allocate the storage space for the norm of the gradient of rhoa and of rhob separatedly
582 IF (needs%norm_drho_spin) THEN
583 NULLIFY (rho_set%norm_drhoa, rho_set%norm_drhob)
584 ALLOCATE (rho_set%norm_drhoa(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
585 ALLOCATE (rho_set%norm_drhob(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
586 rho_set%owns%norm_drho_spin = .true.
587 rho_set%has%norm_drho_spin = .false.
588 END IF
589! Allocate the storage space for the components of the gradient for the total rho
590 IF (needs%drho) THEN
591 DO idir = 1, 3
592 NULLIFY (rho_set%drho(idir)%array)
593 ALLOCATE (rho_set%drho(idir)%array(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
594 END DO
595 rho_set%owns%drho = .true.
596 rho_set%has%drho = .false.
597 END IF
598! Allocate the storage space for the components of the gradient for rhoa and rhob
599 IF (needs%drho_spin) THEN
600 DO idir = 1, 3
601 NULLIFY (rho_set%drhoa(idir)%array, rho_set%drhob(idir)%array)
602 ALLOCATE (rho_set%drhoa(idir)%array(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
603 ALLOCATE (rho_set%drhob(idir)%array(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
604 END DO
605 rho_set%owns%drho_spin = .true.
606 rho_set%has%drho_spin = .false.
607 END IF
608!
609 END SELECT
610
611 ! tau part
612 IF (needs%tau) THEN
613 NULLIFY (rho_set%tau)
614 ALLOCATE (rho_set%tau(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
615 rho_set%owns%tau = .true.
616 END IF
617 IF (needs%tau_spin) THEN
618 NULLIFY (rho_set%tau_a, rho_set%tau_b)
619 ALLOCATE (rho_set%tau_a(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
620 ALLOCATE (rho_set%tau_b(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
621 rho_set%owns%tau_spin = .true.
622 rho_set%has%tau_spin = .false.
623 END IF
624
625 ! Laplace part
626 IF (needs%laplace_rho) THEN
627 NULLIFY (rho_set%laplace_rho)
628 ALLOCATE (rho_set%laplace_rho(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
629 rho_set%owns%laplace_rho = .true.
630 END IF
631 IF (needs%laplace_rho_spin) THEN
632 NULLIFY (rho_set%laplace_rhoa)
633 NULLIFY (rho_set%laplace_rhob)
634 ALLOCATE (rho_set%laplace_rhoa(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
635 ALLOCATE (rho_set%laplace_rhob(bo(1, 1):bo(2, 1), bo(1, 2):bo(2, 2), bo(1, 3):bo(2, 3)))
636 rho_set%owns%laplace_rho_spin = .true.
637 rho_set%has%laplace_rho_spin = .true.
638 END IF
639
640 END SUBROUTINE xc_rho_set_atom_update
641
642! **************************************************************************************************
643!> \brief ...
644!> \param rho_set ...
645!> \param lsd ...
646!> \param nspins ...
647!> \param needs ...
648!> \param rho ...
649!> \param drho ...
650!> \param tau ...
651!> \param na ...
652!> \param ir ...
653! **************************************************************************************************
654 SUBROUTINE fill_rho_set(rho_set, lsd, nspins, needs, rho, drho, tau, na, ir)
655
656 TYPE(xc_rho_set_type), INTENT(INOUT) :: rho_set
657 LOGICAL, INTENT(IN) :: lsd
658 INTEGER, INTENT(IN) :: nspins
659 TYPE(xc_rho_cflags_type), INTENT(IN) :: needs
660 REAL(dp), DIMENSION(:, :, :), INTENT(IN) :: rho
661 REAL(dp), DIMENSION(:, :, :, :), INTENT(IN) :: drho
662 REAL(dp), DIMENSION(:, :, :), INTENT(IN) :: tau
663 INTEGER, INTENT(IN) :: na, ir
664
665 REAL(kind=dp), PARAMETER :: f13 = (1.0_dp/3.0_dp)
666
667 INTEGER :: ia, idir, my_nspins
668 LOGICAL :: gradient_f, tddft_split
669
670 my_nspins = nspins
671 tddft_split = .false.
672 IF (lsd .AND. nspins == 1) THEN
673 my_nspins = 2
674 tddft_split = .true.
675 END IF
676
677 ! some checks
678 IF (lsd) THEN
679 ELSE
680 cpassert(SIZE(rho, 3) == 1)
681 END IF
682 SELECT CASE (my_nspins)
683 CASE (1)
684 cpassert(.NOT. needs%rho_spin)
685 cpassert(.NOT. needs%drho_spin)
686 cpassert(.NOT. needs%norm_drho_spin)
687 cpassert(.NOT. needs%rho_spin_1_3)
688 CASE (2)
689 CASE default
690 cpabort("Unsupported number of spins")
691 END SELECT
692
693 gradient_f = (needs%drho_spin .OR. needs%norm_drho_spin .OR. &
694 needs%drho .OR. needs%norm_drho)
695
696 SELECT CASE (my_nspins)
697 CASE (1)
698 ! Give rho to 1/3
699 IF (needs%rho_1_3) THEN
700 DO ia = 1, na
701 rho_set%rho_1_3(ia, ir, 1) = max(rho(ia, ir, 1), 0.0_dp)**f13
702 END DO
703 rho_set%owns%rho_1_3 = .true.
704 rho_set%has%rho_1_3 = .true.
705 END IF
706 ! Give the density
707 IF (needs%rho) THEN
708 DO ia = 1, na
709 rho_set%rho(ia, ir, 1) = rho(ia, ir, 1)
710 END DO
711 rho_set%owns%rho = .true.
712 rho_set%has%rho = .true.
713 END IF
714 ! Give the norm of the gradient of the density
715 IF (needs%norm_drho) THEN
716 DO ia = 1, na
717 rho_set%norm_drho(ia, ir, 1) = drho(4, ia, ir, 1)
718 END DO
719 rho_set%owns%norm_drho = .true.
720 rho_set%has%norm_drho = .true.
721 END IF
722 ! Give the three components of the gradient of the density
723 IF (needs%drho) THEN
724 DO idir = 1, 3
725 DO ia = 1, na
726 rho_set%drho(idir)%array(ia, ir, 1) = drho(idir, ia, ir, 1)
727 END DO
728 END DO
729 rho_set%owns%drho = .true.
730 rho_set%has%drho = .true.
731 END IF
732 CASE (2)
733 ! Give the total density
734 IF (needs%rho) THEN
735 ! this should never be the case unless you use LDA functionals with LSD
736 IF (.NOT. tddft_split) THEN
737 DO ia = 1, na
738 rho_set%rho(ia, ir, 1) = rho(ia, ir, 1) + rho(ia, ir, 2)
739 END DO
740 ELSE
741 DO ia = 1, na
742 rho_set%rho(ia, ir, 1) = rho(ia, ir, 1)
743 END DO
744 END IF
745 rho_set%owns%rho = .true.
746 rho_set%has%rho = .true.
747 END IF
748 ! Give the total density to 1/3
749 IF (needs%rho_1_3) THEN
750 IF (.NOT. tddft_split) THEN
751 DO ia = 1, na
752 rho_set%rho_1_3(ia, ir, 1) = max(rho(ia, ir, 1) + rho(ia, ir, 2), 0.0_dp)**f13
753 END DO
754 ELSE
755 DO ia = 1, na
756 rho_set%rho_1_3(ia, ir, 1) = max(rho(ia, ir, 1), 0.0_dp)**f13
757 END DO
758 END IF
759 rho_set%owns%rho_1_3 = .true.
760 rho_set%has%rho_1_3 = .true.
761 END IF
762 ! Give the spin densities to 1/3
763 IF (needs%rho_spin_1_3) THEN
764 IF (.NOT. tddft_split) THEN
765 DO ia = 1, na
766 rho_set%rhoa_1_3(ia, ir, 1) = max(rho(ia, ir, 1), 0.0_dp)**f13
767 rho_set%rhob_1_3(ia, ir, 1) = max(rho(ia, ir, 2), 0.0_dp)**f13
768 END DO
769 ELSE
770 DO ia = 1, na
771 rho_set%rhoa_1_3(ia, ir, 1) = max(0.5_dp*rho(ia, ir, 1), 0.0_dp)**f13
772 rho_set%rhob_1_3(ia, ir, 1) = rho_set%rhoa_1_3(ia, ir, 1)
773 END DO
774 END IF
775 rho_set%owns%rho_spin_1_3 = .true.
776 rho_set%has%rho_spin_1_3 = .true.
777 END IF
778 ! Give the spin densities rhoa and rhob
779 IF (needs%rho_spin) THEN
780 IF (.NOT. tddft_split) THEN
781 DO ia = 1, na
782 rho_set%rhoa(ia, ir, 1) = rho(ia, ir, 1)
783 rho_set%rhob(ia, ir, 1) = rho(ia, ir, 2)
784 END DO
785 ELSE
786 DO ia = 1, na
787 rho_set%rhoa(ia, ir, 1) = 0.5_dp*rho(ia, ir, 1)
788 rho_set%rhob(ia, ir, 1) = rho_set%rhoa(ia, ir, 1)
789 END DO
790 END IF
791 rho_set%owns%rho_spin = .true.
792 rho_set%has%rho_spin = .true.
793 END IF
794 ! Give the norm of the gradient of the total density
795 IF (needs%norm_drho) THEN
796 IF (.NOT. tddft_split) THEN
797 DO ia = 1, na
798 rho_set%norm_drho(ia, ir, 1) = sqrt( &
799 (drho(1, ia, ir, 1) + drho(1, ia, ir, 2))**2 + &
800 (drho(2, ia, ir, 1) + drho(2, ia, ir, 2))**2 + &
801 (drho(3, ia, ir, 1) + drho(3, ia, ir, 2))**2)
802 END DO
803 ELSE
804 DO ia = 1, na
805 rho_set%norm_drho(ia, ir, 1) = drho(4, ia, ir, 1)
806 END DO
807 END IF
808 rho_set%owns%norm_drho = .true.
809 rho_set%has%norm_drho = .true.
810 END IF
811 ! Give the norm of the gradient of rhoa and of rhob separatedly
812 IF (needs%norm_drho_spin) THEN
813 IF (.NOT. tddft_split) THEN
814 DO ia = 1, na
815 rho_set%norm_drhoa(ia, ir, 1) = drho(4, ia, ir, 1)
816 rho_set%norm_drhob(ia, ir, 1) = drho(4, ia, ir, 2)
817 END DO
818 ELSE
819 DO ia = 1, na
820 rho_set%norm_drhoa(ia, ir, 1) = 0.5_dp*drho(4, ia, ir, 1)
821 rho_set%norm_drhob(ia, ir, 1) = rho_set%norm_drhoa(ia, ir, 1)
822 END DO
823 END IF
824 rho_set%owns%norm_drho_spin = .true.
825 rho_set%has%norm_drho_spin = .true.
826 END IF
827 ! Give the components of the gradient for the total rho
828 IF (needs%drho) THEN
829 IF (.NOT. tddft_split) THEN
830 DO idir = 1, 3
831 DO ia = 1, na
832 rho_set%drho(idir)%array(ia, ir, 1) = drho(idir, ia, ir, 1) + drho(idir, ia, ir, 2)
833 END DO
834 END DO
835 ELSE
836 DO idir = 1, 3
837 DO ia = 1, na
838 rho_set%drho(idir)%array(ia, ir, 1) = drho(idir, ia, ir, 1)
839 END DO
840 END DO
841 END IF
842 rho_set%owns%drho = .true.
843 rho_set%has%drho = .true.
844 END IF
845 ! Give the components of the gradient for rhoa and rhob
846 IF (needs%drho_spin) THEN
847 IF (.NOT. tddft_split) THEN
848 DO idir = 1, 3
849 DO ia = 1, na
850 rho_set%drhoa(idir)%array(ia, ir, 1) = drho(idir, ia, ir, 1)
851 rho_set%drhob(idir)%array(ia, ir, 1) = drho(idir, ia, ir, 2)
852 END DO
853 END DO
854 ELSE
855 DO idir = 1, 3
856 DO ia = 1, na
857 rho_set%drhoa(idir)%array(ia, ir, 1) = 0.5_dp*drho(idir, ia, ir, 1)
858 rho_set%drhob(idir)%array(ia, ir, 1) = rho_set%drhoa(idir)%array(ia, ir, 1)
859 END DO
860 END DO
861 END IF
862 rho_set%owns%drho_spin = .true.
863 rho_set%has%drho_spin = .true.
864 END IF
865 !
866 END SELECT
867
868 ! tau part
869 IF (needs%tau .OR. needs%tau_spin) THEN
870 cpassert(SIZE(tau, 3) == my_nspins)
871 END IF
872 IF (needs%tau) THEN
873 IF (my_nspins == 2) THEN
874 DO ia = 1, na
875 rho_set%tau(ia, ir, 1) = tau(ia, ir, 1) + tau(ia, ir, 2)
876 END DO
877 rho_set%owns%tau = .true.
878 rho_set%has%tau = .true.
879 ELSE
880 DO ia = 1, na
881 rho_set%tau(ia, ir, 1) = tau(ia, ir, 1)
882 END DO
883 rho_set%owns%tau = .true.
884 rho_set%has%tau = .true.
885 END IF
886 END IF
887 IF (needs%tau_spin) THEN
888 DO ia = 1, na
889 rho_set%tau_a(ia, ir, 1) = tau(ia, ir, 1)
890 rho_set%tau_b(ia, ir, 1) = tau(ia, ir, 2)
891 END DO
892 rho_set%owns%tau_spin = .true.
893 rho_set%has%tau_spin = .true.
894 END IF
895
896 END SUBROUTINE fill_rho_set
897
898END MODULE xc_atom
cp_linked_list_xc_deriv
Definition cp_linked_list_xc_deriv.F:9
cp_linked_list_xc_deriv::cp_sll_xc_deriv_next
logical function, public cp_sll_xc_deriv_next(iterator, el_att)
returns true if the actual element is valid (i.e. iterator ont at end) moves the iterator to the next...
Definition cp_linked_list_xc_deriv.F:395
input_section_types
objects that represent the structure of input sections and the data contained in an input section
Definition input_section_types.F:15
input_section_types::section_vals_get_subs_vals
recursive type(section_vals_type) function, pointer, public section_vals_get_subs_vals(section_vals, subsection_name, i_rep_section, can_return_null)
returns the values of the requested subsection
Definition input_section_types.F:731
kinds
Defines the basic variable types.
Definition kinds.F:23
kinds::dp
integer, parameter, public dp
Definition kinds.F:34
pw_pool_types
Manages a pool of grids (to be used for example as tmp objects), but can also be used to instantiate ...
Definition pw_pool_types.F:24
pw_types
Definition pw_types.F:24
xc_atom
Definition xc_atom.F:9
xc_atom::vxc_of_r_epr
subroutine, public vxc_of_r_epr(xc_fun_section, rho_set, deriv_set, needs, w, lsd, na, nr, exc, vxc, vxg, vtau)
Specific EPR version of vxc_of_r_new.
Definition xc_atom.F:297
xc_atom::vxc_of_r_new
subroutine, public vxc_of_r_new(xc_fun_section, rho_set, deriv_set, deriv_order, needs, w, lsd, na, nr, exc, vxc, vxg, vtau, energy_only, adiabatic_rescale_factor)
...
Definition xc_atom.F:64
xc_atom::xc_rho_set_atom_update
subroutine, public xc_rho_set_atom_update(rho_set, needs, nspins, bo)
...
Definition xc_atom.F:495
xc_atom::fill_rho_set
subroutine, public fill_rho_set(rho_set, lsd, nspins, needs, rho, drho, tau, na, ir)
...
Definition xc_atom.F:655
xc_atom::xc_2nd_deriv_of_r
subroutine, public xc_2nd_deriv_of_r(rho_set, rho1_set, xc_section, deriv_set, w, vxc, vxg, do_triplet, do_sf)
...
Definition xc_atom.F:402
xc_derivative_desc
Module with functions to handle derivative descriptors. derivative description are strings have the f...
Definition xc_derivative_desc.F:26
xc_derivative_desc::deriv_norm_drho
integer, parameter, public deriv_norm_drho
Definition xc_derivative_desc.F:35
xc_derivative_desc::deriv_norm_drhoa
integer, parameter, public deriv_norm_drhoa
Definition xc_derivative_desc.F:35
xc_derivative_desc::deriv_rhob
integer, parameter, public deriv_rhob
Definition xc_derivative_desc.F:35
xc_derivative_desc::deriv_rhoa
integer, parameter, public deriv_rhoa
Definition xc_derivative_desc.F:35
xc_derivative_desc::deriv_tau
integer, parameter, public deriv_tau
Definition xc_derivative_desc.F:35
xc_derivative_desc::deriv_tau_b
integer, parameter, public deriv_tau_b
Definition xc_derivative_desc.F:35
xc_derivative_desc::deriv_tau_a
integer, parameter, public deriv_tau_a
Definition xc_derivative_desc.F:35
xc_derivative_desc::deriv_rho
integer, parameter, public deriv_rho
Definition xc_derivative_desc.F:35
xc_derivative_desc::deriv_norm_drhob
integer, parameter, public deriv_norm_drhob
Definition xc_derivative_desc.F:35
xc_derivative_set_types
represent a group ofunctional derivatives
Definition xc_derivative_set_types.F:14
xc_derivative_set_types::xc_dset_get_derivative
type(xc_derivative_type) function, pointer, public xc_dset_get_derivative(derivative_set, description, allocate_deriv)
returns the requested xc_derivative
Definition xc_derivative_set_types.F:65
xc_derivative_types
Provides types for the management of the xc-functionals and their derivatives.
Definition xc_derivative_types.F:12
xc_derivative_types::xc_derivative_get
subroutine, public xc_derivative_get(deriv, split_desc, order, deriv_data, accept_null_data)
returns various information on the given derivative
Definition xc_derivative_types.F:103
xc_derivatives
Definition xc_derivatives.F:9
xc_derivatives::xc_functionals_eval
subroutine, public xc_functionals_eval(functionals, lsd, rho_set, deriv_set, deriv_order)
...
Definition xc_derivatives.F:565
xc_rho_cflags_types
contains the structure
Definition xc_rho_cflags_types.F:14
xc_rho_set_types
contains the structure
Definition xc_rho_set_types.F:14
xc_rho_set_types::xc_rho_set_get
subroutine, public xc_rho_set_get(rho_set, can_return_null, rho, drho, norm_drho, rhoa, rhob, norm_drhoa, norm_drhob, rho_1_3, rhoa_1_3, rhob_1_3, laplace_rho, laplace_rhoa, laplace_rhob, drhoa, drhob, rho_cutoff, drho_cutoff, tau_cutoff, tau, tau_a, tau_b, local_bounds)
returns the various attributes of rho_set
Definition xc_rho_set_types.F:281
xc
Exchange and Correlation functional calculations.
Definition xc.F:17
xc::divide_by_norm_drho
subroutine, public divide_by_norm_drho(deriv_set, rho_set, lsd)
divides derivatives from deriv_set by norm_drho
Definition xc.F:5647
xc::xc_calc_2nd_deriv_analytical
subroutine, public xc_calc_2nd_deriv_analytical(v_xc, v_xc_tau, deriv_set, rho_set, rho1_set, pw_pool, xc_section, gapw, vxg, tddfpt_fac, compute_virial, virial_xc, spinflip)
Calculates the second derivative of E_xc at rho in the direction rho1 (if you see the second derivati...
Definition xc.F:2051
cp_linked_list_xc_deriv::cp_sll_xc_deriv_type
represent a single linked list that stores pointers to the elements
Definition cp_linked_list_xc_deriv.F:136
input_section_types::section_vals_type
stores the values of a section
Definition input_section_types.F:127
pw_pool_types::pw_pool_type
Manages a pool of grids (to be used for example as tmp objects), but can also be used to instantiate ...
Definition pw_pool_types.F:83
pw_types::pw_r3d_rs_type
Definition pw_types.F:62
xc_derivative_set_types::xc_derivative_set_type
A derivative set contains the different derivatives of a xc-functional in form of a linked list.
Definition xc_derivative_set_types.F:48
xc_derivative_types::xc_derivative_type
represent a derivative of a functional
Definition xc_derivative_types.F:32
xc_rho_cflags_types::xc_rho_cflags_type
contains a flag for each component of xc_rho_set, so that you can use it to tell which components you...
Definition xc_rho_cflags_types.F:48
xc_rho_set_types::xc_rho_set_type
represent a density, with all the representation and data needed to perform a functional evaluation
Definition xc_rho_set_types.F:78