RF calculation (iscf=-3) is very slow...

Phonons, DFPT, electron-phonon, electric-field response, mechanical response…

Moderators: mverstra, joaocarloscabreu

Locked
Hiroki
Posts: 5
Joined: Mon Nov 19, 2018 12:33 am

RF calculation (iscf=-3) is very slow...

Post by Hiroki » Sun Dec 09, 2018 8:55 am

Dear ABINIT developers/users

I tried to calculate the response functions with iscf = -3 and previous SCF result(WFK file).
The calculation terminated successfully, but took a lot of time (3-4 times longer than SCF cycle). I think this is a problem.
Looking the log file to survey the cause, I found below comments:

Code: Select all

--- !WARNING
src_file: mpi_setup.F90
src_line: 672
message: |
    Your number of spins*k-points (=4) and bands (=400) will not distribute      correctly
    with the current number of processors (=144).
    You will leave some empty.
    YOU ARE STRONGLY ADVICED TO ACTIVATE AUTOMATIC PARALLELIZATION!
    PUT "AUTOPARAL=1" IN THE INPUT FILE.
...

 -   nproc =  144   -> not optimal: autoparal keyword recommended in input file

In this run, however, I put "autoparal=1" in my input file.
I can't figure out what this message means...

And I tried to calculate without autoparal, set parameters with respect to parallelization using paral_kgb.
But I got the following comment again:

Code: Select all

--- !WARNING
src_file: mpi_setup.F90
src_line: 672
message: |
    Your number of spins*k-points (=4) and bands (=400) will not distribute      correctly
    with the current number of processors (=144).
    You will leave some empty.
    YOU ARE STRONGLY ADVICED TO ACTIVATE AUTOMATIC PARALLELIZATION!
    PUT "AUTOPARAL=1" IN THE INPUT FILE.

--- !WARNING
src_file: mpi_setup.F90
src_line: 240
message: |
    For non ground state calculation, set bandpp, npfft, npband, npspinor npkpt and nphf to 1
...

So, according to error comments, I need to use autoparal...
What should I do to optimize the parallelization for RF calculation?

The attached files are what I used in this calculation.

Sincerely yours,
Attachments
RF_abi.in
Input file for RF calculation
(4.88 KiB) Downloaded 233 times
SCF-nSCF_abi.in
Input file for SCF and nSCF calculation
(5.1 KiB) Downloaded 227 times

ebousquet
Posts: 469
Joined: Tue Apr 19, 2011 11:13 am
Location: University of Liege, Belgium

Re: RF calculation (iscf=-3) is very slow...

Post by ebousquet » Fri Dec 14, 2018 4:25 pm

Dear Hiroki,
You should indeed adapt the number of CPU w.r.t. the number of k-points and the number of bands, etc.
Please have alook here on advises to do so here:
viewtopic.php?f=8&t=3837

For a convergence problem in this calculation, I also advise you to set nline to a larger value than the default, something like 8 to 12.

Let us know if this fix your problems,
Best wishes,
Eric

Locked