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

« back to all changes in this revision

Viewing changes to src/nwdft/rt_tddft/closedshell/zfock_cs_coul.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
C
 
2
C     zfock_cs_coul.F
 
3
C
 
4
C     Calculates 2 electron part of a closed shell complex Fock matrix.
 
5
C     This routine is used when you *only* need to compute the Coulomb
 
6
C     term (i.e., pure DFT, with no Hartree-Fock exchange).  Note that
 
7
C     input is the real part of the density matrix, and output is real
 
8
C     part of the Fock matrix.
 
9
C
 
10
      subroutine zfock_cs_coul (params, g_densre, Ecoul, g_fc)
 
11
      implicit none
 
12
 
 
13
#include "bas.fh"
 
14
#include "errquit.fh"
 
15
#include "mafdecls.fh"
 
16
#include "stdio.fh"
 
17
#include "global.fh"
 
18
#include "util.fh"
 
19
#include "cdft.fh"
 
20
#include "rt_tddft.fh"
 
21
 
 
22
      
 
23
C     == Inputs ==
 
24
      type(rt_params_t), intent(in) :: params
 
25
      integer, intent(in)           :: g_densre
 
26
 
 
27
      
 
28
C     == Outputs ==
 
29
      double precision, intent(out) :: Ecoul
 
30
      integer, intent(in)           :: g_fc
 
31
 
 
32
      
 
33
C     == Parameters ==
 
34
      character(*), parameter     :: pname = "zfock_cs_coul:"
 
35
      
 
36
 
 
37
C     == Variables ==
 
38
      double precision jfac, kfac
 
39
C      double precision tol2e
 
40
      double precision elapsed
 
41
 
 
42
 
 
43
      call rt_tddft_cs_confirm (params)
 
44
 
 
45
      if (params%prof) call prof_start (elapsed)
 
46
 
 
47
      
 
48
C     == Initializations ==
 
49
C      tol2e=10.d0**(-itol2e)
 
50
      Ecoul = 0.d0
 
51
 
 
52
      call ga_zero (g_fc)
 
53
 
 
54
      jfac = 1d0
 
55
      kfac = 0d0
 
56
      
 
57
 
 
58
C     == Calculate Coulomb part of F and energy ==
 
59
      call fock_2e (geom, ao_bas_han, 1, jfac, kfac,
 
60
     $     params%tol2e_re, oskel, g_densre, g_fc, .false.)
 
61
 
 
62
      Ecoul = 0.5d0 * ga_ddot (g_densre, g_fc)
 
63
 
 
64
 
 
65
      if (params%prof) call prof_end (elapsed,"Fock CS standard coul")
 
66
      end subroutine