[CPMD-list] Different number of iteration optimizing the wavefunction.
Axel Kohlmeyer
akohlmey at vitae.cmm.upenn.edu
Wed Nov 16 00:02:22 CET 2005
On Tue, 15 Nov 2005 rafapa at us.es wrote:
RP> Dear all,
dear rafael,
RP> I have run a calculation in several machines (SGI Altix and Athlon
RP> MP). The final result is the same but the number of iteration s
RP> needed to optimize the wavefunction varies wildly. On the Altix it
RP> takes over 200 cycles whereas on the Athlon takes less than 80. Both
RP> runs where done using 4 CPUs. CPMD was compiled with Intel Fortran.
RP> Is this behaviour normal?
for some systems and using the default DIIS optimizer, unfortunately
yes. DIIS is very fast when it works, but occasionally a bit flaky.
a little numerical noise, different BLAS/LAPACK/FFT implementations
and different (aggressive) compiler optimizations can result
in that. parallel vs. serial and the number of nodes can make
a difference, too. with PCG MINIMIZE it should be more consistent.
e.g. the bulk 63 si example (in the benchmark examples on my
old homepage) is such a notoriously instable case. nowadays,
i prefer using:
ODIIS NO_RESET=100
5
MAXSTEP
50
to get consistent numbers for benchmarking a DIIS step across
platforms (since the example usually does not converge within 50 steps).
best regards,
axel.
RP>
RP> Thanks for the attention.
RP>
RP> Rafael
RP> _______________________________________________
RP> CPMD-list mailing list
RP> CPMD-list at cpmd.org
RP> http://cpmd.org/mailman/listinfo/cpmd-list
RP>
--
=======================================================================
Axel Kohlmeyer akohlmey at cmm.chem.upenn.edu http://www.cmm.upenn.edu
Center for Molecular Modeling -- University of Pennsylvania
Department of Chemistry, 231 S.34th Street, Philadelphia, PA 19104-6323
tel: 1-215-898-1582, fax: 1-215-573-6233, office-tel: 1-215-898-5425
=======================================================================
If you make something idiot-proof, the universe creates a better idiot.
More information about the CPMD-list
mailing list