<div dir="ltr"><div>Dear pw users,</div><div>the patch provided in reply to my previous email by Paolo Giannozzi works fine.</div><div>Thanks !</div><div><br></div><div>Best,<br></div><div>D.</div><div><br></div><div>P.S.: sorry for not directly replying to the email by Paolo,</div><div>but I used the "digest mode" for the mailing list and I did not understand<br></div><div>how to reply to the email from there.<br></div></div><br><div class="gmail_quote"><div dir="ltr">On Wed, 12 Dec 2018 at 10:23, Davide Sangalli <<a href="mailto:davide.sangalli@gmail.com">davide.sangalli@gmail.com</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div class="gmail_quote"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>Dear pw users,</div><div>I recently tried to run projwfc.x on three systems with SOC and systematically run into the following error:</div><div><br></div><div><div>%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>     Error in routine d_matrix_so (7):<br>     D_S (j=1/2) <span class="gmail-m_3654515318314890285m_94754221507729258m_1426920682240473270m_-118291583349210998gmail-il">for</span> <span class="gmail-m_3654515318314890285m_94754221507729258m_1426920682240473270m_-118291583349210998gmail-il">this</span> <span class="gmail-m_3654515318314890285m_94754221507729258m_1426920682240473270m_-118291583349210998gmail-il">symmetry</span> <span class="gmail-m_3654515318314890285m_94754221507729258m_1426920682240473270m_-118291583349210998gmail-il">operation</span> <span class="gmail-m_3654515318314890285m_94754221507729258m_1426920682240473270m_-118291583349210998gmail-il">is</span> <span class="gmail-m_3654515318314890285m_94754221507729258m_1426920682240473270m_-118291583349210998gmail-il">not</span> <span class="gmail-m_3654515318314890285m_94754221507729258m_1426920682240473270m_-118291583349210998gmail-il">unitary</span><br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br></div><div><br></div><div>I reproduced the error onto three different systems:</div><div>MoS2, bulk iron, FeRh<br></div><div>with both the last git revision ( commit 3e920635bff849e891bb3fcec99acc096c3bdcf6)</div><div>and qe-6.2.1<br></div><div><br></div><div>I had a look to the value of "capel" in the subroutine</div><div>"PP/src/d_matrix_so.f90"</div><div>for the case of MoS2 (inputs attached).</div><div><br></div><div>In the loop over symmetries, "capel" becomes 2 for the symmetry number 9, since</div><div>ZDOTU(2, d12_con(1,m), 1, dyj12(1,n,isym), 1)</div><div>is identically zero for all n and m.</div><div><br></div><div>Kind regards,<br></div><div>D.</div><div><br></div><div><br></div><div><br></div><div><br></div></div></div></div></div></div>
</div></div>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>Davide Sangalli, PhD<br></div><div>CNR-ISM, Division of Ultrafast Processes in Materials (FLASHit) and MaX Centre</div><div>Area della Ricerca di Roma 1, 00016 Monterotondo Scalo, Italy<br></div><a href="http://www.ism.cnr.it/en/davide-sangalli-cv/" target="_blank">http://www.ism.cnr.it/en/davide-sangalli-cv/<br></a><a href="http://www.max-center.eu/" target="_blank">http://www.max-centre.eu/</a><br></div></div></div></div></div></div>