summaryrefslogtreecommitdiff
path: root/examples/fortran/ccd.f
blob: 12d834dafc02edd85b6ba67f0a96f52b13e6d7a3 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
c
c These have been separated out from ccsd_t_singles_l.F and ccsd_t_doubles_l.F
c
      subroutine clean_sd_t_s1_1(h3d,h2d,h1d,p6d,p5d,p4d,
     2                     triplesx,t1sub,v2sub)
      IMPLICIT NONE
      integer h3d,h2d,h1d,p6d,p5d,p4d
      integer h3,h2,h1,p6,p5,p4
      integer N
			double precision triplesx(16,16,16,16,16,16)
      double precision t1sub(16,16)
      double precision v2sub(16,16,16,16)
      
      N = 16       

      do p4=1,10
      do p5=1,10
      do p6=1,10
      do h1=1,10
      do h2=1,10
      do h3=1,10
       triplesx(h3,h2,h1,p6,p5,p4)=triplesx(h3,h2,h1,p6,p5,p4)
     1   + t1sub(p4,h1)*v2sub(h3,h2,p6,p5)
      enddo
      enddo
      enddo
      enddo
      enddo
      enddo
      return
      end