<div dir="ltr">Dear Prof. Lorenzo Paulatto,<br>                                    I was trying to run projected density of states from an nscf calculation
 with spin-orbit included and I am getting the following error:<br><br><br clear="all"><div> Calling projwave_nc ....<br><br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>     Error in routine d_matrix_so (5):<br>     D_S (j=1/2) for this symmetry operation is not unitary<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br><br>     stopping ...<br><br></div><div>I looked at a previous discussion with you at pw_forum in <br><a href="http://qe-forge.org/pipermail/pw_forum/2012-June/098763.html">http://qe-forge.org/pipermail/pw_forum/2012-June/098763.html</a> <br></div><div>where you have suggested to increase eps and print out capel in d_matrix_so.f90 . So I have done that and the capel value in my case is 1.3756250103975467, which is quite large. I am running it on espresso-5.0.1 and am attaching the scf, nscf and 
pdos inputs so that one can reproduce the error. It would be very helpful for me if you can suggest anything to resolve this issue.<br><br> Thanks in advance,<br clear="all"></div>-- <br><div class="gmail_signature"><div dir="ltr">Rajdeep Banerjee<br>PhD student<div>JNCASR, Bangalore</div><div>India</div></div></div>
</div>