/*
 *  Copyright (C) 2004-2026 Edward F. Valeev
 *
 *  This file is part of Libint library.
 *
 *  Libint library is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU Lesser General Public License as published by
 *  the Free Software Foundation, either version 3 of the License, or
 *  (at your option) any later version.
 *
 *  Libint library is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU Lesser General Public License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public License
 *  along with Libint library.  If not, see <http://www.gnu.org/licenses/>.
 *
 */

#include <libint2.h>

#ifdef __cplusplus
LIBINT_PRAGMA_CLANG(diagnostic push)
LIBINT_PRAGMA_CLANG(diagnostic ignored "-Wunused-variable")
LIBINT_PRAGMA_GCC(diagnostic push)
LIBINT_PRAGMA_GCC(diagnostic ignored "-Wunused-variable")
extern "C" {
#endif
void CR_aB_p__0___Kinetic_h__0___Ab__up_(const Libint_t* inteval, LIBINT2_REALTYPE* target, const LIBINT2_REALTYPE* src0, const LIBINT2_REALTYPE* src1, const LIBINT2_REALTYPE* src2) {

LIBINT2_REALTYPE*const  stack = target;
{
const int hsi = 0;
{
const int lsi = 0;
{
const int vi = 0;
LIBINT2_REALTYPE fp75;
fp75 = 1.5000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp74;
fp74 = fp75 * src2[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp78;
fp78 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp77;
fp77 = fp78 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp76;
fp76 = fp77 * src2[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp73;
fp73 = fp76 - fp74;
LIBINT2_REALTYPE fp540;
fp540 = src0[((hsi*21+2)*1+lsi)*1] * src1[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp539;
fp539 = fp540 * fp73;
LIBINT2_REALTYPE fp151;
fp151 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp150;
fp150 = fp151 * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp154;
fp154 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp153;
fp153 = fp154 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp152;
fp152 = fp153 * src1[((hsi*21+15)*1+lsi)*1];
LIBINT2_REALTYPE fp149;
fp149 = fp152 - fp150;
LIBINT2_REALTYPE fp543;
fp543 = src0[((hsi*21+2)*1+lsi)*1] * fp149;
LIBINT2_REALTYPE fp542;
fp542 = fp543 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp199;
fp199 = 1.0000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp198;
fp198 = fp199 * src0[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp202;
fp202 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp201;
fp201 = fp202 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp200;
fp200 = fp201 * src0[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp197;
fp197 = fp200 - fp198;
LIBINT2_REALTYPE fp545;
fp545 = fp197 * src1[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp544;
fp544 = fp545 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp541;
fp541 = fp544 + fp542;
LIBINT2_REALTYPE fp538;
fp538 = fp541 + fp539;
target[((hsi*63+30)*1+lsi)*1] = fp538;
LIBINT2_REALTYPE fp81;
fp81 = 1.0000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp80;
fp80 = fp81 * src2[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp84;
fp84 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp83;
fp83 = fp84 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp82;
fp82 = fp83 * src2[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp79;
fp79 = fp82 - fp80;
LIBINT2_REALTYPE fp548;
fp548 = src0[((hsi*21+2)*1+lsi)*1] * src1[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp547;
fp547 = fp548 * fp79;
LIBINT2_REALTYPE fp139;
fp139 = 5.0000000000000000e-01 * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp142;
fp142 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp141;
fp141 = fp142 * src1[((hsi*21+14)*1+lsi)*1];
LIBINT2_REALTYPE fp145;
fp145 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp144;
fp144 = fp145 * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp148;
fp148 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp147;
fp147 = fp148 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp146;
fp146 = fp147 * src1[((hsi*21+16)*1+lsi)*1];
LIBINT2_REALTYPE fp143;
fp143 = fp146 - fp144;
LIBINT2_REALTYPE fp140;
fp140 = fp143 - fp141;
LIBINT2_REALTYPE fp138;
fp138 = fp140 + fp139;
LIBINT2_REALTYPE fp551;
fp551 = src0[((hsi*21+2)*1+lsi)*1] * fp138;
LIBINT2_REALTYPE fp550;
fp550 = fp551 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp553;
fp553 = fp197 * src1[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp552;
fp552 = fp553 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp549;
fp549 = fp552 + fp550;
LIBINT2_REALTYPE fp546;
fp546 = fp549 + fp547;
target[((hsi*63+29)*1+lsi)*1] = fp546;
LIBINT2_REALTYPE fp87;
fp87 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp86;
fp86 = fp87 * src2[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp90;
fp90 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp89;
fp89 = fp90 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp88;
fp88 = fp89 * src2[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp85;
fp85 = fp88 - fp86;
LIBINT2_REALTYPE fp556;
fp556 = src0[((hsi*21+2)*1+lsi)*1] * src1[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp555;
fp555 = fp556 * fp85;
LIBINT2_REALTYPE fp128;
fp128 = 1.0000000000000000e+00 * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp131;
fp131 = 1.0000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp130;
fp130 = fp131 * src1[((hsi*21+15)*1+lsi)*1];
LIBINT2_REALTYPE fp134;
fp134 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp133;
fp133 = fp134 * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp137;
fp137 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp136;
fp136 = fp137 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp135;
fp135 = fp136 * src1[((hsi*21+17)*1+lsi)*1];
LIBINT2_REALTYPE fp132;
fp132 = fp135 - fp133;
LIBINT2_REALTYPE fp129;
fp129 = fp132 - fp130;
LIBINT2_REALTYPE fp127;
fp127 = fp129 + fp128;
LIBINT2_REALTYPE fp559;
fp559 = src0[((hsi*21+2)*1+lsi)*1] * fp127;
LIBINT2_REALTYPE fp558;
fp558 = fp559 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp561;
fp561 = fp197 * src1[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp560;
fp560 = fp561 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp557;
fp557 = fp560 + fp558;
LIBINT2_REALTYPE fp554;
fp554 = fp557 + fp555;
target[((hsi*63+28)*1+lsi)*1] = fp554;
LIBINT2_REALTYPE fp93;
fp93 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp92;
fp92 = fp93 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp91;
fp91 = fp92 * src2[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp564;
fp564 = src0[((hsi*21+2)*1+lsi)*1] * src1[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp563;
fp563 = fp564 * fp91;
LIBINT2_REALTYPE fp117;
fp117 = 1.5000000000000000e+00 * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp120;
fp120 = 1.5000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp119;
fp119 = fp120 * src1[((hsi*21+16)*1+lsi)*1];
LIBINT2_REALTYPE fp123;
fp123 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp122;
fp122 = fp123 * src1[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp126;
fp126 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp125;
fp125 = fp126 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp124;
fp124 = fp125 * src1[((hsi*21+18)*1+lsi)*1];
LIBINT2_REALTYPE fp121;
fp121 = fp124 - fp122;
LIBINT2_REALTYPE fp118;
fp118 = fp121 - fp119;
LIBINT2_REALTYPE fp116;
fp116 = fp118 + fp117;
LIBINT2_REALTYPE fp567;
fp567 = src0[((hsi*21+2)*1+lsi)*1] * fp116;
LIBINT2_REALTYPE fp566;
fp566 = fp567 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp569;
fp569 = fp197 * src1[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp568;
fp568 = fp569 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp565;
fp565 = fp568 + fp566;
LIBINT2_REALTYPE fp562;
fp562 = fp565 + fp563;
target[((hsi*63+27)*1+lsi)*1] = fp562;
LIBINT2_REALTYPE fp572;
fp572 = src0[((hsi*21+3)*1+lsi)*1] * src1[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp571;
fp571 = fp572 * fp79;
LIBINT2_REALTYPE fp575;
fp575 = src0[((hsi*21+3)*1+lsi)*1] * fp149;
LIBINT2_REALTYPE fp574;
fp574 = fp575 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp205;
fp205 = 1.5000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp204;
fp204 = fp205 * src0[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp208;
fp208 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp207;
fp207 = fp208 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp206;
fp206 = fp207 * src0[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp203;
fp203 = fp206 - fp204;
LIBINT2_REALTYPE fp577;
fp577 = fp203 * src1[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp576;
fp576 = fp577 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp573;
fp573 = fp576 + fp574;
LIBINT2_REALTYPE fp570;
fp570 = fp573 + fp571;
target[((hsi*63+26)*1+lsi)*1] = fp570;
LIBINT2_REALTYPE fp580;
fp580 = src0[((hsi*21+3)*1+lsi)*1] * src1[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp579;
fp579 = fp580 * fp85;
LIBINT2_REALTYPE fp583;
fp583 = src0[((hsi*21+3)*1+lsi)*1] * fp138;
LIBINT2_REALTYPE fp582;
fp582 = fp583 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp585;
fp585 = fp203 * src1[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp584;
fp584 = fp585 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp581;
fp581 = fp584 + fp582;
LIBINT2_REALTYPE fp578;
fp578 = fp581 + fp579;
target[((hsi*63+25)*1+lsi)*1] = fp578;
LIBINT2_REALTYPE fp588;
fp588 = src0[((hsi*21+3)*1+lsi)*1] * src1[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp587;
fp587 = fp588 * fp91;
LIBINT2_REALTYPE fp591;
fp591 = src0[((hsi*21+3)*1+lsi)*1] * fp127;
LIBINT2_REALTYPE fp590;
fp590 = fp591 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp593;
fp593 = fp203 * src1[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp592;
fp592 = fp593 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp589;
fp589 = fp592 + fp590;
LIBINT2_REALTYPE fp586;
fp586 = fp589 + fp587;
target[((hsi*63+24)*1+lsi)*1] = fp586;
LIBINT2_REALTYPE fp596;
fp596 = src0[((hsi*21+4)*1+lsi)*1] * src1[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp595;
fp595 = fp596 * fp85;
LIBINT2_REALTYPE fp599;
fp599 = src0[((hsi*21+4)*1+lsi)*1] * fp149;
LIBINT2_REALTYPE fp598;
fp598 = fp599 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp211;
fp211 = 2.0000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp210;
fp210 = fp211 * src0[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp214;
fp214 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp213;
fp213 = fp214 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp212;
fp212 = fp213 * src0[((hsi*21+12)*1+lsi)*1];
LIBINT2_REALTYPE fp209;
fp209 = fp212 - fp210;
LIBINT2_REALTYPE fp601;
fp601 = fp209 * src1[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp600;
fp600 = fp601 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp597;
fp597 = fp600 + fp598;
LIBINT2_REALTYPE fp594;
fp594 = fp597 + fp595;
target[((hsi*63+23)*1+lsi)*1] = fp594;
LIBINT2_REALTYPE fp604;
fp604 = src0[((hsi*21+4)*1+lsi)*1] * src1[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp603;
fp603 = fp604 * fp91;
LIBINT2_REALTYPE fp607;
fp607 = src0[((hsi*21+4)*1+lsi)*1] * fp138;
LIBINT2_REALTYPE fp606;
fp606 = fp607 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp609;
fp609 = fp209 * src1[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp608;
fp608 = fp609 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp605;
fp605 = fp608 + fp606;
LIBINT2_REALTYPE fp602;
fp602 = fp605 + fp603;
target[((hsi*63+22)*1+lsi)*1] = fp602;
LIBINT2_REALTYPE fp612;
fp612 = src0[((hsi*21+5)*1+lsi)*1] * src1[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp611;
fp611 = fp612 * fp91;
LIBINT2_REALTYPE fp615;
fp615 = src0[((hsi*21+5)*1+lsi)*1] * fp149;
LIBINT2_REALTYPE fp614;
fp614 = fp615 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp217;
fp217 = 2.5000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp216;
fp216 = fp217 * src0[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp220;
fp220 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp219;
fp219 = fp220 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp218;
fp218 = fp219 * src0[((hsi*21+13)*1+lsi)*1];
LIBINT2_REALTYPE fp215;
fp215 = fp218 - fp216;
LIBINT2_REALTYPE fp617;
fp617 = fp215 * src1[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp616;
fp616 = fp617 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp613;
fp613 = fp616 + fp614;
LIBINT2_REALTYPE fp610;
fp610 = fp613 + fp611;
target[((hsi*63+21)*1+lsi)*1] = fp610;
LIBINT2_REALTYPE fp63;
fp63 = 2.5000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp62;
fp62 = fp63 * src2[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp66;
fp66 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp65;
fp65 = fp66 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp64;
fp64 = fp65 * src2[((hsi*21+13)*1+lsi)*1];
LIBINT2_REALTYPE fp61;
fp61 = fp64 - fp62;
LIBINT2_REALTYPE fp620;
fp620 = src0[((hsi*21+7)*1+lsi)*1] * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp619;
fp619 = fp620 * fp61;
LIBINT2_REALTYPE fp187;
fp187 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp186;
fp186 = fp187 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp185;
fp185 = fp186 * src1[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp623;
fp623 = src0[((hsi*21+7)*1+lsi)*1] * fp185;
LIBINT2_REALTYPE fp622;
fp622 = fp623 * src2[((hsi*21+5)*1+lsi)*1];
LIBINT2_REALTYPE fp223;
fp223 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp222;
fp222 = fp223 * src0[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp226;
fp226 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp225;
fp225 = fp226 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp224;
fp224 = fp225 * src0[((hsi*21+15)*1+lsi)*1];
LIBINT2_REALTYPE fp221;
fp221 = fp224 - fp222;
LIBINT2_REALTYPE fp625;
fp625 = fp221 * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp624;
fp624 = fp625 * src2[((hsi*21+5)*1+lsi)*1];
LIBINT2_REALTYPE fp621;
fp621 = fp624 + fp622;
LIBINT2_REALTYPE fp618;
fp618 = fp621 + fp619;
target[((hsi*63+20)*1+lsi)*1] = fp618;
LIBINT2_REALTYPE fp69;
fp69 = 2.0000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp68;
fp68 = fp69 * src2[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp72;
fp72 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp71;
fp71 = fp72 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp70;
fp70 = fp71 * src2[((hsi*21+12)*1+lsi)*1];
LIBINT2_REALTYPE fp67;
fp67 = fp70 - fp68;
LIBINT2_REALTYPE fp628;
fp628 = src0[((hsi*21+7)*1+lsi)*1] * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp627;
fp627 = fp628 * fp67;
LIBINT2_REALTYPE fp181;
fp181 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp180;
fp180 = fp181 * src1[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp184;
fp184 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp183;
fp183 = fp184 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp182;
fp182 = fp183 * src1[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp179;
fp179 = fp182 - fp180;
LIBINT2_REALTYPE fp631;
fp631 = src0[((hsi*21+7)*1+lsi)*1] * fp179;
LIBINT2_REALTYPE fp630;
fp630 = fp631 * src2[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp633;
fp633 = fp221 * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp632;
fp632 = fp633 * src2[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp629;
fp629 = fp632 + fp630;
LIBINT2_REALTYPE fp626;
fp626 = fp629 + fp627;
target[((hsi*63+19)*1+lsi)*1] = fp626;
LIBINT2_REALTYPE fp636;
fp636 = src0[((hsi*21+7)*1+lsi)*1] * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp635;
fp635 = fp636 * fp73;
LIBINT2_REALTYPE fp175;
fp175 = 1.0000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp174;
fp174 = fp175 * src1[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp178;
fp178 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp177;
fp177 = fp178 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp176;
fp176 = fp177 * src1[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp173;
fp173 = fp176 - fp174;
LIBINT2_REALTYPE fp639;
fp639 = src0[((hsi*21+7)*1+lsi)*1] * fp173;
LIBINT2_REALTYPE fp638;
fp638 = fp639 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp641;
fp641 = fp221 * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp640;
fp640 = fp641 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp637;
fp637 = fp640 + fp638;
LIBINT2_REALTYPE fp634;
fp634 = fp637 + fp635;
target[((hsi*63+18)*1+lsi)*1] = fp634;
LIBINT2_REALTYPE fp644;
fp644 = src0[((hsi*21+7)*1+lsi)*1] * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp643;
fp643 = fp644 * fp79;
LIBINT2_REALTYPE fp169;
fp169 = 1.5000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp168;
fp168 = fp169 * src1[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp172;
fp172 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp171;
fp171 = fp172 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp170;
fp170 = fp171 * src1[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp167;
fp167 = fp170 - fp168;
LIBINT2_REALTYPE fp647;
fp647 = src0[((hsi*21+7)*1+lsi)*1] * fp167;
LIBINT2_REALTYPE fp646;
fp646 = fp647 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp649;
fp649 = fp221 * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp648;
fp648 = fp649 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp645;
fp645 = fp648 + fp646;
LIBINT2_REALTYPE fp642;
fp642 = fp645 + fp643;
target[((hsi*63+17)*1+lsi)*1] = fp642;
LIBINT2_REALTYPE fp780;
fp780 = src0[((hsi*21+12)*1+lsi)*1] * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp779;
fp779 = fp780 * fp91;
LIBINT2_REALTYPE fp783;
fp783 = src0[((hsi*21+12)*1+lsi)*1] * fp185;
LIBINT2_REALTYPE fp782;
fp782 = fp783 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp272;
fp272 = 2.5000000000000000e+00 * src0[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp275;
fp275 = 2.5000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp274;
fp274 = fp275 * src0[((hsi*21+18)*1+lsi)*1];
LIBINT2_REALTYPE fp278;
fp278 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp277;
fp277 = fp278 * src0[((hsi*21+6)*1+lsi)*1];
LIBINT2_REALTYPE fp281;
fp281 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp280;
fp280 = fp281 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp279;
fp279 = fp280 * src0[((hsi*21+20)*1+lsi)*1];
LIBINT2_REALTYPE fp276;
fp276 = fp279 - fp277;
LIBINT2_REALTYPE fp273;
fp273 = fp276 - fp274;
LIBINT2_REALTYPE fp271;
fp271 = fp273 + fp272;
LIBINT2_REALTYPE fp785;
fp785 = fp271 * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp784;
fp784 = fp785 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp781;
fp781 = fp784 + fp782;
LIBINT2_REALTYPE fp778;
fp778 = fp781 + fp779;
target[((hsi*63+0)*1+lsi)*1] = fp778;
LIBINT2_REALTYPE fp660;
fp660 = src0[((hsi*21+7)*1+lsi)*1] * src1[((hsi*21+5)*1+lsi)*1];
LIBINT2_REALTYPE fp659;
fp659 = fp660 * fp91;
LIBINT2_REALTYPE fp157;
fp157 = 2.5000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp156;
fp156 = fp157 * src1[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp160;
fp160 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp159;
fp159 = fp160 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp158;
fp158 = fp159 * src1[((hsi*21+13)*1+lsi)*1];
LIBINT2_REALTYPE fp155;
fp155 = fp158 - fp156;
LIBINT2_REALTYPE fp663;
fp663 = src0[((hsi*21+7)*1+lsi)*1] * fp155;
LIBINT2_REALTYPE fp662;
fp662 = fp663 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp665;
fp665 = fp221 * src1[((hsi*21+5)*1+lsi)*1];
LIBINT2_REALTYPE fp664;
fp664 = fp665 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp661;
fp661 = fp664 + fp662;
LIBINT2_REALTYPE fp658;
fp658 = fp661 + fp659;
target[((hsi*63+15)*1+lsi)*1] = fp658;
LIBINT2_REALTYPE fp668;
fp668 = src0[((hsi*21+8)*1+lsi)*1] * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp667;
fp667 = fp668 * fp67;
LIBINT2_REALTYPE fp671;
fp671 = src0[((hsi*21+8)*1+lsi)*1] * fp185;
LIBINT2_REALTYPE fp670;
fp670 = fp671 * src2[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp228;
fp228 = 5.0000000000000000e-01 * src0[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp231;
fp231 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp230;
fp230 = fp231 * src0[((hsi*21+14)*1+lsi)*1];
LIBINT2_REALTYPE fp234;
fp234 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp233;
fp233 = fp234 * src0[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp237;
fp237 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp236;
fp236 = fp237 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp235;
fp235 = fp236 * src0[((hsi*21+16)*1+lsi)*1];
LIBINT2_REALTYPE fp232;
fp232 = fp235 - fp233;
LIBINT2_REALTYPE fp229;
fp229 = fp232 - fp230;
LIBINT2_REALTYPE fp227;
fp227 = fp229 + fp228;
LIBINT2_REALTYPE fp673;
fp673 = fp227 * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp672;
fp672 = fp673 * src2[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp669;
fp669 = fp672 + fp670;
LIBINT2_REALTYPE fp666;
fp666 = fp669 + fp667;
target[((hsi*63+14)*1+lsi)*1] = fp666;
LIBINT2_REALTYPE fp676;
fp676 = src0[((hsi*21+8)*1+lsi)*1] * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp675;
fp675 = fp676 * fp73;
LIBINT2_REALTYPE fp679;
fp679 = src0[((hsi*21+8)*1+lsi)*1] * fp179;
LIBINT2_REALTYPE fp678;
fp678 = fp679 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp681;
fp681 = fp227 * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp680;
fp680 = fp681 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp677;
fp677 = fp680 + fp678;
LIBINT2_REALTYPE fp674;
fp674 = fp677 + fp675;
target[((hsi*63+13)*1+lsi)*1] = fp674;
LIBINT2_REALTYPE fp684;
fp684 = src0[((hsi*21+8)*1+lsi)*1] * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp683;
fp683 = fp684 * fp79;
LIBINT2_REALTYPE fp687;
fp687 = src0[((hsi*21+8)*1+lsi)*1] * fp173;
LIBINT2_REALTYPE fp686;
fp686 = fp687 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp689;
fp689 = fp227 * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp688;
fp688 = fp689 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp685;
fp685 = fp688 + fp686;
LIBINT2_REALTYPE fp682;
fp682 = fp685 + fp683;
target[((hsi*63+12)*1+lsi)*1] = fp682;
LIBINT2_REALTYPE fp692;
fp692 = src0[((hsi*21+8)*1+lsi)*1] * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp691;
fp691 = fp692 * fp85;
LIBINT2_REALTYPE fp695;
fp695 = src0[((hsi*21+8)*1+lsi)*1] * fp167;
LIBINT2_REALTYPE fp694;
fp694 = fp695 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp697;
fp697 = fp227 * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp696;
fp696 = fp697 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp693;
fp693 = fp696 + fp694;
LIBINT2_REALTYPE fp690;
fp690 = fp693 + fp691;
target[((hsi*63+11)*1+lsi)*1] = fp690;
LIBINT2_REALTYPE fp700;
fp700 = src0[((hsi*21+8)*1+lsi)*1] * src1[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp699;
fp699 = fp700 * fp91;
LIBINT2_REALTYPE fp163;
fp163 = 2.0000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp162;
fp162 = fp163 * src1[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp166;
fp166 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp165;
fp165 = fp166 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp164;
fp164 = fp165 * src1[((hsi*21+12)*1+lsi)*1];
LIBINT2_REALTYPE fp161;
fp161 = fp164 - fp162;
LIBINT2_REALTYPE fp703;
fp703 = src0[((hsi*21+8)*1+lsi)*1] * fp161;
LIBINT2_REALTYPE fp702;
fp702 = fp703 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp705;
fp705 = fp227 * src1[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp704;
fp704 = fp705 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp701;
fp701 = fp704 + fp702;
LIBINT2_REALTYPE fp698;
fp698 = fp701 + fp699;
target[((hsi*63+10)*1+lsi)*1] = fp698;
LIBINT2_REALTYPE fp708;
fp708 = src0[((hsi*21+9)*1+lsi)*1] * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp707;
fp707 = fp708 * fp73;
LIBINT2_REALTYPE fp711;
fp711 = src0[((hsi*21+9)*1+lsi)*1] * fp185;
LIBINT2_REALTYPE fp710;
fp710 = fp711 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp239;
fp239 = 1.0000000000000000e+00 * src0[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp242;
fp242 = 1.0000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp241;
fp241 = fp242 * src0[((hsi*21+15)*1+lsi)*1];
LIBINT2_REALTYPE fp245;
fp245 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp244;
fp244 = fp245 * src0[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp248;
fp248 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp247;
fp247 = fp248 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp246;
fp246 = fp247 * src0[((hsi*21+17)*1+lsi)*1];
LIBINT2_REALTYPE fp243;
fp243 = fp246 - fp244;
LIBINT2_REALTYPE fp240;
fp240 = fp243 - fp241;
LIBINT2_REALTYPE fp238;
fp238 = fp240 + fp239;
LIBINT2_REALTYPE fp713;
fp713 = fp238 * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp712;
fp712 = fp713 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp709;
fp709 = fp712 + fp710;
LIBINT2_REALTYPE fp706;
fp706 = fp709 + fp707;
target[((hsi*63+9)*1+lsi)*1] = fp706;
LIBINT2_REALTYPE fp716;
fp716 = src0[((hsi*21+9)*1+lsi)*1] * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp715;
fp715 = fp716 * fp79;
LIBINT2_REALTYPE fp719;
fp719 = src0[((hsi*21+9)*1+lsi)*1] * fp179;
LIBINT2_REALTYPE fp718;
fp718 = fp719 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp721;
fp721 = fp238 * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp720;
fp720 = fp721 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp717;
fp717 = fp720 + fp718;
LIBINT2_REALTYPE fp714;
fp714 = fp717 + fp715;
target[((hsi*63+8)*1+lsi)*1] = fp714;
LIBINT2_REALTYPE fp724;
fp724 = src0[((hsi*21+9)*1+lsi)*1] * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp723;
fp723 = fp724 * fp85;
LIBINT2_REALTYPE fp727;
fp727 = src0[((hsi*21+9)*1+lsi)*1] * fp173;
LIBINT2_REALTYPE fp726;
fp726 = fp727 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp729;
fp729 = fp238 * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp728;
fp728 = fp729 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp725;
fp725 = fp728 + fp726;
LIBINT2_REALTYPE fp722;
fp722 = fp725 + fp723;
target[((hsi*63+7)*1+lsi)*1] = fp722;
LIBINT2_REALTYPE fp732;
fp732 = src0[((hsi*21+9)*1+lsi)*1] * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp731;
fp731 = fp732 * fp91;
LIBINT2_REALTYPE fp735;
fp735 = src0[((hsi*21+9)*1+lsi)*1] * fp167;
LIBINT2_REALTYPE fp734;
fp734 = fp735 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp737;
fp737 = fp238 * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp736;
fp736 = fp737 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp733;
fp733 = fp736 + fp734;
LIBINT2_REALTYPE fp730;
fp730 = fp733 + fp731;
target[((hsi*63+6)*1+lsi)*1] = fp730;
LIBINT2_REALTYPE fp740;
fp740 = src0[((hsi*21+10)*1+lsi)*1] * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp739;
fp739 = fp740 * fp79;
LIBINT2_REALTYPE fp743;
fp743 = src0[((hsi*21+10)*1+lsi)*1] * fp185;
LIBINT2_REALTYPE fp742;
fp742 = fp743 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp250;
fp250 = 1.5000000000000000e+00 * src0[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp253;
fp253 = 1.5000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp252;
fp252 = fp253 * src0[((hsi*21+16)*1+lsi)*1];
LIBINT2_REALTYPE fp256;
fp256 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp255;
fp255 = fp256 * src0[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp259;
fp259 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp258;
fp258 = fp259 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp257;
fp257 = fp258 * src0[((hsi*21+18)*1+lsi)*1];
LIBINT2_REALTYPE fp254;
fp254 = fp257 - fp255;
LIBINT2_REALTYPE fp251;
fp251 = fp254 - fp252;
LIBINT2_REALTYPE fp249;
fp249 = fp251 + fp250;
LIBINT2_REALTYPE fp745;
fp745 = fp249 * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp744;
fp744 = fp745 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp741;
fp741 = fp744 + fp742;
LIBINT2_REALTYPE fp738;
fp738 = fp741 + fp739;
target[((hsi*63+5)*1+lsi)*1] = fp738;
LIBINT2_REALTYPE fp748;
fp748 = src0[((hsi*21+10)*1+lsi)*1] * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp747;
fp747 = fp748 * fp85;
LIBINT2_REALTYPE fp751;
fp751 = src0[((hsi*21+10)*1+lsi)*1] * fp179;
LIBINT2_REALTYPE fp750;
fp750 = fp751 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp753;
fp753 = fp249 * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp752;
fp752 = fp753 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp749;
fp749 = fp752 + fp750;
LIBINT2_REALTYPE fp746;
fp746 = fp749 + fp747;
target[((hsi*63+4)*1+lsi)*1] = fp746;
LIBINT2_REALTYPE fp756;
fp756 = src0[((hsi*21+10)*1+lsi)*1] * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp755;
fp755 = fp756 * fp91;
LIBINT2_REALTYPE fp759;
fp759 = src0[((hsi*21+10)*1+lsi)*1] * fp173;
LIBINT2_REALTYPE fp758;
fp758 = fp759 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp761;
fp761 = fp249 * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp760;
fp760 = fp761 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp757;
fp757 = fp760 + fp758;
LIBINT2_REALTYPE fp754;
fp754 = fp757 + fp755;
target[((hsi*63+3)*1+lsi)*1] = fp754;
LIBINT2_REALTYPE fp764;
fp764 = src0[((hsi*21+11)*1+lsi)*1] * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp763;
fp763 = fp764 * fp85;
LIBINT2_REALTYPE fp767;
fp767 = src0[((hsi*21+11)*1+lsi)*1] * fp185;
LIBINT2_REALTYPE fp766;
fp766 = fp767 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp261;
fp261 = 2.0000000000000000e+00 * src0[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp264;
fp264 = 2.0000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp263;
fp263 = fp264 * src0[((hsi*21+17)*1+lsi)*1];
LIBINT2_REALTYPE fp267;
fp267 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp266;
fp266 = fp267 * src0[((hsi*21+5)*1+lsi)*1];
LIBINT2_REALTYPE fp270;
fp270 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp269;
fp269 = fp270 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp268;
fp268 = fp269 * src0[((hsi*21+19)*1+lsi)*1];
LIBINT2_REALTYPE fp265;
fp265 = fp268 - fp266;
LIBINT2_REALTYPE fp262;
fp262 = fp265 - fp263;
LIBINT2_REALTYPE fp260;
fp260 = fp262 + fp261;
LIBINT2_REALTYPE fp769;
fp769 = fp260 * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp768;
fp768 = fp769 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp765;
fp765 = fp768 + fp766;
LIBINT2_REALTYPE fp762;
fp762 = fp765 + fp763;
target[((hsi*63+2)*1+lsi)*1] = fp762;
LIBINT2_REALTYPE fp772;
fp772 = src0[((hsi*21+11)*1+lsi)*1] * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp771;
fp771 = fp772 * fp91;
LIBINT2_REALTYPE fp775;
fp775 = src0[((hsi*21+11)*1+lsi)*1] * fp179;
LIBINT2_REALTYPE fp774;
fp774 = fp775 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp777;
fp777 = fp260 * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp776;
fp776 = fp777 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp773;
fp773 = fp776 + fp774;
LIBINT2_REALTYPE fp770;
fp770 = fp773 + fp771;
target[((hsi*63+1)*1+lsi)*1] = fp770;
LIBINT2_REALTYPE fp652;
fp652 = src0[((hsi*21+7)*1+lsi)*1] * src1[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp651;
fp651 = fp652 * fp85;
LIBINT2_REALTYPE fp655;
fp655 = src0[((hsi*21+7)*1+lsi)*1] * fp161;
LIBINT2_REALTYPE fp654;
fp654 = fp655 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp657;
fp657 = fp221 * src1[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp656;
fp656 = fp657 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp653;
fp653 = fp656 + fp654;
LIBINT2_REALTYPE fp650;
fp650 = fp653 + fp651;
target[((hsi*63+16)*1+lsi)*1] = fp650;
LIBINT2_REALTYPE fp1;
fp1 = 2.5000000000000000e+00 * src2[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp4;
fp4 = 2.5000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp3;
fp3 = fp4 * src2[((hsi*21+18)*1+lsi)*1];
LIBINT2_REALTYPE fp7;
fp7 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp6;
fp6 = fp7 * src2[((hsi*21+6)*1+lsi)*1];
LIBINT2_REALTYPE fp10;
fp10 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp9;
fp9 = fp10 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp8;
fp8 = fp9 * src2[((hsi*21+20)*1+lsi)*1];
LIBINT2_REALTYPE fp5;
fp5 = fp8 - fp6;
LIBINT2_REALTYPE fp2;
fp2 = fp5 - fp3;
LIBINT2_REALTYPE fp0;
fp0 = fp2 + fp1;
LIBINT2_REALTYPE fp284;
fp284 = src0[((hsi*21+0)*1+lsi)*1] * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp283;
fp283 = fp284 * fp0;
LIBINT2_REALTYPE fp287;
fp287 = src0[((hsi*21+0)*1+lsi)*1] * fp185;
LIBINT2_REALTYPE fp286;
fp286 = fp287 * src2[((hsi*21+12)*1+lsi)*1];
LIBINT2_REALTYPE fp190;
fp190 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp189;
fp189 = fp190 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp188;
fp188 = fp189 * src0[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp289;
fp289 = fp188 * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp288;
fp288 = fp289 * src2[((hsi*21+12)*1+lsi)*1];
LIBINT2_REALTYPE fp285;
fp285 = fp288 + fp286;
LIBINT2_REALTYPE fp282;
fp282 = fp285 + fp283;
target[((hsi*63+62)*1+lsi)*1] = fp282;
LIBINT2_REALTYPE fp12;
fp12 = 2.0000000000000000e+00 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp15;
fp15 = 2.0000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp14;
fp14 = fp15 * src2[((hsi*21+17)*1+lsi)*1];
LIBINT2_REALTYPE fp18;
fp18 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp17;
fp17 = fp18 * src2[((hsi*21+5)*1+lsi)*1];
LIBINT2_REALTYPE fp21;
fp21 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp20;
fp20 = fp21 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp19;
fp19 = fp20 * src2[((hsi*21+19)*1+lsi)*1];
LIBINT2_REALTYPE fp16;
fp16 = fp19 - fp17;
LIBINT2_REALTYPE fp13;
fp13 = fp16 - fp14;
LIBINT2_REALTYPE fp11;
fp11 = fp13 + fp12;
LIBINT2_REALTYPE fp292;
fp292 = src0[((hsi*21+0)*1+lsi)*1] * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp291;
fp291 = fp292 * fp11;
LIBINT2_REALTYPE fp295;
fp295 = src0[((hsi*21+0)*1+lsi)*1] * fp179;
LIBINT2_REALTYPE fp294;
fp294 = fp295 * src2[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp297;
fp297 = fp188 * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp296;
fp296 = fp297 * src2[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp293;
fp293 = fp296 + fp294;
LIBINT2_REALTYPE fp290;
fp290 = fp293 + fp291;
target[((hsi*63+61)*1+lsi)*1] = fp290;
LIBINT2_REALTYPE fp23;
fp23 = 1.5000000000000000e+00 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp26;
fp26 = 1.5000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp25;
fp25 = fp26 * src2[((hsi*21+16)*1+lsi)*1];
LIBINT2_REALTYPE fp29;
fp29 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp28;
fp28 = fp29 * src2[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp32;
fp32 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp31;
fp31 = fp32 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp30;
fp30 = fp31 * src2[((hsi*21+18)*1+lsi)*1];
LIBINT2_REALTYPE fp27;
fp27 = fp30 - fp28;
LIBINT2_REALTYPE fp24;
fp24 = fp27 - fp25;
LIBINT2_REALTYPE fp22;
fp22 = fp24 + fp23;
LIBINT2_REALTYPE fp300;
fp300 = src0[((hsi*21+0)*1+lsi)*1] * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp299;
fp299 = fp300 * fp22;
LIBINT2_REALTYPE fp303;
fp303 = src0[((hsi*21+0)*1+lsi)*1] * fp173;
LIBINT2_REALTYPE fp302;
fp302 = fp303 * src2[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp305;
fp305 = fp188 * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp304;
fp304 = fp305 * src2[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp301;
fp301 = fp304 + fp302;
LIBINT2_REALTYPE fp298;
fp298 = fp301 + fp299;
target[((hsi*63+60)*1+lsi)*1] = fp298;
LIBINT2_REALTYPE fp34;
fp34 = 1.0000000000000000e+00 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp37;
fp37 = 1.0000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp36;
fp36 = fp37 * src2[((hsi*21+15)*1+lsi)*1];
LIBINT2_REALTYPE fp40;
fp40 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp39;
fp39 = fp40 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp43;
fp43 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp42;
fp42 = fp43 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp41;
fp41 = fp42 * src2[((hsi*21+17)*1+lsi)*1];
LIBINT2_REALTYPE fp38;
fp38 = fp41 - fp39;
LIBINT2_REALTYPE fp35;
fp35 = fp38 - fp36;
LIBINT2_REALTYPE fp33;
fp33 = fp35 + fp34;
LIBINT2_REALTYPE fp308;
fp308 = src0[((hsi*21+0)*1+lsi)*1] * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp307;
fp307 = fp308 * fp33;
LIBINT2_REALTYPE fp311;
fp311 = src0[((hsi*21+0)*1+lsi)*1] * fp167;
LIBINT2_REALTYPE fp310;
fp310 = fp311 * src2[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp313;
fp313 = fp188 * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp312;
fp312 = fp313 * src2[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp309;
fp309 = fp312 + fp310;
LIBINT2_REALTYPE fp306;
fp306 = fp309 + fp307;
target[((hsi*63+59)*1+lsi)*1] = fp306;
LIBINT2_REALTYPE fp45;
fp45 = 5.0000000000000000e-01 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp48;
fp48 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp47;
fp47 = fp48 * src2[((hsi*21+14)*1+lsi)*1];
LIBINT2_REALTYPE fp51;
fp51 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp50;
fp50 = fp51 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp54;
fp54 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp53;
fp53 = fp54 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp52;
fp52 = fp53 * src2[((hsi*21+16)*1+lsi)*1];
LIBINT2_REALTYPE fp49;
fp49 = fp52 - fp50;
LIBINT2_REALTYPE fp46;
fp46 = fp49 - fp47;
LIBINT2_REALTYPE fp44;
fp44 = fp46 + fp45;
LIBINT2_REALTYPE fp316;
fp316 = src0[((hsi*21+0)*1+lsi)*1] * src1[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp315;
fp315 = fp316 * fp44;
LIBINT2_REALTYPE fp319;
fp319 = src0[((hsi*21+0)*1+lsi)*1] * fp161;
LIBINT2_REALTYPE fp318;
fp318 = fp319 * src2[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp321;
fp321 = fp188 * src1[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp320;
fp320 = fp321 * src2[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp317;
fp317 = fp320 + fp318;
LIBINT2_REALTYPE fp314;
fp314 = fp317 + fp315;
target[((hsi*63+58)*1+lsi)*1] = fp314;
LIBINT2_REALTYPE fp57;
fp57 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp56;
fp56 = fp57 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp60;
fp60 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp59;
fp59 = fp60 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp58;
fp58 = fp59 * src2[((hsi*21+15)*1+lsi)*1];
LIBINT2_REALTYPE fp55;
fp55 = fp58 - fp56;
LIBINT2_REALTYPE fp324;
fp324 = src0[((hsi*21+0)*1+lsi)*1] * src1[((hsi*21+5)*1+lsi)*1];
LIBINT2_REALTYPE fp323;
fp323 = fp324 * fp55;
LIBINT2_REALTYPE fp327;
fp327 = src0[((hsi*21+0)*1+lsi)*1] * fp155;
LIBINT2_REALTYPE fp326;
fp326 = fp327 * src2[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp329;
fp329 = fp188 * src1[((hsi*21+5)*1+lsi)*1];
LIBINT2_REALTYPE fp328;
fp328 = fp329 * src2[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp325;
fp325 = fp328 + fp326;
LIBINT2_REALTYPE fp322;
fp322 = fp325 + fp323;
target[((hsi*63+57)*1+lsi)*1] = fp322;
LIBINT2_REALTYPE fp332;
fp332 = src0[((hsi*21+1)*1+lsi)*1] * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp331;
fp331 = fp332 * fp11;
LIBINT2_REALTYPE fp335;
fp335 = src0[((hsi*21+1)*1+lsi)*1] * fp185;
LIBINT2_REALTYPE fp334;
fp334 = fp335 * src2[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp193;
fp193 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp192;
fp192 = fp193 * src0[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp196;
fp196 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp195;
fp195 = fp196 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp194;
fp194 = fp195 * src0[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp191;
fp191 = fp194 - fp192;
LIBINT2_REALTYPE fp337;
fp337 = fp191 * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp336;
fp336 = fp337 * src2[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp333;
fp333 = fp336 + fp334;
LIBINT2_REALTYPE fp330;
fp330 = fp333 + fp331;
target[((hsi*63+56)*1+lsi)*1] = fp330;
LIBINT2_REALTYPE fp340;
fp340 = src0[((hsi*21+1)*1+lsi)*1] * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp339;
fp339 = fp340 * fp22;
LIBINT2_REALTYPE fp343;
fp343 = src0[((hsi*21+1)*1+lsi)*1] * fp179;
LIBINT2_REALTYPE fp342;
fp342 = fp343 * src2[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp345;
fp345 = fp191 * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp344;
fp344 = fp345 * src2[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp341;
fp341 = fp344 + fp342;
LIBINT2_REALTYPE fp338;
fp338 = fp341 + fp339;
target[((hsi*63+55)*1+lsi)*1] = fp338;
LIBINT2_REALTYPE fp348;
fp348 = src0[((hsi*21+1)*1+lsi)*1] * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp347;
fp347 = fp348 * fp33;
LIBINT2_REALTYPE fp351;
fp351 = src0[((hsi*21+1)*1+lsi)*1] * fp173;
LIBINT2_REALTYPE fp350;
fp350 = fp351 * src2[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp353;
fp353 = fp191 * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp352;
fp352 = fp353 * src2[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp349;
fp349 = fp352 + fp350;
LIBINT2_REALTYPE fp346;
fp346 = fp349 + fp347;
target[((hsi*63+54)*1+lsi)*1] = fp346;
LIBINT2_REALTYPE fp356;
fp356 = src0[((hsi*21+1)*1+lsi)*1] * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp355;
fp355 = fp356 * fp44;
LIBINT2_REALTYPE fp359;
fp359 = src0[((hsi*21+1)*1+lsi)*1] * fp167;
LIBINT2_REALTYPE fp358;
fp358 = fp359 * src2[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp361;
fp361 = fp191 * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp360;
fp360 = fp361 * src2[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp357;
fp357 = fp360 + fp358;
LIBINT2_REALTYPE fp354;
fp354 = fp357 + fp355;
target[((hsi*63+53)*1+lsi)*1] = fp354;
LIBINT2_REALTYPE fp364;
fp364 = src0[((hsi*21+1)*1+lsi)*1] * src1[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp363;
fp363 = fp364 * fp55;
LIBINT2_REALTYPE fp367;
fp367 = src0[((hsi*21+1)*1+lsi)*1] * fp161;
LIBINT2_REALTYPE fp366;
fp366 = fp367 * src2[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp369;
fp369 = fp191 * src1[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp368;
fp368 = fp369 * src2[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp365;
fp365 = fp368 + fp366;
LIBINT2_REALTYPE fp362;
fp362 = fp365 + fp363;
target[((hsi*63+52)*1+lsi)*1] = fp362;
LIBINT2_REALTYPE fp372;
fp372 = src0[((hsi*21+2)*1+lsi)*1] * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp371;
fp371 = fp372 * fp22;
LIBINT2_REALTYPE fp375;
fp375 = src0[((hsi*21+2)*1+lsi)*1] * fp185;
LIBINT2_REALTYPE fp374;
fp374 = fp375 * src2[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp377;
fp377 = fp197 * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp376;
fp376 = fp377 * src2[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp373;
fp373 = fp376 + fp374;
LIBINT2_REALTYPE fp370;
fp370 = fp373 + fp371;
target[((hsi*63+51)*1+lsi)*1] = fp370;
LIBINT2_REALTYPE fp380;
fp380 = src0[((hsi*21+2)*1+lsi)*1] * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp379;
fp379 = fp380 * fp33;
LIBINT2_REALTYPE fp383;
fp383 = src0[((hsi*21+2)*1+lsi)*1] * fp179;
LIBINT2_REALTYPE fp382;
fp382 = fp383 * src2[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp385;
fp385 = fp197 * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp384;
fp384 = fp385 * src2[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp381;
fp381 = fp384 + fp382;
LIBINT2_REALTYPE fp378;
fp378 = fp381 + fp379;
target[((hsi*63+50)*1+lsi)*1] = fp378;
LIBINT2_REALTYPE fp388;
fp388 = src0[((hsi*21+2)*1+lsi)*1] * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp387;
fp387 = fp388 * fp44;
LIBINT2_REALTYPE fp391;
fp391 = src0[((hsi*21+2)*1+lsi)*1] * fp173;
LIBINT2_REALTYPE fp390;
fp390 = fp391 * src2[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp393;
fp393 = fp197 * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp392;
fp392 = fp393 * src2[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp389;
fp389 = fp392 + fp390;
LIBINT2_REALTYPE fp386;
fp386 = fp389 + fp387;
target[((hsi*63+49)*1+lsi)*1] = fp386;
LIBINT2_REALTYPE fp396;
fp396 = src0[((hsi*21+2)*1+lsi)*1] * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp395;
fp395 = fp396 * fp55;
LIBINT2_REALTYPE fp399;
fp399 = src0[((hsi*21+2)*1+lsi)*1] * fp167;
LIBINT2_REALTYPE fp398;
fp398 = fp399 * src2[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp401;
fp401 = fp197 * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp400;
fp400 = fp401 * src2[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp397;
fp397 = fp400 + fp398;
LIBINT2_REALTYPE fp394;
fp394 = fp397 + fp395;
target[((hsi*63+48)*1+lsi)*1] = fp394;
LIBINT2_REALTYPE fp532;
fp532 = src0[((hsi*21+1)*1+lsi)*1] * src1[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp531;
fp531 = fp532 * fp91;
LIBINT2_REALTYPE fp106;
fp106 = 2.0000000000000000e+00 * src1[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp109;
fp109 = 2.0000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp108;
fp108 = fp109 * src1[((hsi*21+17)*1+lsi)*1];
LIBINT2_REALTYPE fp112;
fp112 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp111;
fp111 = fp112 * src1[((hsi*21+5)*1+lsi)*1];
LIBINT2_REALTYPE fp115;
fp115 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp114;
fp114 = fp115 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp113;
fp113 = fp114 * src1[((hsi*21+19)*1+lsi)*1];
LIBINT2_REALTYPE fp110;
fp110 = fp113 - fp111;
LIBINT2_REALTYPE fp107;
fp107 = fp110 - fp108;
LIBINT2_REALTYPE fp105;
fp105 = fp107 + fp106;
LIBINT2_REALTYPE fp535;
fp535 = src0[((hsi*21+1)*1+lsi)*1] * fp105;
LIBINT2_REALTYPE fp534;
fp534 = fp535 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp537;
fp537 = fp191 * src1[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp536;
fp536 = fp537 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp533;
fp533 = fp536 + fp534;
LIBINT2_REALTYPE fp530;
fp530 = fp533 + fp531;
target[((hsi*63+31)*1+lsi)*1] = fp530;
LIBINT2_REALTYPE fp412;
fp412 = src0[((hsi*21+3)*1+lsi)*1] * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp411;
fp411 = fp412 * fp44;
LIBINT2_REALTYPE fp415;
fp415 = src0[((hsi*21+3)*1+lsi)*1] * fp179;
LIBINT2_REALTYPE fp414;
fp414 = fp415 * src2[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp417;
fp417 = fp203 * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp416;
fp416 = fp417 * src2[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp413;
fp413 = fp416 + fp414;
LIBINT2_REALTYPE fp410;
fp410 = fp413 + fp411;
target[((hsi*63+46)*1+lsi)*1] = fp410;
LIBINT2_REALTYPE fp420;
fp420 = src0[((hsi*21+3)*1+lsi)*1] * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp419;
fp419 = fp420 * fp55;
LIBINT2_REALTYPE fp423;
fp423 = src0[((hsi*21+3)*1+lsi)*1] * fp173;
LIBINT2_REALTYPE fp422;
fp422 = fp423 * src2[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp425;
fp425 = fp203 * src1[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp424;
fp424 = fp425 * src2[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp421;
fp421 = fp424 + fp422;
LIBINT2_REALTYPE fp418;
fp418 = fp421 + fp419;
target[((hsi*63+45)*1+lsi)*1] = fp418;
LIBINT2_REALTYPE fp428;
fp428 = src0[((hsi*21+4)*1+lsi)*1] * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp427;
fp427 = fp428 * fp44;
LIBINT2_REALTYPE fp431;
fp431 = src0[((hsi*21+4)*1+lsi)*1] * fp185;
LIBINT2_REALTYPE fp430;
fp430 = fp431 * src2[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp433;
fp433 = fp209 * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp432;
fp432 = fp433 * src2[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp429;
fp429 = fp432 + fp430;
LIBINT2_REALTYPE fp426;
fp426 = fp429 + fp427;
target[((hsi*63+44)*1+lsi)*1] = fp426;
LIBINT2_REALTYPE fp436;
fp436 = src0[((hsi*21+4)*1+lsi)*1] * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp435;
fp435 = fp436 * fp55;
LIBINT2_REALTYPE fp439;
fp439 = src0[((hsi*21+4)*1+lsi)*1] * fp179;
LIBINT2_REALTYPE fp438;
fp438 = fp439 * src2[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp441;
fp441 = fp209 * src1[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp440;
fp440 = fp441 * src2[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp437;
fp437 = fp440 + fp438;
LIBINT2_REALTYPE fp434;
fp434 = fp437 + fp435;
target[((hsi*63+43)*1+lsi)*1] = fp434;
LIBINT2_REALTYPE fp444;
fp444 = src0[((hsi*21+5)*1+lsi)*1] * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp443;
fp443 = fp444 * fp55;
LIBINT2_REALTYPE fp447;
fp447 = src0[((hsi*21+5)*1+lsi)*1] * fp185;
LIBINT2_REALTYPE fp446;
fp446 = fp447 * src2[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp449;
fp449 = fp215 * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp448;
fp448 = fp449 * src2[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp445;
fp445 = fp448 + fp446;
LIBINT2_REALTYPE fp442;
fp442 = fp445 + fp443;
target[((hsi*63+42)*1+lsi)*1] = fp442;
LIBINT2_REALTYPE fp452;
fp452 = src0[((hsi*21+0)*1+lsi)*1] * src1[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp451;
fp451 = fp452 * fp61;
LIBINT2_REALTYPE fp455;
fp455 = src0[((hsi*21+0)*1+lsi)*1] * fp149;
LIBINT2_REALTYPE fp454;
fp454 = fp455 * src2[((hsi*21+5)*1+lsi)*1];
LIBINT2_REALTYPE fp457;
fp457 = fp188 * src1[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp456;
fp456 = fp457 * src2[((hsi*21+5)*1+lsi)*1];
LIBINT2_REALTYPE fp453;
fp453 = fp456 + fp454;
LIBINT2_REALTYPE fp450;
fp450 = fp453 + fp451;
target[((hsi*63+41)*1+lsi)*1] = fp450;
LIBINT2_REALTYPE fp460;
fp460 = src0[((hsi*21+0)*1+lsi)*1] * src1[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp459;
fp459 = fp460 * fp67;
LIBINT2_REALTYPE fp463;
fp463 = src0[((hsi*21+0)*1+lsi)*1] * fp138;
LIBINT2_REALTYPE fp462;
fp462 = fp463 * src2[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp465;
fp465 = fp188 * src1[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp464;
fp464 = fp465 * src2[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp461;
fp461 = fp464 + fp462;
LIBINT2_REALTYPE fp458;
fp458 = fp461 + fp459;
target[((hsi*63+40)*1+lsi)*1] = fp458;
LIBINT2_REALTYPE fp468;
fp468 = src0[((hsi*21+0)*1+lsi)*1] * src1[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp467;
fp467 = fp468 * fp73;
LIBINT2_REALTYPE fp471;
fp471 = src0[((hsi*21+0)*1+lsi)*1] * fp127;
LIBINT2_REALTYPE fp470;
fp470 = fp471 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp473;
fp473 = fp188 * src1[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp472;
fp472 = fp473 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp469;
fp469 = fp472 + fp470;
LIBINT2_REALTYPE fp466;
fp466 = fp469 + fp467;
target[((hsi*63+39)*1+lsi)*1] = fp466;
LIBINT2_REALTYPE fp476;
fp476 = src0[((hsi*21+0)*1+lsi)*1] * src1[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp475;
fp475 = fp476 * fp79;
LIBINT2_REALTYPE fp479;
fp479 = src0[((hsi*21+0)*1+lsi)*1] * fp116;
LIBINT2_REALTYPE fp478;
fp478 = fp479 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp481;
fp481 = fp188 * src1[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp480;
fp480 = fp481 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp477;
fp477 = fp480 + fp478;
LIBINT2_REALTYPE fp474;
fp474 = fp477 + fp475;
target[((hsi*63+38)*1+lsi)*1] = fp474;
LIBINT2_REALTYPE fp484;
fp484 = src0[((hsi*21+0)*1+lsi)*1] * src1[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp483;
fp483 = fp484 * fp85;
LIBINT2_REALTYPE fp487;
fp487 = src0[((hsi*21+0)*1+lsi)*1] * fp105;
LIBINT2_REALTYPE fp486;
fp486 = fp487 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp489;
fp489 = fp188 * src1[((hsi*21+11)*1+lsi)*1];
LIBINT2_REALTYPE fp488;
fp488 = fp489 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp485;
fp485 = fp488 + fp486;
LIBINT2_REALTYPE fp482;
fp482 = fp485 + fp483;
target[((hsi*63+37)*1+lsi)*1] = fp482;
LIBINT2_REALTYPE fp492;
fp492 = src0[((hsi*21+0)*1+lsi)*1] * src1[((hsi*21+12)*1+lsi)*1];
LIBINT2_REALTYPE fp491;
fp491 = fp492 * fp91;
LIBINT2_REALTYPE fp95;
fp95 = 2.5000000000000000e+00 * src1[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp98;
fp98 = 2.5000000000000000e+00 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp97;
fp97 = fp98 * src1[((hsi*21+18)*1+lsi)*1];
LIBINT2_REALTYPE fp101;
fp101 = 5.0000000000000000e-01 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp100;
fp100 = fp101 * src1[((hsi*21+6)*1+lsi)*1];
LIBINT2_REALTYPE fp104;
fp104 = 5.0000000000000000e-01 * inteval->two_alpha0_bra[vi];
LIBINT2_REALTYPE fp103;
fp103 = fp104 * inteval->two_alpha0_ket[vi];
LIBINT2_REALTYPE fp102;
fp102 = fp103 * src1[((hsi*21+20)*1+lsi)*1];
LIBINT2_REALTYPE fp99;
fp99 = fp102 - fp100;
LIBINT2_REALTYPE fp96;
fp96 = fp99 - fp97;
LIBINT2_REALTYPE fp94;
fp94 = fp96 + fp95;
LIBINT2_REALTYPE fp495;
fp495 = src0[((hsi*21+0)*1+lsi)*1] * fp94;
LIBINT2_REALTYPE fp494;
fp494 = fp495 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp497;
fp497 = fp188 * src1[((hsi*21+12)*1+lsi)*1];
LIBINT2_REALTYPE fp496;
fp496 = fp497 * src2[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp493;
fp493 = fp496 + fp494;
LIBINT2_REALTYPE fp490;
fp490 = fp493 + fp491;
target[((hsi*63+36)*1+lsi)*1] = fp490;
LIBINT2_REALTYPE fp500;
fp500 = src0[((hsi*21+1)*1+lsi)*1] * src1[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp499;
fp499 = fp500 * fp67;
LIBINT2_REALTYPE fp503;
fp503 = src0[((hsi*21+1)*1+lsi)*1] * fp149;
LIBINT2_REALTYPE fp502;
fp502 = fp503 * src2[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp505;
fp505 = fp191 * src1[((hsi*21+7)*1+lsi)*1];
LIBINT2_REALTYPE fp504;
fp504 = fp505 * src2[((hsi*21+4)*1+lsi)*1];
LIBINT2_REALTYPE fp501;
fp501 = fp504 + fp502;
LIBINT2_REALTYPE fp498;
fp498 = fp501 + fp499;
target[((hsi*63+35)*1+lsi)*1] = fp498;
LIBINT2_REALTYPE fp508;
fp508 = src0[((hsi*21+1)*1+lsi)*1] * src1[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp507;
fp507 = fp508 * fp73;
LIBINT2_REALTYPE fp511;
fp511 = src0[((hsi*21+1)*1+lsi)*1] * fp138;
LIBINT2_REALTYPE fp510;
fp510 = fp511 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp513;
fp513 = fp191 * src1[((hsi*21+8)*1+lsi)*1];
LIBINT2_REALTYPE fp512;
fp512 = fp513 * src2[((hsi*21+3)*1+lsi)*1];
LIBINT2_REALTYPE fp509;
fp509 = fp512 + fp510;
LIBINT2_REALTYPE fp506;
fp506 = fp509 + fp507;
target[((hsi*63+34)*1+lsi)*1] = fp506;
LIBINT2_REALTYPE fp516;
fp516 = src0[((hsi*21+1)*1+lsi)*1] * src1[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp515;
fp515 = fp516 * fp79;
LIBINT2_REALTYPE fp519;
fp519 = src0[((hsi*21+1)*1+lsi)*1] * fp127;
LIBINT2_REALTYPE fp518;
fp518 = fp519 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp521;
fp521 = fp191 * src1[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp520;
fp520 = fp521 * src2[((hsi*21+2)*1+lsi)*1];
LIBINT2_REALTYPE fp517;
fp517 = fp520 + fp518;
LIBINT2_REALTYPE fp514;
fp514 = fp517 + fp515;
target[((hsi*63+33)*1+lsi)*1] = fp514;
LIBINT2_REALTYPE fp524;
fp524 = src0[((hsi*21+1)*1+lsi)*1] * src1[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp523;
fp523 = fp524 * fp85;
LIBINT2_REALTYPE fp527;
fp527 = src0[((hsi*21+1)*1+lsi)*1] * fp116;
LIBINT2_REALTYPE fp526;
fp526 = fp527 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp529;
fp529 = fp191 * src1[((hsi*21+10)*1+lsi)*1];
LIBINT2_REALTYPE fp528;
fp528 = fp529 * src2[((hsi*21+1)*1+lsi)*1];
LIBINT2_REALTYPE fp525;
fp525 = fp528 + fp526;
LIBINT2_REALTYPE fp522;
fp522 = fp525 + fp523;
target[((hsi*63+32)*1+lsi)*1] = fp522;
LIBINT2_REALTYPE fp404;
fp404 = src0[((hsi*21+3)*1+lsi)*1] * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp403;
fp403 = fp404 * fp33;
LIBINT2_REALTYPE fp407;
fp407 = src0[((hsi*21+3)*1+lsi)*1] * fp185;
LIBINT2_REALTYPE fp406;
fp406 = fp407 * src2[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp409;
fp409 = fp203 * src1[((hsi*21+0)*1+lsi)*1];
LIBINT2_REALTYPE fp408;
fp408 = fp409 * src2[((hsi*21+9)*1+lsi)*1];
LIBINT2_REALTYPE fp405;
fp405 = fp408 + fp406;
LIBINT2_REALTYPE fp402;
fp402 = fp405 + fp403;
target[((hsi*63+47)*1+lsi)*1] = fp402;
}
}
}
const int hsi = 0;
const int lsi = 0;
const int vi = 0;
/** Number of flops = 786 */
}

#ifdef __cplusplus
};
LIBINT_PRAGMA_CLANG(diagnostic pop)
LIBINT_PRAGMA_GCC(diagnostic pop)
#endif
