[Wannier] a little problem with pwscf 3.2.1?

Tadashi Ogitsu ogitsu@llnl.gov
Mon, 7 May 2007 11:13:44 -0700


Hi Arash,

Thanks for the response. Well, then, I got a non-official version  
probably...? I thought, I downloaded from the webpage, but I cannot  
be 100% sure.

I'll download the official distribution and replace with the one I  
have now just in case. Sorry about the confusion.

Tadashi
On May 4, 2007, at 7:18 PM, Arash A Mostofi wrote:

> Dear Tadashi
>
> I'm slightly confused: the eigenvalues have always (ie, since the  
> original incarnation of pw2wannier90 in December 2005) been written  
> by pw2wannier90 as:
>
>   write (iun_band,'(2i5,f18.12)') ibnd1, ikevc, et(ibnd,ik)*rytoev
>
> (ie, without any reference to weights, wg). This is easily verified  
> by checking the espresso cvs repository.
>
> Consequently, Wannier90 has always read the eigenvalues as:
>
>    read(eig_unit,*,err=106,end=106) i,j,eigval(n,k)
>
> as I've verified by checking our Wannier90 cvs repository.
>
> Can you remember from where you got your copy of the Wannier90  
> source? Was it from the www.wannier.org or from someone who may  
> have made some modifications for his/her own purposes? In any case,  
> the change you've made is correct.
>
> Best regards
>
> Arash
>
> -----------------------------------------
> Dr Arash A Mostofi
> Dept of Materials Science and Engineering
> 13-4084 MIT
> 77 Massachusetts Avenue
> Cambridge MA 02139-4307 USA
> Tel: 617.2536026
> -----------------------------------------
>
> On Fri, 4 May 2007, Tadashi Ogitsu wrote:
>
>> Hi Wannier90s,
>>
>> I've found a small problem with, I think, pw2wannier90 came with  
>> pwscf3.2.1 (3.2 plus the distributed patch on their wabpage).
>>
>> Somewhat they have changed the output format of .eig file so that  
>> the last data, the weight of each eigenvalues, does not exist  
>> anymore, while the current distribution of wannier90 try to read  
>> and complains that the data is not consistent.
>>
>> I just removed the last read as follows,
>>
>> old:
>> read(eig_unit,*,err=106,end=106) i,j,eigval(n,k),wg(n,k)
>> new:
>> read(eig_unit,*,err=106,end=106) i,j,eigval(n,k)
>>
>> First of all, is my understanding correct? (i.e., this is a  
>> problem caused by the change in pw2wannier90 in 3.2.1 distribution)
>>
>> Second, did I do the right correction? (I have no idea if wg is  
>> used anywhere, although it does not seem to be...)
>>
>> Thanking in advance!
>>
>> Tadashi
>> PS. Jonathan, I've started to work on the gold transport and as  
>> the first try, I'm playing with the example6.
>> _______________________________________________
>> Wannier mailing list
>> Wannier@quantum-espresso.org
>> http://www.democritos.it/mailman/listinfo/wannier
> _______________________________________________
> Wannier mailing list
> Wannier@quantum-espresso.org
> http://www.democritos.it/mailman/listinfo/wannier