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
|