<div dir="ltr"><div><div>Unfortunately, the calculation is still going and statistics are not yet printed.<br><br></div>Best<br></div>Lorenzo Lagasco <br><div><div><br></div></div></div><br><div class="gmail_quote gmail_quote_container"><div dir="ltr" class="gmail_attr">Il giorno mar 19 mag 2026 alle ore 13:46 Johann Pototschnig <<a href="mailto:pototschnig.johann@gmail.com">pototschnig.johann@gmail.com</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">If the computation gets slower with more OMP threads, most of the time they are run on the same CPU. <br><div>You can add to mpirun "--bind-to none" or explicit binding. </div><div>Since you have 52 CPUs you can 13 MPI processes with 4 threads each. Depending on the method it can be faster than 52 MPI processes. </div><div><br></div><div>At the end of the output file statistics are provided (showing which routine/ library used how much time). </div><div>Can you ether provide the output or this information? This should show which library/routine is costly. </div><div><br></div><div><br></div><div>best</div><br><div class="gmail_quote"><div dir="auto" class="gmail_attr">On Tuesday, May 19, 2026 at 11:46:57 AM UTC+2 Lorenzo Lagasco wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><div>Good morning everyone,<br></div>I'm launching a HSEO6  calculation on a 395 atom large system for computing  PDOS and inter-eletronic coupling using Kondov diabatization scheme and Im trying to increase the velocity calculation. The convergence criteria are not particularly strict for scf convergence (as you can check from the input here copied). <br>&GLOBAL<br>  PROJECT CuAlO2-HSE06<br>  RUN_TYPE ENERGY<br>  PRINT_LEVEL LOW<br>&END GLOBAL<br>&FORCE_EVAL<br><br>  METHOD QS<br>  &DFT<br>    WFN_RESTART_FILE_NAME /auto/tms7/llagasco/DELAFOSSITI-proj/INTERFACE/DOS-interface-HSE06/1.PBE/CuAlO2-PBE-RESTART.wfn<br>    LSD<br>    MULTIPLICITY 2<br>    BASIS_SET_FILE_NAME ./BASIS_MOLOPT<br>    BASIS_SET_FILE_NAME ./BASIS_ADMM<br>    BASIS_SET_FILE_NAME ./BASIS_ADMM_MOLOPT<br>    POTENTIAL_FILE_NAME ./POTENTIAL<br>    &MGRID<br>      NGRIDS 5<br>      CUTOFF 500<br>      REL_CUTOFF 50<br>    &END MGRID<br>    &AUXILIARY_DENSITY_MATRIX_METHOD<br>      METHOD BASIS_PROJECTION<br>      ADMM_PURIFICATION_METHOD NONE<br>    &END<br>    &QS<br>      EPS_DEFAULT          1.0E-10<br>      EPS_PGF_ORB          1.0E-6<br>      EXTRAPOLATION PS<br>      EXTRAPOLATION_ORDER 3<br>     #&DDAPC_RESTRAINT<br>     #  STRENGTH   1.0<br>     #  TARGET    -1.0<br>     #  ATOMS       43<br>     #  TYPE_OF_DENSITY SPIN<br>     #  FUNCTIONAL_FORM RESTRAINT<br>     #&END<br>    &END QS<br>    &SCF<br>      EPS_SCF     1.0E-6<br>     #SCF_GUESS ATOMIC<br>      SCF_GUESS RESTART<br>      MAX_SCF 50<br>      &OUTER_SCF<br>        EPS_SCF 1.0E-6<br>        MAX_SCF 50<br>      &END<br>     &OT ON<br>        MINIMIZER DIIS<br>        PRECONDITIONER FULL_SINGLE_INVERSE<br>        ENERGY_GAP 0.25<br>        LINESEARCH 2PNT<br>      &END OT<br>    &END SCF<br>    &XC<br>      &XC_FUNCTIONAL<br>       &PBE<br>         SCALE_X 0.88<br>         SCALE_C 1.0<br>       &END PBE<br>       &PBE_HOLE_T_C_LR<br>           SCALE_X 0.12       ! + 25% of truncated PBE0 functional - that includes exact hfx<br>           CUTOFF_RADIUS 4.0  ! that has interaction truncated at 4.0 A from the atomic core<br>        &END<br>      &END XC_FUNCTIONAL<br>      &HF<br>        &SCREENING<br>          EPS_SCHWARZ 1.0E-6<br>          SCREEN_ON_INITIAL_P TRUE<br>        &END<br>        &INTERACTION_POTENTIAL<br>          POTENTIAL_TYPE TRUNCATED<br>          CUTOFF_RADIUS 4.0<br>          T_C_G_DATA /soft_rocky8/prod/cp2k/12may_2025/gnu14.2.0-openmpi4.1.6/cp2k/data/t_c_g.dat<br>        &END<br>        &MEMORY<br>          MAX_MEMORY 2000<br>          EPS_STORAGE_SCALING 0.1<br>        &END<br>        &END<br>        FRACTION 0.12<br>       &END<br>      &END XC<br>      &DENSITY_FITTING<br>      &END<br>    &PRINT<br>      &PDOS<br>        NLUMO 100<br>        COMPONENTS .TRUE.<br>      &END PDOS<br>    &END PRINT<br>    &END DFT<br>  &SUBSYS<br>    &CELL<br>    A     1.7166299819900004E+01    0.0000000000000000E+00    0.0000000000000000E+00<br>    B     0.0000000000000000E+00    1.7678800582899999E+01    0.0000000000000000E+00<br>    C     0.0000000000000000E+00    0.0000000000000000E+00    3.3554000854500003E+01<br>    PERIODIC  XYZ<br>    MULTIPLE_UNIT_CELL  1 1 1<br>    &END CELL<br> &TOPOLOGY<br>    COORD_FILE_NAME <a href="http://final-geom-opt.xyz" rel="nofollow" target="_blank">final-geom-opt.xyz</a><br>    COORD_FILE_FORMAT XYZ<br>  &END TOPOLOGY<br><br>    &KIND Cu<br>      ELEMENT Cu<br>      BASIS_SET DZVP-MOLOPT-SR-GTH<br>      BASIS_SET AUX_FIT cFIT10<br>      POTENTIAL GTH-PBE-q11<br>       &BS  T<br>         &ALPHA<br>           NEL  2 2<br>           L  2 0<br>           N  3 4<br>         &END ALPHA<br>         &BETA<br>           NEL  -2 -2<br>           L  2 0<br>           N  3 4<br>         &END BETA<br>       &END BS<br>    &END KIND<br><br>    &KIND Al<br>      ELEMENT Al<br>      BASIS_SET DZVP-MOLOPT-SR-GTH<br>      BASIS_SET AUX_FIT cFIT3<br>      POTENTIAL GTH-PBE-q3<br>    &END KIND<br><br>    &KIND O_D<br>      ELEMENT O<br>      BASIS_SET DZVP-MOLOPT-SR-GTH<br>      BASIS_SET AUX_FIT cFIT3<br>      POTENTIAL GTH-PBE-q6<br>    &END KIND<br>   &KIND O_COOH<br>      ELEMENT O<br>      BASIS_SET DZVP-MOLOPT-SR-GTH<br>      BASIS_SET AUX_FIT cFIT3<br>      POTENTIAL GTH-PBE-q6<br>    &END KIND<br><br>     &KIND H<br>       ELEMENT H<br>       BASIS_SET DZVP-MOLOPT-SR-GTH<br>       BASIS_SET AUX_FIT cFIT3<br>       POTENTIAL GTH-PBE-q1<br>     &END KIND<br><br>     &KIND H_COOH<br>       ELEMENT H<br>       BASIS_SET DZVP-MOLOPT-SR-GTH<br>       BASIS_SET AUX_FIT cFIT3<br>       POTENTIAL GTH-PBE-q1<br>     &END KIND<br><br>    &KIND O<br>      ELEMENT O<br>      BASIS_SET DZVP-MOLOPT-SR-GTH<br>      BASIS_SET AUX_FIT cFIT3<br>      POTENTIAL GTH-PBE-q6<br>    &END KIND<br><br>    &KIND N<br>       ELEMENT N<br>       BASIS_SET DZVP-MOLOPT-SR-GTH<br>       BASIS_SET AUX_FIT cFIT3<br>       POTENTIAL GTH-PBE-q5<br>    &END KIND<br><br>    &KIND C<br>       ELEMENT C<br>       BASIS_SET DZVP-MOLOPT-SR-GTH<br>       BASIS_SET AUX_FIT cFIT3<br>       POTENTIAL GTH-PBE-q4<br>    &END KIND<br><br>    &KIND C_COOH<br>       ELEMENT C<br>       BASIS_SET DZVP-MOLOPT-SR-GTH<br>       BASIS_SET AUX_FIT cFIT3<br>       POTENTIAL GTH-PBE-q4<br>    &END KIND<br><br>  &END SUBSYS<br><br> Moreover, I am running the calculation on 4 nodes with 52 processors each. The most surprising thing is that increasing the number of OMP threads beyond 1 actually makes the calculation slower. It is already quite slow — a single SCF step takes about 11 minutes.This is my SLUM submission script:<br><br>#!/bin/bash<br>#SBATCH --job-name=HSE06<br>#SBATCH --nodes=4<br>#SBATCH --ntasks-per-node=52        <br>#SBATCH --partition=taras2-6230r<br><br>module purge<br>module load cp2k/may2025-gnu14.2.0-openmpi4.1.6-psm211.2.230<br><br># OpenMP settings<br>export OMP_NUM_THREADS=1<br><br># Run<br>mpirun cp2k.popt -i SLAB+DYE.inp -o SLAB+DYE.out<br><br></div>I was wondering if it is possible to increase the calculation speed. I hope this email is clear<div>~                                                       <br><br><br></div></div>
</blockquote></div>

<p></p>

-- <br>
You received this message because you are subscribed to the Google Groups "cp2k" group.<br>
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com" target="_blank">cp2k+unsubscribe@googlegroups.com</a>.<br>
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/1be59005-2f09-41c5-a952-d71f5f2526c9n%40googlegroups.com?utm_medium=email&utm_source=footer" target="_blank">https://groups.google.com/d/msgid/cp2k/1be59005-2f09-41c5-a952-d71f5f2526c9n%40googlegroups.com</a>.<br>
</blockquote></div>

<p></p>

-- <br />
You received this message because you are subscribed to the Google Groups "cp2k" group.<br />
To unsubscribe from this group and stop receiving emails from it, send an email to <a href="mailto:cp2k+unsubscribe@googlegroups.com">cp2k+unsubscribe@googlegroups.com</a>.<br />
To view this discussion visit <a href="https://groups.google.com/d/msgid/cp2k/CAKRVyRsrAdzoZaMisQZ-8kwN8_kS_qZxUST0c%3DtB6wCyv3OZUQ%40mail.gmail.com?utm_medium=email&utm_source=footer">https://groups.google.com/d/msgid/cp2k/CAKRVyRsrAdzoZaMisQZ-8kwN8_kS_qZxUST0c%3DtB6wCyv3OZUQ%40mail.gmail.com</a>.<br />