<div dir="ltr"><div><div><div><div><div><div><div><div><div>Hi all,<br><br></div>When I performed phonon calculation in QE 5.3.0 on BGQ,<br></div>I got set_irr_sym_new<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br>     Error in routine set_irr_sym_new (211):<br>     wrong representation<br> %%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br></div>This error doesn't appear all the time but on some specific q points in a grid.<br><br></div>By comparing with the same calculation performed on x86,<br></div>I noticed that some of the eigenvalues (autoval) computed for the q points emitting errors are quite different. So I went to check how those values are computed. I found that the cdiagh (in PW/src/cdiagh.f90) are branched using ZHPEV in ESSL or ZHEEV in LAPACK. So I changed the route to force the LAPACK one and my issue disappears. Though this solution seems solve my issue but I was not sure it address the key problem.<br><br></div><div>My library linking order is SCALAPACK ESSL LAPACK ESSL (for FFT and BLAS) <br></div>I think this is the right order for BGQ.<br></div>Since, the cdiagh is not reported at the end of pw or ph output, I guess it takes little time during the run and I won't get any performance penalty.<br><br></div>Any comments or suggestions? particularly from developers.<br><br></div><div>PS: I send an email a few weeks ago and got no reply. I hope someone can also answer these questions.<br></div><div>1) Any update for releasing DFT+U phonon? I saw the paper was out years ago.<br></div><div>2) Any ways to reduce the phonon disk I/O to 'none' similar to pw.x? Reduce_io seems help very little during the calculation.<br><br></div><div>Best regards,<br></div>Ye <br><div><div><div><br clear="all"><div><div><div><div><div><div><div><div><div><div><div><div><div><div dir="ltr">===================<br>
Ye Luo, Ph.D.<br>
Leadership Computing Facility<br>
Argonne National Laboratory</div></div></div>
</div></div></div></div></div></div></div></div></div></div></div></div></div></div></div>