bug fix help in paral_kgb case  [SOLVED]

Documentation, Web site and code modifications

Moderators: baguetl, routerov

Locked
User avatar
jzwanzig
Posts: 504
Joined: Mon Aug 17, 2009 9:25 am

bug fix help in paral_kgb case

Post by jzwanzig » Thu Jun 04, 2015 8:09 pm

Hi,
my routine make_efg_el.F90 fails in the fft parallel case. The idea in this routine is to FFT the converged ground state density rhor to reciprocal space, and do some sums involving the G vectors and the density. In the non-parallel case I knew the size of rhor and its transform based on elements of nfft(18) but in the parallel case apparently I don't anymore. Can someone direct me to the right documentation to understand how the parallel fft case works?

thanks,
Joe
Josef W. Zwanziger
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com

User avatar
torrent
Posts: 127
Joined: Fri Aug 14, 2009 7:40 pm

Re: bug fix help in paral_kgb case  [SOLVED]

Post by torrent » Tue Jun 09, 2015 5:21 pm

Hi Joe,

You should have a look in routines doing a similar job, for instance hartre.F90.

I see several modifications to perform:
- eliminate mpi_enreg_seq and replace it by the true mpi_enreg (passed as argument)
- call the ptabs_fourdp routine to inialize (if necessary) the parallel fft tabs
- add a line "if (fftn2_distrib(i2) == me_fft) then" in the igfft2 loop
- call xmpi_reduce on efg(:,:,:) after the loops to merge the computation of the procs

This done, you'll see if debug is necessary or not...

A bientôt,
Marc
Marc Torrent
CEA - Bruyères-le-Chatel
France

User avatar
jzwanzig
Posts: 504
Joined: Mon Aug 17, 2009 9:25 am

Re: bug fix help in paral_kgb case

Post by jzwanzig » Wed Jun 10, 2015 5:49 pm

thanks, Marc, this is really helpful--I figured I could understand what's going on if I had an example to compare to.

By way did you see my question posted in the response function section about DDK eigenvalues? I'd be grateful to get your take on that--

thanks
Joe
Josef W. Zwanziger
Professor, Department of Chemistry
Canada Research Chair in NMR Studies of Materials
Dalhousie University
Halifax, NS B3H 4J3 Canada
jzwanzig@gmail.com

Locked