#include <libint2.h>

#ifdef __cplusplus
extern "C" {
#endif
void CR_aB_s__0___CartesianMultipole_sB_1_c_0_c_0_Sb__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 fp0;
fp0 = 0 + src2[((hsi*7+5)*1+lsi)*1];
LIBINT2_REALTYPE fp11;
fp11 = 0 + src1[((hsi*7+0)*1+lsi)*1];
LIBINT2_REALTYPE fp15;
fp15 = 0 + src0[((hsi*7+1)*1+lsi)*1];
LIBINT2_REALTYPE fp12;
fp12 = 0 + src0[((hsi*7+0)*1+lsi)*1];
LIBINT2_REALTYPE fp14;
fp14 = inteval->BO_x[vi] * fp12;
LIBINT2_REALTYPE fp13;
fp13 = fp15 + fp14;
LIBINT2_REALTYPE fp32;
fp32 = fp13 * fp11;
LIBINT2_REALTYPE fp31;
fp31 = fp32 * fp0;
target[((hsi*21+20)*1+lsi)*1] = fp31;
LIBINT2_REALTYPE fp1;
fp1 = 0 + src2[((hsi*7+4)*1+lsi)*1];
LIBINT2_REALTYPE fp10;
fp10 = 0 + src1[((hsi*7+1)*1+lsi)*1];
LIBINT2_REALTYPE fp34;
fp34 = fp13 * fp10;
LIBINT2_REALTYPE fp33;
fp33 = fp34 * fp1;
target[((hsi*21+19)*1+lsi)*1] = fp33;
LIBINT2_REALTYPE fp2;
fp2 = 0 + src2[((hsi*7+3)*1+lsi)*1];
LIBINT2_REALTYPE fp9;
fp9 = 0 + src1[((hsi*7+2)*1+lsi)*1];
LIBINT2_REALTYPE fp36;
fp36 = fp13 * fp9;
LIBINT2_REALTYPE fp35;
fp35 = fp36 * fp2;
target[((hsi*21+18)*1+lsi)*1] = fp35;
LIBINT2_REALTYPE fp3;
fp3 = 0 + src2[((hsi*7+2)*1+lsi)*1];
LIBINT2_REALTYPE fp8;
fp8 = 0 + src1[((hsi*7+3)*1+lsi)*1];
LIBINT2_REALTYPE fp38;
fp38 = fp13 * fp8;
LIBINT2_REALTYPE fp37;
fp37 = fp38 * fp3;
target[((hsi*21+17)*1+lsi)*1] = fp37;
LIBINT2_REALTYPE fp4;
fp4 = 0 + src2[((hsi*7+1)*1+lsi)*1];
LIBINT2_REALTYPE fp7;
fp7 = 0 + src1[((hsi*7+4)*1+lsi)*1];
LIBINT2_REALTYPE fp40;
fp40 = fp13 * fp7;
LIBINT2_REALTYPE fp39;
fp39 = fp40 * fp4;
target[((hsi*21+16)*1+lsi)*1] = fp39;
LIBINT2_REALTYPE fp5;
fp5 = 0 + src2[((hsi*7+0)*1+lsi)*1];
LIBINT2_REALTYPE fp6;
fp6 = 0 + src1[((hsi*7+5)*1+lsi)*1];
LIBINT2_REALTYPE fp42;
fp42 = fp13 * fp6;
LIBINT2_REALTYPE fp41;
fp41 = fp42 * fp5;
target[((hsi*21+15)*1+lsi)*1] = fp41;
LIBINT2_REALTYPE fp18;
fp18 = 0 + src0[((hsi*7+2)*1+lsi)*1];
LIBINT2_REALTYPE fp17;
fp17 = inteval->BO_x[vi] * fp15;
LIBINT2_REALTYPE fp16;
fp16 = fp18 + fp17;
LIBINT2_REALTYPE fp44;
fp44 = fp16 * fp11;
LIBINT2_REALTYPE fp43;
fp43 = fp44 * fp1;
target[((hsi*21+14)*1+lsi)*1] = fp43;
LIBINT2_REALTYPE fp46;
fp46 = fp16 * fp10;
LIBINT2_REALTYPE fp45;
fp45 = fp46 * fp2;
target[((hsi*21+13)*1+lsi)*1] = fp45;
LIBINT2_REALTYPE fp48;
fp48 = fp16 * fp9;
LIBINT2_REALTYPE fp47;
fp47 = fp48 * fp3;
target[((hsi*21+12)*1+lsi)*1] = fp47;
LIBINT2_REALTYPE fp50;
fp50 = fp16 * fp8;
LIBINT2_REALTYPE fp49;
fp49 = fp50 * fp4;
target[((hsi*21+11)*1+lsi)*1] = fp49;
LIBINT2_REALTYPE fp52;
fp52 = fp16 * fp7;
LIBINT2_REALTYPE fp51;
fp51 = fp52 * fp5;
target[((hsi*21+10)*1+lsi)*1] = fp51;
LIBINT2_REALTYPE fp21;
fp21 = 0 + src0[((hsi*7+3)*1+lsi)*1];
LIBINT2_REALTYPE fp20;
fp20 = inteval->BO_x[vi] * fp18;
LIBINT2_REALTYPE fp19;
fp19 = fp21 + fp20;
LIBINT2_REALTYPE fp54;
fp54 = fp19 * fp11;
LIBINT2_REALTYPE fp53;
fp53 = fp54 * fp2;
target[((hsi*21+9)*1+lsi)*1] = fp53;
LIBINT2_REALTYPE fp56;
fp56 = fp19 * fp10;
LIBINT2_REALTYPE fp55;
fp55 = fp56 * fp3;
target[((hsi*21+8)*1+lsi)*1] = fp55;
LIBINT2_REALTYPE fp58;
fp58 = fp19 * fp9;
LIBINT2_REALTYPE fp57;
fp57 = fp58 * fp4;
target[((hsi*21+7)*1+lsi)*1] = fp57;
LIBINT2_REALTYPE fp60;
fp60 = fp19 * fp8;
LIBINT2_REALTYPE fp59;
fp59 = fp60 * fp5;
target[((hsi*21+6)*1+lsi)*1] = fp59;
LIBINT2_REALTYPE fp24;
fp24 = 0 + src0[((hsi*7+4)*1+lsi)*1];
LIBINT2_REALTYPE fp23;
fp23 = inteval->BO_x[vi] * fp21;
LIBINT2_REALTYPE fp22;
fp22 = fp24 + fp23;
LIBINT2_REALTYPE fp62;
fp62 = fp22 * fp11;
LIBINT2_REALTYPE fp61;
fp61 = fp62 * fp3;
target[((hsi*21+5)*1+lsi)*1] = fp61;
LIBINT2_REALTYPE fp64;
fp64 = fp22 * fp10;
LIBINT2_REALTYPE fp63;
fp63 = fp64 * fp4;
target[((hsi*21+4)*1+lsi)*1] = fp63;
LIBINT2_REALTYPE fp66;
fp66 = fp22 * fp9;
LIBINT2_REALTYPE fp65;
fp65 = fp66 * fp5;
target[((hsi*21+3)*1+lsi)*1] = fp65;
LIBINT2_REALTYPE fp27;
fp27 = 0 + src0[((hsi*7+5)*1+lsi)*1];
LIBINT2_REALTYPE fp26;
fp26 = inteval->BO_x[vi] * fp24;
LIBINT2_REALTYPE fp25;
fp25 = fp27 + fp26;
LIBINT2_REALTYPE fp68;
fp68 = fp25 * fp11;
LIBINT2_REALTYPE fp67;
fp67 = fp68 * fp4;
target[((hsi*21+2)*1+lsi)*1] = fp67;
LIBINT2_REALTYPE fp70;
fp70 = fp25 * fp10;
LIBINT2_REALTYPE fp69;
fp69 = fp70 * fp5;
target[((hsi*21+1)*1+lsi)*1] = fp69;
LIBINT2_REALTYPE fp30;
fp30 = inteval->BO_x[vi] * fp27;
LIBINT2_REALTYPE fp28;
fp28 = 0 + src0[((hsi*7+6)*1+lsi)*1];
LIBINT2_REALTYPE fp29;
fp29 = fp28 + fp30;
LIBINT2_REALTYPE fp72;
fp72 = fp29 * fp11;
LIBINT2_REALTYPE fp71;
fp71 = fp72 * fp5;
target[((hsi*21+0)*1+lsi)*1] = fp71;
}
}
}
const int hsi = 0;
const int lsi = 0;
const int vi = 0;
/** Number of flops = 73 */
}

#ifdef __cplusplus
};
#endif
