1 SUBROUTINE ccsdt_lambda_o1(d_i0,d_o1,d_t1,d_t2,d_t3,d_y1,d_y2,d_y3 2 &,k_i0_offset,k_o1_offset,k_t1_offset,k_t2_offset,k_t3_offset,k_y1_ 3 &offset,k_y2_offset,k_y3_offset) 4C $Id$ 5C This is a Fortran77 program generated by Tensor Contraction Engine v.1.0 6C Copyright (c) Battelle & Pacific Northwest National Laboratory (2002) 7C i0 ( h2 p1 )_yo + = -1 * Sum ( h5 ) * y ( h5 p1 )_y * i1 ( h2 h5 )_o 8C i1 ( h2 h5 )_o + = 1 * o ( h2 h5 )_o 9C i1 ( h2 h5 )_ot + = 1 * Sum ( p3 ) * o ( h2 p3 )_o * t ( p3 h5 )_t 10C i0 ( h2 p1 )_yo + = 1 * Sum ( p3 ) * o ( p3 p1 )_o * y ( h2 p3 )_y 11C i0 ( h2 p1 )_yo + = 1 * Sum ( h4 p3 ) * o ( p3 h4 )_o * y ( h2 h4 p1 p3 )_y 12C i0 ( h2 p1 )_yto + = -1 * Sum ( h3 ) * o ( h3 p1 )_o * i1 ( h2 h3 )_yt 13C i1 ( h2 h3 )_yt + = 1 * Sum ( p4 ) * t ( p4 h3 )_t * y ( h2 p4 )_y 14C i1 ( h2 h3 )_yt + = 1/2 * Sum ( h6 p5 p4 ) * t ( p4 p5 h3 h6 )_t * y ( h2 h6 p4 p5 )_y 15C i1 ( h2 h3 )_yt + = 1/12 * Sum ( h8 h7 p6 p5 p4 ) * t ( p4 p5 p6 h3 h7 h8 )_t * y ( h2 h7 h8 p4 p5 p6 )_y 16C i0 ( h2 p1 )_oty + = 1/2 * Sum ( p4 h5 h6 ) * y ( h5 h6 p1 p4 )_y * i1 ( h2 p4 h5 h6 )_ot 17C i1 ( h2 p4 h5 h6 )_ot + = -1 * Sum ( p3 ) * o ( h2 p3 )_o * t ( p3 p4 h5 h6 )_t 18C i0 ( h2 p1 )_oty + = -1/12 * Sum ( p4 p5 h6 h7 h8 ) * y ( h6 h7 h8 p1 p4 p5 )_y * i1 ( h2 p4 p5 h6 h7 h8 )_ot 19C i1 ( h2 p4 p5 h6 h7 h8 )_ot + = 1 * Sum ( p3 ) * o ( h2 p3 )_o * t ( p3 p4 p5 h6 h7 h8 )_t 20 IMPLICIT NONE 21#include "global.fh" 22#include "mafdecls.fh" 23#include "util.fh" 24#include "errquit.fh" 25#include "tce.fh" 26 INTEGER d_i0 27 INTEGER k_i0_offset 28 INTEGER d_y1 29 INTEGER k_y1_offset 30 INTEGER d_i1 31 INTEGER k_i1_offset 32 INTEGER d_o1 33 INTEGER k_o1_offset 34 INTEGER d_y2 35 INTEGER k_y2_offset 36 INTEGER d_y3 37 INTEGER k_y3_offset 38 INTEGER l_i1_offset 39 INTEGER size_i1 40 INTEGER d_t1 41 INTEGER k_t1_offset 42 INTEGER d_t2 43 INTEGER k_t2_offset 44 INTEGER d_t3 45 INTEGER k_t3_offset 46 CHARACTER*255 filename 47#if !defined(CCSDTLR) 48 CALL ERRQUIT('CCSDTLR not compiled (ccsdt_lambda_o1.F)', 49 & 0,CAPMIS_ERR) 50 RETURN 51 END 52#else 53 CALL OFFSET_ccsdt_lambda_o1_1_1(l_i1_offset,k_i1_offset,size_i1) 54 CALL TCE_FILENAME('ccsdt_lambda_o1_1_1_i1',filename) 55 CALL CREATEFILE(filename,d_i1,size_i1) 56 CALL ccsdt_lambda_o1_1_1(d_o1,k_o1_offset,d_i1,k_i1_offset) 57 CALL ccsdt_lambda_o1_1_2(d_o1,k_o1_offset,d_t1,k_t1_offset,d_i1,k_ 58 &i1_offset) 59 CALL RECONCILEFILE(d_i1,size_i1) 60 CALL ccsdt_lambda_o1_1(d_y1,k_y1_offset,d_i1,k_i1_offset,d_i0,k_i0 61 &_offset) 62 CALL DELETEFILE(d_i1) 63 IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('ccsdt_lambda_o1' 64 &,-1,MA_ERR) 65 CALL ccsdt_lambda_o1_2(d_o1,k_o1_offset,d_y1,k_y1_offset,d_i0,k_i0 66 &_offset) 67 CALL ccsdt_lambda_o1_3(d_o1,k_o1_offset,d_y2,k_y2_offset,d_i0,k_i0 68 &_offset) 69 CALL OFFSET_ccsdt_lambda_o1_4_1(l_i1_offset,k_i1_offset,size_i1) 70 CALL TCE_FILENAME('ccsdt_lambda_o1_4_1_i1',filename) 71 CALL CREATEFILE(filename,d_i1,size_i1) 72 CALL ccsdt_lambda_o1_4_1(d_t1,k_t1_offset,d_y1,k_y1_offset,d_i1,k_ 73 &i1_offset) 74 CALL ccsdt_lambda_o1_4_2(d_t2,k_t2_offset,d_y2,k_y2_offset,d_i1,k_ 75 &i1_offset) 76 CALL ccsdt_lambda_o1_4_3(d_t3,k_t3_offset,d_y3,k_y3_offset,d_i1,k_ 77 &i1_offset) 78 CALL RECONCILEFILE(d_i1,size_i1) 79 CALL ccsdt_lambda_o1_4(d_o1,k_o1_offset,d_i1,k_i1_offset,d_i0,k_i0 80 &_offset) 81 CALL DELETEFILE(d_i1) 82 IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('ccsdt_lambda_o1' 83 &,-1,MA_ERR) 84 CALL OFFSET_ccsdt_lambda_o1_5_1(l_i1_offset,k_i1_offset,size_i1) 85 CALL TCE_FILENAME('ccsdt_lambda_o1_5_1_i1',filename) 86 CALL CREATEFILE(filename,d_i1,size_i1) 87 CALL ccsdt_lambda_o1_5_1(d_o1,k_o1_offset,d_t2,k_t2_offset,d_i1,k_ 88 &i1_offset) 89 CALL RECONCILEFILE(d_i1,size_i1) 90 CALL ccsdt_lambda_o1_5(d_y2,k_y2_offset,d_i1,k_i1_offset,d_i0,k_i0 91 &_offset) 92 CALL DELETEFILE(d_i1) 93 IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('ccsdt_lambda_o1' 94 &,-1,MA_ERR) 95 CALL OFFSET_ccsdt_lambda_o1_6_1(l_i1_offset,k_i1_offset,size_i1) 96 CALL TCE_FILENAME('ccsdt_lambda_o1_6_1_i1',filename) 97 CALL CREATEFILE(filename,d_i1,size_i1) 98 CALL ccsdt_lambda_o1_6_1(d_o1,k_o1_offset,d_t3,k_t3_offset,d_i1,k_ 99 &i1_offset) 100 CALL RECONCILEFILE(d_i1,size_i1) 101 CALL ccsdt_lambda_o1_6(d_y3,k_y3_offset,d_i1,k_i1_offset,d_i0,k_i0 102 &_offset) 103 CALL DELETEFILE(d_i1) 104 IF (.not.MA_POP_STACK(l_i1_offset)) CALL ERRQUIT('ccsdt_lambda_o1' 105 &,-1,MA_ERR) 106 RETURN 107 END 108#endif 109