~roberto-robles/siesta/trunk-RRR

« back to all changes in this revision

Viewing changes to Src/siesta_options.F90

  • Committer: Nick Papior
  • Date: 2016-08-18 05:56:45 UTC
  • Revision ID: nickpapior@gmail.com-20160818055645-judotnqdl75x5y9b
Run all tests, solved merge problems and updated convergence criteria

- Re-run all tests.

  This showed a couple of interesting things.

  1. There was a mistake in the compute_dm code after
  the PEXSI merge when spin-orbit coupling was introduced
  2. I have removed Hprev as it essentially is the same as
  Hold. Either way it shouldn't produce a huge difference
  in the tracking of the dEbs, etc. (after all they are
  not physically used other than for convergence criteria)
  3. The change to SCF.Mix Hamiltonian resulted in a huge
  number of changes in the output. This is because the first
  step prints out the energies at INIT. However, the Hamiltonian
  is different because it is initialized after the compute_dm step.

- Changed the logic in convergence criteria.
  Now the convergence criterias are additive and may be fully controlled.
  However, at least one convergence criteria must be used.

  Now the default convergence criteria is both the Hamiltonian and the density
  matrix.

  This is updated in the manual and the compatibility note.

- Initially I thought the above differences in the energies was
  due to inconsistencies after r538. Hence I have created some
  simple routines in the m_energies.f90 code which updates a selected
  few of the energies.
  I think this should be adopted in the future to ensure that all
  calls to update energies are consistent.
  This will make changes to energy calculations less error-prone.

- Implemented the spin-type in the following routines:
   compute_dm
   final_H_f_stress
   state_init

- Changed the m_compute_max_diff to an interface code with appropriate
  size calculations. There was no reason for explicitly using the
  sparse pattern.

- Fixed a bug in the molecularmechanics code (introduced by Nick r542)

- Added a cyclediffs.sh script which loops on OUT.diffs files and
  it lets one easily cycle the diffs, simply do:

  cd Tests
  make check
  ./cyclediffs.sh

  and answer all the questions. Basically it makes deletes OUT.diffs
  which you have agreed isn't really a change.

- The tests may now be runned via:

   make MPI="mpirun -np 4"
   
  which then uses the default SIESTA location.

  Currently the script checks whether mpirun/mpiexec
  is in SIESTA variable, and if so, does not use MPI variable.
  This makes it easier to decide on the number of cores without
  writing the full path.

Show diffs side-by-side

added added

removed removed

Lines of Context:
128
128
  logical :: atmonly       ! Set up pseudoatom information only?
129
129
  logical :: harrisfun     ! Use Harris functional?
130
130
  logical :: muldeb        ! Write Mulliken polpulations at every SCF step?
131
 
  logical :: require_energy_convergence ! free Energy conv. to finish SCF iteration?
132
 
  logical :: require_harris_convergence ! to finish SCF iteration?
133
 
  logical :: require_hamiltonian_convergence ! to finish SCF iteration?
 
131
  logical :: converge_FreeE    ! free Energy conv. to finish SCF iteration?
 
132
  real(dp) :: tolerance_FreeE  ! Free-energy tolerance
 
133
  logical :: converge_Eharr    ! to finish SCF iteration?
 
134
  real(dp) :: tolerance_Eharr  ! Harris tolerance
 
135
  logical :: converge_DM       ! to finish SCF iteration?
 
136
  real(dp) :: dDtol            ! Tolerance in change of DM elements to finish SCF iteration
 
137
  logical :: converge_H        ! to finish SCF iteration?
 
138
  real(dp) :: dHtol            ! Tolerance in change of H elements to finish SCF iteration
134
139
  logical :: broyden_optim ! Use Broyden method to optimize geometry?
135
140
  logical :: fire_optim    ! Use FIRE method to optimize geometry?
136
141
  logical :: struct_only   ! Output initial structure only?
174
179
  real(dp) :: beta          ! Inverse temperature for Chebishev expansion.
175
180
  real(dp) :: bulkm         ! Bulk modulus
176
181
  real(dp) :: charnet       ! Net electric charge
177
 
  real(dp) :: Energy_tolerance   ! Free-energy tolerance
178
 
  real(dp) :: Harris_tolerance
179
182
  real(dp) :: rijmin        ! Min. permited interatomic distance without warning
180
183
  real(dp) :: dm_normalization_tol    ! Threshold for DM normalization mismatch error
181
184
  logical  :: normalize_dm_during_scf ! Whether we normalize the DM 
182
 
  real(dp) :: dDtol         ! Tolerance in change of DM elements to finish SCF iteration
183
 
  real(dp) :: dHtol         ! Tolerance in change of H elements to finish SCF iteration
184
185
  real(dp) :: dt            ! Time step in dynamics
185
186
  real(dp) :: dx            ! Atomic displacement used to calculate Hessian matrix
186
187
  real(dp) :: dxmax         ! Max. atomic displacement allowed during geom. relaxation