[Pw_forum] El-Ph funny error
Miguel Martínez Canales
wmbmacam at lg.ehu.es
Tue Oct 4 18:20:15 CEST 2005
Dear Dr. Gianozzi
First, thanks to you and Dr.Isaev for your answers.
>
> IF ( ABS( xq(1) - ( xk(1,ikq) - xk(1,ikk) ) ) > eps8 .OR. &
> ABS( xq(2) - ( xk(2,ikq) - xk(2,ikk) ) ) > eps8 .OR. &
> ABS( xq(3) - ( xk(3,ikq) - xk(3,ikk) ) ) > eps8 ) THEN
> WRITE( stdout, * ) ikk, ikq, nksq
> WRITE( stdout, * ) ( xq(ipol), ipol = 1, 3)
> WRITE( stdout, * ) ( xk(ipol,ikq), ipol = 1, 3)
> WRITE( stdout, * ) ( xk(ipol,ikk), ipol = 1, 3)
> CALL errore( 'phq_init', 'wrong order of k points', 1 )
> END IF
>
> what does it print immediately before that?
>
The outpur (from the last k point to the error is):
k(16250) = ( 0.4200000 1.4200000 -0.2500000), wk = 0.0000000
pseudo 1 is Al (US) zval = 3.0 lmax= 1 lloc= 0
Version 0 0 0 of US pseudo code
Using log mesh of 893 points
The pseudopotential has 2 beta functions with:
l(1) = 0
l(2) = 1
Q(r) pseudized with 8 coefficients, rinner = 0.900 0.900 0.900
Atomic displacements:
There are 3 irreducible representations
Representation 1 1 modes - To be done
Representation 2 1 modes - To be done
Representation 3 1 modes - To be done
1 2 8125
-0.125000000000000 -0.125000000000000 0.375000000000000
-2.00000000000000D-002 2.00000000000000D-002 0.270000000000000
-2.00000000000000D-002 2.00000000000000D-002 2.00000000000000D-002
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
from phq_init : error # 1
wrong order of k points
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
stopping ...
If it helps, this tipically occurs in low symmetry points (i.e. outside 00x
or xxx lines). Thanks again, and sorry for the late answer,
Miguel
--
----------------------------------------
Miguel Martínez Canales
Dto. Física de la Materia Condensada
UPV/EHU
Facultad de Ciencia y Tecnología
Apdo. 644
48080 Bilbao (Spain)
Fax: +34 94 601 3500
Tlf: +34 94 601 5437
----------------------------------------
2005 International Year of Physics
More information about the Pw_forum
mailing list