Skip to Content.
Sympa Menu

forum - Re: [abinit-forum] frustrating RF calculation

forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

Re: [abinit-forum] frustrating RF calculation


Chronological Thread 
  • From: "matthieu verstraete" <matthieu.jean.verstraete@gmail.com>
  • To: sarayazdani@yahoo.com
  • Cc: forum@abinit.org
  • Subject: Re: [abinit-forum] frustrating RF calculation
  • Date: Tue, 9 Dec 2008 20:54:21 +0100
  • Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:to:subject:cc:in-reply-to:mime-version :content-type:references; b=GRlEc6PVbULOinWRBmwG10CEWCev3i31FZwwjjiFdxHSqEpGKdVFon6ku+yARMiShS X4ctpEpBdasmPUUPHgmwWd7BCmVFj19DhnMw1In9gvfebeWKlBBtE2KAD5M8pM9B/gCF khQldGKTtW9+Wk4T6AvDO/UE4VD2FGX7EGxPw=


Hello Sara,

a dirty fix is fairly simple: in src/18seqpar/loper3.F90, around line 1220, comment the allocation of gh1_rbz and replace it by

!      allocate(gh1_rbz(nkpt_rbz,dtset%mband,2,mpw1*nspinor))
      allocate(gh1_rbz(1,1,1,1))

Recompile and run. This will of course crash if you try to use Paul's developments, but should be fine for you (keep ieig2rf to its default value of 0).

Cheers

Matthieu


On Tue, Dec 9, 2008 at 11:13 AM, sara yazdani <sarayazdani@yahoo.com> wrote:
Hello Paul and abinit users,

Thank you for your reply. Since I am in short of time to present my results and I have done some other good calculations with abinit, please let me know what I should do instead to get better phonon frequencies. Since I used "rfasr=1 and 2" for all the calculations as well and the result was worse.

Regards,
Sara Yazdani



--- On Tue, 9/12/08, Paul Boulanger <paul.boulanger@umontreal.ca> wrote:
From: Paul Boulanger <paul.boulanger@umontreal.ca>
Subject: Re: [abinit-forum] frustrating RF calculation
To: forum@abinit.org, sarayazdani@yahoo.com
Date: Tuesday, 9 December, 2008, 1:58 AM


Hello Sara,

For the 8x8x8 grid, the memory allocation error is actually my fault, I
included in past versions an array that scales badly with the number of
K-points. This is fixed in my 5.6.5 public branch. I am terribly sorry if this
issue as bothered/hindered anybody else. I will discuss with Xavier to see what
is the best way to make it accessible for you and others before the release of
5.6.5.

Paul




--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Dr. Matthieu Verstraete

European Theoretical Spectroscopy Facility (ETSF) 
Dpto. Fisica de Materiales, 
U. del Pais Vasco,            
Centro Joxe Mari Korta, Av. de Tolosa, 72,   Phone: +34-943018393
E-20018 Donostia-San Sebastian, Spain        Fax  : +34-943018390

Mail : matthieu.jean.verstraete@gmail.com
http://www-users.york.ac.uk/~mjv500



Archive powered by MHonArc 2.6.15.

Top of Page