<DIV>Dear Dr. Lorenzo Paulatto and QE experts,</DIV>
<DIV> </DIV>
<DIV>    I have a question about phonon calculation with qe-6.0 for large unit cell. Recently, I encountered</DIV>
<DIV>an error like "Error in routine set_irr_sym_new (1322): wrong representation" in the phonon calculations </DIV>
<DIV>for two compounds, in which there are 64 and 70 atoms in the unit cells, respectively. I searched the </DIV>
<DIV>PW_forum and found you have given a reply to this problem. As a temporary workaround, you wrote </DIV>
<DIV>"you can fiddle a bit with the threshold at line 258 of PH/set_irr.f90", i.e. (ABS(eigen(imode)) + ABS (eigen (jmode) )) < 1.d-4) ), </DIV>
<DIV>which is now at line 240 in qe-6.0. Beside above modification, do I need to increase the threshold </DIV>
<DIV>value in 156-159 line in set_irr_sym.f90? Since this error comes from the routine set_irr_sym.f90.</DIV>
<DIV> </DIV>
<DIV>[ line 156 to line 159 in set_irr_sym.f90:<BR>              if (jpert.ne.ipert .and. abs(wrk)> 1.d-6 ) &<BR>                     call errore('set_irr_sym_new','wrong representation',100*irr+10*jpert+ipert)<BR>              if (jpert.eq.ipert .and. abs(wrk-1.d0)> 1.d-6 ) &<BR>                     call errore('set_irr_sym_new','wrong representation',100*irr+10*jpert+ipert)  ]</DIV>
<DIV> </DIV>
<DIV>Beside above temporary solution, is there a better one for this problem? I also noticed that you mentioned</DIV>
<DIV>that you have a personalised version of the random_matrix routine, could you please send a copy to me? </DIV>
<DIV>Your reply is very valuable to me, any help will be highly appreciated! Thanks in advance!</DIV>
<DIV><BR>Best wishes,<BR>Miao Gao</DIV>
<DIV> </DIV>
<DIV>Ningbo University, Zhejiang, China.<BR><BR></DIV><br>