~ubuntu-branches/ubuntu/saucy/nwchem/saucy

« back to all changes in this revision

Viewing changes to src/nwdft/scf_dft_cg/dft_roks_hessv.F

  • Committer: Package Import Robot
  • Author(s): Michael Banck, Michael Banck, Daniel Leidert
  • Date: 2012-02-09 20:02:41 UTC
  • mfrom: (1.1.1)
  • Revision ID: package-import@ubuntu.com-20120209200241-jgk03qfsphal4ug2
Tags: 6.1-1
* New upstream release.

[ Michael Banck ]
* debian/patches/02_makefile_flags.patch: Updated.
* debian/patches/02_makefile_flags.patch: Use internal blas and lapack code.
* debian/patches/02_makefile_flags.patch: Define GCC4 for LINUX and LINUX64
  (Closes: #632611 and LP: #791308).
* debian/control (Build-Depends): Added openssh-client.
* debian/rules (USE_SCALAPACK, SCALAPACK): Removed variables (Closes:
  #654658).
* debian/rules (LIBDIR, USE_MPIF4, ARMCI_NETWORK): New variables.
* debian/TODO: New file.
* debian/control (Build-Depends): Removed libblas-dev, liblapack-dev and
  libscalapack-mpi-dev.
* debian/patches/04_show_testsuite_diff_output.patch: New patch, shows the
  diff output for failed tests.
* debian/patches/series: Adjusted.
* debian/testsuite: Optionally run all tests if "all" is passed as option.
* debian/rules: Run debian/testsuite with "all" if DEB_BUILD_OPTIONS
  contains "checkall".

[ Daniel Leidert ]
* debian/control: Used wrap-and-sort. Added Vcs-Svn and Vcs-Browser fields.
  (Priority): Moved to extra according to policy section 2.5.
  (Standards-Version): Bumped to 3.9.2.
  (Description): Fixed a typo.
* debian/watch: Added.
* debian/patches/03_hurd-i386_define_path_max.patch: Added.
  - Define MAX_PATH if not defines to fix FTBFS on hurd.
* debian/patches/series: Adjusted.

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
      subroutine dft_roks_hessv( acc, g_x, g_ax )
 
2
      implicit none
 
3
#include "errquit.fh"
 
4
#include "crohf.fh"
 
5
#include "cscf.fh"
 
6
c     
 
7
c     $Id: dft_roks_hessv.F 21176 2011-10-10 06:35:49Z d3y133 $
 
8
c
 
9
      integer g_x, g_ax
 
10
      double precision acc
 
11
c     
 
12
      integer gtype,grow,gcol,growp,gcolp
 
13
c     
 
14
c     Check
 
15
c     
 
16
      if (crohf_init_flag.ne.1)
 
17
     $  call errquit('dft_roks_hessv: ROKS internal block invalid',0,
 
18
     &       UNKNOWN_ERR)
 
19
      call ga_inquire(g_x,gtype,grow,gcol)
 
20
      if (grow.ne.crohf_vlen)
 
21
     $  call errquit('dft_roks_hessv: invalid vector length',0,
 
22
     &       UNKNOWN_ERR)
 
23
      call ga_inquire(g_ax,gtype,growp,gcolp)
 
24
      if (growp.ne.crohf_vlen)
 
25
     $  call errquit('dft_roks_hessv: invalid vector length',0,
 
26
     &       UNKNOWN_ERR)
 
27
      if (gcol.ne.gcolp)
 
28
     $  call errquit('dft_roks_hessv: invalid no. of vectors',0,
 
29
     &       UNKNOWN_ERR)
 
30
c     
 
31
c     Call internal routine
 
32
c     
 
33
      call dft_roks_hessv_xx( basis, geom, nbf, nmo,
 
34
     $     nclosed, nopen,
 
35
     $     pflg, g_movecs, oskel, 
 
36
     $     crohf_g_fcv, crohf_g_fpv, crohf_g_fcp,
 
37
     $     acc, lshift, g_x, g_ax )
 
38
c
 
39
c     Zap numbers much smaller than acc to ensure hard zeroes 
 
40
c     remain unpolluted ... cannot use a threshold larger than the
 
41
c     integral accuracy since can break symmetry in non-abelian groups
 
42
c     Also must ensure that the threshold tends to zero to permit
 
43
c     tight convergence.
 
44
c
 
45
      call ga_screen(g_ax, min(acc*acc,acc*0.01d0,1d-12))
 
46
c
 
47
      end