/*
 *  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 OSVRRP1InBra_aB_s001__0__h__1___TwoPRep_unit__0__unit__1___Ab__up_0(const Libint_t* inteval, LIBINT2_REALTYPE* target, const LIBINT2_REALTYPE* src0, const LIBINT2_REALTYPE* src1, const LIBINT2_REALTYPE* src2, const LIBINT2_REALTYPE* src3) {

LIBINT2_REALTYPE*const  stack = target;
{
const int hsi = 0;
{
const int lsi = 0;
{
const int vi = 0;
LIBINT2_REALTYPE fp134;
fp134 = inteval->WQ_y[vi] * src0[((hsi*15+0)*1+lsi)*1];
target[((hsi*21+1)*1+lsi)*1] = fp134;
LIBINT2_REALTYPE fp132;
fp132 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp131;
fp131 = fp132 * src3[((hsi*15+0)*1+lsi)*1];
LIBINT2_REALTYPE fp133;
fp133 = inteval->WQ_z[vi] * src0[((hsi*15+0)*1+lsi)*1];
LIBINT2_REALTYPE fp130;
fp130 = fp133 + fp131;
target[((hsi*21+2)*1+lsi)*1] = fp130;
LIBINT2_REALTYPE fp127;
fp127 = inteval->roe[vi] * src2[((hsi*10+0)*1+lsi)*1];
LIBINT2_REALTYPE fp126;
fp126 = src1[((hsi*10+0)*1+lsi)*1] - fp127;
LIBINT2_REALTYPE fp128;
fp128 = 1.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp125;
fp125 = fp128 * fp126;
LIBINT2_REALTYPE fp129;
fp129 = inteval->WQ_y[vi] * src0[((hsi*15+1)*1+lsi)*1];
LIBINT2_REALTYPE fp124;
fp124 = fp129 + fp125;
target[((hsi*21+3)*1+lsi)*1] = fp124;
LIBINT2_REALTYPE fp122;
fp122 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp121;
fp121 = fp122 * src3[((hsi*15+1)*1+lsi)*1];
LIBINT2_REALTYPE fp123;
fp123 = inteval->WQ_z[vi] * src0[((hsi*15+1)*1+lsi)*1];
LIBINT2_REALTYPE fp120;
fp120 = fp123 + fp121;
target[((hsi*21+4)*1+lsi)*1] = fp120;
LIBINT2_REALTYPE fp113;
fp113 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp112;
fp112 = fp113 * src3[((hsi*15+2)*1+lsi)*1];
LIBINT2_REALTYPE fp117;
fp117 = inteval->roe[vi] * src2[((hsi*10+0)*1+lsi)*1];
LIBINT2_REALTYPE fp116;
fp116 = src1[((hsi*10+0)*1+lsi)*1] - fp117;
LIBINT2_REALTYPE fp118;
fp118 = 1.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp115;
fp115 = fp118 * fp116;
LIBINT2_REALTYPE fp119;
fp119 = inteval->WQ_z[vi] * src0[((hsi*15+2)*1+lsi)*1];
LIBINT2_REALTYPE fp114;
fp114 = fp119 + fp115;
LIBINT2_REALTYPE fp111;
fp111 = fp114 + fp112;
target[((hsi*21+5)*1+lsi)*1] = fp111;
LIBINT2_REALTYPE fp108;
fp108 = inteval->roe[vi] * src2[((hsi*10+1)*1+lsi)*1];
LIBINT2_REALTYPE fp107;
fp107 = src1[((hsi*10+1)*1+lsi)*1] - fp108;
LIBINT2_REALTYPE fp109;
fp109 = 2.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp106;
fp106 = fp109 * fp107;
LIBINT2_REALTYPE fp110;
fp110 = inteval->WQ_y[vi] * src0[((hsi*15+3)*1+lsi)*1];
LIBINT2_REALTYPE fp105;
fp105 = fp110 + fp106;
target[((hsi*21+6)*1+lsi)*1] = fp105;
LIBINT2_REALTYPE fp103;
fp103 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp102;
fp102 = fp103 * src3[((hsi*15+3)*1+lsi)*1];
LIBINT2_REALTYPE fp104;
fp104 = inteval->WQ_z[vi] * src0[((hsi*15+3)*1+lsi)*1];
LIBINT2_REALTYPE fp101;
fp101 = fp104 + fp102;
target[((hsi*21+7)*1+lsi)*1] = fp101;
LIBINT2_REALTYPE fp94;
fp94 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp93;
fp93 = fp94 * src3[((hsi*15+4)*1+lsi)*1];
LIBINT2_REALTYPE fp98;
fp98 = inteval->roe[vi] * src2[((hsi*10+1)*1+lsi)*1];
LIBINT2_REALTYPE fp97;
fp97 = src1[((hsi*10+1)*1+lsi)*1] - fp98;
LIBINT2_REALTYPE fp99;
fp99 = 1.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp96;
fp96 = fp99 * fp97;
LIBINT2_REALTYPE fp100;
fp100 = inteval->WQ_z[vi] * src0[((hsi*15+4)*1+lsi)*1];
LIBINT2_REALTYPE fp95;
fp95 = fp100 + fp96;
LIBINT2_REALTYPE fp92;
fp92 = fp95 + fp93;
target[((hsi*21+8)*1+lsi)*1] = fp92;
LIBINT2_REALTYPE fp85;
fp85 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp84;
fp84 = fp85 * src3[((hsi*15+5)*1+lsi)*1];
LIBINT2_REALTYPE fp89;
fp89 = inteval->roe[vi] * src2[((hsi*10+2)*1+lsi)*1];
LIBINT2_REALTYPE fp88;
fp88 = src1[((hsi*10+2)*1+lsi)*1] - fp89;
LIBINT2_REALTYPE fp90;
fp90 = 2.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp87;
fp87 = fp90 * fp88;
LIBINT2_REALTYPE fp91;
fp91 = inteval->WQ_z[vi] * src0[((hsi*15+5)*1+lsi)*1];
LIBINT2_REALTYPE fp86;
fp86 = fp91 + fp87;
LIBINT2_REALTYPE fp83;
fp83 = fp86 + fp84;
target[((hsi*21+9)*1+lsi)*1] = fp83;
LIBINT2_REALTYPE fp138;
fp138 = inteval->roe[vi] * src2[((hsi*10+0)*1+lsi)*1];
LIBINT2_REALTYPE fp137;
fp137 = src1[((hsi*10+0)*1+lsi)*1] - fp138;
LIBINT2_REALTYPE fp139;
fp139 = 4.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp136;
fp136 = fp139 * fp137;
LIBINT2_REALTYPE fp140;
fp140 = inteval->WQ_x[vi] * src0[((hsi*15+0)*1+lsi)*1];
LIBINT2_REALTYPE fp135;
fp135 = fp140 + fp136;
target[((hsi*21+0)*1+lsi)*1] = fp135;
LIBINT2_REALTYPE fp75;
fp75 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp74;
fp74 = fp75 * src3[((hsi*15+6)*1+lsi)*1];
LIBINT2_REALTYPE fp76;
fp76 = inteval->WQ_z[vi] * src0[((hsi*15+6)*1+lsi)*1];
LIBINT2_REALTYPE fp73;
fp73 = fp76 + fp74;
target[((hsi*21+11)*1+lsi)*1] = fp73;
LIBINT2_REALTYPE fp66;
fp66 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp65;
fp65 = fp66 * src3[((hsi*15+7)*1+lsi)*1];
LIBINT2_REALTYPE fp70;
fp70 = inteval->roe[vi] * src2[((hsi*10+3)*1+lsi)*1];
LIBINT2_REALTYPE fp69;
fp69 = src1[((hsi*10+3)*1+lsi)*1] - fp70;
LIBINT2_REALTYPE fp71;
fp71 = 1.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp68;
fp68 = fp71 * fp69;
LIBINT2_REALTYPE fp72;
fp72 = inteval->WQ_z[vi] * src0[((hsi*15+7)*1+lsi)*1];
LIBINT2_REALTYPE fp67;
fp67 = fp72 + fp68;
LIBINT2_REALTYPE fp64;
fp64 = fp67 + fp65;
target[((hsi*21+12)*1+lsi)*1] = fp64;
LIBINT2_REALTYPE fp57;
fp57 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp56;
fp56 = fp57 * src3[((hsi*15+8)*1+lsi)*1];
LIBINT2_REALTYPE fp61;
fp61 = inteval->roe[vi] * src2[((hsi*10+4)*1+lsi)*1];
LIBINT2_REALTYPE fp60;
fp60 = src1[((hsi*10+4)*1+lsi)*1] - fp61;
LIBINT2_REALTYPE fp62;
fp62 = 2.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp59;
fp59 = fp62 * fp60;
LIBINT2_REALTYPE fp63;
fp63 = inteval->WQ_z[vi] * src0[((hsi*15+8)*1+lsi)*1];
LIBINT2_REALTYPE fp58;
fp58 = fp63 + fp59;
LIBINT2_REALTYPE fp55;
fp55 = fp58 + fp56;
target[((hsi*21+13)*1+lsi)*1] = fp55;
LIBINT2_REALTYPE fp48;
fp48 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp47;
fp47 = fp48 * src3[((hsi*15+9)*1+lsi)*1];
LIBINT2_REALTYPE fp52;
fp52 = inteval->roe[vi] * src2[((hsi*10+5)*1+lsi)*1];
LIBINT2_REALTYPE fp51;
fp51 = src1[((hsi*10+5)*1+lsi)*1] - fp52;
LIBINT2_REALTYPE fp53;
fp53 = 3.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp50;
fp50 = fp53 * fp51;
LIBINT2_REALTYPE fp54;
fp54 = inteval->WQ_z[vi] * src0[((hsi*15+9)*1+lsi)*1];
LIBINT2_REALTYPE fp49;
fp49 = fp54 + fp50;
LIBINT2_REALTYPE fp46;
fp46 = fp49 + fp47;
target[((hsi*21+14)*1+lsi)*1] = fp46;
LIBINT2_REALTYPE fp43;
fp43 = inteval->roe[vi] * src2[((hsi*10+6)*1+lsi)*1];
LIBINT2_REALTYPE fp42;
fp42 = src1[((hsi*10+6)*1+lsi)*1] - fp43;
LIBINT2_REALTYPE fp44;
fp44 = 4.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp41;
fp41 = fp44 * fp42;
LIBINT2_REALTYPE fp45;
fp45 = inteval->WQ_y[vi] * src0[((hsi*15+10)*1+lsi)*1];
LIBINT2_REALTYPE fp40;
fp40 = fp45 + fp41;
target[((hsi*21+15)*1+lsi)*1] = fp40;
LIBINT2_REALTYPE fp38;
fp38 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp37;
fp37 = fp38 * src3[((hsi*15+10)*1+lsi)*1];
LIBINT2_REALTYPE fp39;
fp39 = inteval->WQ_z[vi] * src0[((hsi*15+10)*1+lsi)*1];
LIBINT2_REALTYPE fp36;
fp36 = fp39 + fp37;
target[((hsi*21+16)*1+lsi)*1] = fp36;
LIBINT2_REALTYPE fp29;
fp29 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp28;
fp28 = fp29 * src3[((hsi*15+11)*1+lsi)*1];
LIBINT2_REALTYPE fp33;
fp33 = inteval->roe[vi] * src2[((hsi*10+6)*1+lsi)*1];
LIBINT2_REALTYPE fp32;
fp32 = src1[((hsi*10+6)*1+lsi)*1] - fp33;
LIBINT2_REALTYPE fp34;
fp34 = 1.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp31;
fp31 = fp34 * fp32;
LIBINT2_REALTYPE fp35;
fp35 = inteval->WQ_z[vi] * src0[((hsi*15+11)*1+lsi)*1];
LIBINT2_REALTYPE fp30;
fp30 = fp35 + fp31;
LIBINT2_REALTYPE fp27;
fp27 = fp30 + fp28;
target[((hsi*21+17)*1+lsi)*1] = fp27;
LIBINT2_REALTYPE fp20;
fp20 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp19;
fp19 = fp20 * src3[((hsi*15+12)*1+lsi)*1];
LIBINT2_REALTYPE fp24;
fp24 = inteval->roe[vi] * src2[((hsi*10+7)*1+lsi)*1];
LIBINT2_REALTYPE fp23;
fp23 = src1[((hsi*10+7)*1+lsi)*1] - fp24;
LIBINT2_REALTYPE fp25;
fp25 = 2.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp22;
fp22 = fp25 * fp23;
LIBINT2_REALTYPE fp26;
fp26 = inteval->WQ_z[vi] * src0[((hsi*15+12)*1+lsi)*1];
LIBINT2_REALTYPE fp21;
fp21 = fp26 + fp22;
LIBINT2_REALTYPE fp18;
fp18 = fp21 + fp19;
target[((hsi*21+18)*1+lsi)*1] = fp18;
LIBINT2_REALTYPE fp11;
fp11 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp10;
fp10 = fp11 * src3[((hsi*15+13)*1+lsi)*1];
LIBINT2_REALTYPE fp15;
fp15 = inteval->roe[vi] * src2[((hsi*10+8)*1+lsi)*1];
LIBINT2_REALTYPE fp14;
fp14 = src1[((hsi*10+8)*1+lsi)*1] - fp15;
LIBINT2_REALTYPE fp16;
fp16 = 3.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp13;
fp13 = fp16 * fp14;
LIBINT2_REALTYPE fp17;
fp17 = inteval->WQ_z[vi] * src0[((hsi*15+13)*1+lsi)*1];
LIBINT2_REALTYPE fp12;
fp12 = fp17 + fp13;
LIBINT2_REALTYPE fp9;
fp9 = fp12 + fp10;
target[((hsi*21+19)*1+lsi)*1] = fp9;
LIBINT2_REALTYPE fp2;
fp2 = 1.0000000000000000e+00 * inteval->alpha1_over_zetapluseta[vi];
LIBINT2_REALTYPE fp1;
fp1 = fp2 * src3[((hsi*15+14)*1+lsi)*1];
LIBINT2_REALTYPE fp6;
fp6 = inteval->roe[vi] * src2[((hsi*10+9)*1+lsi)*1];
LIBINT2_REALTYPE fp5;
fp5 = src1[((hsi*10+9)*1+lsi)*1] - fp6;
LIBINT2_REALTYPE fp7;
fp7 = 4.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp4;
fp4 = fp7 * fp5;
LIBINT2_REALTYPE fp8;
fp8 = inteval->WQ_z[vi] * src0[((hsi*15+14)*1+lsi)*1];
LIBINT2_REALTYPE fp3;
fp3 = fp8 + fp4;
LIBINT2_REALTYPE fp0;
fp0 = fp3 + fp1;
target[((hsi*21+20)*1+lsi)*1] = fp0;
LIBINT2_REALTYPE fp80;
fp80 = inteval->roe[vi] * src2[((hsi*10+3)*1+lsi)*1];
LIBINT2_REALTYPE fp79;
fp79 = src1[((hsi*10+3)*1+lsi)*1] - fp80;
LIBINT2_REALTYPE fp81;
fp81 = 3.0000000000000000e+00 * inteval->oo2e[vi];
LIBINT2_REALTYPE fp78;
fp78 = fp81 * fp79;
LIBINT2_REALTYPE fp82;
fp82 = inteval->WQ_y[vi] * src0[((hsi*15+6)*1+lsi)*1];
LIBINT2_REALTYPE fp77;
fp77 = fp82 + fp78;
target[((hsi*21+10)*1+lsi)*1] = fp77;
}
}
}
const int hsi = 0;
const int lsi = 0;
const int vi = 0;
/** Number of flops = 141 */
}

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