<div dir="ltr"><div><div>Dear all,</div><div>I find some problems in running qs calculations with cp2K 4.1 compiled with intelmpi-2017 (the same run goes fine with the executable obtained with openmpi-gnu compiler).</div><div>In particular in output I obtain</div></div><div><br></div><div><div>----------------------------------- OT ---------------------------------------</div><div><br></div><div>  Step     Update method      Time    Convergence         Total energy    Change</div><div>  ------------------------------------------------------------------------------</div><div><br></div><div>  Trace(PS):                                 1200.0000000051</div><div>  Electronic density on regular grids:                   NaN                 NaN</div><div>  Core density on regular grids:             1200.0000000000       -0.0000000000</div><div>  Total charge density on r-space grids:                 NaN</div><div>  Total charge density g-space grids:          -5.8357006210</div></div><div><br></div><div>Unlike the code compiled with openmpi-gnu gives:</div><div><br></div><div><div> ----------------------------------- OT ---------------------------------------</div><div><br></div><div>  Step     Update method      Time    Convergence         Total energy    Change</div><div>  ------------------------------------------------------------------------------</div><div><br></div><div>  Trace(PS):                                 1199.9999998902</div><div>  Electronic density on regular grids:      -1199.9999998901        0.0000001099</div><div>  Core density on regular grids:             1199.9999999999       -0.0000000001</div><div>  Total charge density on r-space grids:        0.0000001098</div><div>  Total charge density g-space grids:           0.0000001099</div><div><br></div></div><div>Clearly there is something wrong with the quantities</div><div>Electronic density on regular grids<br></div><div>Total charge density on r-space grids<br></div><div><br></div><div>Looking at the source code I find that the problem may come from the quantities tot_tho_r and tot_rho_r_arr, in qs_ks_utils.F</div><div><br></div><div>Line 855 in qs_ks_utils.F<br><div>CALL qs_rho_get(rho, tot_rho_r=tot_rho_r_arr, rho_ao_kp=rho_ao)<br></div></div><div><br></div><div>If I print tot_rho_r_arr after this call I obtain NaN for both its components</div><div>and as consequences also</div><div>tot_rho_r = accurate_sum(tot_rho_r_arr) <br></div><div>is NaN </div><div> while if I run the executable gnu it gives the right value</div><div>tot_rho_r = accurate_sum(tot_rho_r_arr) =- 1199.99999989<br></div><div><br></div><div>I attach the restart file used for the calculations.</div><div><br></div><div>Can you help me to fix this problem?<br></div><div><br></div><div>Best regards,</div><div>Mariella</div></div>