~ubuntu-branches/ubuntu/trusty/nwchem/trusty-proposed

« back to all changes in this revision

Viewing changes to src/tools/ga-5-1/global/examples/scf/cscf.h

  • Committer: Package Import Robot
  • Author(s): Michael Banck, Daniel Leidert, Andreas Tille, Michael Banck
  • Date: 2013-07-04 12:14:55 UTC
  • mfrom: (1.1.2)
  • Revision ID: package-import@ubuntu.com-20130704121455-5tvsx2qabor3nrui
Tags: 6.3-1
* New upstream release.
* Fixes anisotropic properties (Closes: #696361).
* New features include:
  + Multi-reference coupled cluster (MRCC) approaches
  + Hybrid DFT calculations with short-range HF 
  + New density-functionals including Minnesota (M08, M11) and HSE hybrid
    functionals
  + X-ray absorption spectroscopy (XAS) with TDDFT
  + Analytical gradients for the COSMO solvation model
  + Transition densities from TDDFT 
  + DFT+U and Electron-Transfer (ET) methods for plane wave calculations
  + Exploitation of space group symmetry in plane wave geometry optimizations
  + Local density of states (LDOS) collective variable added to Metadynamics
  + Various new XC functionals added for plane wave calculations, including
    hybrid and range-corrected ones
  + Electric field gradients with relativistic corrections 
  + Nudged Elastic Band optimization method
  + Updated basis sets and ECPs 

[ Daniel Leidert ]
* debian/watch: Fixed.

[ Andreas Tille ]
* debian/upstream: References

[ Michael Banck ]
* debian/upstream (Name): New field.
* debian/patches/02_makefile_flags.patch: Refreshed.
* debian/patches/06_statfs_kfreebsd.patch: Likewise.
* debian/patches/07_ga_target_force_linux.patch: Likewise.
* debian/patches/05_avoid_inline_assembler.patch: Removed, no longer needed.
* debian/patches/09_backported_6.1.1_fixes.patch: Likewise.
* debian/control (Build-Depends): Added gfortran-4.7 and gcc-4.7.
* debian/patches/10_force_gcc-4.7.patch: New patch, explicitly sets
  gfortran-4.7 and gcc-4.7, fixes test suite hang with gcc-4.8 (Closes:
  #701328, #713262).
* debian/testsuite: Added tests for COSMO analytical gradients and MRCC.
* debian/rules (MRCC_METHODS): New variable, required to enable MRCC methods.

Show diffs side-by-side

added added

removed removed

Lines of Context:
1
 
c
2
 
c     include file defining common /cscf/
3
 
c
4
 
c     constant parameters are set in cscf, parameters are set
5
 
c     after reading input.
6
 
c
7
 
c     maxatom  = maximum no. of atoms     (constant parameter)
8
 
c     maxnbfn  = maximum no. of bas. fnct.(constant parameter)
9
 
c     maxnnbfn = maxnbfn*(maxnbfn+1)/2    (constant parameter)
10
 
c     natom    = no. of atoms             (parameter)
11
 
c     nbfn     = no. of basis functions   (parameter)
12
 
c     nnbfn    = nbfn*(nbfn+1)/2          (parameter)
13
 
c     nocc     = no. of occupied orbitals (parameter)
14
 
c     mxiter   = maximim no. of iterations(constant parameter)
15
 
c     tol      = convergence criterion    (constant parameter)
16
 
c     pi       = a familiar constant      (constant parameter)
17
 
c     tol2e    = 2-e integral screening   (constant parameter)
18
 
c
19
 
c     the remainder is initialized in block data or in the
20
 
c     routine ininrm (rnorm and iky)
21
 
c
22
 
c     enrep        = nuclear repulsion energy
23
 
c     q(1:natom)   = nuclear charge of atom
24
 
c     ax(1:natom)  = x co-ordinate of atom
25
 
c     ay(1:natom)  = y ...
26
 
c     az(1:natom)  = z ...
27
 
c     x(1:nbfn)    = x co-ordinate of basis function
28
 
c     y(1:nbfn)    = y ...
29
 
c     z(1:nbfn)    = z ...
30
 
c     expnt(1:nbfn)= exponent of gaussian
31
 
c     rnorm(1:nbfn)= normalization constant of gaussian
32
 
c     iky(1:nbfn)  = iky(i) = i*(i-1)/2 to speed up fock build
33
 
c     icut1        = no. of successful ij   2-e screening tests
34
 
c     icut2        = no. of successful ijkl 2-e screening tests
35
 
c     icut3        = no. of 2-e integrals computed
36
 
c
37
 
      parameter (maxatom = 286)        !cste original value 50
38
 
      parameter (maxnbfn =15*maxatom, mxiter = 30)
39
 
      parameter (maxnnbfn = maxnbfn*(maxnbfn+1)/2)
40
 
      parameter (pi = 3.141592653589793d0)
41
 
      parameter (tol= 0.5d-3)
42
 
      parameter (tol2e=1.0d-6)
43
 
c
44
 
      common /cscf/
45
 
     $     enrep, q(maxatom), ax(maxatom), ay(maxatom), az(maxatom),
46
 
     $     x(maxnbfn), y(maxnbfn), z(maxnbfn), expnt(maxnbfn),
47
 
     $     rnorm(maxnbfn),iky(maxnbfn), icut1, icut2, icut3, icut4,
48
 
     $     natom, nocc, nbfn, nnbfn
49
 
      double precision enrep, q, ax, ay, az, x, y, z, expnt, rnorm
50
 
      integer*8 iky, icut1, icut2, icut3, icut4, natom, nocc, nbfn,
51
 
     $     nnbfn
52
 
c
53
 
c    Global array parameters used in calculations:
54
 
c
55
 
c    ichunk:    chunk size for distributing workload
56
 
c
57
 
c    g_counter: global array used to assign next task
58
 
c    g_dens:    global array used to store density matrix
59
 
c    g_fock:    global array used to store fock matrix
60
 
c    g_tfock:   global array used to store transformed fock matrix
61
 
c    g_schwarz: global array used to store schwarz matrix
62
 
c    g_work:    global array used to store work matrix
63
 
c    g_ident:   global array used to store identity matrix
64
 
c    g_orbs:    global array used to store orbital vectors
65
 
c
66
 
      parameter (ichunk = 20)        !cste original value 10
67
 
      common /g_arrays/ eigv(maxnbfn),
68
 
     $     g_counter, g_dens, g_fock, g_tfock, g_schwarz, g_work,
69
 
     $     g_ident, g_orbs
70
 
      double precision eigv
71
 
      integer g_counter, g_dens, g_fock, g_tfock, g_schwarz, g_work,
72
 
     $     g_ident, g_orbs