Branch data Line data Source code
1 : :
2 : : #if defined HAVE_CONFIG_H
3 : : #include "config.h"
4 : : #endif
5 : : !!****m* ABINIT/m_CtqmcoffdiagInterfaceComplex
6 : : !! NAME
7 : : !! m_CtqmcoffdiagInterfaceComplex
8 : : !!
9 : : !! FUNCTION
10 : : !! Manage a ctqmc simulation.
11 : : !! friendly interface for the user
12 : : !!
13 : : !! COPYRIGHT
14 : : !! Copyright (C) 2013-2025 ABINIT group (J. Bieder, B. Amadon, J. Denier)
15 : : !! This file is distributed under the terms of the
16 : : !! GNU General Public License, see ~abinit/COPYING
17 : : !! or http://www.gnu.org/copyleft/gpl.txt .
18 : : !!
19 : : !! NOTES
20 : : !!
21 : : !! SOURCE
22 : :
23 : : #include "defs.h"
24 : : MODULE m_CtqmcoffdiagInterfaceComplex
25 : : USE m_CtqmcoffdiagComplex
26 : : use defs_basis
27 : :
28 : : IMPLICIT NONE
29 : :
30 : : !!***
31 : :
32 : : !!****t* m_CtqmcoffdiagInterfaceComplex/CtqmcoffdiagInterfaceComplex
33 : : !! NAME
34 : : !! CtqmcoffdiagInterfaceComplex
35 : : !!
36 : : !! FUNCTION
37 : : !! This structured datatype contains the necessary data
38 : : !!
39 : : !! COPYRIGHT
40 : : !! Copyright (C) 2013-2025 ABINIT group (J. Bieder)
41 : : !! This file is distributed under the terms of the
42 : : !! GNU General Public License, see ~abinit/COPYING
43 : : !! or http://www.gnu.org/copyleft/gpl.txt .
44 : : !!
45 : : !! SOURCE
46 : :
47 : : TYPE CtqmcoffdiagInterfaceComplex
48 : : TYPE(CtqmcoffdiagComplex) :: Hybrid
49 : : INTEGER :: opt_fk = 0
50 : : INTEGER :: opt_order = 0
51 : : INTEGER :: opt_histo = 0
52 : : INTEGER :: opt_movie = 0
53 : : INTEGER :: opt_analysis = 0
54 : : INTEGER :: opt_check = 0
55 : : INTEGER :: opt_spectra = 0
56 : : INTEGER :: opt_noise = 0
57 : : INTEGER :: opt_gMove = 0
58 : : END TYPE CtqmcoffdiagInterfaceComplex
59 : : !!***
60 : :
61 : : CONTAINS
62 : : !!***
63 : :
64 : : !!****f* ABINIT/m_CtqmcoffdiagInterfaceComplex/CtqmcoffdiagInterfaceComplex_init
65 : : !! NAME
66 : : !! CtqmcoffdiagInterfaceComplex_init
67 : : !!
68 : : !! FUNCTION
69 : : !! Initialize with permanent parameters
70 : : !!
71 : : !! COPYRIGHT
72 : : !! Copyright (C) 2013-2025 ABINIT group (J. Bieder)
73 : : !! This file is distributed under the terms of the
74 : : !! GNU General Public License, see ~abinit/COPYING
75 : : !! or http://www.gnu.org/copyleft/gpl.txt .
76 : : !!
77 : : !! INPUTS
78 : : !! op=ctqmcinterface
79 : : !! iseed=seed for rng
80 : : !! sweeps=number of sweeps for the run
81 : : !! thermalization=number of sweeps to thermalize
82 : : !! measurements=how often we measure (modulo)
83 : : !! flavors=number of orbitals (spin degenerated)
84 : : !! samples=imaginary time slices
85 : : !! beta=inverse temperature
86 : : !! U=interaction parameter
87 : : !! ostream=where to write output
88 : : !! MPI_COMM=mpi communicator for the run
89 : : !! nspinor=number of spinor
90 : : !!
91 : : !! OUTPUT
92 : : !!
93 : : !! SIDE EFFECTS
94 : : !!
95 : : !! NOTES
96 : : !!
97 : : !! SOURCE
98 : :
99 : 4 : SUBROUTINE CtqmcoffdiagInterfaceComplex_init(op,iseed,sweeps,thermalization,&
100 : : &measurements,flavors,samples,beta,U,ostream,MPI_COMM,opt_nondiag,nspinor)
101 : :
102 : : !Arguments ------------------------------------
103 : : TYPE(CtqmcoffdiagInterfaceComplex), INTENT(INOUT) :: op
104 : : INTEGER, OPTIONAL, INTENT(IN) :: MPI_COMM
105 : : INTEGER, INTENT(IN) :: iseed
106 : : DOUBLE PRECISION, INTENT(IN) :: sweeps
107 : : INTEGER, INTENT(IN) :: thermalization
108 : : INTEGER, INTENT(IN) :: measurements
109 : : INTEGER, INTENT(IN) :: flavors
110 : : INTEGER, INTENT(IN) :: samples
111 : : !INTEGER, INTENT(IN) :: Wmax
112 : : INTEGER, INTENT(IN) :: ostream
113 : : INTEGER, INTENT(IN) :: opt_nondiag
114 : : DOUBLE PRECISION, INTENT(IN) :: beta
115 : : DOUBLE PRECISION, INTENT(IN) :: u
116 : : INTEGER, INTENT(IN) :: nspinor
117 : : !DOUBLE PRECISION, INTENT(IN) :: mu
118 : : !Local arguements -----------------------------
119 : : INTEGER :: ifstream
120 : : DOUBLE PRECISION, DIMENSION(1:11) :: buffer
121 : :
122 : 4 : ifstream = 42
123 : :
124 : 4 : buffer(1)=DBLE(iseed)
125 : 4 : buffer(2)=sweeps
126 : 4 : buffer(3)=DBLE(thermalization)
127 : 4 : buffer(4)=DBLE(measurements)
128 : 4 : buffer(5)=DBLE(flavors)
129 : 4 : buffer(6)=DBLE(samples)
130 : 4 : buffer(7)=beta
131 : 4 : buffer(8)=U
132 : 4 : buffer(9)=GREENHYB_TAU
133 : 4 : buffer(10)=DBLE(opt_nondiag)
134 : 4 : buffer(11)=nspinor
135 : : !buffer(9)=0.d0!mu
136 : : !buffer(9)=DBLE(Wmax)
137 : :
138 [ + - ]: 4 : IF ( PRESENT( MPI_COMM ) ) THEN
139 : 4 : CALL CtqmcoffdiagComplex_init(op%Hybrid, ostream, ifstream, .FALSE., MPI_COMM,buffer)
140 : : ELSE
141 : 0 : CALL CtqmcoffdiagComplex_init(op%Hybrid, ostream, ifstream, .FALSE.,iBuffer=buffer)
142 : : END IF
143 : 4 : op%opt_fk = 0
144 : 4 : op%opt_order = 0
145 : 4 : op%opt_histo = 0
146 : 4 : op%opt_movie = 0
147 : 4 : op%opt_analysis = 0
148 : 4 : op%opt_check = 0
149 : 4 : op%opt_noise = 0
150 : 4 : op%opt_spectra = 0
151 : 4 : END SUBROUTINE CtqmcoffdiagInterfaceComplex_init
152 : : !!***
153 : :
154 : : !!****f* ABINIT/m_CtqmcoffdiagInterfaceComplex/CtqmcoffdiagInterfaceComplex_setOpts
155 : : !! NAME
156 : : !! CtqmcoffdiagInterfaceComplex_setOpts
157 : : !!
158 : : !! FUNCTION
159 : : !! Set and save options for many runs
160 : : !!
161 : : !! COPYRIGHT
162 : : !! Copyright (C) 2013-2025 ABINIT group (J. Bieder)
163 : : !! This file is distributed under the terms of the
164 : : !! GNU General Public License, see ~abinit/COPYING
165 : : !! or http://www.gnu.org/copyleft/gpl.txt .
166 : : !!
167 : : !! INPUTS
168 : : !! op=ctqmcinterface
169 : : !! opt_Fk=0 if we give us Gw0 and 1 if 1/Gw0+iwn
170 : : !! opt_order=maximal perturbation order to scope(>0)
171 : : !! opt_movie=print a latex file (0 or 1)
172 : : !! opt_analysis=measure correlations (0 or 1)
173 : : !! opt_check=check fast calculation :0 nothing
174 : : !! 1 Impurity
175 : : !! 2 Bath
176 : : !! 3 Both
177 : : !! opt_noise=calculate noise ofr green functions(0 ro 1)
178 : : !! opt_spectra=fourier transform of time evolution of number of electrons
179 : : !! (0 or 1)
180 : : !! opt_gMove=number of global moves (>0)
181 : : !!
182 : : !! OUTPUT
183 : : !!
184 : : !! SIDE EFFECTS
185 : : !!
186 : : !! NOTES
187 : : !!
188 : : !! SOURCE
189 : :
190 : 4 : SUBROUTINE CtqmcoffdiagInterfaceComplex_setOpts(op,opt_Fk,opt_order,opt_histo,opt_movie,&
191 : : & opt_analysis,opt_check, opt_noise, opt_spectra, opt_gMove)
192 : :
193 : : !Arguments ------------------------------------
194 : : TYPE(CtqmcoffdiagInterfaceComplex), INTENT(INOUT) :: op
195 : : INTEGER , OPTIONAL , INTENT(IN ) :: opt_Fk
196 : : INTEGER , OPTIONAL , INTENT(IN ) :: opt_order
197 : : INTEGER , OPTIONAL , INTENT(IN ) :: opt_histo
198 : : INTEGER , OPTIONAL , INTENT(IN ) :: opt_movie
199 : : INTEGER , OPTIONAL , INTENT(IN ) :: opt_analysis
200 : : INTEGER , OPTIONAL , INTENT(IN ) :: opt_check
201 : : INTEGER , OPTIONAL , INTENT(IN ) :: opt_noise
202 : : INTEGER , OPTIONAL , INTENT(IN ) :: opt_spectra
203 : : INTEGER , OPTIONAL , INTENT(IN ) :: opt_gMove
204 : :
205 [ + - ]: 4 : IF ( PRESENT(opt_Fk) ) &
206 : 4 : op%opt_Fk = opt_fk
207 [ + - ]: 4 : IF ( PRESENT(opt_order) ) &
208 : 4 : op%opt_order = opt_order
209 [ + - ]: 4 : IF ( PRESENT(opt_histo) ) &
210 : 4 : op%opt_histo = opt_histo
211 [ + - ]: 4 : IF ( PRESENT(opt_analysis) ) &
212 : 4 : op%opt_analysis = opt_analysis
213 [ + - ]: 4 : IF ( PRESENT(opt_check) ) &
214 : 4 : op%opt_check = opt_check
215 [ + - ]: 4 : IF ( PRESENT(opt_movie) ) &
216 : 4 : op%opt_movie = opt_movie
217 [ + - ]: 4 : IF ( PRESENT(opt_noise) ) &
218 : 4 : op%opt_noise = opt_noise
219 [ + - ]: 4 : IF ( PRESENT(opt_spectra) ) &
220 : 4 : op%opt_spectra = opt_spectra
221 [ + - ]: 4 : IF ( PRESENT(opt_gMove) ) &
222 : 4 : op%opt_gMove = opt_gMove
223 : :
224 : 4 : END SUBROUTINE CtqmcoffdiagInterfaceComplex_setOpts
225 : : !!***
226 : :
227 : : !!****f* ABINIT/m_CtqmcoffdiagInterfaceComplex/CtqmcoffdiagInterfaceComplex_run
228 : : !! NAME
229 : : !! CtqmcoffdiagInterfaceComplex_run
230 : : !!
231 : : !! FUNCTION
232 : : !! run a ctqmc simu and get results
233 : : !!
234 : : !! COPYRIGHT
235 : : !! Copyright (C) 2013-2025 ABINIT group (J. Bieder)
236 : : !! This file is distributed under the terms of the
237 : : !! GNU General Public License, see ~abinit/COPYING
238 : : !! or http://www.gnu.org/copyleft/gpl.txt .
239 : : !!
240 : : !! INPUTS
241 : : !! op=ctqmcinterface
242 : : !! G0omega=Gw0 (according to opt_Fk)
243 : : !! matU=interaction matrice
244 : : !! opt_sym=weight factors to symmetrise G
245 : : !! opt_levels=energy for each level (with respect to fermi level)
246 : : !!
247 : : !! OUTPUT
248 : : !! Gtau=G(tau)
249 : : !! Gw=fourier transform of Gtau
250 : : !! D=full double occupancy
251 : : !! E=Interaction energy
252 : : !! Noise=Noise on E
253 : : !!
254 : : !! SIDE EFFECTS
255 : : !!
256 : : !! NOTES
257 : : !!
258 : : !! SOURCE
259 : :
260 [ + - + - : 8 : SUBROUTINE CtqmcoffdiagInterfaceComplex_run(op,G0omega, Gtau, Gw, D,E,Noise,matU,Docc,opt_sym,opt_levels,hybri_limit,Magmom_orb,&
- + + - +
- - + - +
+ - + - -
+ + - + -
- + - + -
- - - + -
+ - - + +
- + - - +
+ - + - -
+ ]
261 [ + - + - : 4 : &Magmom_spin,Magmom_tot,Iatom,fname,opthybri)
- + + - +
- - + ]
262 : :
263 : : !Arguments ------------------------------------
264 : : TYPE(CtqmcoffdiagInterfaceComplex), INTENT(INOUT) :: op
265 : : COMPLEX(KIND=8) , DIMENSION(:,:,:), INTENT(IN ) :: G0omega
266 : : COMPLEX(KIND=8), DIMENSION(:,:,:), OPTIONAL, INTENT(OUT) :: Gtau
267 : : COMPLEX(KIND=8) , DIMENSION(:,:,:), OPTIONAL, INTENT(INOUT) :: Gw
268 : : DOUBLE PRECISION, OPTIONAL , INTENT(OUT) :: D
269 : : DOUBLE PRECISION, OPTIONAL , INTENT(OUT) :: E
270 : : DOUBLE PRECISION, OPTIONAL , INTENT(OUT) :: Noise
271 : : COMPLEX(KIND=8), DIMENSION(:,:), OPTIONAL, INTENT(IN ) :: matU
272 : : DOUBLE PRECISION, DIMENSION(:,:),OPTIONAL, INTENT(OUT ) :: Docc
273 : : DOUBLE PRECISION, DIMENSION(:,:),OPTIONAL, INTENT(IN ) :: opt_sym
274 : : COMPLEX(KIND=8), DIMENSION(:), OPTIONAL, INTENT(IN ) :: opt_levels
275 : : COMPLEX(KIND=8) , DIMENSION(:,:),OPTIONAL, INTENT(IN ) :: hybri_limit
276 : : DOUBLE PRECISION, DIMENSION(:,:),OPTIONAL, INTENT(IN ) :: Magmom_orb
277 : : DOUBLE PRECISION, DIMENSION(:,:),OPTIONAL, INTENT(IN ) :: Magmom_spin
278 : : DOUBLE PRECISION, DIMENSION(:,:),OPTIONAL, INTENT(IN ) :: Magmom_tot
279 : : INTEGER, INTENT(IN ) :: Iatom
280 : : INTEGER, INTENT(IN ) :: opthybri
281 : : character(len=fnlen), INTENT(INOUT) :: fname
282 : : !local variables--------------------------------
283 : : ! INTEGER :: iflavor1,iflavor2
284 : :
285 : : ! do iflavor1=1,10
286 : : ! do iflavor2=1,10
287 : : ! if(iflavor1==iflavor2) THEN
288 : : ! write(6,*) iflavor1, iflavor2, Magmom(iflavor1,iflavor2)
289 : : ! end if
290 : : ! end do
291 : : ! end do
292 : :
293 : 4 : CALL CtqmcoffdiagComplex_reset(op%Hybrid)
294 : :
295 [ + - ]: 4 : IF ( PRESENT(opt_levels)) &
296 : 4 : CALL CtqmcoffdiagComplex_setMu(op%Hybrid, opt_levels)
297 : :
298 [ + - ]: 4 : IF ( PRESENT(hybri_limit)) &
299 : 4 : CALL CtqmcoffdiagComplex_sethybri_limit(op%Hybrid, hybri_limit,opthybri)
300 : :
301 : : !call xmpi_barrier(op%Hybrid%MY_COMM)
302 : 4 : CALL CtqmcoffdiagComplex_setG0wTab(op%Hybrid, G0omega,op%opt_fk,Iatom,fname)
303 : : !call xmpi_barrier(op%Hybrid%MY_COMM)
304 : :
305 [ + - ]: 4 : IF ( PRESENT(matU) ) &
306 : 4 : CALL CtqmcoffdiagComplex_setU(op%Hybrid, matU)
307 : : !call xmpi_barrier(op%Hybrid%MY_COMM)
308 : :
309 [ + - ]: 4 : IF ( PRESENT(Magmom_orb) ) &
310 : 4 : CALL CtqmcoffdiagComplex_setMagmom(op%Hybrid, Magmom_orb, Magmom_spin, Magmom_tot)
311 : : !call xmpi_barrier(op%Hybrid%MY_COMM)
312 : :
313 : :
314 : : CALL CtqmcoffdiagComplex_run(op%Hybrid,opt_order=op%opt_order, &
315 : : opt_histo=op%opt_histo, &
316 : : opt_movie=op%opt_movie, &
317 : : opt_analysis=op%opt_analysis, &
318 : : opt_check=op%opt_check, &
319 : : opt_noise=op%opt_noise, &
320 : : opt_spectra=op%opt_spectra, &
321 : 4 : opt_gMove=op%opt_gMove)
322 : : !call xmpi_barrier(op%Hybrid%MY_COMM)
323 : :
324 : : ! write(6,*) "op%Hybrid%stats",op%Hybrid%stats
325 : : ! write(6,*) "opt_gMove",op%opt_gMove
326 : :
327 : 4 : CALL CtqmcoffdiagComplex_getResult(op%Hybrid,Iatom,fname)
328 : :
329 [ - + ]: 4 : IF ( PRESENT(opt_sym) ) THEN
330 : 0 : CALL CtqmcoffdiagComplex_symmetrizeGreen(op%Hybrid,opt_sym)
331 : : END IF
332 : :
333 : : ! write(6,*) "op%Hybrid%stats",op%Hybrid%stats
334 : :
335 [ + - ]: 4 : IF ( PRESENT(Gtau) .AND. PRESENT(Gw) ) THEN
336 : 4 : CALL CtqmcoffdiagComplex_getGreen(op%Hybrid, Gtau=Gtau, Gw=Gw)
337 : : ! write(6,*) "size",size(Gw,dim=1),size(gw,dim=2)
338 : : ! IF ( hybrid%rank .EQ. 0 ) write(389,*) Gw(:,hybrid%flavors+1)
339 : : ! call flush(389)
340 [ # # ]: 0 : ELSE IF ( PRESENT(Gtau) .AND. .NOT. PRESENT(Gw) ) THEN
341 [ # # ]: 0 : CALL CtqmcoffdiagComplex_getGreen(op%Hybrid, Gtau=Gtau)
342 [ # # ]: 0 : ELSE IF ( .NOT. PRESENT(Gtau) .AND. PRESENT(Gw) ) THEN
343 [ # # ]: 0 : CALL CtqmcoffdiagComplex_getGreen(op%Hybrid, Gw=Gw)
344 : : END IF
345 : :
346 : : !write(6,*) "op%Hybrid%stats",op%Hybrid%stats
347 : :
348 [ + - ]: 4 : IF ( PRESENT(D) ) &
349 : 4 : CALL CtqmcoffdiagComplex_getD(op%Hybrid, D)
350 [ + + + + ]: 844 : Docc=op%Hybrid%measDE
351 : :
352 : : !write(6,*) "op%Hybrid%stats",op%Hybrid%stats
353 : :
354 [ + - ]: 4 : IF ( PRESENT(E) .AND. PRESENT(Noise) ) &
355 : 4 : CALL CtqmcoffdiagComplex_getE(op%Hybrid, E, Noise)
356 : :
357 : 4 : CALL CtqmcoffdiagComplex_printAll(op%Hybrid)
358 : : !CALL CtqmcoffdiagComplex_printQMC(op%Hybrid)
359 : :
360 : 4 : END SUBROUTINE CtqmcoffdiagInterfaceComplex_run
361 : : !!***
362 : :
363 : : !!****f* ABINIT/m_CtqmcoffdiagInterfaceComplex/CtqmcoffdiagInterfaceComplex_setSweeps
364 : : !! NAME
365 : : !! CtqmcoffdiagInterfaceComplex_setSweeps
366 : : !!
367 : : !! FUNCTION
368 : : !! change sweeps on the fly
369 : : !!
370 : : !! COPYRIGHT
371 : : !! Copyright (C) 2013-2025 ABINIT group (J. Bieder)
372 : : !! This file is distributed under the terms of the
373 : : !! GNU General Public License, see ~abinit/COPYING
374 : : !! or http://www.gnu.org/copyleft/gpl.txt .
375 : : !!
376 : : !! INPUTS
377 : : !! op=ctqmcinterface
378 : : !! sweeps=new number of sweeps
379 : : !!
380 : : !! OUTPUT
381 : : !!
382 : : !! SIDE EFFECTS
383 : : !!
384 : : !! NOTES
385 : : !!
386 : : !! SOURCE
387 : :
388 : 0 : SUBROUTINE CtqmcoffdiagInterfaceComplex_setSweeps(op, sweeps)
389 : :
390 : : !Arguments ------------------------------------
391 : : TYPE(CtqmcoffdiagInterfaceComplex), INTENT(INOUT) :: op
392 : : DOUBLE PRECISION, INTENT(IN) :: sweeps
393 : :
394 : 0 : CALL CtqmcoffdiagComplex_setSweeps(op%Hybrid,sweeps)
395 : 0 : END SUBROUTINE CtqmcoffdiagInterfaceComplex_setSweeps
396 : : !!***
397 : :
398 : : !!****f* ABINIT/m_CtqmcoffdiagInterfaceComplex/CtqmcoffdiagInterfaceComplex_finalize
399 : : !! NAME
400 : : !! CtqmcoffdiagInterfaceComplex_finalize
401 : : !!
402 : : !! FUNCTION
403 : : !! Destroy simulation
404 : : !!
405 : : !! COPYRIGHT
406 : : !! Copyright (C) 2013-2025 ABINIT group (J. Bieder)
407 : : !! This file is distributed under the terms of the
408 : : !! GNU General Public License, see ~abinit/COPYING
409 : : !! or http://www.gnu.org/copyleft/gpl.txt .
410 : : !!
411 : : !! INPUTS
412 : : !! op=ctqmcinterface
413 : : !!
414 : : !! OUTPUT
415 : : !!
416 : : !! SIDE EFFECTS
417 : : !!
418 : : !! NOTES
419 : : !!
420 : : !! SOURCE
421 : :
422 : 4 : SUBROUTINE CtqmcoffdiagInterfaceComplex_finalize(op)
423 : :
424 : : !Arguments ------------------------------------
425 : : TYPE(CtqmcoffdiagInterfaceComplex), INTENT(INOUT) :: op
426 : :
427 : : !IF ( op%Hybrid%init .EQV. .TRUE. ) THEN
428 : : ! CALL CtqmcoffdiagComplex_printAll(op%Hybrid)
429 : : !write(6,*) "before ctqmc_destroy in CtqmcoffdiagComplexinterface_finalize"
430 : 4 : CALL CtqmcoffdiagComplex_destroy(op%Hybrid)
431 : : !END IF
432 : :
433 : 4 : END SUBROUTINE CtqmcoffdiagInterfaceComplex_finalize
434 : : !!***
435 : :
436 [ # # # # : 0 : END MODULE m_CtqmcoffdiagInterfaceComplex
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # # #
# # # #
# ]
437 : : !!***
|
