← Back to branch summary

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

« back to all changes in this revision

Viewing changes to src/tce/mrcc/create_mr_evl_sorted.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
http://bugs.debian.org/696361
* 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:
src/tce/mrcc/create_mr_evl_sorted.F
 
1
       subroutine create_mr_evl_sorted(k_f1_offsetm)
 
2
        implicit none
 
3
#include "tce.fh"
 
4
#include "mafdecls.fh"
 
5
#include "stdio.fh"
 
6
#include "rtdb.fh"
 
7
#include "errquit.fh"
 
8
#include "sym.fh"
 
9
#include "tce_mrcc.fh"
 
10
#include "tce_main.fh"
 
11
#include "global.fh"
 
12
 
 
13
 
 
14
 
 
15
      integer k_f1_offsetm(maxref)
 
16
c      integer nref,iref
 
17
      integer iref
 
18
c      logical nodezero 
 
19
      integer g6b,g5b
 
20
      integer i,j,k,l
 
21
      integer l_f1,k_f1
 
22
      integer size
 
23
      integer offsetf
 
24
 
 
25
c      nodezero = (ga_nodeid().eq.0)
 
26
 
 
27
      do iref=1,nref
 
28
 
 
29
 
 
30
         k = 0
 
31
         k_sym = k_symm(iref)
 
32
         k_offset = k_offsetm(iref)
 
33
         k_range = k_rangem(iref)
 
34
         k_spin = k_spinm(iref)
 
35
         k_movecs_sorted = k_movecs_sortedm(iref)
 
36
 
 
37
         noa = nblcks(1,iref)
 
38
         nob = nblcks(2,iref)
 
39
         nva = nblcks(3,iref)
 
40
         nvb = nblcks(4,iref)
 
41
 
 
42
         noab = noa+nob
 
43
         nvab = nva+nvb
 
44
 
 
45
      DO g6b = 1,noab+nvab
 
46
      g5b = g6b
 
47
      IF (int_mb(k_spin+g6b-1) .eq. int_mb(k_spin+g5b-1)) THEN
 
48
      IF (ieor(int_mb(k_sym+g6b-1),int_mb(k_sym+g5b-1)) .eq. irrep_f) TH
 
49
     &EN
 
50
      IF ((.not.restricted).or.(int_mb(k_spin+g6b-1)+int_mb(k_spin+g5b-1
 
51
     &).ne.4)) THEN
 
52
 
 
53
        size = int_mb(k_range+g6b-1) * int_mb(k_range+g5b-1)
 
54
 
 
55
        if (.not.ma_push_get(mt_dbl,size,'f1mi',l_f1,k_f1))
 
56
     1   call errquit('tce_mrcc_iface_f1: MA problem get',0,MA_ERR)
 
57
 
 
58
 
 
59
        call get_hash_block(d_f1m(iref),dbl_mb(k_f1),size,
 
60
     1   int_mb(k_f1_offsetm(iref)),g6b-1+(noab+nvab)*(g5b-1))
 
61
        l = 0
 
62
        do i=1,int_mb(k_range+g6b-1)
 
63
        do j=1,int_mb(k_range+g5b-1)
 
64
        l = l + 1
 
65
           if(i.eq.j) then
 
66
             k = k + 1
 
67
             dbl_mb(k_evl_sortedm(iref)+k-1) = dbl_mb(k_f1-1+l)
 
68
           endif
 
69
        enddo
 
70
        enddo
 
71
 
 
72
        if (.not.ma_pop_stack(l_f1))
 
73
     1   call errquit('tce_mrcc_iface_f1: MA problem pop',1,MA_ERR)
 
74
 
 
75
      
 
76
      END IF
 
77
      END IF
 
78
      END IF
 
79
      END DO
 
80
 
 
81
      enddo
 
82
 
 
83
      return
 
84
      end
 
85