Commit Graph

72 Commits

Author SHA1 Message Date
dalcorso a8b7dcea58 Added the possibility to make a spin-orbit calculation with
zero magnetization. (Still to be cleaned)
Removed the routine to calculate k - -k needed just in
this case.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@1034 c92efa57-630b-4861-b058-cf58834340f0
2004-07-01 15:42:40 +00:00
deinzer b36b5527d0 Cleanup: 'dispersion' removed from PW, is now entirely in PH
module disp moved from pwcom to phcom
         variable nqx changed to nqmax


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@1032 c92efa57-630b-4861-b058-cf58834340f0
2004-06-30 17:56:05 +00:00
sbraccia 463137ad88 Removed references (now everywhere) to module basis for those variables already in module ions_base.
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@958 c92efa57-630b-4861-b058-cf58834340f0
2004-06-14 07:30:47 +00:00
sbraccia 2058195fda Cleanup of the molecular dynamics routines. Some debug messages removed. Some duplicated variables (fixatom and if_pos) removed from pwcom (already present in ions_base). A logical variable lfixatom added in control_flags. The output of a molecular dynamics run has been slightlychanged to be (hopefully) clearer.
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@951 c92efa57-630b-4861-b058-cf58834340f0
2004-06-11 16:53:16 +00:00
giannozz 25ff6785e3 More uspp variables rearrangement (beware undesired side effects!)
qq_spinorb moved into Modules/uspp.f90 as qq_so
dvan => dvan + dvan_so (for spin-orbit only)
PW/clean_pw.f90 replaces PWNC/clean_pw.f90 (PG)


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@945 c92efa57-630b-4861-b058-cf58834340f0
2004-06-08 20:00:04 +00:00
deinzer 595bda9fe8 First step to calculate all dyn. matrices needed for the entire dispersion
in one step with pw.
A new calculation is possible called 'dispersion'
in the phonon card new input parameters are:
nq1, nq2, nq3 : number of q points in the direction
tr2_ph : threshold for the phonon calculation (see also INPUT_PH)
Therefore tr2_ph has been moved from phcom.f90 to Module disp in pwcom.f90
UP TO NOW ONLY THE Q POINTS FOR THE DYN. MATRICES ARE CALCULATED AND PRINTED.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@944 c92efa57-630b-4861-b058-cf58834340f0
2004-06-08 16:20:41 +00:00
giannozz df51dc0b2a More USPP variables moving around
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@937 c92efa57-630b-4861-b058-cf58834340f0
2004-06-01 21:22:19 +00:00
giannozz c16e3106a9 More tweaking with USPP variables
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@936 c92efa57-630b-4861-b058-cf58834340f0
2004-06-01 16:50:03 +00:00
giannozz 47c5de62f8 More USPP_related variables moved to Modules/uspp.f90
Note that lqx => lmaxq for consistency with other names
(those ending in x are static dimensioning)
Beware unexpected side effects (PG)


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@933 c92efa57-630b-4861-b058-cf58834340f0
2004-05-31 17:55:33 +00:00
giannozz 4a5732790a Part of the variables for ultrasoft pseudopotentials have been moved to a
new module uspp_param, shared between PW and CP (in file Modules/uspp.f90)
The flag newpseudo is no longer used to decide whether states with oc=0 are
bound or not: for those states, and only for the old RRKJ format, oc is set
to a negative value, and this signals that they are not bound. Not elegant
but simpler than it used to be. The converter to UPF format did not include
those states anyway. Beware unexpected side effects!


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@914 c92efa57-630b-4861-b058-cf58834340f0
2004-05-26 11:04:07 +00:00
deinzer f1874de924 Changes in pw to calculate the Raman tensor via the 2n+1 theorem.
A new calculation called 'raman' is now available.
For each k point additional points at k+b will be generated and a nscf
calculation will be performed.
The b vectors are small vectors either in cartesian coordinates or along the
direction along to the nearest neighbours.
A new input card named raman is introduced.
New variables:
b_length:  the length of the b vectors
lcart:     if .true. cartesian direction (default)
           if .false. in direction to the nearest neighbours


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@911 c92efa57-630b-4861-b058-cf58834340f0
2004-05-25 14:22:43 +00:00
sbraccia 12079048da A new check has been added in move_ions in the case of BFGS+LSDA: if when the convergence is achieved the absolute magnetization (computed in electrons and now saved in the module lsda_mod) is less than 1.D-8 a new scf cycle is performed starting from atomic potential and wcfs and the BFGS convergence is checked again. To be tested.
Removed some unused variables from pwscf.
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@905 c92efa57-630b-4861-b058-cf58834340f0
2004-05-21 11:15:39 +00:00
sbraccia 1ab31b8fdc Yet another version of the DIIS algorithm (it substitutes the previous one). This version is slower (but not dramatically) than davidson for small systems, faster for large systems (the cross-over is around 300 bands), but there is space for further improvements. A gamma-point version is also implemented (previously missing). An "holes-sniffer" is implemented to prevent holes in the eigenvalues spectrum.
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@876 c92efa57-630b-4861-b058-cf58834340f0
2004-05-11 08:11:52 +00:00
giannozz 68b075c1ae Cleanup: pwcom no longer used in PP/
qgm is now a local variable as it should be


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@867 c92efa57-630b-4861-b058-cf58834340f0
2004-05-08 10:24:24 +00:00
dalcorso aef6ec462c Introduced variables for spin-orbit support. Changes on these
routines are due to the fact that dvan becomes spin dependent and has
an additional index.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@839 c92efa57-630b-4861-b058-cf58834340f0
2004-04-30 08:51:31 +00:00
giannozz d69e112706 More PP cleanup
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@828 c92efa57-630b-4861-b058-cf58834340f0
2004-04-28 10:25:36 +00:00
giannozz c8c16f1995 ...oops
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@825 c92efa57-630b-4861-b058-cf58834340f0
2004-04-27 15:42:14 +00:00
giannozz 150ae509f3 more pseudopotential cleanup
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@824 c92efa57-630b-4861-b058-cf58834340f0
2004-04-27 13:44:55 +00:00
giannozz 6c76b02482 Pseudopotential cleanup: vnl => vloc_at
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@818 c92efa57-630b-4861-b058-cf58834340f0
2004-04-26 07:25:01 +00:00
giannozz d198ebad0c ndm, mmaxx => ndmx ; 0:ndm => ndmx (starting from 0 no longer needed)
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@815 c92efa57-630b-4861-b058-cf58834340f0
2004-04-23 17:10:44 +00:00
giannozz 18658b544f Cleanup: zv moved to module "basis"
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@809 c92efa57-630b-4861-b058-cf58834340f0
2004-04-21 17:14:38 +00:00
giannozz e277741dc8 First steps in USPP merge between CP and PW:
PW: nqfm => nqfx, some us variables and aainit moved to Modules/uspp.f90 (new)
CP: lx, lqx => lqmax, lix => lmaxx+1, variables moved, aainit not yet merged
invmat moved to flib/ and merged with invmat3 of CP/FPMD,
rndm and ylmr2 moved to flib/ in rranf.f90


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@793 c92efa57-630b-4861-b058-cf58834340f0
2004-04-13 21:49:51 +00:00
degironc 3deee25b0d 2 Apr 2004 A new (experimental) feature is added to LDA+U method.
The +U part of the functional can be built using
              for the projection on localized orbitals some simple
              Wannier functions generated by the post-processing
              code poormanwannier.f90 (pmw.x) from the LDA band
              structure. (SdG)


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@764 c92efa57-630b-4861-b058-cf58834340f0
2004-04-02 16:05:17 +00:00
sbraccia 5f4676dda1 Cleanup: all references to varie replaced by references to control_flags
(also when varie is indirectly called by pwcom).
Fake module varie has been removed from pwcom.
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@761 c92efa57-630b-4861-b058-cf58834340f0
2004-04-02 13:01:15 +00:00
sbraccia 764ce7f663 References to brilz replaced everywhere by cell_base.
There are still references to pwcom: please do not use pwcom in new code !!!
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@760 c92efa57-630b-4861-b058-cf58834340f0
2004-04-02 10:30:15 +00:00
cavazzon c880a0a848 - Almost all neb routines moved to Modules
- New module check_stop used by all codes
  to check for exit conditions ( maximum time
  or EXIT file )


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@718 c92efa57-630b-4861-b058-cf58834340f0
2004-03-15 00:08:26 +00:00
cavazzon 744fb1cbb2 - neb basic routines moved to Modules, to be shared among codes
- basic cell and ionic variables moved frpm pwcom to
  cell_base and ions_base .
  cell_base modules and ions_base modules should be initialized
  through the use of cell_base_init and ions_base_init, this after
  having read the namelists and the cards .
  IN THE FUTURE NONE OF THE VARIABLES IN THE BASIC MODULES SHOULD
  BE SET IN THE MAIN CODES, BUT THROUGH A CALL TO SOME INIT
  FUNCTIONS CONTAINED IN THE BASIC MODULE ITSELF.
  THE IDEA IS THAT THE BASIC VARIABLES SHOULD BE USED AS
  READ-ONLY. UNFORTUNATELY IN F90 THERE IS NO POSSIBILITY
  TO SET A READ ONLY ATTRIBUTE.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@702 c92efa57-630b-4861-b058-cf58834340f0
2004-03-08 08:43:25 +00:00
cavazzon 2b9008ea8a - More merging between CP and FPMD
- some common NEB modules moved to Modules directory
  others NEB modules will be moved soon


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@687 c92efa57-630b-4861-b058-cf58834340f0
2004-02-29 23:30:07 +00:00
profeta ed16640a70 start adding a paw module outside pwcom
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@680 c92efa57-630b-4861-b058-cf58834340f0
2004-02-26 15:33:58 +00:00
profeta 3f541d0d8e Add electric field gradient calculation with PAW reconstruction
move setv and symz to PW
MP


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@672 c92efa57-630b-4861-b058-cf58834340f0
2004-02-25 15:20:28 +00:00
degironc ce9d9804c9 14 feb 2004 LDA+U: variable starting_ns_eigenvalue has been added in
order to allow to specify how to fill starting ns occupation
             matrices. Useful when the default choice take the wrong
             (or not desired) path.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@598 c92efa57-630b-4861-b058-cf58834340f0
2004-02-14 08:39:34 +00:00
sbraccia 4367d42d15 The new bfgs minimization scheme is now the default.
To use the old scheme the keyword ion_dynamics must be set to "old-bfgs".
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@576 c92efa57-630b-4861-b058-cf58834340f0
2004-02-09 11:15:33 +00:00
giannozz 5c6bba2639 Workaround in PW/h_psi.f90 for ifc v.8 (works, only pwcond doesn't)
Manual updated, example6/run_example modified to prevent SP4 weirdness
Empty workspace module removed. Misc cleanup.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@538 c92efa57-630b-4861-b058-cf58834340f0
2004-01-28 18:14:17 +00:00
degironc f10d5b1866 LDA+U forces and stress are now much faster (yet not optimal)
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@517 c92efa57-630b-4861-b058-cf58834340f0
2004-01-24 14:44:20 +00:00
giannozz 0b1a5c9fd8 Extensive module cleanup: DP moved from wrong place (parameters) to the
correct place (kinds); module "varie" replaced by "control_flags" (not
yet in pwcom, though) - many many files changed.
64-bit cpus (Opteron, maybe Itanium) should now work if __LINUX64 is defined


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@513 c92efa57-630b-4861-b058-cf58834340f0
2004-01-23 15:08:03 +00:00
sbraccia b59ad69456 Again some problems with ethr: the diagonalization threshold can now be read by the input file (diago_thr_init).
This value of ethr is overwritten (with a warning message) in the case of a phonon calculation (lphonon flag added). If diago_thr_init is not specified in the input file it is automatically set in setup.f90 (for example in the case of a nscf calculation ethr = tr2 / nelec).
I hope this policy is appropriate for all kind of calculations.
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@511 c92efa57-630b-4861-b058-cf58834340f0
2004-01-23 09:50:00 +00:00
cavazzon bb4db10003 - FPMD prepared for NEB
- variables in module pwcom/varie moved to Modules/control_flags
- module PW/format moved to Modules/format
- in module PW/supercell added workaround for AUTOMATIC variable and
  AIX compilers
- PH Makefile added rbecmod.o


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@468 c92efa57-630b-4861-b058-cf58834340f0
2004-01-07 22:54:02 +00:00
giannozz 5a8c4a05dc Tabs removed (again)
NEB/ removed from makedeps.sh
plotrho: added possibility to have contour lines of various type
A few undefined variables defined (again)


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@465 c92efa57-630b-4861-b058-cf58834340f0
2004-01-05 18:11:01 +00:00
sbraccia 2f4d7baded PW, Gamma and NEB are now unified.
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@449 c92efa57-630b-4861-b058-cf58834340f0
2003-12-10 14:57:07 +00:00
sbraccia 506600d510 Yet another bug in CI_scheme = 'manual' (neb code).
Cleanup.
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@446 c92efa57-630b-4861-b058-cf58834340f0
2003-12-09 08:16:59 +00:00
cavazzon 9ce7263953 - bug (alpha compiler) reported by Sergey has been fixed
- module filnam in pwcom moved to io_files module
- module units in pwcom moved to io_file

From now on, all file names and their unit and attributes (i.e. reclen)
should be defined in Modules/io_files.f90


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@382 c92efa57-630b-4861-b058-cf58834340f0
2003-11-09 18:30:08 +00:00
cavazzon 40e943b6f4 ibm xlf bug workaround related to the "use only" sintax construct
carlo


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@377 c92efa57-630b-4861-b058-cf58834340f0
2003-11-06 15:35:23 +00:00
sbraccia 33921e1f1e iunsdtout removed by PW/pwcom.f90. The unit for standar output (stdout)
is contained in Modules/io_global.f90.
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@362 c92efa57-630b-4861-b058-cf58834340f0
2003-11-03 18:13:21 +00:00
sbraccia 03fe3ac840 Two units have been added (iunstdout and iunneb) to rationalize the output.
iunstdout will substitute 6 everywhere.
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@361 c92efa57-630b-4861-b058-cf58834340f0
2003-11-03 16:00:45 +00:00
sbraccia f25926898e LOGICAL variables lbfgs, lmd, lneb added in MODULE varie. They are
intented to replace iswitch, at least when new code is written.
C.S.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@350 c92efa57-630b-4861-b058-cf58834340f0
2003-10-29 10:46:22 +00:00
giannozz 57d733a8b0 Added SAVE to all modules
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@336 c92efa57-630b-4861-b058-cf58834340f0
2003-10-17 12:03:47 +00:00
giannozz bba38fd761 bfgs modified to work with arbitrary components of the atomic positions
kept fixed (Carlo Sbraccia)
"Total force" now prints a more sensible quantity:
\sqrt(\sum_{i,n} f_i(n)^2) instead of \sum_{i,n} |f_i(n)|
Cleanup in several routines and in modules.


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@334 c92efa57-630b-4861-b058-cf58834340f0
2003-10-16 16:30:12 +00:00
giannozz 9e3fedec2c phonon: added comments to remind why dbecsum is summed
over processors; various calls to setv removed
Second occurrence of wsinit and wsweight removed
Unit number cleanup (solves problem with conflicting units)


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@329 c92efa57-630b-4861-b058-cf58834340f0
2003-10-15 11:40:07 +00:00
cavazzon 1d34b5c0a8 non collinear spin version of PW,
contributed by Adriano Mosca Conte,
added to the repository.
The non collinear code has required
a new module "wavefunctions" containing
"psic" and "evc", that no longer belong
to pwcom.

use: cvs update -d

there's a new example (example13), that is not
yet complete, I'm waiting for the appropriate
pseudopotential from Adriano


git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@322 c92efa57-630b-4861-b058-cf58834340f0
2003-10-03 14:01:11 +00:00
fabris 033a080ddd parallel and more robust version of DIIS diagonalization
git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@297 c92efa57-630b-4861-b058-cf58834340f0
2003-09-02 12:43:13 +00:00