- BGL porting of CP almost complete. New command
line argument "-ntask_groups" to specify the
number of task group used in parallelization.
to generate an executable with BGL extensions,
in the file make.sys the macro __BGL should
be specified.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@3069 c92efa57-630b-4861-b058-cf58834340f0
This allowes for the possibility of setting constraints that are not collective
variables. Example and documentation updated.
C.S.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@3059 c92efa57-630b-4861-b058-cf58834340f0
the bug was introduced with the reorganization of the fft.
- NEB working with CP/FPMD (new calculation fpmd-neb added to check
all the scf engines)
- Clean-ups
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2993 c92efa57-630b-4861-b058-cf58834340f0
density file is not written (perhaps the default could be reverted)
2) added an option ion_velocities='change_step'. Here one restarts from
a file with CG while changing the timestep. The old timestep should be supplied
in tolp. The main goal is to make the [t-1] ion positions to be proportionally
smaller. This way one can go through CP->CG->CP sequence while changing the
timestep and keeping all the ionic and Nose positions & velocities fully
consistent as to keep previous equilibration. Perhaps variable cell variables
should also be rescaled in the same way.
Option 2) is somewhat ugly in how things are specified, but practically useful
if one does not want to re-equilibrate with the time step change.
I could not extract the timestep value from the restart file, is it not
being written?
Kostya
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2930 c92efa57-630b-4861-b058-cf58834340f0
- CP/FPMD post processing updated, charge and states written in xsf/grd format
- CP rhoofr moved to file chargedensity.f90
- merged from_scratch subroutine
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2839 c92efa57-630b-4861-b058-cf58834340f0
an index). Removed the two cases constrained_verlet and constrained_damp: a constrained
dynamics is automatically done if the CONSTRAINTS card is present. Added a keyword to
monitor the value of the constraints during the simulation. In this case the constraints are
not imposed.
C.S.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2740 c92efa57-630b-4861-b058-cf58834340f0
- Many ortho auxiliary functions (tauset, rhoset, sigset, calphi, updatc)
are now in common between FPMD/CP, and moved to module ortho_base.f90
- In FPMD, three index vectors, related to real space like charge and potential
have been substituted with single index vector like in CP, for compatibility
and efficiency.
- Bug fix in pwtools/matdyn.f90 a logical variable was used in place of a
character variable
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2694 c92efa57-630b-4861-b058-cf58834340f0
principle since the same information can be specified with nelup and neldw.
However in systems containing many atoms of various kinds it can help reduce
human errors. If unspecified, this parameter is ignored. [silviu]
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2578 c92efa57-630b-4861-b058-cf58834340f0
could also be specified with a proper choise of nelec. The later option however
is more sensitive to human errors, especially when working on large systems
containing diverse elements (or when working in the middle of the night...) [silviu]
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2577 c92efa57-630b-4861-b058-cf58834340f0
implemented is the one described in PRL 92, 17061 (2004) and, at present, only
works at zero temperature.
The restart files are now written in xml format and are saved ina sub-directory
of the main "save" directory.
C.S.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2562 c92efa57-630b-4861-b058-cf58834340f0
ensemble DFT and NOT conjugate gradient
berry phase electric field and nspin==2
conjugate gradient and nspin==2
P.U.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2465 c92efa57-630b-4861-b058-cf58834340f0
The idea is that module funct contains data defining the DFT functional in use
and a number of functions and subroutines to manage them.
Data are PRIVATE and are (shuold be) accessed and set only by function calls.
Basic drivers to compute XC quantities are also included.
It should be the main place where definitions for and calculation of xc
functionals are centralized.
In CPV printout [A.U.] replaced by [HARTREE A.U.] in many places.
stefano
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2411 c92efa57-630b-4861-b058-cf58834340f0
NEB: removed unused variables; removed allocatable input arrays from input_parameters (to be
consistent with the rest of the input stuff) and replaced by static arrays.
XML: unit 99 used to read/write wfc's replaced by a call to iotk_free_unit.
constraints: the static array "fion" is passed as an array-section to constraint subroutines.
metadyn: fe_step must be defined as an array of "time-steps" (one for each deg-of-freedom).
C.S.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2322 c92efa57-630b-4861-b058-cf58834340f0
vofrho overwrites rhor, what was written in the restart was the
potential and not the charge density.
- clean-up in input_parameters and INPUT_CP doc
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2248 c92efa57-630b-4861-b058-cf58834340f0
accurancy of the pseudopotential table for CP/FPMD (default values = 0.05Ryd)
mmx now is calculated depending on refg and ecutrho
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2155 c92efa57-630b-4861-b058-cf58834340f0
dion, beta, bec everyware.
- subroutines formfn, compute_beta, nlsm1, nlsm2, ecc ... now are common
between FPMD and CPV, a lot of clean ups!
- Changes in stdout: relevant physical quantities ( positions velocities an cell )
are now printed with the seme format of the corresponding input card,
like in PW, as was suggested by SdG.
- exemple23 updated to reflect the new input namelist "wannier"
- Subroutine init_run now is used in FPMD too.
- WARNING in the stress computed with CP, for a pseudo with core-corrections,
a contribution is missing! Not yet fixed, I need to talk with PG for the
box staff.
- WARNING the examples reference are not updated, I'm on the IBM sp, and
I prefer to update them from a linux machine.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2110 c92efa57-630b-4861-b058-cf58834340f0
not specified in the input, like in PW
- FPMD/CP directory for restart file now could be specified with
the input parameter "scradir" in the control namelist
- FPMD: NEB is working again!
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2051 c92efa57-630b-4861-b058-cf58834340f0
Convert electron_dynamics='cg' to tcg=.true. in CP code
In principle, input variable 'tcg' should be removed at some point
Kostya
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2026 c92efa57-630b-4861-b058-cf58834340f0
This patch adds "massive" Nose-Hoover chains for ions (i.e. each ion
can have a separate NH chain attached to it)
Some fixes are still needed in different places:
-the information on the number of NH chains [nhpdim] needs to be
saved and read from the restart file (not done now)
-the NH velocities also need to be all saved [nhpdim*nhpcl]
-an input option needs to be added to zero out the NH velocities during
a restart in order to permit "on the fly" thermostat changes
-deallocation of the module variables gives glibc error with IFC 8.0 & 8.1 in
cpr.f90
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@1998 c92efa57-630b-4861-b058-cf58834340f0
- core charge, same module between CPV and FPMD
- Makefile: added dependency between "mods" and "libiotk"
- Makefile: fix for dependencies in PH and Gamma when __NEWPUNCH is defined
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@1995 c92efa57-630b-4861-b058-cf58834340f0