/* Model step function */
void model2_simulink_step(void)
{
  __m128d tmp_1;
  real_T rtb_Add_b[50];
  real_T rtb_Max[50];
  real_T tmp[5];
  real_T tmp_0;
  int32_T i;
  int32_T i_0;

  /* FC Layer 1: input(5) -> hidden(50) */
  for (i_0 = 0; i_0 < 50; i_0++) {
    tmp_0 = 0.0;
    for (i = 0; i < 5; i++) {
      tmp_0 += model2_simulink_ConstP.Weights_Value[50 * i + i_0] *
        model2_simulink_U.input[i];
    }
    rtb_Add_b[i_0] = tmp_0 + model2_simulink_ConstP.Bias_Value[i_0];
  }

  /* ReLU 1 */
  model2_simulink_relu1(rtb_Add_b, rtb_Max);

  /* FC Layers 2-6: hidden(50) -> hidden(50) with ReLU */
  /* ... (repeated pattern for each layer) ... */

  /* FC Layer 7: hidden(50) -> output(5) */
  /* Uses SSE2 intrinsics for vectorized MAC operations */
}

/* Shared ReLU implementation for all 6 hidden layers */
void model2_simulink_relu1(const real_T rtu_In1[50], real_T rty_Out1[50])
{
  int32_T i;
  for (i = 0; i < 50; i++) {
    rty_Out1[i] = fmax(rtu_In1[i], 0.0);
  }
}

/* Model step function (INT8) */
void model2_simulink_int8_step(void)
{
  int32_T i;
  int32_T i_0;
  int32_T tmp_0;
  int8_T rtb_DataTypeConversion_m[50];
  int8_T rtb_Max_p[50];
  int8_T tmp[50];

  /* FC Layer 1 with INT8 weights and activations */
  for (i = 0; i < 50; i++) {
    tmp_0 = 0;
    for (i_0 = 0; i_0 < 5; i_0++) {
      tmp_0 += model2_simulink_int8_ConstP.Weights_Value_f[50 * i_0 + i] *
        model2_simulink_int8_U.input[i_0];
    }
    /* Saturate to INT8 range [-128, 127] */
    /* ... */
  }

  /* INT8 ReLU: simple comparison with 0 */
  model2_simulink_int8_relu1(rtb_DataTypeConversion_m, rtb_Max_p);
}

/* INT8 ReLU - branch instead of fmax */
void model2_simulink_int8_relu1(const int8_T rtu_In1[50], int8_T rty_Out1[50])
{
  int32_T i;
  for (i = 0; i < 50; i++) {
    int8_T u0;
    u0 = rtu_In1[i];
    if (u0 >= 0) {
      rty_Out1[i] = u0;
    } else {
      rty_Out1[i] = 0;
    }
  }
}