- PW: do not try to remove the entire prefix.save directory, it doesn't work
for many compilers. Remove the prefix.save/data-file.xml file instead when
starting a calculation from_scratch to prevent restarting from inconsistent
data. Notice that the original content of the prefix.save directory is
destroyed (except for the charge density) in the case of a non-scf or band
calculation as well.
- PW+CP : write rho in binary format; in the LSDA case write rhoup+rhodw in
one file and rhoup-rhodw in 'polarization' file (or should it be called
'magnetization' like in the noncolinear case?). Beware: you cannot read
LSDA files produced by early versions.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@3065 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
states. This bug (introduced at the end of May 2005) was responsible for
somehow wrong results in the case of calculations with many k-points.
Cleanup of diagonalization thresholds: wg_set variable removed and btype array
is now set in sum_band.
Added a logical input variable (diago_full_acc) to decide if the empty states
(defined as states with occupation < 1%) have to be converged at full accuracy
or not. Defaul is .FALSE. (I am not sure this is the best choice, comments
appreciated).
C.S.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@3040 c92efa57-630b-4861-b058-cf58834340f0
Mariginal improvement of the broyden optimisation algorithm. The reference
output of example17 not yet regenerated.
C.S.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@3033 c92efa57-630b-4861-b058-cf58834340f0
LDA+U yields normalized values for projections even if the
atomic wavefunctions are not correctly normalized (Matteo C)
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2942 c92efa57-630b-4861-b058-cf58834340f0
(he claims it is more tested); documentation updated (sort of)
- routines find_free_unit, delete_if_present, moved from 'parser'
to 'io_files' (it is a more natural place)
- routines int_to_char moved from 'parser' to flib/ (it is quite
unrelated with the rest of the module)
- routine trimcheck moved from flib/ to 'io_files' (more natural place)
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2753 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
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
If ibrav=0 ph.x write also at(3,3) and symm_type and so do q2.x and matdyn.x
It could actually be simpler to always write/read at(3,3) and symm_type and
forget about ibrav and celldm ...
latgen modified so that it deals properly with ibrav=0 case.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2543 c92efa57-630b-4861-b058-cf58834340f0
charges are computed in real space only in spheres around the atoms and
these patches are added to rho on the dense fft grid.
Affected routines: addusdens and newd.
Should scale linearly with system size, as opposed to the quadratic
scaling of the G-space implentation.
Experimental version:
- no forces nor stress yet.
- parallel version to be tested
- not to mention phonon and the other codes ...
DISCLAIMER
There are accuracy issues to be understood better, we may need properly
modifies pseudopotentials in order to avoid (or reduce) errors when
ecutrho is still insufficient to describe all Fourier components of the
augmentation charge.
In this case (that is nearly always) also the usal G-space has errors
(negative charges for molecules and surfaces, for instance) but they are
usually small and we are used to live with them.
Some experience is needed to understand the relevance of the error related
to R-space augmentation charges.
CP has R-space augmentation with small-boxes.
They also produce some different, usually small, errors.
I think (I may be wrong) this is the origin of the fact that sometime PWscf
and CP total energies are not exactly the same.
SdG and Antonio Suriano
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2531 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
If absent, it defaults to 'outdir'
Whatever files opened in openfil.f90 will be opened in 'nodedir' (wfc, igk, ...)
Everything else is unaffected
The idea is to make it easier to store large *.wfc files on nodes while keeping
the other files in a central location
The patch is minimally intrusive and only replaces temporarily a variable
in openfil.f90 ( this is an ugly but "clean" hack for now )
Please test if it does what is claimed
Kostya
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2305 c92efa57-630b-4861-b058-cf58834340f0
updated with the BFGS algorithm. The performance is similar to the quasi-newton method,
but the algorithm is much simpler. Moreover it works also with constrained dynamics.
C.S.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2264 c92efa57-630b-4861-b058-cf58834340f0
support in pw.
in call ro routines ccgdiagg cegterg cinitcgg
added ik (integer) k-point considerd
h_epsi_her calculates V_el|Psi>
c_phase_field.f90 calculates berry phase
during a scf calculation
P.U.
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@2098 c92efa57-630b-4861-b058-cf58834340f0