diff options
Diffstat (limited to 'examples/fortran/ccd.f')
| -rw-r--r-- | examples/fortran/ccd.f | 32 |
1 files changed, 32 insertions, 0 deletions
diff --git a/examples/fortran/ccd.f b/examples/fortran/ccd.f new file mode 100644 index 0000000..12d834d --- /dev/null +++ b/examples/fortran/ccd.f @@ -0,0 +1,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 + |