forum@abinit.org

Subject: The ABINIT Users Mailing List ( CLOSED )

List archive

[abinit-forum] Reading wavefunction output file: how to obtain GS or RF from header?

From: Rob <spamrefuse@yahoo.com>
To: abinit <forum@abinit.org>
Subject: [abinit-forum] Reading wavefunction output file: how to obtain GS or RF from header?
Date: Sun, 21 Jun 2009 19:43:51 -0700 (PDT)
Domainkey-signature: a=rsa-sha1; q=dns; c=nofws; s=s1024; d=yahoo.com; h=Message-ID:X-YMail-OSG:Received:X-Mailer:Date:From:Subject:To:MIME-Version:Content-Type; b=XJjtHIOY08wt8B+vig9jkmCtItkOZGBVx68Sli9Lt7OBzS4WNe0K57flrg+QS2zc9mNpMGrA+KR63lo+Awnl4eb6suq8PohiPOK3UfE2BLrqJhaZiY52YMQ3N/DS8jakrmPQn2IXdHGOXslGxOA1tFT/aVC2DBR7p97jEg1Mnkc=;

Hello,

I'm writing a routine to read the binary wavefunction output file
(foobar_WFK).
The format depends on whether it was a Ground State (GS) or Response
Function (RF) calculation.

How can I determine the GS or RF from the binary file?
I can't see which info in the header tells me this.

(I know how to determine cplex and for cplex = 2 it must be RF;
however, both GS and RF can have cplex = 1).

================================

Another question:
The description in the abinit-help website says for the ground state wf file:

bantot=0
do isppol=1,nsppol
do ikpt=1,nkpt
write(unit) npw,nspinor,nband
write(unit) kg(1:3,1:npw)
write(unit) eigen(1+bantot:nband+bantot),occ(1+bantot:nband+bantot)
do iband=1,nband
write(unit) (cg(ii+index),ii=1,2*npw*nspinor)
enddo
bantot=bantot+nband
enddo
enddo

What is 'index' here? Should that be 'iband' or 'ikpt' instead?
Thanks.
Rob.

[abinit-forum] Reading wavefunction output file: how to obtain GS or RF from header?, Rob, 06/22/2009
- Re: [abinit-forum] Reading wavefunction output file: how to obtain GS or RF from header?, Xavier Gonze, 06/24/2009