(git:9111030)
Loading...
Searching...
No Matches
admm_types.F
Go to the documentation of this file.
1!--------------------------------------------------------------------------------------------------!
2! CP2K: A general program to perform molecular dynamics simulations !
3! Copyright 2000-2026 CP2K developers group <https://cp2k.org> !
4! !
5! SPDX-License-Identifier: GPL-2.0-or-later !
6!--------------------------------------------------------------------------------------------------!
7
8! **************************************************************************************************
9!> \brief Types and set/get functions for auxiliary density matrix methods
10!> \par History
11!> 05.2008 created [Manuel Guidon]
12!> 12.2019 Made GAPW compatiblae [Augustin Bussy]
13!> \author Manuel Guidon
14! **************************************************************************************************
15MODULE admm_types
16 USE admm_dm_types, ONLY: admm_dm_release,&
17 admm_dm_type
18 USE bibliography, ONLY: guidon2010,&
19 cite_reference
20 USE cp_blacs_env, ONLY: cp_blacs_env_type
21 USE cp_control_types, ONLY: admm_control_type
22 USE cp_dbcsr_api, ONLY: dbcsr_p_type
23 USE cp_dbcsr_operations, ONLY: dbcsr_deallocate_matrix_set
24 USE cp_fm_struct, ONLY: cp_fm_struct_create,&
25 cp_fm_struct_release,&
26 cp_fm_struct_type
27 USE cp_fm_types, ONLY: cp_fm_create,&
28 cp_fm_release,&
29 cp_fm_type
30 USE input_constants, ONLY: do_admm_aux_exch_func_none,&
31 do_admm_blocked_projection,&
32 do_admm_blocking_purify_full,&
33 do_admm_charge_constrained_projection,&
34 do_admm_exch_scaling_merlot,&
35 do_admm_exch_scaling_none,&
36 do_admm_purify_none
37 USE input_section_types, ONLY: section_vals_release,&
38 section_vals_type
39 USE kinds, ONLY: dp
40 USE kpoint_transitional, ONLY: get_1d_pointer,&
41 get_2d_pointer,&
42 kpoint_transitional_release,&
43 kpoint_transitional_type,&
44 set_1d_pointer,&
45 set_2d_pointer
46 USE message_passing, ONLY: mp_para_env_type
47 USE qs_kind_types, ONLY: deallocate_qs_kind_set,&
48 qs_kind_type
49 USE qs_local_rho_types, ONLY: local_rho_set_release,&
50 local_rho_type
51 USE qs_mo_types, ONLY: deallocate_mo_set,&
52 get_mo_set,&
53 mo_set_type
54 USE qs_neighbor_list_types, ONLY: neighbor_list_set_p_type,&
55 release_neighbor_list_sets
56 USE qs_oce_types, ONLY: deallocate_oce_set,&
57 oce_matrix_type
58 USE qs_rho_types, ONLY: qs_rho_release,&
59 qs_rho_type
60 USE task_list_types, ONLY: deallocate_task_list,&
61 task_list_type
62#include "./base/base_uses.f90"
63
64 IMPLICIT NONE
65 PRIVATE
66 PUBLIC :: admm_env_create, admm_env_release, admm_type, admm_gapw_r3d_rs_type, set_admm_env, get_admm_env
67
68 CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'admm_types'
69
70 TYPE eigvals_type
71 REAL(dp), DIMENSION(:), POINTER :: DATA => null()
72 END TYPE eigvals_type
73
74 TYPE eigvals_p_type
75 TYPE(eigvals_type), POINTER :: eigvals => null()
76 END TYPE eigvals_p_type
77
78! **************************************************************************************************
79!> \brief A subtype of the admm_env that contains the extra data needed for an ADMM GAPW calculation
80!> \param admm_kind_set gets its own qs_kind set to store all relevant basis/grid/etc info
81!> \param local_rho_set caontains soft and hard AUX_FIT atomoc densities
82!> \param task_list the task list used for all soft density pw operations
83!> \param oce stores the precomputed oce integrals
84! **************************************************************************************************
85 TYPE admm_gapw_r3d_rs_type
86 TYPE(qs_kind_type), DIMENSION(:), POINTER :: admm_kind_set => null()
87 TYPE(local_rho_type), POINTER :: local_rho_set => null()
88 TYPE(task_list_type), POINTER :: task_list => null()
89 TYPE(oce_matrix_type), POINTER :: oce => null()
90 END TYPE admm_gapw_r3d_rs_type
91
92! **************************************************************************************************
93!> \brief stores some data used in wavefunction fitting
94!> \param S overlap matrix for auxiliary fit basis set
95!> \param P overlap matrix for mixed aux_fit/orb basis set
96!> \param A contains inv(S)*P
97!> \param B contains transpose(P)*inv(S)*P = transpose(P)*A
98!> \param lambda contains transpose(mo_coeff_aux_fit)*B*mo_coeff_aux_fit
99!> \param lambda_inv_sqrt contains inv(SQRT(lambda))
100!> \param R contains eigenvectors of lambda
101!> \param work_aux_aux temporary matrix
102!> \param work_orb_nmo temporary matrix
103!> \param work_nmo_nmo1 temporary matrix
104!> \param work_nmo_nmo2 temporary matrix
105!> \param work_aux_nmo temporary matrix
106!> \param H contains KS_matrix * mo_coeff for auxiliary basis set
107!> \param K contains KS matrix for auxiliary basis set
108!> \param M contains matrix holding the 2nd order residues
109!> \param nao_orb number of atomic orbitals in orb basis set
110!> \param nao_aux_fit number of atomic orbitals in aux basis set
111!> \param nmo number of molecular orbitals per spin
112!> \param eigvals_lamda eigenvalues of lambda matrix
113!> \param gsi contains ratio N_dens_m/N_aux_dens_m
114!> \param admm_gapw_env the type containing ADMM GAPW specific data
115!> \param do_gapw an internal logical switch for GAPW
116!> \par History
117!> 05.2008 created [Manuel Guidon]
118!> \author Manuel Guidon
119! **************************************************************************************************
120 TYPE admm_type
121 TYPE(cp_fm_type), POINTER :: s_inv => null(), &
122 s => null(), &
123 q => null(), &
124 a => null(), &
125 b => null(), &
126 work_orb_orb => null(), &
127 work_orb_orb2 => null(), &
128 work_orb_orb3 => null(), &
129 work_aux_orb => null(), &
130 work_aux_orb2 => null(), &
131 work_aux_orb3 => null(), &
132 work_aux_aux => null(), &
133 work_aux_aux2 => null(), &
134 work_aux_aux3 => null(), &
135 work_aux_aux4 => null(), &
136 work_aux_aux5 => null()
137
138 TYPE(cp_fm_type), DIMENSION(:), &
139 POINTER :: lambda => null(), &
140 lambda_inv => null(), &
141 lambda_inv_sqrt => null(), &
142 r => null(), &
143 r_purify => null(), &
144 work_orb_nmo => null(), &
145 work_nmo_nmo1 => null(), &
146 r_schur_r_t => null(), &
147 work_nmo_nmo2 => null(), &
148 work_aux_nmo => null(), &
149 work_aux_nmo2 => null(), &
150 h => null(), &
151 h_corr => null(), &
152 mo_derivs_tmp => null(), &
153 k => null(), &
154 m => null(), &
155 m_purify => null(), &
156 p_to_be_purified => null(), &
157 lambda_inv2 => null(), &
158 c_hat => null(), &
159 p_tilde => null(), &
160 ks_to_be_merged => null(), &
161 scf_work_aux_fit => null()
162 TYPE(eigvals_p_type), DIMENSION(:), &
163 POINTER :: eigvals_lambda => null(), &
164 eigvals_p_to_be_purified => null()
165 TYPE(section_vals_type), POINTER :: xc_section_primary => null(), &
166 xc_section_aux => null()
167 REAL(kind=dp) :: gsi(3) = 0.0_dp, &
168 lambda_merlot(2) = 0.0_dp, &
169 n_large_basis(3) = 0.0_dp
170 INTEGER :: nao_orb = 0, nao_aux_fit = 0, nmo(2) = 0
171 INTEGER :: purification_method = do_admm_purify_none
172 LOGICAL :: charge_constrain = .false., do_admmp = .false., &
173 do_admmq = .false., do_admms = .false.
174 INTEGER :: scaling_model = do_admm_exch_scaling_none, &
175 aux_exch_func = do_admm_aux_exch_func_none
176 LOGICAL :: aux_exch_func_param = .false.
177 REAL(kind=dp), DIMENSION(3) :: aux_x_param = 0.0_dp
178 LOGICAL :: block_dm = .false.
179 LOGICAL :: block_fit = .false.
180 INTEGER, DIMENSION(:, :), POINTER :: block_map => null()
181 TYPE(admm_gapw_r3d_rs_type), POINTER :: admm_gapw_env => null()
182 LOGICAL :: do_gapw = .false.
183 TYPE(admm_dm_type), POINTER :: admm_dm => null()
184
185 TYPE(mo_set_type), DIMENSION(:), &
186 POINTER :: mos_aux_fit => null()
187 TYPE(neighbor_list_set_p_type), &
188 DIMENSION(:), POINTER :: sab_aux_fit => null(), sab_aux_fit_asymm => null(), sab_aux_fit_vs_orb => null()
189 TYPE(dbcsr_p_type), DIMENSION(:), &
190 POINTER :: matrix_ks_aux_fit_im => null()
191 TYPE(kpoint_transitional_type) :: matrix_ks_aux_fit, &
192 matrix_ks_aux_fit_dft, &
193 matrix_ks_aux_fit_hfx, &
194 matrix_s_aux_fit, &
195 matrix_s_aux_fit_vs_orb
196 TYPE(qs_rho_type), POINTER :: rho_aux_fit => null(), rho_aux_fit_buffer => null()
197 TYPE(task_list_type), POINTER :: task_list_aux_fit => null()
198 TYPE(cp_fm_type), DIMENSION(:), &
199 POINTER :: mo_derivs_aux_fit => null()
200
201 END TYPE admm_type
202
203CONTAINS
204
205! **************************************************************************************************
206!> \brief creates ADMM environment, initializes the basic types
207!>
208!> \param admm_env The ADMM env
209!> \param admm_control ...
210!> \param mos the MO's of the orbital basis set
211!> \param para_env The parallel env
212!> \param natoms ...
213!> \param nao_aux_fit ...
214!> \param blacs_env_ext ...
215!> \par History
216!> 05.2008 created [Manuel Guidon]
217!> \author Manuel Guidon
218! **************************************************************************************************
219 SUBROUTINE admm_env_create(admm_env, admm_control, mos, para_env, natoms, nao_aux_fit, blacs_env_ext)
220
221 TYPE(admm_type), POINTER :: admm_env
222 TYPE(admm_control_type), POINTER :: admm_control
223 TYPE(mo_set_type), DIMENSION(:), INTENT(IN) :: mos
224 TYPE(mp_para_env_type), POINTER :: para_env
225 INTEGER, INTENT(IN) :: natoms, nao_aux_fit
226 TYPE(cp_blacs_env_type), OPTIONAL, POINTER :: blacs_env_ext
227
228 INTEGER :: i, iatom, iblock, ispin, j, jatom, &
229 nao_orb, nmo, nspins
230 TYPE(cp_blacs_env_type), POINTER :: blacs_env
231 TYPE(cp_fm_struct_type), POINTER :: fm_struct_aux_aux, fm_struct_aux_nmo, fm_struct_aux_orb, &
232 fm_struct_nmo_nmo, fm_struct_orb_nmo, fm_struct_orb_orb
233 TYPE(cp_fm_type), POINTER :: mo_coeff
234
235 CALL cite_reference(guidon2010)
236
237 ALLOCATE (admm_env)
238
239 nspins = SIZE(mos)
240 CALL get_mo_set(mos(1), mo_coeff=mo_coeff, nmo=nmo, nao=nao_orb)
241 blacs_env => mo_coeff%matrix_struct%context
242 IF (PRESENT(blacs_env_ext)) blacs_env => blacs_env_ext
243
244 admm_env%nmo = 0
245 admm_env%nao_aux_fit = nao_aux_fit
246 admm_env%nao_orb = nao_orb
247 CALL cp_fm_struct_create(fm_struct_aux_aux, &
248 context=blacs_env, &
249 nrow_global=nao_aux_fit, &
250 ncol_global=nao_aux_fit, &
251 para_env=para_env)
252 CALL cp_fm_struct_create(fm_struct_aux_orb, &
253 context=blacs_env, &
254 nrow_global=nao_aux_fit, &
255 ncol_global=nao_orb, &
256 para_env=para_env)
257 CALL cp_fm_struct_create(fm_struct_orb_orb, &
258 context=blacs_env, &
259 nrow_global=nao_orb, &
260 ncol_global=nao_orb, &
261 para_env=para_env)
262
263 NULLIFY (admm_env%S, admm_env%S_inv, admm_env%Q, admm_env%A, admm_env%B, &
264 admm_env%work_orb_orb, admm_env%work_orb_orb2, admm_env%work_orb_orb3, &
265 admm_env%work_aux_orb, admm_env%work_aux_orb2, admm_env%work_aux_orb3, &
266 admm_env%work_aux_aux, admm_env%work_aux_aux2, admm_env%work_aux_aux3, &
267 admm_env%work_aux_aux4, admm_env%work_aux_aux5)
268 ALLOCATE (admm_env%S, admm_env%S_inv, admm_env%Q, admm_env%A, admm_env%B, &
269 admm_env%work_orb_orb, admm_env%work_orb_orb2, admm_env%work_orb_orb3, &
270 admm_env%work_aux_orb, admm_env%work_aux_orb2, admm_env%work_aux_orb3, &
271 admm_env%work_aux_aux, admm_env%work_aux_aux2, admm_env%work_aux_aux3, &
272 admm_env%work_aux_aux4, admm_env%work_aux_aux5)
273 CALL cp_fm_create(admm_env%S, fm_struct_aux_aux, name="aux_fit_overlap")
274 CALL cp_fm_create(admm_env%S_inv, fm_struct_aux_aux, name="aux_fit_overlap_inv")
275 CALL cp_fm_create(admm_env%Q, fm_struct_aux_orb, name="mixed_overlap")
276 CALL cp_fm_create(admm_env%A, fm_struct_aux_orb, name="work_A")
277 CALL cp_fm_create(admm_env%B, fm_struct_orb_orb, name="work_B")
278 CALL cp_fm_create(admm_env%work_orb_orb, fm_struct_orb_orb, name="work_orb_orb")
279 CALL cp_fm_create(admm_env%work_orb_orb2, fm_struct_orb_orb, name="work_orb_orb")
280 CALL cp_fm_create(admm_env%work_orb_orb3, fm_struct_orb_orb, name="work_orb_orb3")
281 CALL cp_fm_create(admm_env%work_aux_orb, fm_struct_aux_orb, name="work_aux_orb")
282 CALL cp_fm_create(admm_env%work_aux_orb2, fm_struct_aux_orb, name="work_aux_orb2")
283 CALL cp_fm_create(admm_env%work_aux_orb3, fm_struct_aux_orb, name="work_aux_orb3")
284 CALL cp_fm_create(admm_env%work_aux_aux, fm_struct_aux_aux, name="work_aux_aux")
285 CALL cp_fm_create(admm_env%work_aux_aux2, fm_struct_aux_aux, name="work_aux_aux2")
286 CALL cp_fm_create(admm_env%work_aux_aux3, fm_struct_aux_aux, name="work_aux_aux3")
287 CALL cp_fm_create(admm_env%work_aux_aux4, fm_struct_aux_aux, name="work_aux_aux4")
288 CALL cp_fm_create(admm_env%work_aux_aux5, fm_struct_aux_aux, name="work_aux_aux5")
289
290 ALLOCATE (admm_env%lambda_inv(nspins))
291 ALLOCATE (admm_env%lambda(nspins))
292 ALLOCATE (admm_env%lambda_inv_sqrt(nspins))
293 ALLOCATE (admm_env%R(nspins))
294 ALLOCATE (admm_env%R_purify(nspins))
295 ALLOCATE (admm_env%work_orb_nmo(nspins))
296 ALLOCATE (admm_env%work_nmo_nmo1(nspins))
297 ALLOCATE (admm_env%R_schur_R_t(nspins))
298 ALLOCATE (admm_env%work_nmo_nmo2(nspins))
299 ALLOCATE (admm_env%eigvals_lambda(nspins))
300 ALLOCATE (admm_env%eigvals_P_to_be_purified(nspins))
301 ALLOCATE (admm_env%H(nspins))
302 ALLOCATE (admm_env%K(nspins))
303 ALLOCATE (admm_env%M(nspins))
304 ALLOCATE (admm_env%M_purify(nspins))
305 ALLOCATE (admm_env%P_to_be_purified(nspins))
306 ALLOCATE (admm_env%work_aux_nmo(nspins))
307 ALLOCATE (admm_env%work_aux_nmo2(nspins))
308 ALLOCATE (admm_env%mo_derivs_tmp(nspins))
309 ALLOCATE (admm_env%H_corr(nspins))
310 ALLOCATE (admm_env%ks_to_be_merged(nspins))
311 ALLOCATE (admm_env%lambda_inv2(nspins))
312 ALLOCATE (admm_env%C_hat(nspins))
313 ALLOCATE (admm_env%P_tilde(nspins))
314
315 DO ispin = 1, nspins
316 CALL get_mo_set(mos(ispin), mo_coeff=mo_coeff, nmo=nmo)
317 admm_env%nmo(ispin) = nmo
318 CALL cp_fm_struct_create(fm_struct_aux_nmo, &
319 context=blacs_env, &
320 nrow_global=nao_aux_fit, &
321 ncol_global=nmo, &
322 para_env=para_env)
323 CALL cp_fm_struct_create(fm_struct_orb_nmo, &
324 context=blacs_env, &
325 nrow_global=nao_orb, &
326 ncol_global=nmo, &
327 para_env=para_env)
328 CALL cp_fm_struct_create(fm_struct_nmo_nmo, &
329 context=blacs_env, &
330 nrow_global=nmo, &
331 ncol_global=nmo, &
332 para_env=para_env)
333
334 CALL cp_fm_create(admm_env%work_orb_nmo(ispin), fm_struct_orb_nmo, name="work_orb_nmo")
335 CALL cp_fm_create(admm_env%work_nmo_nmo1(ispin), fm_struct_nmo_nmo, name="work_nmo_nmo1")
336 CALL cp_fm_create(admm_env%R_schur_R_t(ispin), fm_struct_nmo_nmo, name="R_schur_R_t")
337 CALL cp_fm_create(admm_env%work_nmo_nmo2(ispin), fm_struct_nmo_nmo, name="work_nmo_nmo2")
338 CALL cp_fm_create(admm_env%lambda(ispin), fm_struct_nmo_nmo, name="lambda")
339 CALL cp_fm_create(admm_env%lambda_inv(ispin), fm_struct_nmo_nmo, name="lambda_inv")
340 CALL cp_fm_create(admm_env%lambda_inv_sqrt(ispin), fm_struct_nmo_nmo, name="lambda_inv_sqrt")
341 CALL cp_fm_create(admm_env%R(ispin), fm_struct_nmo_nmo, name="R")
342 CALL cp_fm_create(admm_env%R_purify(ispin), fm_struct_aux_aux, name="R_purify")
343 CALL cp_fm_create(admm_env%K(ispin), fm_struct_aux_aux, name="K")
344 CALL cp_fm_create(admm_env%H(ispin), fm_struct_aux_nmo, name="H")
345 CALL cp_fm_create(admm_env%H_corr(ispin), fm_struct_orb_orb, name="H_corr")
346 CALL cp_fm_create(admm_env%M(ispin), fm_struct_nmo_nmo, name="M")
347 CALL cp_fm_create(admm_env%M_purify(ispin), fm_struct_aux_aux, name="M aux")
348 CALL cp_fm_create(admm_env%P_to_be_purified(ispin), fm_struct_aux_aux, name="P_to_be_purified")
349 CALL cp_fm_create(admm_env%work_aux_nmo(ispin), fm_struct_aux_nmo, name="work_aux_nmo")
350 CALL cp_fm_create(admm_env%work_aux_nmo2(ispin), fm_struct_aux_nmo, name="work_aux_nmo2")
351 CALL cp_fm_create(admm_env%mo_derivs_tmp(ispin), fm_struct_orb_nmo, name="mo_derivs_tmp")
352 CALL cp_fm_create(admm_env%lambda_inv2(ispin), fm_struct_nmo_nmo, name="lambda_inv2")
353 CALL cp_fm_create(admm_env%C_hat(ispin), fm_struct_aux_nmo, name="C_hat")
354 CALL cp_fm_create(admm_env%P_tilde(ispin), fm_struct_aux_aux, name="P_tilde")
355 CALL cp_fm_create(admm_env%ks_to_be_merged(ispin), fm_struct_orb_orb, name="KS_to_be_merged")
356
357 ALLOCATE (admm_env%eigvals_lambda(ispin)%eigvals)
358 ALLOCATE (admm_env%eigvals_P_to_be_purified(ispin)%eigvals)
359 ALLOCATE (admm_env%eigvals_lambda(ispin)%eigvals%data(nmo))
360 ALLOCATE (admm_env%eigvals_P_to_be_purified(ispin)%eigvals%data(nao_aux_fit))
361 admm_env%eigvals_lambda(ispin)%eigvals%data = 0.0_dp
362 admm_env%eigvals_P_to_be_purified(ispin)%eigvals%data = 0.0_dp
363 CALL cp_fm_struct_release(fm_struct_aux_nmo)
364 CALL cp_fm_struct_release(fm_struct_orb_nmo)
365 CALL cp_fm_struct_release(fm_struct_nmo_nmo)
366 END DO
367
368 CALL cp_fm_struct_release(fm_struct_aux_aux)
369 CALL cp_fm_struct_release(fm_struct_aux_orb)
370 CALL cp_fm_struct_release(fm_struct_orb_orb)
371
372 ! Copy settings from admm_control
373 cpassert(ASSOCIATED(admm_control))
374 admm_env%purification_method = admm_control%purification_method
375 admm_env%scaling_model = admm_control%scaling_model
376 admm_env%aux_exch_func = admm_control%aux_exch_func
377 admm_env%charge_constrain = (admm_control%method == do_admm_charge_constrained_projection)
378 admm_env%block_dm = ((admm_control%method == do_admm_blocking_purify_full) .OR. &
379 (admm_control%method == do_admm_blocked_projection))
380 admm_env%block_fit = admm_control%method == do_admm_blocked_projection
381 admm_env%aux_exch_func_param = admm_control%aux_exch_func_param
382 admm_env%aux_x_param(:) = admm_control%aux_x_param(:)
383
384 !ADMMP, ADMMQ, ADMMS
385 IF ((.NOT. admm_env%charge_constrain) .AND. (admm_env%scaling_model == do_admm_exch_scaling_merlot)) THEN
386 admm_env%do_admmp = .true.
387 END IF
388 IF (admm_env%charge_constrain .AND. (admm_env%scaling_model == do_admm_exch_scaling_none)) THEN
389 admm_env%do_admmq = .true.
390 END IF
391 IF (admm_env%charge_constrain .AND. (admm_env%scaling_model == do_admm_exch_scaling_merlot)) THEN
392 admm_env%do_admms = .true.
393 END IF
394
395 IF ((admm_control%method == do_admm_blocking_purify_full) .OR. &
396 (admm_control%method == do_admm_blocked_projection)) THEN
397 ! Create block map
398 ALLOCATE (admm_env%block_map(natoms, natoms))
399 admm_env%block_map(:, :) = 0
400 DO iblock = 1, SIZE(admm_control%blocks)
401 DO i = 1, SIZE(admm_control%blocks(iblock)%list)
402 iatom = admm_control%blocks(iblock)%list(i)
403 DO j = 1, SIZE(admm_control%blocks(iblock)%list)
404 jatom = admm_control%blocks(iblock)%list(j)
405 admm_env%block_map(iatom, jatom) = 1
406 END DO
407 END DO
408 END DO
409 END IF
410
411 NULLIFY (admm_env%admm_gapw_env)
412 admm_env%do_gapw = .false.
413
414 NULLIFY (admm_env%mos_aux_fit, admm_env%sab_aux_fit, admm_env%sab_aux_fit_asymm, admm_env%sab_aux_fit_vs_orb)
415 NULLIFY (admm_env%matrix_ks_aux_fit_im)
416 NULLIFY (admm_env%rho_aux_fit, admm_env%rho_aux_fit_buffer, admm_env%task_list_aux_fit, admm_env%mo_derivs_aux_fit)
417
418 END SUBROUTINE admm_env_create
419
420! **************************************************************************************************
421!> \brief releases the ADMM environment, cleans up all types
422!>
423!> \param admm_env The ADMM env
424!> \par History
425!> 05.2008 created [Manuel Guidon]
426!> \author Manuel Guidon
427! **************************************************************************************************
428 SUBROUTINE admm_env_release(admm_env)
429
430 TYPE(admm_type), POINTER :: admm_env
431
432 INTEGER :: ispin
433
434 cpassert(ASSOCIATED(admm_env))
435
436 CALL cp_fm_release(admm_env%S)
437 CALL cp_fm_release(admm_env%S_inv)
438 CALL cp_fm_release(admm_env%Q)
439 CALL cp_fm_release(admm_env%A)
440 CALL cp_fm_release(admm_env%B)
441 CALL cp_fm_release(admm_env%work_orb_orb)
442 CALL cp_fm_release(admm_env%work_orb_orb2)
443 CALL cp_fm_release(admm_env%work_orb_orb3)
444 CALL cp_fm_release(admm_env%work_aux_aux)
445 CALL cp_fm_release(admm_env%work_aux_aux2)
446 CALL cp_fm_release(admm_env%work_aux_aux3)
447 CALL cp_fm_release(admm_env%work_aux_aux4)
448 CALL cp_fm_release(admm_env%work_aux_aux5)
449 CALL cp_fm_release(admm_env%work_aux_orb)
450 CALL cp_fm_release(admm_env%work_aux_orb2)
451 CALL cp_fm_release(admm_env%work_aux_orb3)
452 DEALLOCATE (admm_env%S, admm_env%S_inv, admm_env%Q, admm_env%A, admm_env%B, &
453 admm_env%work_orb_orb, admm_env%work_orb_orb2, admm_env%work_orb_orb3, &
454 admm_env%work_aux_orb, admm_env%work_aux_orb2, admm_env%work_aux_orb3, &
455 admm_env%work_aux_aux, admm_env%work_aux_aux2, admm_env%work_aux_aux3, &
456 admm_env%work_aux_aux4, admm_env%work_aux_aux5)
457 CALL cp_fm_release(admm_env%lambda)
458 CALL cp_fm_release(admm_env%lambda_inv)
459 CALL cp_fm_release(admm_env%lambda_inv_sqrt)
460 CALL cp_fm_release(admm_env%lambda_inv2)
461 CALL cp_fm_release(admm_env%C_hat)
462 CALL cp_fm_release(admm_env%P_tilde)
463 CALL cp_fm_release(admm_env%R)
464 CALL cp_fm_release(admm_env%R_purify)
465 CALL cp_fm_release(admm_env%H)
466 CALL cp_fm_release(admm_env%H_corr)
467 CALL cp_fm_release(admm_env%K)
468 CALL cp_fm_release(admm_env%M)
469 CALL cp_fm_release(admm_env%M_purify)
470 CALL cp_fm_release(admm_env%P_to_be_purified)
471 CALL cp_fm_release(admm_env%work_orb_nmo)
472 CALL cp_fm_release(admm_env%work_nmo_nmo1)
473 CALL cp_fm_release(admm_env%R_schur_R_t)
474 CALL cp_fm_release(admm_env%work_nmo_nmo2)
475 CALL cp_fm_release(admm_env%work_aux_nmo)
476 CALL cp_fm_release(admm_env%work_aux_nmo2)
477 CALL cp_fm_release(admm_env%mo_derivs_tmp)
478 CALL cp_fm_release(admm_env%ks_to_be_merged)
479 CALL cp_fm_release(admm_env%lambda_inv2)
480 DO ispin = 1, SIZE(admm_env%eigvals_lambda)
481 DEALLOCATE (admm_env%eigvals_lambda(ispin)%eigvals%data)
482 DEALLOCATE (admm_env%eigvals_P_to_be_purified(ispin)%eigvals%data)
483 DEALLOCATE (admm_env%eigvals_lambda(ispin)%eigvals)
484 DEALLOCATE (admm_env%eigvals_P_to_be_purified(ispin)%eigvals)
485 END DO
486 DEALLOCATE (admm_env%eigvals_lambda)
487 DEALLOCATE (admm_env%eigvals_P_to_be_purified)
488
489 IF (ASSOCIATED(admm_env%block_map)) THEN
490 DEALLOCATE (admm_env%block_map)
491 END IF
492
493 IF (ASSOCIATED(admm_env%xc_section_primary)) THEN
494 CALL section_vals_release(admm_env%xc_section_primary)
495 END IF
496 IF (ASSOCIATED(admm_env%xc_section_aux)) THEN
497 CALL section_vals_release(admm_env%xc_section_aux)
498 END IF
499
500 IF (ASSOCIATED(admm_env%admm_gapw_env)) CALL admm_gapw_env_release(admm_env%admm_gapw_env)
501 IF (ASSOCIATED(admm_env%admm_dm)) CALL admm_dm_release(admm_env%admm_dm)
502
503 IF (ASSOCIATED(admm_env%mos_aux_fit)) THEN
504 DO ispin = 1, SIZE(admm_env%mos_aux_fit)
505 CALL deallocate_mo_set(admm_env%mos_aux_fit(ispin))
506 END DO
507 DEALLOCATE (admm_env%mos_aux_fit)
508 END IF
509 CALL cp_fm_release(admm_env%mo_derivs_aux_fit)
510
511 IF (ASSOCIATED(admm_env%scf_work_aux_fit)) THEN
512 CALL cp_fm_release(admm_env%scf_work_aux_fit)
513 END IF
514
515 IF (ASSOCIATED(admm_env%sab_aux_fit)) CALL release_neighbor_list_sets(admm_env%sab_aux_fit)
516 IF (ASSOCIATED(admm_env%sab_aux_fit_vs_orb)) CALL release_neighbor_list_sets(admm_env%sab_aux_fit_vs_orb)
517 IF (ASSOCIATED(admm_env%sab_aux_fit_asymm)) CALL release_neighbor_list_sets(admm_env%sab_aux_fit_asymm)
518
519 CALL kpoint_transitional_release(admm_env%matrix_ks_aux_fit)
520 CALL kpoint_transitional_release(admm_env%matrix_ks_aux_fit_dft)
521 CALL kpoint_transitional_release(admm_env%matrix_ks_aux_fit_hfx)
522 CALL kpoint_transitional_release(admm_env%matrix_s_aux_fit)
523 CALL kpoint_transitional_release(admm_env%matrix_s_aux_fit_vs_orb)
524
525 IF (ASSOCIATED(admm_env%matrix_ks_aux_fit_im)) CALL dbcsr_deallocate_matrix_set(admm_env%matrix_ks_aux_fit_im)
526
527 IF (ASSOCIATED(admm_env%rho_aux_fit)) THEN
528 CALL qs_rho_release(admm_env%rho_aux_fit)
529 DEALLOCATE (admm_env%rho_aux_fit)
530 END IF
531 IF (ASSOCIATED(admm_env%rho_aux_fit_buffer)) THEN
532 CALL qs_rho_release(admm_env%rho_aux_fit_buffer)
533 DEALLOCATE (admm_env%rho_aux_fit_buffer)
534 END IF
535
536 IF (ASSOCIATED(admm_env%task_list_aux_fit)) CALL deallocate_task_list(admm_env%task_list_aux_fit)
537
538 DEALLOCATE (admm_env)
539
540 END SUBROUTINE admm_env_release
541
542! **************************************************************************************************
543!> \brief Release the ADMM GAPW stuff
544!> \param admm_gapw_env ...
545! **************************************************************************************************
546 SUBROUTINE admm_gapw_env_release(admm_gapw_env)
547
548 TYPE(admm_gapw_r3d_rs_type), POINTER :: admm_gapw_env
549
550 IF (ASSOCIATED(admm_gapw_env%admm_kind_set)) THEN
551 CALL deallocate_qs_kind_set(admm_gapw_env%admm_kind_set)
552 END IF
553
554 IF (ASSOCIATED(admm_gapw_env%local_rho_set)) THEN
555 CALL local_rho_set_release(admm_gapw_env%local_rho_set)
556 END IF
557
558 IF (ASSOCIATED(admm_gapw_env%task_list)) THEN
559 CALL deallocate_task_list(admm_gapw_env%task_list)
560 END IF
561
562 IF (ASSOCIATED(admm_gapw_env%oce)) THEN
563 CALL deallocate_oce_set(admm_gapw_env%oce)
564 END IF
565
566 DEALLOCATE (admm_gapw_env)
567
568 END SUBROUTINE admm_gapw_env_release
569
570! **************************************************************************************************
571!> \brief Get routine for the ADMM env
572!> \param admm_env ...
573!> \param mo_derivs_aux_fit ...
574!> \param mos_aux_fit ...
575!> \param sab_aux_fit ...
576!> \param sab_aux_fit_asymm ...
577!> \param sab_aux_fit_vs_orb ...
578!> \param matrix_s_aux_fit ...
579!> \param matrix_s_aux_fit_kp ...
580!> \param matrix_s_aux_fit_vs_orb ...
581!> \param matrix_s_aux_fit_vs_orb_kp ...
582!> \param task_list_aux_fit ...
583!> \param matrix_ks_aux_fit ...
584!> \param matrix_ks_aux_fit_kp ...
585!> \param matrix_ks_aux_fit_im ...
586!> \param matrix_ks_aux_fit_dft ...
587!> \param matrix_ks_aux_fit_hfx ...
588!> \param matrix_ks_aux_fit_dft_kp ...
589!> \param matrix_ks_aux_fit_hfx_kp ...
590!> \param rho_aux_fit ...
591!> \param rho_aux_fit_buffer ...
592!> \param admm_dm ...
593! **************************************************************************************************
594 SUBROUTINE get_admm_env(admm_env, mo_derivs_aux_fit, mos_aux_fit, sab_aux_fit, sab_aux_fit_asymm, &
595 sab_aux_fit_vs_orb, matrix_s_aux_fit, matrix_s_aux_fit_kp, matrix_s_aux_fit_vs_orb, matrix_s_aux_fit_vs_orb_kp, &
596 task_list_aux_fit, matrix_ks_aux_fit, matrix_ks_aux_fit_kp, matrix_ks_aux_fit_im, &
597 matrix_ks_aux_fit_dft, matrix_ks_aux_fit_hfx, matrix_ks_aux_fit_dft_kp, matrix_ks_aux_fit_hfx_kp, rho_aux_fit, &
598 rho_aux_fit_buffer, admm_dm)
599
600 TYPE(admm_type), INTENT(IN), POINTER :: admm_env
601 TYPE(cp_fm_type), DIMENSION(:), OPTIONAL, POINTER :: mo_derivs_aux_fit
602 TYPE(mo_set_type), DIMENSION(:), OPTIONAL, POINTER :: mos_aux_fit
603 TYPE(neighbor_list_set_p_type), DIMENSION(:), &
604 OPTIONAL, POINTER :: sab_aux_fit, sab_aux_fit_asymm, &
605 sab_aux_fit_vs_orb
606 TYPE(dbcsr_p_type), DIMENSION(:), OPTIONAL, &
607 POINTER :: matrix_s_aux_fit
608 TYPE(dbcsr_p_type), DIMENSION(:, :), OPTIONAL, &
609 POINTER :: matrix_s_aux_fit_kp
610 TYPE(dbcsr_p_type), DIMENSION(:), OPTIONAL, &
611 POINTER :: matrix_s_aux_fit_vs_orb
612 TYPE(dbcsr_p_type), DIMENSION(:, :), OPTIONAL, &
613 POINTER :: matrix_s_aux_fit_vs_orb_kp
614 TYPE(task_list_type), OPTIONAL, POINTER :: task_list_aux_fit
615 TYPE(dbcsr_p_type), DIMENSION(:), OPTIONAL, &
616 POINTER :: matrix_ks_aux_fit
617 TYPE(dbcsr_p_type), DIMENSION(:, :), OPTIONAL, &
618 POINTER :: matrix_ks_aux_fit_kp
619 TYPE(dbcsr_p_type), DIMENSION(:), OPTIONAL, &
620 POINTER :: matrix_ks_aux_fit_im, &
621 matrix_ks_aux_fit_dft, &
622 matrix_ks_aux_fit_hfx
623 TYPE(dbcsr_p_type), DIMENSION(:, :), OPTIONAL, &
624 POINTER :: matrix_ks_aux_fit_dft_kp, &
625 matrix_ks_aux_fit_hfx_kp
626 TYPE(qs_rho_type), OPTIONAL, POINTER :: rho_aux_fit, rho_aux_fit_buffer
627 TYPE(admm_dm_type), OPTIONAL, POINTER :: admm_dm
628
629 cpassert(ASSOCIATED(admm_env))
630
631 IF (PRESENT(mo_derivs_aux_fit)) mo_derivs_aux_fit => admm_env%mo_derivs_aux_fit
632 IF (PRESENT(mos_aux_fit)) mos_aux_fit => admm_env%mos_aux_fit
633 IF (PRESENT(sab_aux_fit)) sab_aux_fit => admm_env%sab_aux_fit
634 IF (PRESENT(sab_aux_fit_asymm)) sab_aux_fit_asymm => admm_env%sab_aux_fit_asymm
635 IF (PRESENT(sab_aux_fit_vs_orb)) sab_aux_fit_vs_orb => admm_env%sab_aux_fit_vs_orb
636 IF (PRESENT(task_list_aux_fit)) task_list_aux_fit => admm_env%task_list_aux_fit
637 IF (PRESENT(matrix_ks_aux_fit_im)) matrix_ks_aux_fit_im => admm_env%matrix_ks_aux_fit_im
638 IF (PRESENT(rho_aux_fit)) rho_aux_fit => admm_env%rho_aux_fit
639 IF (PRESENT(rho_aux_fit_buffer)) rho_aux_fit_buffer => admm_env%rho_aux_fit_buffer
640 IF (PRESENT(admm_dm)) admm_dm => admm_env%admm_dm
641
642 IF (PRESENT(matrix_ks_aux_fit)) matrix_ks_aux_fit => get_1d_pointer(admm_env%matrix_ks_aux_fit)
643 IF (PRESENT(matrix_ks_aux_fit_kp)) matrix_ks_aux_fit_kp => get_2d_pointer(admm_env%matrix_ks_aux_fit)
644 IF (PRESENT(matrix_ks_aux_fit_dft)) matrix_ks_aux_fit_dft => get_1d_pointer(admm_env%matrix_ks_aux_fit_dft)
645 IF (PRESENT(matrix_ks_aux_fit_dft_kp)) matrix_ks_aux_fit_dft_kp => get_2d_pointer(admm_env%matrix_ks_aux_fit_dft)
646 IF (PRESENT(matrix_ks_aux_fit_hfx)) matrix_ks_aux_fit_hfx => get_1d_pointer(admm_env%matrix_ks_aux_fit_hfx)
647 IF (PRESENT(matrix_ks_aux_fit_hfx_kp)) matrix_ks_aux_fit_hfx_kp => get_2d_pointer(admm_env%matrix_ks_aux_fit_hfx)
648 IF (PRESENT(matrix_s_aux_fit)) matrix_s_aux_fit => get_1d_pointer(admm_env%matrix_s_aux_fit)
649 IF (PRESENT(matrix_s_aux_fit_kp)) matrix_s_aux_fit_kp => get_2d_pointer(admm_env%matrix_s_aux_fit)
650 IF (PRESENT(matrix_s_aux_fit_vs_orb)) matrix_s_aux_fit_vs_orb => get_1d_pointer(admm_env%matrix_s_aux_fit_vs_orb)
651 IF (PRESENT(matrix_s_aux_fit_vs_orb_kp)) matrix_s_aux_fit_vs_orb_kp => get_2d_pointer(admm_env%matrix_s_aux_fit_vs_orb)
652 END SUBROUTINE get_admm_env
653
654! **************************************************************************************************
655!> \brief Set routine for the ADMM env
656!> \param admm_env ...
657!> \param mo_derivs_aux_fit ...
658!> \param mos_aux_fit ...
659!> \param sab_aux_fit ...
660!> \param sab_aux_fit_asymm ...
661!> \param sab_aux_fit_vs_orb ...
662!> \param matrix_s_aux_fit ...
663!> \param matrix_s_aux_fit_kp ...
664!> \param matrix_s_aux_fit_vs_orb ...
665!> \param matrix_s_aux_fit_vs_orb_kp ...
666!> \param task_list_aux_fit ...
667!> \param matrix_ks_aux_fit ...
668!> \param matrix_ks_aux_fit_kp ...
669!> \param matrix_ks_aux_fit_im ...
670!> \param matrix_ks_aux_fit_dft ...
671!> \param matrix_ks_aux_fit_hfx ...
672!> \param matrix_ks_aux_fit_dft_kp ...
673!> \param matrix_ks_aux_fit_hfx_kp ...
674!> \param rho_aux_fit ...
675!> \param rho_aux_fit_buffer ...
676!> \param admm_dm ...
677! **************************************************************************************************
678 SUBROUTINE set_admm_env(admm_env, mo_derivs_aux_fit, mos_aux_fit, sab_aux_fit, sab_aux_fit_asymm, &
679 sab_aux_fit_vs_orb, matrix_s_aux_fit, matrix_s_aux_fit_kp, matrix_s_aux_fit_vs_orb, matrix_s_aux_fit_vs_orb_kp, &
680 task_list_aux_fit, matrix_ks_aux_fit, matrix_ks_aux_fit_kp, matrix_ks_aux_fit_im, &
681 matrix_ks_aux_fit_dft, matrix_ks_aux_fit_hfx, matrix_ks_aux_fit_dft_kp, matrix_ks_aux_fit_hfx_kp, rho_aux_fit, &
682 rho_aux_fit_buffer, admm_dm)
683
684 TYPE(admm_type), INTENT(INOUT), POINTER :: admm_env
685 TYPE(cp_fm_type), DIMENSION(:), OPTIONAL, POINTER :: mo_derivs_aux_fit
686 TYPE(mo_set_type), DIMENSION(:), OPTIONAL, POINTER :: mos_aux_fit
687 TYPE(neighbor_list_set_p_type), DIMENSION(:), &
688 OPTIONAL, POINTER :: sab_aux_fit, sab_aux_fit_asymm, &
689 sab_aux_fit_vs_orb
690 TYPE(dbcsr_p_type), DIMENSION(:), OPTIONAL, &
691 POINTER :: matrix_s_aux_fit
692 TYPE(dbcsr_p_type), DIMENSION(:, :), OPTIONAL, &
693 POINTER :: matrix_s_aux_fit_kp
694 TYPE(dbcsr_p_type), DIMENSION(:), OPTIONAL, &
695 POINTER :: matrix_s_aux_fit_vs_orb
696 TYPE(dbcsr_p_type), DIMENSION(:, :), OPTIONAL, &
697 POINTER :: matrix_s_aux_fit_vs_orb_kp
698 TYPE(task_list_type), OPTIONAL, POINTER :: task_list_aux_fit
699 TYPE(dbcsr_p_type), DIMENSION(:), OPTIONAL, &
700 POINTER :: matrix_ks_aux_fit
701 TYPE(dbcsr_p_type), DIMENSION(:, :), OPTIONAL, &
702 POINTER :: matrix_ks_aux_fit_kp
703 TYPE(dbcsr_p_type), DIMENSION(:), OPTIONAL, &
704 POINTER :: matrix_ks_aux_fit_im, &
705 matrix_ks_aux_fit_dft, &
706 matrix_ks_aux_fit_hfx
707 TYPE(dbcsr_p_type), DIMENSION(:, :), OPTIONAL, &
708 POINTER :: matrix_ks_aux_fit_dft_kp, &
709 matrix_ks_aux_fit_hfx_kp
710 TYPE(qs_rho_type), OPTIONAL, POINTER :: rho_aux_fit, rho_aux_fit_buffer
711 TYPE(admm_dm_type), OPTIONAL, POINTER :: admm_dm
712
713 cpassert(ASSOCIATED(admm_env))
714
715 IF (PRESENT(mo_derivs_aux_fit)) admm_env%mo_derivs_aux_fit => mo_derivs_aux_fit
716 IF (PRESENT(mos_aux_fit)) admm_env%mos_aux_fit => mos_aux_fit
717 IF (PRESENT(sab_aux_fit)) admm_env%sab_aux_fit => sab_aux_fit
718 IF (PRESENT(sab_aux_fit_asymm)) admm_env%sab_aux_fit_asymm => sab_aux_fit_asymm
719 IF (PRESENT(sab_aux_fit_vs_orb)) admm_env%sab_aux_fit_vs_orb => sab_aux_fit_vs_orb
720 IF (PRESENT(task_list_aux_fit)) admm_env%task_list_aux_fit => task_list_aux_fit
721 IF (PRESENT(matrix_ks_aux_fit_im)) admm_env%matrix_ks_aux_fit_im => matrix_ks_aux_fit_im
722 IF (PRESENT(rho_aux_fit)) admm_env%rho_aux_fit => rho_aux_fit
723 IF (PRESENT(rho_aux_fit_buffer)) admm_env%rho_aux_fit_buffer => rho_aux_fit_buffer
724 IF (PRESENT(admm_dm)) admm_env%admm_dm => admm_dm
725
726 IF (PRESENT(matrix_ks_aux_fit)) CALL set_1d_pointer(admm_env%matrix_ks_aux_fit, matrix_ks_aux_fit)
727 IF (PRESENT(matrix_ks_aux_fit_kp)) CALL set_2d_pointer(admm_env%matrix_ks_aux_fit, matrix_ks_aux_fit_kp)
728 IF (PRESENT(matrix_ks_aux_fit_dft)) CALL set_1d_pointer(admm_env%matrix_ks_aux_fit_dft, matrix_ks_aux_fit_dft)
729 IF (PRESENT(matrix_ks_aux_fit_dft_kp)) CALL set_2d_pointer(admm_env%matrix_ks_aux_fit_dft, matrix_ks_aux_fit_dft_kp)
730 IF (PRESENT(matrix_ks_aux_fit_hfx)) CALL set_1d_pointer(admm_env%matrix_ks_aux_fit_hfx, matrix_ks_aux_fit)
731 IF (PRESENT(matrix_ks_aux_fit_hfx_kp)) CALL set_2d_pointer(admm_env%matrix_ks_aux_fit_hfx, matrix_ks_aux_fit_hfx_kp)
732 IF (PRESENT(matrix_s_aux_fit)) CALL set_1d_pointer(admm_env%matrix_s_aux_fit, matrix_s_aux_fit)
733 IF (PRESENT(matrix_s_aux_fit_kp)) CALL set_2d_pointer(admm_env%matrix_s_aux_fit, matrix_s_aux_fit_kp)
734 IF (PRESENT(matrix_s_aux_fit_vs_orb)) CALL set_1d_pointer(admm_env%matrix_s_aux_fit_vs_orb, matrix_s_aux_fit_vs_orb)
735 IF (PRESENT(matrix_s_aux_fit_vs_orb_kp)) CALL set_2d_pointer(admm_env%matrix_s_aux_fit_vs_orb, matrix_s_aux_fit_vs_orb_kp)
736
737 END SUBROUTINE set_admm_env
738
739END MODULE admm_types
740
cp_dbcsr_operations::dbcsr_deallocate_matrix_set
Definition cp_dbcsr_operations.F:85
cp_fm_types::cp_fm_release
Definition cp_fm_types.F:88
admm_dm_types
Types and set/get functions for auxiliary density matrix methods.
Definition admm_dm_types.F:14
admm_dm_types::admm_dm_release
subroutine, public admm_dm_release(admm_dm)
Release a admm_dm type.
Definition admm_dm_types.F:99
admm_types
Types and set/get functions for auxiliary density matrix methods.
Definition admm_types.F:15
admm_types::get_admm_env
subroutine, public get_admm_env(admm_env, mo_derivs_aux_fit, mos_aux_fit, sab_aux_fit, sab_aux_fit_asymm, sab_aux_fit_vs_orb, matrix_s_aux_fit, matrix_s_aux_fit_kp, matrix_s_aux_fit_vs_orb, matrix_s_aux_fit_vs_orb_kp, task_list_aux_fit, matrix_ks_aux_fit, matrix_ks_aux_fit_kp, matrix_ks_aux_fit_im, matrix_ks_aux_fit_dft, matrix_ks_aux_fit_hfx, matrix_ks_aux_fit_dft_kp, matrix_ks_aux_fit_hfx_kp, rho_aux_fit, rho_aux_fit_buffer, admm_dm)
Get routine for the ADMM env.
Definition admm_types.F:599
admm_types::set_admm_env
subroutine, public set_admm_env(admm_env, mo_derivs_aux_fit, mos_aux_fit, sab_aux_fit, sab_aux_fit_asymm, sab_aux_fit_vs_orb, matrix_s_aux_fit, matrix_s_aux_fit_kp, matrix_s_aux_fit_vs_orb, matrix_s_aux_fit_vs_orb_kp, task_list_aux_fit, matrix_ks_aux_fit, matrix_ks_aux_fit_kp, matrix_ks_aux_fit_im, matrix_ks_aux_fit_dft, matrix_ks_aux_fit_hfx, matrix_ks_aux_fit_dft_kp, matrix_ks_aux_fit_hfx_kp, rho_aux_fit, rho_aux_fit_buffer, admm_dm)
Set routine for the ADMM env.
Definition admm_types.F:683
admm_types::admm_env_release
subroutine, public admm_env_release(admm_env)
releases the ADMM environment, cleans up all types
Definition admm_types.F:429
admm_types::admm_env_create
subroutine, public admm_env_create(admm_env, admm_control, mos, para_env, natoms, nao_aux_fit, blacs_env_ext)
creates ADMM environment, initializes the basic types
Definition admm_types.F:220
bibliography
collects all references to literature in CP2K as new algorithms / method are included from literature...
Definition bibliography.F:21
bibliography::guidon2010
integer, save, public guidon2010
Definition bibliography.F:36
cp_blacs_env
methods related to the blacs parallel environment
Definition cp_blacs_env.F:15
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_dbcsr_api
Definition cp_dbcsr_api.F:8
cp_dbcsr_operations
DBCSR operations in CP2K.
Definition cp_dbcsr_operations.F:18
cp_fm_struct
represent the structure of a full matrix
Definition cp_fm_struct.F:14
cp_fm_struct::cp_fm_struct_create
subroutine, public cp_fm_struct_create(fmstruct, para_env, context, nrow_global, ncol_global, nrow_block, ncol_block, descriptor, first_p_pos, local_leading_dimension, template_fmstruct, square_blocks, force_block)
allocates and initializes a full matrix structure
Definition cp_fm_struct.F:132
cp_fm_struct::cp_fm_struct_release
subroutine, public cp_fm_struct_release(fmstruct)
releases a full matrix structure
Definition cp_fm_struct.F:355
cp_fm_types
represent a full matrix distributed on many processors
Definition cp_fm_types.F:15
cp_fm_types::cp_fm_create
subroutine, public cp_fm_create(matrix, matrix_struct, name, nrow, ncol, set_zero)
creates a new full matrix with the given structure
Definition cp_fm_types.F:165
input_constants
collects all constants needed in input so that they can be used without circular dependencies
Definition input_constants.F:18
input_constants::do_admm_purify_none
integer, parameter, public do_admm_purify_none
Definition input_constants.F:885
input_constants::do_admm_blocking_purify_full
integer, parameter, public do_admm_blocking_purify_full
Definition input_constants.F:894
input_constants::do_admm_blocked_projection
integer, parameter, public do_admm_blocked_projection
Definition input_constants.F:894
input_constants::do_admm_aux_exch_func_none
integer, parameter, public do_admm_aux_exch_func_none
Definition input_constants.F:900
input_constants::do_admm_charge_constrained_projection
integer, parameter, public do_admm_charge_constrained_projection
Definition input_constants.F:894
input_constants::do_admm_exch_scaling_none
integer, parameter, public do_admm_exch_scaling_none
Definition input_constants.F:912
input_constants::do_admm_exch_scaling_merlot
integer, parameter, public do_admm_exch_scaling_merlot
Definition input_constants.F:912
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_release
recursive subroutine, public section_vals_release(section_vals)
releases the given object
Definition input_section_types.F:657
kinds
Defines the basic variable types.
Definition kinds.F:23
kinds::dp
integer, parameter, public dp
Definition kinds.F:34
kpoint_transitional
Datatype to translate between k-points (2d) and gamma-point (1d) code.
Definition kpoint_transitional.F:19
kpoint_transitional::get_1d_pointer
type(dbcsr_p_type) function, dimension(:), pointer, public get_1d_pointer(this)
Smart getter, raises an error when called during a k-point calculation.
Definition kpoint_transitional.F:46
kpoint_transitional::get_2d_pointer
type(dbcsr_p_type) function, dimension(:, :), pointer, public get_2d_pointer(this)
Simple getter, needed because of PRIVATE.
Definition kpoint_transitional.F:65
kpoint_transitional::kpoint_transitional_release
subroutine, public kpoint_transitional_release(this)
Release the matrix set, using the right pointer.
Definition kpoint_transitional.F:120
kpoint_transitional::set_1d_pointer
subroutine, public set_1d_pointer(this, ptr_1d)
Assigns a 1D pointer.
Definition kpoint_transitional.F:78
kpoint_transitional::set_2d_pointer
subroutine, public set_2d_pointer(this, ptr_2d)
Assigns a 2D pointer.
Definition kpoint_transitional.F:101
message_passing
Interface to the message passing library MPI.
Definition message_passing.F:23
qs_kind_types
Define the quickstep kind type and their sub types.
Definition qs_kind_types.F:23
qs_kind_types::deallocate_qs_kind_set
subroutine, public deallocate_qs_kind_set(qs_kind_set)
Destructor routine for a set of qs kinds.
Definition qs_kind_types.F:273
qs_local_rho_types
Definition qs_local_rho_types.F:8
qs_local_rho_types::local_rho_set_release
subroutine, public local_rho_set_release(local_rho_set)
...
Definition qs_local_rho_types.F:224
qs_mo_types
Definition and initialisation of the mo data type.
Definition qs_mo_types.F:22
qs_mo_types::deallocate_mo_set
subroutine, public deallocate_mo_set(mo_set)
Deallocate a wavefunction data structure.
Definition qs_mo_types.F:356
qs_mo_types::get_mo_set
subroutine, public get_mo_set(mo_set, maxocc, homo, lfomo, nao, nelectron, n_el_f, nmo, eigenvalues, occupation_numbers, mo_coeff, mo_coeff_b, uniform_occupation, kts, mu, flexible_electron_count)
Get the components of a MO set data structure.
Definition qs_mo_types.F:401
qs_neighbor_list_types
Define the neighbor list data types and the corresponding functionality.
Definition qs_neighbor_list_types.F:17
qs_neighbor_list_types::release_neighbor_list_sets
subroutine, public release_neighbor_list_sets(nlists)
releases an array of neighbor_list_sets
Definition qs_neighbor_list_types.F:1096
qs_oce_types
Definition qs_oce_types.F:8
qs_oce_types::deallocate_oce_set
subroutine, public deallocate_oce_set(oce_set)
Deallocate the matrix set of oce coefficients.
Definition qs_oce_types.F:98
qs_rho_types
superstucture that hold various representations of the density and keeps track of which ones are vali...
Definition qs_rho_types.F:18
qs_rho_types::qs_rho_release
subroutine, public qs_rho_release(rho_struct)
releases a rho_struct by decreasing the reference count by one and deallocating if it reaches 0 (to b...
Definition qs_rho_types.F:115
task_list_types
types for task lists
Definition task_list_types.F:14
task_list_types::deallocate_task_list
subroutine, public deallocate_task_list(task_list)
deallocates the components and the object itself
Definition task_list_types.F:145
admm_dm_types::admm_dm_type
Definition admm_dm_types.F:41
admm_types::admm_gapw_r3d_rs_type
A subtype of the admm_env that contains the extra data needed for an ADMM GAPW calculation.
Definition admm_types.F:85
admm_types::admm_type
stores some data used in wavefunction fitting
Definition admm_types.F:120
cp_blacs_env::cp_blacs_env_type
represent a blacs multidimensional parallel environment (for the mpi corrispective see cp_paratypes/m...
Definition cp_blacs_env.F:53
cp_control_types::admm_control_type
Definition cp_control_types.F:455
cp_dbcsr_api::dbcsr_p_type
Definition cp_dbcsr_api.F:172
cp_fm_struct::cp_fm_struct_type
keeps the information about the structure of a full matrix
Definition cp_fm_struct.F:82
cp_fm_types::cp_fm_type
represent a full matrix
Definition cp_fm_types.F:114
input_section_types::section_vals_type
stores the values of a section
Definition input_section_types.F:127
kpoint_transitional::kpoint_transitional_type
Definition kpoint_transitional.F:30
message_passing::mp_para_env_type
stores all the informations relevant to an mpi environment
Definition message_passing.F:746
qs_kind_types::qs_kind_type
Provides all information about a quickstep kind.
Definition qs_kind_types.F:177
qs_local_rho_types::local_rho_type
Definition qs_local_rho_types.F:45
qs_mo_types::mo_set_type
Definition qs_mo_types.F:47
qs_neighbor_list_types::neighbor_list_set_p_type
Definition qs_neighbor_list_types.F:67
qs_oce_types::oce_matrix_type
Definition qs_oce_types.F:37
qs_rho_types::qs_rho_type
keeps the density in various representations, keeping track of which ones are valid.
Definition qs_rho_types.F:62
task_list_types::task_list_type
Definition task_list_types.F:58