← Back to branch summary

~ubuntu-branches/ubuntu/utopic/nwchem/utopic

« back to all changes in this revision

Viewing changes to src/tce/eomccsd/eomccsd_1prdm_ph.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/eomccsd/eomccsd_1prdm_ph.F
1
1
      SUBROUTINE eomccsd_1prdm_ph(d_i0,d_t1,d_t2,d_x0,d_x1,d_x2,d_y1,d_y
2
2
     &2,k_i0_offset,k_t1_offset,k_t2_offset,k_x0_offset,k_x1_offset,k_x2
3
3
     &_offset,k_y1_offset,k_y2_offset)
4
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
4
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
5
5
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
6
6
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
7
7
C     i0 ( p2 h1 )_yxt + = 1 * Sum ( h4 p3 ) * t ( p2 p3 h1 h4 )_t * i1 ( h4 p3 )_yx
187
187
      END
188
188
      SUBROUTINE eomccsd_1prdm_ph_1(d_a,k_a_offset,d_b,k_b_offset,d_c,k_
189
189
     &c_offset)
190
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
190
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
191
191
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
192
192
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
193
193
C     i0 ( p2 h1 )_yxt + = 1 * Sum ( h4 p3 ) * t ( p2 p3 h1 h4 )_t * i1 ( h4 p3 )_yx
203
203
      INTEGER k_b_offset
204
204
      INTEGER d_c
205
205
      INTEGER k_c_offset
206
 
      INTEGER NXTVAL
 
206
      INTEGER NXTASK
207
207
      INTEGER next
208
208
      INTEGER nprocs
209
209
      INTEGER count
235
235
      INTEGER k_b
236
236
      INTEGER l_c
237
237
      INTEGER k_c
238
 
      EXTERNAL NXTVAL
 
238
      EXTERNAL NXTASK
239
239
      nprocs = GA_NNODES()
240
240
      count = 0
241
 
      next = NXTVAL(nprocs)
 
241
      next = NXTASK(nprocs,1)
242
242
      DO p2b = noab+1,noab+nvab
243
243
      DO h1b = 1,noab
244
244
      IF (next.eq.count) THEN
338
338
      END IF
339
339
      END IF
340
340
      END IF
341
 
      next = NXTVAL(nprocs)
 
341
      next = NXTASK(nprocs,1)
342
342
      END IF
343
343
      count = count + 1
344
344
      END DO
345
345
      END DO
346
 
      next = NXTVAL(-nprocs)
 
346
      next = NXTASK(-nprocs,1)
347
347
      call GA_SYNC()
348
348
      RETURN
349
349
      END
350
350
      SUBROUTINE eomccsd_1prdm_ph_1_1(d_a,k_a_offset,d_b,k_b_offset,d_c,
351
351
     &k_c_offset)
352
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
352
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
353
353
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
354
354
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
355
355
C     i1 ( h4 p3 )_yx + = 1 * x ( )_x * y ( h4 p3 )_y
365
365
      INTEGER k_b_offset
366
366
      INTEGER d_c
367
367
      INTEGER k_c_offset
368
 
      INTEGER NXTVAL
 
368
      INTEGER NXTASK
369
369
      INTEGER next
370
370
      INTEGER nprocs
371
371
      INTEGER count
391
391
      INTEGER k_b
392
392
      INTEGER l_c
393
393
      INTEGER k_c
394
 
      EXTERNAL NXTVAL
 
394
      EXTERNAL NXTASK
395
395
      nprocs = GA_NNODES()
396
396
      count = 0
397
 
      next = NXTVAL(nprocs)
 
397
      next = NXTASK(nprocs,1)
398
398
      DO h4b = 1,noab
399
399
      DO p3b = noab+1,noab+nvab
400
400
      IF (next.eq.count) THEN
455
455
      END IF
456
456
      END IF
457
457
      END IF
458
 
      next = NXTVAL(nprocs)
 
458
      next = NXTASK(nprocs,1)
459
459
      END IF
460
460
      count = count + 1
461
461
      END DO
462
462
      END DO
463
 
      next = NXTVAL(-nprocs)
 
463
      next = NXTASK(-nprocs,1)
464
464
      call GA_SYNC()
465
465
      RETURN
466
466
      END
467
467
      SUBROUTINE OFFSET_eomccsd_1prdm_ph_1_1(l_a_offset,k_a_offset,size)
468
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
468
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
469
469
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
470
470
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
471
471
C     i1 ( h4 p3 )_yx
521
521
      END
522
522
      SUBROUTINE eomccsd_1prdm_ph_1_2(d_a,k_a_offset,d_b,k_b_offset,d_c,
523
523
     &k_c_offset)
524
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
524
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
525
525
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
526
526
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
527
527
C     i1 ( h4 p3 )_yx + = 1 * Sum ( h6 p5 ) * x ( p5 h6 )_x * y ( h4 h6 p3 p5 )_y
537
537
      INTEGER k_b_offset
538
538
      INTEGER d_c
539
539
      INTEGER k_c_offset
540
 
      INTEGER NXTVAL
 
540
      INTEGER NXTASK
541
541
      INTEGER next
542
542
      INTEGER nprocs
543
543
      INTEGER count
569
569
      INTEGER k_b
570
570
      INTEGER l_c
571
571
      INTEGER k_c
572
 
      EXTERNAL NXTVAL
 
572
      EXTERNAL NXTASK
573
573
      nprocs = GA_NNODES()
574
574
      count = 0
575
 
      next = NXTVAL(nprocs)
 
575
      next = NXTASK(nprocs,1)
576
576
      DO h4b = 1,noab
577
577
      DO p3b = noab+1,noab+nvab
578
578
      IF (next.eq.count) THEN
671
671
      END IF
672
672
      END IF
673
673
      END IF
674
 
      next = NXTVAL(nprocs)
 
674
      next = NXTASK(nprocs,1)
675
675
      END IF
676
676
      count = count + 1
677
677
      END DO
678
678
      END DO
679
 
      next = NXTVAL(-nprocs)
 
679
      next = NXTASK(-nprocs,1)
680
680
      call GA_SYNC()
681
681
      RETURN
682
682
      END
683
683
      SUBROUTINE eomccsd_1prdm_ph_2(d_a,k_a_offset,d_b,k_b_offset,d_c,k_
684
684
     &c_offset)
685
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
685
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
686
686
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
687
687
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
688
688
C     i0 ( p2 h1 )_ytxt + = -1 * Sum ( h8 ) * t ( p2 h8 )_t * i1 ( h8 h1 )_ytx
698
698
      INTEGER k_b_offset
699
699
      INTEGER d_c
700
700
      INTEGER k_c_offset
701
 
      INTEGER NXTVAL
 
701
      INTEGER NXTASK
702
702
      INTEGER next
703
703
      INTEGER nprocs
704
704
      INTEGER count
727
727
      INTEGER k_b
728
728
      INTEGER l_c
729
729
      INTEGER k_c
730
 
      EXTERNAL NXTVAL
 
730
      EXTERNAL NXTASK
731
731
      nprocs = GA_NNODES()
732
732
      count = 0
733
 
      next = NXTVAL(nprocs)
 
733
      next = NXTASK(nprocs,1)
734
734
      DO p2b = noab+1,noab+nvab
735
735
      DO h1b = 1,noab
736
736
      IF (next.eq.count) THEN
799
799
      END IF
800
800
      END IF
801
801
      END IF
802
 
      next = NXTVAL(nprocs)
 
802
      next = NXTASK(nprocs,1)
803
803
      END IF
804
804
      count = count + 1
805
805
      END DO
806
806
      END DO
807
 
      next = NXTVAL(-nprocs)
 
807
      next = NXTASK(-nprocs,1)
808
808
      call GA_SYNC()
809
809
      RETURN
810
810
      END
811
811
      SUBROUTINE eomccsd_1prdm_ph_2_1(d_a,k_a_offset,d_b,k_b_offset,d_c,
812
812
     &k_c_offset)
813
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
813
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
814
814
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
815
815
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
816
816
C     i1 ( h8 h1 )_ytx + = 1 * x ( )_x * i2 ( h8 h1 )_yt
826
826
      INTEGER k_b_offset
827
827
      INTEGER d_c
828
828
      INTEGER k_c_offset
829
 
      INTEGER NXTVAL
 
829
      INTEGER NXTASK
830
830
      INTEGER next
831
831
      INTEGER nprocs
832
832
      INTEGER count
852
852
      INTEGER k_b
853
853
      INTEGER l_c
854
854
      INTEGER k_c
855
 
      EXTERNAL NXTVAL
 
855
      EXTERNAL NXTASK
856
856
      nprocs = GA_NNODES()
857
857
      count = 0
858
 
      next = NXTVAL(nprocs)
 
858
      next = NXTASK(nprocs,1)
859
859
      DO h8b = 1,noab
860
860
      DO h1b = 1,noab
861
861
      IF (next.eq.count) THEN
916
916
      END IF
917
917
      END IF
918
918
      END IF
919
 
      next = NXTVAL(nprocs)
 
919
      next = NXTASK(nprocs,1)
920
920
      END IF
921
921
      count = count + 1
922
922
      END DO
923
923
      END DO
924
 
      next = NXTVAL(-nprocs)
 
924
      next = NXTASK(-nprocs,1)
925
925
      call GA_SYNC()
926
926
      RETURN
927
927
      END
928
928
      SUBROUTINE OFFSET_eomccsd_1prdm_ph_2_1(l_a_offset,k_a_offset,size)
929
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
929
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
930
930
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
931
931
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
932
932
C     i1 ( h8 h1 )_ytx
982
982
      END
983
983
      SUBROUTINE eomccsd_1prdm_ph_2_1_1(d_a,k_a_offset,d_b,k_b_offset,d_
984
984
     &c,k_c_offset)
985
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
985
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
986
986
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
987
987
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
988
988
C     i2 ( h8 h1 )_yt + = 1 * Sum ( p3 ) * t ( p3 h1 )_t * y ( h8 p3 )_y
998
998
      INTEGER k_b_offset
999
999
      INTEGER d_c
1000
1000
      INTEGER k_c_offset
1001
 
      INTEGER NXTVAL
 
1001
      INTEGER NXTASK
1002
1002
      INTEGER next
1003
1003
      INTEGER nprocs
1004
1004
      INTEGER count
1027
1027
      INTEGER k_b
1028
1028
      INTEGER l_c
1029
1029
      INTEGER k_c
1030
 
      EXTERNAL NXTVAL
 
1030
      EXTERNAL NXTASK
1031
1031
      nprocs = GA_NNODES()
1032
1032
      count = 0
1033
 
      next = NXTVAL(nprocs)
 
1033
      next = NXTASK(nprocs,1)
1034
1034
      DO h8b = 1,noab
1035
1035
      DO h1b = 1,noab
1036
1036
      IF (next.eq.count) THEN
1099
1099
      END IF
1100
1100
      END IF
1101
1101
      END IF
1102
 
      next = NXTVAL(nprocs)
 
1102
      next = NXTASK(nprocs,1)
1103
1103
      END IF
1104
1104
      count = count + 1
1105
1105
      END DO
1106
1106
      END DO
1107
 
      next = NXTVAL(-nprocs)
 
1107
      next = NXTASK(-nprocs,1)
1108
1108
      call GA_SYNC()
1109
1109
      RETURN
1110
1110
      END
1111
1111
      SUBROUTINE OFFSET_eomccsd_1prdm_ph_2_1_1(l_a_offset,k_a_offset,siz
1112
1112
     &e)
1113
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
1113
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
1114
1114
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
1115
1115
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
1116
1116
C     i2 ( h8 h1 )_yt
1166
1166
      END
1167
1167
      SUBROUTINE eomccsd_1prdm_ph_2_1_2(d_a,k_a_offset,d_b,k_b_offset,d_
1168
1168
     &c,k_c_offset)
1169
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
1169
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
1170
1170
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
1171
1171
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
1172
1172
C     i2 ( h8 h1 )_yt + = -1/2 * Sum ( h5 p3 p4 ) * t ( p3 p4 h1 h5 )_t * y ( h5 h8 p3 p4 )_y
1182
1182
      INTEGER k_b_offset
1183
1183
      INTEGER d_c
1184
1184
      INTEGER k_c_offset
1185
 
      INTEGER NXTVAL
 
1185
      INTEGER NXTASK
1186
1186
      INTEGER next
1187
1187
      INTEGER nprocs
1188
1188
      INTEGER count
1220
1220
      INTEGER l_c
1221
1221
      INTEGER k_c
1222
1222
      DOUBLE PRECISION FACTORIAL
1223
 
      EXTERNAL NXTVAL
 
1223
      EXTERNAL NXTASK
1224
1224
      EXTERNAL FACTORIAL
1225
1225
      nprocs = GA_NNODES()
1226
1226
      count = 0
1227
 
      next = NXTVAL(nprocs)
 
1227
      next = NXTASK(nprocs,1)
1228
1228
      DO h8b = 1,noab
1229
1229
      DO h1b = 1,noab
1230
1230
      IF (next.eq.count) THEN
1331
1331
      END IF
1332
1332
      END IF
1333
1333
      END IF
1334
 
      next = NXTVAL(nprocs)
 
1334
      next = NXTASK(nprocs,1)
1335
1335
      END IF
1336
1336
      count = count + 1
1337
1337
      END DO
1338
1338
      END DO
1339
 
      next = NXTVAL(-nprocs)
 
1339
      next = NXTASK(-nprocs,1)
1340
1340
      call GA_SYNC()
1341
1341
      RETURN
1342
1342
      END
1343
1343
      SUBROUTINE eomccsd_1prdm_ph_2_2(d_a,k_a_offset,d_b,k_b_offset,d_c,
1344
1344
     &k_c_offset)
1345
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
1345
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
1346
1346
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
1347
1347
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
1348
1348
C     i1 ( h8 h1 )_yx + = 1 * Sum ( p4 ) * x ( p4 h1 )_x * y ( h8 p4 )_y
1358
1358
      INTEGER k_b_offset
1359
1359
      INTEGER d_c
1360
1360
      INTEGER k_c_offset
1361
 
      INTEGER NXTVAL
 
1361
      INTEGER NXTASK
1362
1362
      INTEGER next
1363
1363
      INTEGER nprocs
1364
1364
      INTEGER count
1387
1387
      INTEGER k_b
1388
1388
      INTEGER l_c
1389
1389
      INTEGER k_c
1390
 
      EXTERNAL NXTVAL
 
1390
      EXTERNAL NXTASK
1391
1391
      nprocs = GA_NNODES()
1392
1392
      count = 0
1393
 
      next = NXTVAL(nprocs)
 
1393
      next = NXTASK(nprocs,1)
1394
1394
      DO h8b = 1,noab
1395
1395
      DO h1b = 1,noab
1396
1396
      IF (next.eq.count) THEN
1459
1459
      END IF
1460
1460
      END IF
1461
1461
      END IF
1462
 
      next = NXTVAL(nprocs)
 
1462
      next = NXTASK(nprocs,1)
1463
1463
      END IF
1464
1464
      count = count + 1
1465
1465
      END DO
1466
1466
      END DO
1467
 
      next = NXTVAL(-nprocs)
 
1467
      next = NXTASK(-nprocs,1)
1468
1468
      call GA_SYNC()
1469
1469
      RETURN
1470
1470
      END
1471
1471
      SUBROUTINE eomccsd_1prdm_ph_2_3(d_a,k_a_offset,d_b,k_b_offset,d_c,
1472
1472
     &k_c_offset)
1473
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
1473
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
1474
1474
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
1475
1475
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
1476
1476
C     i1 ( h8 h1 )_yxt + = 1 * Sum ( p3 ) * t ( p3 h1 )_t * i2 ( h8 p3 )_yx
1486
1486
      INTEGER k_b_offset
1487
1487
      INTEGER d_c
1488
1488
      INTEGER k_c_offset
1489
 
      INTEGER NXTVAL
 
1489
      INTEGER NXTASK
1490
1490
      INTEGER next
1491
1491
      INTEGER nprocs
1492
1492
      INTEGER count
1515
1515
      INTEGER k_b
1516
1516
      INTEGER l_c
1517
1517
      INTEGER k_c
1518
 
      EXTERNAL NXTVAL
 
1518
      EXTERNAL NXTASK
1519
1519
      nprocs = GA_NNODES()
1520
1520
      count = 0
1521
 
      next = NXTVAL(nprocs)
 
1521
      next = NXTASK(nprocs,1)
1522
1522
      DO h8b = 1,noab
1523
1523
      DO h1b = 1,noab
1524
1524
      IF (next.eq.count) THEN
1587
1587
      END IF
1588
1588
      END IF
1589
1589
      END IF
1590
 
      next = NXTVAL(nprocs)
 
1590
      next = NXTASK(nprocs,1)
1591
1591
      END IF
1592
1592
      count = count + 1
1593
1593
      END DO
1594
1594
      END DO
1595
 
      next = NXTVAL(-nprocs)
 
1595
      next = NXTASK(-nprocs,1)
1596
1596
      call GA_SYNC()
1597
1597
      RETURN
1598
1598
      END
1599
1599
      SUBROUTINE eomccsd_1prdm_ph_2_3_1(d_a,k_a_offset,d_b,k_b_offset,d_
1600
1600
     &c,k_c_offset)
1601
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
1601
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
1602
1602
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
1603
1603
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
1604
1604
C     i2 ( h8 p3 )_yx + = -1 * Sum ( h6 p5 ) * x ( p5 h6 )_x * y ( h6 h8 p3 p5 )_y
1614
1614
      INTEGER k_b_offset
1615
1615
      INTEGER d_c
1616
1616
      INTEGER k_c_offset
1617
 
      INTEGER NXTVAL
 
1617
      INTEGER NXTASK
1618
1618
      INTEGER next
1619
1619
      INTEGER nprocs
1620
1620
      INTEGER count
1646
1646
      INTEGER k_b
1647
1647
      INTEGER l_c
1648
1648
      INTEGER k_c
1649
 
      EXTERNAL NXTVAL
 
1649
      EXTERNAL NXTASK
1650
1650
      nprocs = GA_NNODES()
1651
1651
      count = 0
1652
 
      next = NXTVAL(nprocs)
 
1652
      next = NXTASK(nprocs,1)
1653
1653
      DO h8b = 1,noab
1654
1654
      DO p3b = noab+1,noab+nvab
1655
1655
      IF (next.eq.count) THEN
1748
1748
      END IF
1749
1749
      END IF
1750
1750
      END IF
1751
 
      next = NXTVAL(nprocs)
 
1751
      next = NXTASK(nprocs,1)
1752
1752
      END IF
1753
1753
      count = count + 1
1754
1754
      END DO
1755
1755
      END DO
1756
 
      next = NXTVAL(-nprocs)
 
1756
      next = NXTASK(-nprocs,1)
1757
1757
      call GA_SYNC()
1758
1758
      RETURN
1759
1759
      END
1760
1760
      SUBROUTINE OFFSET_eomccsd_1prdm_ph_2_3_1(l_a_offset,k_a_offset,siz
1761
1761
     &e)
1762
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
1762
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
1763
1763
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
1764
1764
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
1765
1765
C     i2 ( h8 p3 )_yx
1815
1815
      END
1816
1816
      SUBROUTINE eomccsd_1prdm_ph_2_4(d_a,k_a_offset,d_b,k_b_offset,d_c,
1817
1817
     &k_c_offset)
1818
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
1818
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
1819
1819
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
1820
1820
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
1821
1821
C     i1 ( h8 h1 )_yx + = -1/2 * Sum ( h6 p4 p5 ) * x ( p4 p5 h1 h6 )_x * y ( h6 h8 p4 p5 )_y
1831
1831
      INTEGER k_b_offset
1832
1832
      INTEGER d_c
1833
1833
      INTEGER k_c_offset
1834
 
      INTEGER NXTVAL
 
1834
      INTEGER NXTASK
1835
1835
      INTEGER next
1836
1836
      INTEGER nprocs
1837
1837
      INTEGER count
1869
1869
      INTEGER l_c
1870
1870
      INTEGER k_c
1871
1871
      DOUBLE PRECISION FACTORIAL
1872
 
      EXTERNAL NXTVAL
 
1872
      EXTERNAL NXTASK
1873
1873
      EXTERNAL FACTORIAL
1874
1874
      nprocs = GA_NNODES()
1875
1875
      count = 0
1876
 
      next = NXTVAL(nprocs)
 
1876
      next = NXTASK(nprocs,1)
1877
1877
      DO h8b = 1,noab
1878
1878
      DO h1b = 1,noab
1879
1879
      IF (next.eq.count) THEN
1980
1980
      END IF
1981
1981
      END IF
1982
1982
      END IF
1983
 
      next = NXTVAL(nprocs)
 
1983
      next = NXTASK(nprocs,1)
1984
1984
      END IF
1985
1985
      count = count + 1
1986
1986
      END DO
1987
1987
      END DO
1988
 
      next = NXTVAL(-nprocs)
 
1988
      next = NXTASK(-nprocs,1)
1989
1989
      call GA_SYNC()
1990
1990
      RETURN
1991
1991
      END
1992
1992
      SUBROUTINE eomccsd_1prdm_ph_3(d_a,k_a_offset,d_b,k_b_offset,d_c,k_
1993
1993
     &c_offset)
1994
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
1994
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
1995
1995
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
1996
1996
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
1997
1997
C     i0 ( p2 h1 )_yttx + = 1/2 * x ( )_x * i1 ( p2 h1 )_ytt
2007
2007
      INTEGER k_b_offset
2008
2008
      INTEGER d_c
2009
2009
      INTEGER k_c_offset
2010
 
      INTEGER NXTVAL
 
2010
      INTEGER NXTASK
2011
2011
      INTEGER next
2012
2012
      INTEGER nprocs
2013
2013
      INTEGER count
2033
2033
      INTEGER k_b
2034
2034
      INTEGER l_c
2035
2035
      INTEGER k_c
2036
 
      EXTERNAL NXTVAL
 
2036
      EXTERNAL NXTASK
2037
2037
      nprocs = GA_NNODES()
2038
2038
      count = 0
2039
 
      next = NXTVAL(nprocs)
 
2039
      next = NXTASK(nprocs,1)
2040
2040
      DO p2b = noab+1,noab+nvab
2041
2041
      DO h1b = 1,noab
2042
2042
      IF (next.eq.count) THEN
2097
2097
      END IF
2098
2098
      END IF
2099
2099
      END IF
2100
 
      next = NXTVAL(nprocs)
 
2100
      next = NXTASK(nprocs,1)
2101
2101
      END IF
2102
2102
      count = count + 1
2103
2103
      END DO
2104
2104
      END DO
2105
 
      next = NXTVAL(-nprocs)
 
2105
      next = NXTASK(-nprocs,1)
2106
2106
      call GA_SYNC()
2107
2107
      RETURN
2108
2108
      END
2109
2109
      SUBROUTINE eomccsd_1prdm_ph_3_1(d_a,k_a_offset,d_b,k_b_offset,d_c,
2110
2110
     &k_c_offset)
2111
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
2111
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
2112
2112
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
2113
2113
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
2114
2114
C     i1 ( p2 h1 )_ytt + = -1 * Sum ( h4 h5 p3 ) * t ( p2 p3 h4 h5 )_t * i2 ( h4 h5 h1 p3 )_yt
2124
2124
      INTEGER k_b_offset
2125
2125
      INTEGER d_c
2126
2126
      INTEGER k_c_offset
2127
 
      INTEGER NXTVAL
 
2127
      INTEGER NXTASK
2128
2128
      INTEGER next
2129
2129
      INTEGER nprocs
2130
2130
      INTEGER count
2162
2162
      INTEGER l_c
2163
2163
      INTEGER k_c
2164
2164
      DOUBLE PRECISION FACTORIAL
2165
 
      EXTERNAL NXTVAL
 
2165
      EXTERNAL NXTASK
2166
2166
      EXTERNAL FACTORIAL
2167
2167
      nprocs = GA_NNODES()
2168
2168
      count = 0
2169
 
      next = NXTVAL(nprocs)
 
2169
      next = NXTASK(nprocs,1)
2170
2170
      DO p2b = noab+1,noab+nvab
2171
2171
      DO h1b = 1,noab
2172
2172
      IF (next.eq.count) THEN
2263
2263
      END IF
2264
2264
      END IF
2265
2265
      END IF
2266
 
      next = NXTVAL(nprocs)
 
2266
      next = NXTASK(nprocs,1)
2267
2267
      END IF
2268
2268
      count = count + 1
2269
2269
      END DO
2270
2270
      END DO
2271
 
      next = NXTVAL(-nprocs)
 
2271
      next = NXTASK(-nprocs,1)
2272
2272
      call GA_SYNC()
2273
2273
      RETURN
2274
2274
      END
2275
2275
      SUBROUTINE OFFSET_eomccsd_1prdm_ph_3_1(l_a_offset,k_a_offset,size)
2276
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
2276
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
2277
2277
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
2278
2278
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
2279
2279
C     i1 ( p2 h1 )_ytt
2329
2329
      END
2330
2330
      SUBROUTINE eomccsd_1prdm_ph_3_1_1(d_a,k_a_offset,d_b,k_b_offset,d_
2331
2331
     &c,k_c_offset)
2332
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
2332
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
2333
2333
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
2334
2334
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
2335
2335
C     i2 ( h4 h5 h1 p3 )_yt + = -1 * Sum ( p6 ) * t ( p6 h1 )_t * y ( h4 h5 p3 p6 )_y
2345
2345
      INTEGER k_b_offset
2346
2346
      INTEGER d_c
2347
2347
      INTEGER k_c_offset
2348
 
      INTEGER NXTVAL
 
2348
      INTEGER NXTASK
2349
2349
      INTEGER next
2350
2350
      INTEGER nprocs
2351
2351
      INTEGER count
2378
2378
      INTEGER k_b
2379
2379
      INTEGER l_c
2380
2380
      INTEGER k_c
2381
 
      EXTERNAL NXTVAL
 
2381
      EXTERNAL NXTASK
2382
2382
      nprocs = GA_NNODES()
2383
2383
      count = 0
2384
 
      next = NXTVAL(nprocs)
 
2384
      next = NXTASK(nprocs,1)
2385
2385
      DO h4b = 1,noab
2386
2386
      DO h5b = h4b,noab
2387
2387
      DO h1b = 1,noab
2470
2470
      END IF
2471
2471
      END IF
2472
2472
      END IF
2473
 
      next = NXTVAL(nprocs)
 
2473
      next = NXTASK(nprocs,1)
2474
2474
      END IF
2475
2475
      count = count + 1
2476
2476
      END DO
2477
2477
      END DO
2478
2478
      END DO
2479
2479
      END DO
2480
 
      next = NXTVAL(-nprocs)
 
2480
      next = NXTASK(-nprocs,1)
2481
2481
      call GA_SYNC()
2482
2482
      RETURN
2483
2483
      END
2484
2484
      SUBROUTINE OFFSET_eomccsd_1prdm_ph_3_1_1(l_a_offset,k_a_offset,siz
2485
2485
     &e)
2486
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
2486
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
2487
2487
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
2488
2488
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
2489
2489
C     i2 ( h4 h5 h1 p3 )_yt
2555
2555
      END
2556
2556
      SUBROUTINE eomccsd_1prdm_ph_4(d_a,k_a_offset,d_b,k_b_offset,d_c,k_
2557
2557
     &c_offset)
2558
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
2558
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
2559
2559
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
2560
2560
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
2561
2561
C     i0 ( p2 h1 )_ytx + = -1 * Sum ( h7 ) * x ( p2 h7 )_x * i1 ( h7 h1 )_yt
2571
2571
      INTEGER k_b_offset
2572
2572
      INTEGER d_c
2573
2573
      INTEGER k_c_offset
2574
 
      INTEGER NXTVAL
 
2574
      INTEGER NXTASK
2575
2575
      INTEGER next
2576
2576
      INTEGER nprocs
2577
2577
      INTEGER count
2600
2600
      INTEGER k_b
2601
2601
      INTEGER l_c
2602
2602
      INTEGER k_c
2603
 
      EXTERNAL NXTVAL
 
2603
      EXTERNAL NXTASK
2604
2604
      nprocs = GA_NNODES()
2605
2605
      count = 0
2606
 
      next = NXTVAL(nprocs)
 
2606
      next = NXTASK(nprocs,1)
2607
2607
      DO p2b = noab+1,noab+nvab
2608
2608
      DO h1b = 1,noab
2609
2609
      IF (next.eq.count) THEN
2672
2672
      END IF
2673
2673
      END IF
2674
2674
      END IF
2675
 
      next = NXTVAL(nprocs)
 
2675
      next = NXTASK(nprocs,1)
2676
2676
      END IF
2677
2677
      count = count + 1
2678
2678
      END DO
2679
2679
      END DO
2680
 
      next = NXTVAL(-nprocs)
 
2680
      next = NXTASK(-nprocs,1)
2681
2681
      call GA_SYNC()
2682
2682
      RETURN
2683
2683
      END
2684
2684
      SUBROUTINE eomccsd_1prdm_ph_4_1(d_a,k_a_offset,d_b,k_b_offset,d_c,
2685
2685
     &k_c_offset)
2686
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
2686
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
2687
2687
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
2688
2688
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
2689
2689
C     i1 ( h7 h1 )_yt + = 1 * Sum ( p3 ) * t ( p3 h1 )_t * y ( h7 p3 )_y
2699
2699
      INTEGER k_b_offset
2700
2700
      INTEGER d_c
2701
2701
      INTEGER k_c_offset
2702
 
      INTEGER NXTVAL
 
2702
      INTEGER NXTASK
2703
2703
      INTEGER next
2704
2704
      INTEGER nprocs
2705
2705
      INTEGER count
2728
2728
      INTEGER k_b
2729
2729
      INTEGER l_c
2730
2730
      INTEGER k_c
2731
 
      EXTERNAL NXTVAL
 
2731
      EXTERNAL NXTASK
2732
2732
      nprocs = GA_NNODES()
2733
2733
      count = 0
2734
 
      next = NXTVAL(nprocs)
 
2734
      next = NXTASK(nprocs,1)
2735
2735
      DO h7b = 1,noab
2736
2736
      DO h1b = 1,noab
2737
2737
      IF (next.eq.count) THEN
2800
2800
      END IF
2801
2801
      END IF
2802
2802
      END IF
2803
 
      next = NXTVAL(nprocs)
 
2803
      next = NXTASK(nprocs,1)
2804
2804
      END IF
2805
2805
      count = count + 1
2806
2806
      END DO
2807
2807
      END DO
2808
 
      next = NXTVAL(-nprocs)
 
2808
      next = NXTASK(-nprocs,1)
2809
2809
      call GA_SYNC()
2810
2810
      RETURN
2811
2811
      END
2812
2812
      SUBROUTINE OFFSET_eomccsd_1prdm_ph_4_1(l_a_offset,k_a_offset,size)
2813
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
2813
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
2814
2814
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
2815
2815
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
2816
2816
C     i1 ( h7 h1 )_yt
2866
2866
      END
2867
2867
      SUBROUTINE eomccsd_1prdm_ph_4_2(d_a,k_a_offset,d_b,k_b_offset,d_c,
2868
2868
     &k_c_offset)
2869
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
2869
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
2870
2870
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
2871
2871
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
2872
2872
C     i1 ( h7 h1 )_yt + = -1/2 * Sum ( h5 p3 p4 ) * t ( p3 p4 h1 h5 )_t * y ( h5 h7 p3 p4 )_y
2882
2882
      INTEGER k_b_offset
2883
2883
      INTEGER d_c
2884
2884
      INTEGER k_c_offset
2885
 
      INTEGER NXTVAL
 
2885
      INTEGER NXTASK
2886
2886
      INTEGER next
2887
2887
      INTEGER nprocs
2888
2888
      INTEGER count
2920
2920
      INTEGER l_c
2921
2921
      INTEGER k_c
2922
2922
      DOUBLE PRECISION FACTORIAL
2923
 
      EXTERNAL NXTVAL
 
2923
      EXTERNAL NXTASK
2924
2924
      EXTERNAL FACTORIAL
2925
2925
      nprocs = GA_NNODES()
2926
2926
      count = 0
2927
 
      next = NXTVAL(nprocs)
 
2927
      next = NXTASK(nprocs,1)
2928
2928
      DO h7b = 1,noab
2929
2929
      DO h1b = 1,noab
2930
2930
      IF (next.eq.count) THEN
3031
3031
      END IF
3032
3032
      END IF
3033
3033
      END IF
3034
 
      next = NXTVAL(nprocs)
 
3034
      next = NXTASK(nprocs,1)
3035
3035
      END IF
3036
3036
      count = count + 1
3037
3037
      END DO
3038
3038
      END DO
3039
 
      next = NXTVAL(-nprocs)
 
3039
      next = NXTASK(-nprocs,1)
3040
3040
      call GA_SYNC()
3041
3041
      RETURN
3042
3042
      END
3043
3043
      SUBROUTINE eomccsd_1prdm_ph_5(d_a,k_a_offset,d_b,k_b_offset,d_c,k_
3044
3044
     &c_offset)
3045
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
3045
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
3046
3046
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
3047
3047
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
3048
3048
C     i0 ( p2 h1 )_yxt + = 1 * t ( p2 h1 )_t * i1 ( )_yx
3058
3058
      INTEGER k_b_offset
3059
3059
      INTEGER d_c
3060
3060
      INTEGER k_c_offset
3061
 
      INTEGER NXTVAL
 
3061
      INTEGER NXTASK
3062
3062
      INTEGER next
3063
3063
      INTEGER nprocs
3064
3064
      INTEGER count
3084
3084
      INTEGER k_b
3085
3085
      INTEGER l_c
3086
3086
      INTEGER k_c
3087
 
      EXTERNAL NXTVAL
 
3087
      EXTERNAL NXTASK
3088
3088
      nprocs = GA_NNODES()
3089
3089
      count = 0
3090
 
      next = NXTVAL(nprocs)
 
3090
      next = NXTASK(nprocs,1)
3091
3091
      DO p2b = noab+1,noab+nvab
3092
3092
      DO h1b = 1,noab
3093
3093
      IF (next.eq.count) THEN
3151
3151
      END IF
3152
3152
      END IF
3153
3153
      END IF
3154
 
      next = NXTVAL(nprocs)
 
3154
      next = NXTASK(nprocs,1)
3155
3155
      END IF
3156
3156
      count = count + 1
3157
3157
      END DO
3158
3158
      END DO
3159
 
      next = NXTVAL(-nprocs)
 
3159
      next = NXTASK(-nprocs,1)
3160
3160
      call GA_SYNC()
3161
3161
      RETURN
3162
3162
      END
3163
3163
      SUBROUTINE eomccsd_1prdm_ph_5_1(d_a,k_a_offset,d_b,k_b_offset,d_c,
3164
3164
     &k_c_offset)
3165
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
3165
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
3166
3166
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
3167
3167
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
3168
3168
C     i1 ( )_yx + = 1 * Sum ( h4 p3 ) * x ( p3 h4 )_x * y ( h4 p3 )_y
3178
3178
      INTEGER k_b_offset
3179
3179
      INTEGER d_c
3180
3180
      INTEGER k_c_offset
3181
 
      INTEGER NXTVAL
 
3181
      INTEGER NXTASK
3182
3182
      INTEGER next
3183
3183
      INTEGER nprocs
3184
3184
      INTEGER count
3206
3206
      INTEGER k_b
3207
3207
      INTEGER l_c
3208
3208
      INTEGER k_c
3209
 
      EXTERNAL NXTVAL
 
3209
      EXTERNAL NXTASK
3210
3210
      nprocs = GA_NNODES()
3211
3211
      count = 0
3212
 
      next = NXTVAL(nprocs)
 
3212
      next = NXTASK(nprocs,1)
3213
3213
      IF (next.eq.count) THEN
3214
3214
      IF (0 .eq. ieor(irrep_y,irrep_x)) THEN
3215
3215
      dimc = 1
3270
3270
      IF (.not.MA_POP_STACK(l_c_sort)) CALL ERRQUIT('eomccsd_1prdm_ph_5_
3271
3271
     &1',11,MA_ERR)
3272
3272
      END IF
3273
 
      next = NXTVAL(nprocs)
 
3273
      next = NXTASK(nprocs,1)
3274
3274
      END IF
3275
3275
      count = count + 1
3276
 
      next = NXTVAL(-nprocs)
 
3276
      next = NXTASK(-nprocs,1)
3277
3277
      call GA_SYNC()
3278
3278
      RETURN
3279
3279
      END
3280
3280
      SUBROUTINE OFFSET_eomccsd_1prdm_ph_5_1(l_a_offset,k_a_offset,size)
3281
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
3281
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
3282
3282
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
3283
3283
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
3284
3284
C     i1 ( )_yx
3312
3312
      END
3313
3313
      SUBROUTINE eomccsd_1prdm_ph_5_2(d_a,k_a_offset,d_b,k_b_offset,d_c,
3314
3314
     &k_c_offset)
3315
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
3315
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
3316
3316
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
3317
3317
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
3318
3318
C     i1 ( )_yx + = 1/4 * Sum ( h5 h6 p3 p4 ) * x ( p3 p4 h5 h6 )_x * y ( h5 h6 p3 p4 )_y
3328
3328
      INTEGER k_b_offset
3329
3329
      INTEGER d_c
3330
3330
      INTEGER k_c_offset
3331
 
      INTEGER NXTVAL
 
3331
      INTEGER NXTASK
3332
3332
      INTEGER next
3333
3333
      INTEGER nprocs
3334
3334
      INTEGER count
3367
3367
      INTEGER l_c
3368
3368
      INTEGER k_c
3369
3369
      DOUBLE PRECISION FACTORIAL
3370
 
      EXTERNAL NXTVAL
 
3370
      EXTERNAL NXTASK
3371
3371
      EXTERNAL FACTORIAL
3372
3372
      nprocs = GA_NNODES()
3373
3373
      count = 0
3374
 
      next = NXTVAL(nprocs)
 
3374
      next = NXTASK(nprocs,1)
3375
3375
      IF (next.eq.count) THEN
3376
3376
      IF (0 .eq. ieor(irrep_y,irrep_x)) THEN
3377
3377
      dimc = 1
3459
3459
      IF (.not.MA_POP_STACK(l_c_sort)) CALL ERRQUIT('eomccsd_1prdm_ph_5_
3460
3460
     &2',11,MA_ERR)
3461
3461
      END IF
3462
 
      next = NXTVAL(nprocs)
 
3462
      next = NXTASK(nprocs,1)
3463
3463
      END IF
3464
3464
      count = count + 1
3465
 
      next = NXTVAL(-nprocs)
 
3465
      next = NXTASK(-nprocs,1)
3466
3466
      call GA_SYNC()
3467
3467
      RETURN
3468
3468
      END
3469
3469
      SUBROUTINE eomccsd_1prdm_ph_6(d_a,k_a_offset,d_b,k_b_offset,d_c,k_
3470
3470
     &c_offset)
3471
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
3471
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
3472
3472
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
3473
3473
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
3474
3474
C     i0 ( p2 h1 )_yx + = 1 * Sum ( h4 p3 ) * x ( p2 p3 h1 h4 )_x * y ( h4 p3 )_y
3484
3484
      INTEGER k_b_offset
3485
3485
      INTEGER d_c
3486
3486
      INTEGER k_c_offset
3487
 
      INTEGER NXTVAL
 
3487
      INTEGER NXTASK
3488
3488
      INTEGER next
3489
3489
      INTEGER nprocs
3490
3490
      INTEGER count
3516
3516
      INTEGER k_b
3517
3517
      INTEGER l_c
3518
3518
      INTEGER k_c
3519
 
      EXTERNAL NXTVAL
 
3519
      EXTERNAL NXTASK
3520
3520
      nprocs = GA_NNODES()
3521
3521
      count = 0
3522
 
      next = NXTVAL(nprocs)
 
3522
      next = NXTASK(nprocs,1)
3523
3523
      DO p2b = noab+1,noab+nvab
3524
3524
      DO h1b = 1,noab
3525
3525
      IF (next.eq.count) THEN
3619
3619
      END IF
3620
3620
      END IF
3621
3621
      END IF
3622
 
      next = NXTVAL(nprocs)
 
3622
      next = NXTASK(nprocs,1)
3623
3623
      END IF
3624
3624
      count = count + 1
3625
3625
      END DO
3626
3626
      END DO
3627
 
      next = NXTVAL(-nprocs)
 
3627
      next = NXTASK(-nprocs,1)
3628
3628
      call GA_SYNC()
3629
3629
      RETURN
3630
3630
      END
3631
3631
      SUBROUTINE eomccsd_1prdm_ph_7(d_a,k_a_offset,d_b,k_b_offset,d_c,k_
3632
3632
     &c_offset)
3633
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
3633
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
3634
3634
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
3635
3635
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
3636
3636
C     i0 ( p2 h1 )_yxt + = -1/2 * Sum ( h4 h5 p3 ) * t ( p2 p3 h4 h5 )_t * i1 ( h4 h5 h1 p3 )_yx
3646
3646
      INTEGER k_b_offset
3647
3647
      INTEGER d_c
3648
3648
      INTEGER k_c_offset
3649
 
      INTEGER NXTVAL
 
3649
      INTEGER NXTASK
3650
3650
      INTEGER next
3651
3651
      INTEGER nprocs
3652
3652
      INTEGER count
3684
3684
      INTEGER l_c
3685
3685
      INTEGER k_c
3686
3686
      DOUBLE PRECISION FACTORIAL
3687
 
      EXTERNAL NXTVAL
 
3687
      EXTERNAL NXTASK
3688
3688
      EXTERNAL FACTORIAL
3689
3689
      nprocs = GA_NNODES()
3690
3690
      count = 0
3691
 
      next = NXTVAL(nprocs)
 
3691
      next = NXTASK(nprocs,1)
3692
3692
      DO p2b = noab+1,noab+nvab
3693
3693
      DO h1b = 1,noab
3694
3694
      IF (next.eq.count) THEN
3785
3785
      END IF
3786
3786
      END IF
3787
3787
      END IF
3788
 
      next = NXTVAL(nprocs)
 
3788
      next = NXTASK(nprocs,1)
3789
3789
      END IF
3790
3790
      count = count + 1
3791
3791
      END DO
3792
3792
      END DO
3793
 
      next = NXTVAL(-nprocs)
 
3793
      next = NXTASK(-nprocs,1)
3794
3794
      call GA_SYNC()
3795
3795
      RETURN
3796
3796
      END
3797
3797
      SUBROUTINE eomccsd_1prdm_ph_7_1(d_a,k_a_offset,d_b,k_b_offset,d_c,
3798
3798
     &k_c_offset)
3799
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
3799
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
3800
3800
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
3801
3801
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
3802
3802
C     i1 ( h4 h5 h1 p3 )_yx + = -1 * Sum ( p6 ) * x ( p6 h1 )_x * y ( h4 h5 p3 p6 )_y
3812
3812
      INTEGER k_b_offset
3813
3813
      INTEGER d_c
3814
3814
      INTEGER k_c_offset
3815
 
      INTEGER NXTVAL
 
3815
      INTEGER NXTASK
3816
3816
      INTEGER next
3817
3817
      INTEGER nprocs
3818
3818
      INTEGER count
3845
3845
      INTEGER k_b
3846
3846
      INTEGER l_c
3847
3847
      INTEGER k_c
3848
 
      EXTERNAL NXTVAL
 
3848
      EXTERNAL NXTASK
3849
3849
      nprocs = GA_NNODES()
3850
3850
      count = 0
3851
 
      next = NXTVAL(nprocs)
 
3851
      next = NXTASK(nprocs,1)
3852
3852
      DO h4b = 1,noab
3853
3853
      DO h5b = h4b,noab
3854
3854
      DO h1b = 1,noab
3937
3937
      END IF
3938
3938
      END IF
3939
3939
      END IF
3940
 
      next = NXTVAL(nprocs)
 
3940
      next = NXTASK(nprocs,1)
3941
3941
      END IF
3942
3942
      count = count + 1
3943
3943
      END DO
3944
3944
      END DO
3945
3945
      END DO
3946
3946
      END DO
3947
 
      next = NXTVAL(-nprocs)
 
3947
      next = NXTASK(-nprocs,1)
3948
3948
      call GA_SYNC()
3949
3949
      RETURN
3950
3950
      END
3951
3951
      SUBROUTINE OFFSET_eomccsd_1prdm_ph_7_1(l_a_offset,k_a_offset,size)
3952
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
3952
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
3953
3953
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
3954
3954
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
3955
3955
C     i1 ( h4 h5 h1 p3 )_yx
4021
4021
      END
4022
4022
      SUBROUTINE eomccsd_1prdm_ph_8(d_a,k_a_offset,d_b,k_b_offset,d_c,k_
4023
4023
     &c_offset)
4024
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
4024
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
4025
4025
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
4026
4026
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
4027
4027
C     i0 ( p2 h1 )_ytx + = 1/2 * Sum ( h5 h6 p4 ) * x ( p2 p4 h5 h6 )_x * i1 ( h5 h6 h1 p4 )_yt
4037
4037
      INTEGER k_b_offset
4038
4038
      INTEGER d_c
4039
4039
      INTEGER k_c_offset
4040
 
      INTEGER NXTVAL
 
4040
      INTEGER NXTASK
4041
4041
      INTEGER next
4042
4042
      INTEGER nprocs
4043
4043
      INTEGER count
4075
4075
      INTEGER l_c
4076
4076
      INTEGER k_c
4077
4077
      DOUBLE PRECISION FACTORIAL
4078
 
      EXTERNAL NXTVAL
 
4078
      EXTERNAL NXTASK
4079
4079
      EXTERNAL FACTORIAL
4080
4080
      nprocs = GA_NNODES()
4081
4081
      count = 0
4082
 
      next = NXTVAL(nprocs)
 
4082
      next = NXTASK(nprocs,1)
4083
4083
      DO p2b = noab+1,noab+nvab
4084
4084
      DO h1b = 1,noab
4085
4085
      IF (next.eq.count) THEN
4176
4176
      END IF
4177
4177
      END IF
4178
4178
      END IF
4179
 
      next = NXTVAL(nprocs)
 
4179
      next = NXTASK(nprocs,1)
4180
4180
      END IF
4181
4181
      count = count + 1
4182
4182
      END DO
4183
4183
      END DO
4184
 
      next = NXTVAL(-nprocs)
 
4184
      next = NXTASK(-nprocs,1)
4185
4185
      call GA_SYNC()
4186
4186
      RETURN
4187
4187
      END
4188
4188
      SUBROUTINE eomccsd_1prdm_ph_8_1(d_a,k_a_offset,d_b,k_b_offset,d_c,
4189
4189
     &k_c_offset)
4190
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
4190
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
4191
4191
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
4192
4192
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
4193
4193
C     i1 ( h5 h6 h1 p4 )_yt + = -1 * Sum ( p3 ) * t ( p3 h1 )_t * y ( h5 h6 p3 p4 )_y
4203
4203
      INTEGER k_b_offset
4204
4204
      INTEGER d_c
4205
4205
      INTEGER k_c_offset
4206
 
      INTEGER NXTVAL
 
4206
      INTEGER NXTASK
4207
4207
      INTEGER next
4208
4208
      INTEGER nprocs
4209
4209
      INTEGER count
4236
4236
      INTEGER k_b
4237
4237
      INTEGER l_c
4238
4238
      INTEGER k_c
4239
 
      EXTERNAL NXTVAL
 
4239
      EXTERNAL NXTASK
4240
4240
      nprocs = GA_NNODES()
4241
4241
      count = 0
4242
 
      next = NXTVAL(nprocs)
 
4242
      next = NXTASK(nprocs,1)
4243
4243
      DO h5b = 1,noab
4244
4244
      DO h6b = h5b,noab
4245
4245
      DO h1b = 1,noab
4328
4328
      END IF
4329
4329
      END IF
4330
4330
      END IF
4331
 
      next = NXTVAL(nprocs)
 
4331
      next = NXTASK(nprocs,1)
4332
4332
      END IF
4333
4333
      count = count + 1
4334
4334
      END DO
4335
4335
      END DO
4336
4336
      END DO
4337
4337
      END DO
4338
 
      next = NXTVAL(-nprocs)
 
4338
      next = NXTASK(-nprocs,1)
4339
4339
      call GA_SYNC()
4340
4340
      RETURN
4341
4341
      END
4342
4342
      SUBROUTINE OFFSET_eomccsd_1prdm_ph_8_1(l_a_offset,k_a_offset,size)
4343
 
C     $Id: eomccsd_1prdm_ph.F 19706 2010-10-29 17:52:31Z d3y133 $
 
4343
C     $Id: eomccsd_1prdm_ph.F 24183 2013-05-03 23:47:08Z edo $
4344
4344
C     This is a Fortran77 program generated by Tensor Contraction Engine v.1.0
4345
4345
C     Copyright (c) Battelle & Pacific Northwest National Laboratory (2002)
4346
4346
C     i1 ( h5 h6 h1 p4 )_yt