<div dir="ltr">Dear Juerg, <div><br></div><div><span style="text-align: -webkit-right;">Thanks a lot for your help. I removed the numbers in the atom labels, and the total time of the calculation went down from </span><span style="text-align: -webkit-right;">880.44s to </span>483.034s. </div><div><br></div><div>Dear Matthias, </div><div><br></div><div>There is no specific reason for having an atomic kind for each atom. I generated the cif file with an ASE script, and the numbers were added to the atom lables. Now that I am aware of the issue, I will delete the numbers from the atom label before I use the cif file in a calculation. <br>My input files contain a lot of if clauses and @include statements. I posted the .restart file because I thought that it might be easier to understand what keywords I had used. Now I realize that it might have caused unnecessary confusion. </div><div><br></div><div>Regards,</div><div>Torstein</div><div><br></div><div><br></div><div><div><br></div><div><br><br>onsdag 8. januar 2020 17.59.54 UTC+1 skrev Torstein Fjermestad følgende:<blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div dir="ltr">Dear all, <div><br></div><div>In a geometry optimization, the subroutine "compute_max_radius" is taking about 1/4 of the total time (see the first lines of the timing information below. "compute_max_radius" is shown in bold). </div><div><br></div><div>What is "compute_max_radius" doing?</div><div>Is it normal that it takes that long time?</div><div>Is there anything one can do to reduce the time it takes?</div><div><br></div><div>Thanks in advance for your help.</div><div><br></div><div>Regards,</div><div>Torstein Fjermestad</div><div><br></div><div><br></div><div><table border="0" cellpadding="0" cellspacing="0" width="647" style="border-collapse:collapse;width:487pt">
<colgroup><col width="199" style="width:150pt">
<col width="90" style="width:68pt">
<col width="64" style="width:48pt">
<col width="70" style="width:53pt">
<col width="77" style="width:58pt">
<col width="67" style="width:50pt">
<col width="80" style="width:60pt">
</colgroup><tbody><tr height="19" style="height:14.5pt">
<td height="19" width="199" style="height:14.5pt;width:150pt">
SUBROUTINE
</td>
<td width="90" style="width:68pt"> CALLS</td>
<td width="64" style="width:48pt"> ASD</td>
<td width="70" style="width:53pt"> </td>
<td width="77" style="width:58pt">SELF TIME</td>
<td width="67" style="width:50pt"></td>
<td width="80" style="width:60pt">TOTAL TIME</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div>
</div>
</td>
<td> MAXIMUM </td>
<td> </td>
<td> AVERAGE </td>
<td>MAXIMUM </td>
<td>AVERAGE </td>
<td>MAXIMUM</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> CP2K
</div>
</td>
<td> 1 </td>
<td>1.0 </td>
<td> 0.560 </td>
<td> 0.562 </td>
<td>880.431 </td>
<td align="right">880.44</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> cp_geo_opt
</div>
</td>
<td> 1 </td>
<td>2.0 </td>
<td> 0.060 </td>
<td> 0.060 </td>
<td>632.220 </td>
<td align="right">632.22</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> geoopt_bfgs
</div>
</td>
<td> 1 </td>
<td>3.0 </td>
<td> 0.174 </td>
<td> 0.180 </td>
<td>632.160 </td>
<td align="right">632.16</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> cp_eval_at
</div>
</td>
<td> 35 </td>
<td>4.0 </td>
<td> 0.036 </td>
<td> 0.078 </td>
<td>627.687 </td>
<td align="right">627.71</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> qs_forces
</div>
</td>
<td> 34 </td>
<td>5.0 </td>
<td> 0.113 </td>
<td> 0.118 </td>
<td>614.752 </td>
<td align="right">614.76</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> qs_energies
</div>
</td>
<td> 35 </td>
<td>6.0 </td>
<td> 0.012 </td>
<td> 0.015 </td>
<td>565.518 </td>
<td align="right">565.55</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> scf_env_do_scf
</div>
</td>
<td> 35 </td>
<td>7.0 </td>
<td> 0.002 </td>
<td> 0.005 </td>
<td>381.893 </td>
<td align="right">381.90</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> scf_env_do_scf_inner_loop </div>
</td>
<td> 164 </td>
<td>8.0 </td>
<td> 0.028 </td>
<td> 0.032 </td>
<td>260.120 </td>
<td align="right">261.02</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> qs_init_subsys
</div>
</td>
<td> 1 </td>
<td>2.0 </td>
<td> 0.406 </td>
<td> 0.408 </td>
<td>245.615 </td>
<td align="right">245.62</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> qs_env_setup
</div>
</td>
<td> 1 </td>
<td>3.0 </td>
<td> 0.043 </td>
<td> 0.043 </td>
<td>244.914 </td>
<td align="right">244.92</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> qs_env_rebuild_pw_env
</div>
</td>
<td> 70 </td>
<td>6.5 </td>
<td> 0.069 </td>
<td> 0.070 </td>
<td>244.794 </td>
<td align="right">244.80</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> pw_env_rebuild
</div>
</td>
<td> 1 </td>
<td>5.0 </td>
<td> 0.092 </td>
<td> 0.129 </td>
<td>244.709 </td>
<td align="right">244.71</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> <b>compute_max_radius
</b></div>
</td>
<td><b> 1 </b></td>
<td><b>6.0 </b></td>
<td><b>232.268 </b></td>
<td><b>244.451 </b></td>
<td><b>232.268 </b></td>
<td align="right"><b>244.45</b></td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> rebuild_ks_matrix
</div>
</td>
<td>198 </td>
<td>9.5 </td>
<td> 0.001 </td>
<td> 0.001 </td>
<td>198.294 </td>
<td align="right">198.45</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> qs_ks_build_kohn_sham_matrix</div>
</td>
<td> 198</td>
<td>10.5 </td>
<td> 0.078 </td>
<td> 0.082 </td>
<td>198.293 </td>
<td align="right">198.45</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> qs_ks_update_qs_env
</div>
</td>
<td>199 </td>
<td>9.0 </td>
<td> 0.002 </td>
<td> 0.002 </td>
<td>159.038 </td>
<td align="right">159.18</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> qs_energies_init_hamiltonians</div>
</td>
<td>35 </td>
<td>7.0 </td>
<td> 0.033 </td>
<td> 0.034 </td>
<td>151.114 </td>
<td align="right">151.13</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> build_qs_neighbor_lists
</div>
</td>
<td> 35 </td>
<td>8.0 </td>
<td> 0.187 </td>
<td> 0.212 </td>
<td>100.828 </td>
<td align="right">136.07</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> init_scf_loop
</div>
</td>
<td> 35 </td>
<td>8.0 </td>
<td> 0.001 </td>
<td> 0.002 </td>
<td>120.711 </td>
<td align="right">120.72</td>
</tr>
<tr height="19" style="height:14.5pt">
<td height="19" style="height:14.5pt">
<div> pw_transfer
<div></div>
</div>
</td>
<td> 4199</td>
<td>13.3 </td>
<td> 0.347 </td>
<td> 0.410 </td>
<td> 93.912 </td>
<td> 95.354</td>
</tr></tbody></table></div><div> </div></div></blockquote></div></div></div>