<div dir="ltr">So do I. I just corrected it. Paolo<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Nov 6, 2015 at 4:38 PM, stefano de gironcoli <span dir="ltr"><<a href="mailto:degironc@sissa.it" target="_blank">degironc@sissa.it</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On 06/11/2015 16:21, Lorenzo Paulatto wrote:<br>
> Dear All,<br>
> there is a small bug or a problem with the documentation of ibrav= -9<br>
><br>
><br>
> The Doc says:<br>
><br>
> 9          Orthorhombic base-centered(bco) celldm(2)=b/a<br>
>                                             celldm(3)=c/a<br>
>      v1 = (a/2, b/2,0),  v2 = (-a/2,b/2,0),  v3 = (0,0,c)<br>
> -9          as 9, alternate description<br>
>      v1 = (a/2,-b/2,0),  v2 = (a/2,-b/2,0),  v3 = (0,0,c)<br>
</span>v1 = v2 ?<br>
i think the code is right and the doc is wrong... or both are wrong :-)<br>
<span class="HOEnZb"><font color="#888888">stefano<br>
</font></span><span class="im HOEnZb">><br>
> However the code says:<br>
><br>
>       IF ( ibrav == 9 ) THEN<br>
>          !   old PWscf description<br>
>          a1(1) = 0.5d0 * celldm(1)<br>
>          a1(2) = a1(1) * celldm(2)<br>
>          a2(1) = - a1(1)<br>
>          a2(2) = a1(2)<br>
>       ELSE<br>
>          !   alternate description<br>
>          a1(1) = 0.5d0 * celldm(1)<br>
>          a1(2) =-a1(1) * celldm(2)<br>
>          a2(1) = a1(1)<br>
>          a2(2) =-a1(2)<br>
>       END IF<br>
><br>
> Which means:<br>
> -9          as 9, alternate description<br>
>      v1 = (a/2,-b/2,0),  v2 = (a/2,b/2,0),  v3 = (0,0,c)<br>
><br>
><br>
> I.e the sign of the y-component of v2 is the opposite.<br>
><br>
> I'm writing here because I do not know which of the two is correct. Thanks to<br>
> G Radtke for bringing this to my attention (He thinks that the Doc is correct<br>
> and the code is wrong)<br>
><br>
> LP<br>
><br>
><br>
<br>
</span><div class="HOEnZb"><div class="h5">_______________________________________________<br>
Q-e-developers mailing list<br>
<a href="mailto:Q-e-developers@qe-forge.org">Q-e-developers@qe-forge.org</a><br>
<a href="http://qe-forge.org/mailman/listinfo/q-e-developers" rel="noreferrer" target="_blank">http://qe-forge.org/mailman/listinfo/q-e-developers</a><br>
</div></div></blockquote></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><span><span><font color="#888888">Paolo Giannozzi, Dept. Chemistry&Physics&Environment,<br>
Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>
Phone <a href="tel:%2B39-0432-558216" value="+390432558216" target="_blank">+39-0432-558216</a>, fax <a href="tel:%2B39-0432-558222" value="+390432558222" target="_blank">+39-0432-558222</a></font></span></span></div></div></div></div>
</div>