<div dir="ltr"><div><div>Hello all,</div><div><br></div><div>I am having this trouble when running matdyn.x with fd=.true. and na_ifc=.true. when doing the phDOS with QE 5.4.0.</div><div><font color="#9900ff"><br></font></div><div><font color="#9900ff">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</font></div><div><font color="#9900ff"><br></font></div><div><font color="#9900ff">     Error in routine cdiagh2 (404):</font></div><div><font color="#9900ff">     info =/= 0</font></div><div><font color="#9900ff">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</font></div><div><font color="#9900ff"><br></font></div><div><font color="#9900ff">     stopping ...</font></div><div><br></div><div>It will work without error for some q point grid densities: 5 5 5 or 7 7 7</div><div><br></div><div>but not others: 6 6 6 or 2 2 2 or 30 30 30</div><div><br></div><div>The q point grids that work seem to vary depending on structure and chemistry. I think the problem may be occurring during the call to ZHEEV in lapack (I'm using intel MKL) but I am not sure.</div><div><br></div><div>The FD phonon calcs seem to be fine and matdyn.x for phonon dispersion doesn't seem to have this problem. I appreciate any thoughts on this.</div><div><br></div><div>inputs/outputs as well as details in included libraries are attached as a tar to conserve space. (the ifc file for matdyn.x is over 200kb alone)</div><div><br></div><div><br></div><div><br></div><div>Thanks,</div><div><br></div><div>Andrew Supka</div><div><br></div><div>Graduate Research Asst.</div><div><br></div><div>Central Michigan University</div></div></div>