[Wien] problems in lapw3

B. Yanchitsky yan at im.imag.kiev.ua
Thu Nov 16 15:41:34 CET 2006


Peter Blaha wrote:
> I think interchanging the indices in rotdef.f of SRC_lapw3 
> should fix the problem.

Yes, this works, but also the line
REAL(8),INTENT(IN)  :: pos(3,ndif),tau(3,iord)
must be changed to
REAL(8),INTENT(IN)  :: pos(nat*48,3),tau(3,iord)

I placed changed rotdef.f in attachment, it looks that all is correct now.
I have another question, is it possible to change the code in such a way to print atomic
factors for some predefined mesh of k-points? Total structure factor is nonzero only
for translation vectors in k-space, but atomic factors are nonzero everywhere.

Bogdan Yanchitsky
Institute of Magnetism
Vernadsky Blvd., 36-b
03142  Kiev
UKRAINE

Tel. (+380-44) 444 34 20
Fax. (+380-44) 444 10 20
-------------- next part --------------
A non-text attachment was scrubbed...
Name: rotdef.f.gz
Type: application/x-gzip
Size: 1118 bytes
Desc: not available
Url : http://zeus.theochem.tuwien.ac.at/pipermail/wien/attachments/20061116/31ea70a1/rotdef.f.gz


More information about the Wien mailing list