<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
<div class="">Thanks Paolo. I will generate sr locally and that will solve my problem of transforming the electronic velocity vectors from the wedge to the whole BZ.</div>
<div class=""><br class="">
</div>
<div class="">Best wishes,</div>
<div class="">Vahid</div>
<div class=""><br class="">
</div>
<br class="">
<div>
<blockquote type="cite" class="">
<div class="">On Sep 19, 2016, at 3:12 PM, Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com" class="">p.giannozzi@gmail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="ltr" class="">
<div class="">Maybe rotation matrices in cartesian axis (sr) are not properly initialized, because not used, in EPW. The "historical" way of storing rotation matrices is in crystal axis and as integer matrices. By the way, I consider a very bad idea to have
two distinct representations of the same matrices as global variables.<br class="">
<br class="">
</div>
Paolo<br class="">
</div>
<div class="gmail_extra"><br class="">
<div class="gmail_quote">On Mon, Sep 19, 2016 at 2:57 PM, Vahid Askarpour <span dir="ltr" class="">
<<a href="mailto:vh261281@dal.ca" target="_blank" class="">vh261281@dal.ca</a>></span> wrote:<br class="">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div style="word-wrap:break-word" class="">
<div class="">Hi Paolo,</div>
<div class=""><br class="">
</div>
<div class="">I printed the other s(:,:,isym) and got the following for some of them:</div>
<div class=""><br class="">
</div>
<div class="">for s(3) </div>
<div class=""><span class=""> 1.00000000000000 0.000000000000000E+000 0.<wbr class="">000000000000000E+000<br class="">
0.000000000000000E+000 -1.<wbr class="">00000000000000 0.000000000000000E+000<br class="">
</span> 0.000000000000000E+000 0.<wbr class="">000000000000000E+000 -1.<wbr class="">00000000000000 <br class="">
instead of </div>
<div class=""> cart. s( 3) = ( -1.0000000 0.0000000 0.<wbr class="">0000000 )<br class="">
( 0.0000000 1.0000000 0.<wbr class="">0000000 )<br class="">
( 0.0000000 0.0000000 -1.0000000 )<br class="">
<br class="">
</div>
<div class="">for s(4)</div>
<div class=""> 1.00000000000000 0.000000000000000E+000 0.<wbr class="">000000000000000E+000<br class="">
0.000000000000000E+000 1.<wbr class="">00000000000000 0.000000000000000E+000<br class="">
0.000000000000000E+000 0.<wbr class="">000000000000000E+000 -1.<wbr class="">00000000000000 <br class="">
instead of </div>
<div class=""> cart. s( 4) = ( 1.0000000 0.0000000 0.<wbr class="">0000000 )<span class=""><br class="">
( 0.0000000 -1.0000000 0.0000000 )<br class="">
</span> ( 0.0000000 0.0000000 -1.0000000 )<br class="">
for s(5)</div>
<div class=""> 1.00000000000000 0.000000000000000E+000 0.<wbr class="">000000000000000E+000<br class="">
0.000000000000000E+000 1.<wbr class="">00000000000000 0.000000000000000E+000<br class="">
0.000000000000000E+000 0.<wbr class="">000000000000000E+000 -1.<wbr class="">00000000000000 <br class="">
instead of </div>
<div class=""> cart. s( 5) = ( 0.0000000 1.0000000 0.<wbr class="">0000000 ) f =( -0.2500000 )<br class="">
( 1.0000000 0.0000000 0.<wbr class="">0000000 ) ( 0.2500000 )<br class="">
( 0.0000000 0.0000000 -1.0000000 ) ( 0.2500000 )<br class="">
<br class="">
</div>
<div class="">The following two lines are included at the top of the routine:</div>
<div class=""><br class="">
</div>
<div class=""> USE symm_base, ONLY : s, sr,nsym,t_rev, time_reversal, set_sym_bl, nrot<br class="">
USE cell_base, ONLY : alat, at, omega,bg</div>
<div class=""><br class="">
</div>
<div class="">Interestingly, I get the correct symmetries if in the same EPW routine, I do the following (from the symm_base.f90):</div>
<div class=""><br class="">
</div>
<div class=""> DO isym = 1,nsym<br class="">
sa (:,:) = dble ( s(:,:,isym) )<br class="">
sb = matmul ( bg, sa )<br class="">
sr (:,:, isym) = matmul ( at, transpose (sb) )<br class="">
write(*,*) sr(:,:,isym)<br class="">
ENDDO</div>
<div class=""><br class="">
</div>
<div class="">Somehow, sr needs to be recalculated and it is not enough to include USE symm_base and cell_base at the top of the routine.</div>
<div class=""><br class="">
</div>
<div class="">Best,</div>
<div class=""><br class="">
</div>
<div class="">Vahid<br class="">
<br class="">
</div>
<div class="">
<div class="h5"><br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On Sep 19, 2016, at 8:37 AM, Vahid Askarpour <<a href="mailto:vh261281@dal.ca" target="_blank" class="">vh261281@dal.ca</a>> wrote:</div>
<br class="">
<div class="">
<div style="word-wrap:break-word" class="">
<div class="">Hi Paolo,</div>
<div class=""><br class="">
</div>
<div class="">I was attempting to extract the symmetries in an EPW routine in order to simplify the calculations. I used a loop ns=1,nrot and inside this loop, I printed sr(:,:,2). Perhaps I should have looped over nsym instead of nrot. I will check to see
if this makes a difference.</div>
<div class=""><br class="">
</div>
<div class="">Thanks,</div>
<div class=""><br class="">
</div>
<div class="">Vahid</div>
<div class=""><br class="">
</div>
<br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On Sep 19, 2016, at 2:57 AM, Paolo Giannozzi <<a href="mailto:p.giannozzi@gmail.com" target="_blank" class="">p.giannozzi@gmail.com</a>> wrote:</div>
<br class="">
<div class="">
<div dir="ltr" class="">
<div class="">What is printed on output is exactly sr(:,:,isym), with isym running over symmetries. Where and how did you print sr(:,:,2)? and what about all other symmetries?<br class="">
<br class="">
</div>
Paolo</div>
<div class="gmail_extra"><br class="">
<div class="gmail_quote">On Mon, Sep 19, 2016 at 12:28 AM, Vahid Askarpour <span dir="ltr" class="">
<<a href="mailto:vh261281@dal.ca" target="_blank" class="">vh261281@dal.ca</a>></span> wrote:<br class="">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear QE Users,<br class="">
<br class="">
When doing relaxation on diamond, one of the symmetry operations identified by QE-5.4.0 is the following:<br class="">
<br class="">
isym = 2 180 deg rotation - cart. axis [0,0,1]<br class="">
<br class="">
cryst. s( 2) = ( 0 1 -1 )<br class="">
( 1 0 -1 )<br class="">
( 0 0 -1 )<br class="">
<br class="">
cart. s( 2) = ( -1.0000000 0.0000000 0.0000000 )<br class="">
( 0.0000000 -1.0000000 0.0000000 )<br class="">
( 0.0000000 0.0000000 1.0000000 )<br class="">
<br class="">
However, when I print out sr(:,:,2) for this structure (isym=2), I get :<br class="">
<br class="">
1.00000000000000 0.000000000000000E+000 0.000000000000000E+000<br class="">
0.000000000000000E+000 -1.00000000000000 0.000000000000000E+000<br class="">
0.000000000000000E+000 0.000000000000000E+000 1.00000000000000<br class="">
<br class="">
Since the sr matrix prints the cartesian rotation matrix, there seems to be a mismatch between the sr(1,1) elements. The other elements are consistent.<br class="">
<br class="">
Is this a bug or am I missing another operation embedded in the sr matrix?<br class="">
<br class="">
Thank you,<br class="">
<br class="">
Vahid<br class="">
<br class="">
Vahid Askarpour<br class="">
Department of Physics and Atmospheric Science<br class="">
Dalhousie University,<br class="">
Halifax, NS, Canada<br class="">
<br class="">
<br class="">
______________________________<wbr class="">_________________<br class="">
Pw_forum mailing list<br class="">
<a href="mailto:Pw_forum@pwscf.org" target="_blank" class="">Pw_forum@pwscf.org</a><br class="">
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank" class="">http://pwscf.org/mailman/listi<wbr class="">nfo/pw_forum</a><br class="">
</blockquote>
</div>
<br class="">
<br clear="all" class="">
<br class="">
-- <br class="">
<div data-smartmail="gmail_signature" class="">
<div dir="ltr" class="">
<div class="">
<div dir="ltr" class="">
<div class="">Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br class="">
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br class="">
Phone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank" class="">
+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank" class="">
+39-0432-558222</a><br class="">
<br class="">
</div>
</div>
</div>
</div>
</div>
</div>
______________________________<wbr class="">_________________<br class="">
Pw_forum mailing list<br class="">
<a href="mailto:Pw_forum@pwscf.org" target="_blank" class="">Pw_forum@pwscf.org</a><br class="">
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank" class="">http://pwscf.org/mailman/<wbr class="">listinfo/pw_forum</a></div>
</blockquote>
</div>
<br class="">
</div>
______________________________<wbr class="">_________________<br class="">
Pw_forum mailing list<br class="">
<a href="mailto:Pw_forum@pwscf.org" target="_blank" class="">Pw_forum@pwscf.org</a><br class="">
<a href="http://pwscf.org/mailman/listinfo/pw_forum" target="_blank" class="">http://pwscf.org/mailman/<wbr class="">listinfo/pw_forum</a></div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
<br class="">
______________________________<wbr class="">_________________<br class="">
Pw_forum mailing list<br class="">
<a href="mailto:Pw_forum@pwscf.org" class="">Pw_forum@pwscf.org</a><br class="">
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank" class="">http://pwscf.org/mailman/<wbr class="">listinfo/pw_forum</a><br class="">
</blockquote>
</div>
<br class="">
<br clear="all" class="">
<br class="">
-- <br class="">
<div class="gmail_signature" data-smartmail="gmail_signature">
<div dir="ltr" class="">
<div class="">
<div dir="ltr" class="">
<div class="">Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br class="">
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br class="">
Phone +39-0432-558216, fax +39-0432-558222<br class="">
<br class="">
</div>
</div>
</div>
</div>
</div>
</div>
_______________________________________________<br class="">
Pw_forum mailing list<br class="">
<a href="mailto:Pw_forum@pwscf.org" class="">Pw_forum@pwscf.org</a><br class="">
http://pwscf.org/mailman/listinfo/pw_forum</div>
</blockquote>
</div>
<br class="">
</body>
</html>