<div dir="ltr"><div>Hi Sergey,</div><div><br></div><div>The code works well and can indeed reduce the number of required SCF iterations considerably. There is one caveat though, I had to add the following lines to qs_initial_guess.F before the call to</div><div><br></div><div>CALL calculate_atomic_orbitals(atomic_kind, qs_kind, iunit=ounit, &<br>                                        pmat=pmat(ikind)%mat, nocc=nocc)</div><div><br></div><div>in SUBROUTINE calculate_atomic_block_dm</div><div><br></div><div>         write(*,*) ikind, nkind<br>         write(*,*) nocc<br>         write(*,*) TRIM(atomic_kind%name)</div><div><br></div><div>Without them I always got an undefined reference error and cp2k crashed, with them everything works as intended. This happens only with the ATOMIC_RESTART. Using the ATOMIC guess cp2k runs without the print statements.</div><div><br></div><div>To compile cp2k I use the toolchain with gfortran both on Piz Daint and a local cluster.</div><div><br></div><div>Thanks again for the feature,</div><div>Fabian<br></div></div>