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

« back to all changes in this revision

Viewing changes to src/tools/ga-5-2/global/examples/boltzmann/initpar.F

  • 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
#if HAVE_CONFIG_H
 
2
#   include "config.fh"
 
3
#endif
 
4
      subroutine initpar
 
5
#include "common"
 
6
c
 
7
      double precision q8_ei(18,3)
 
8
      data q8_ei
 
9
     + /1,-1, 0, 0, 1,-1, 1,-1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
 
10
     +  0, 0, 1,-1, 1,-1,-1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
 
11
     +  0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0/
 
12
      integer i,j
 
13
      double precision rdim,b,c2
 
14
c
 
15
      do i=1, 8
 
16
        do j = 1, NDIM
 
17
          ei(i,j) = q8_ei(i,j)
 
18
        end do
 
19
      end do
 
20
c
 
21
      rdim = 4.0d00
 
22
      c2 = cspd**2
 
23
      b = 24.0d00
 
24
c
 
25
      do i =1, 9
 
26
        if (i.eq.1) then
 
27
          ffb(i) = 0.0d00
 
28
          ffc(i) = 0.0d00
 
29
          ffd(i) = -1.0d00/c2-4.0d00*rdim/(2.0d00*b*c2)
 
30
        elseif (i.le.5) then
 
31
          ffb(i) = 4.0d00*rdim/(b*c2)
 
32
          ffc(i) = 4.0d00*rdim*(rdim+2.0d00)/(2.0d00*b*c2**2)
 
33
          ffd(i) = -4.0d00*rdim/(2.0d00*b*c2)
 
34
        else
 
35
          ffb(i) = rdim/(b*c2)
 
36
          ffc(i) = rdim*(rdim+2.0d00)/(2.0d00*b*c2**2)
 
37
          ffd(i) = -rdim/(2.0d00*b*c2)
 
38
        endif
 
39
      end do
 
40
c
 
41
c
 
42
c   Initialize hash table
 
43
c
 
44
      hash(0,0)   = 1
 
45
      hash(1,0)   = 2
 
46
      hash(-1,0)  = 3
 
47
      hash(0,1)   = 4
 
48
      hash(0,-1)  = 5
 
49
      hash(1,1)   = 6
 
50
      hash(-1,-1) = 7
 
51
      hash(1,-1)  = 8
 
52
      hash(-1,1)  = 9
 
53
c
 
54
c   Initialise inverse hash table
 
55
c
 
56
      ihash(0,0)   = 1
 
57
      ihash(1,0)   = 3
 
58
      ihash(-1,0)  = 2
 
59
      ihash(0,1)   = 5
 
60
      ihash(0,-1)  = 4
 
61
      ihash(1,1)   = 7
 
62
      ihash(-1,-1) = 6
 
63
      ihash(1,-1)  = 9
 
64
      ihash(-1,1)  = 8
 
65
c
 
66
      do i = 1, 50
 
67
        tstats(i) = 0.0d00
 
68
      end do
 
69
c
 
70
      return
 
71
      end