// local
static const int N = 6;
int32_t M[N][N];
int32_t tmp[N][N];
int32_t R[N][N];

int32_t sin2t_q27(uint32_t p) {
  int32_t a = (p & ((1 << 20) - 1)) << 8; // q28
  uint32_t i = p >> 20;
  return ___SMMLA(a, sine2t[i + 1] - sine2t[i], sine2t[i] >> 4);
}
int32_t cos2t_q27(uint32_t p) { return sin2t_q27(p + (1UL << 30)); }
// M q27;
void setR(int axis1, int axis2, uint32_t ia) {
  for (int j = 0; j < N; j++) {
    for (int i = 0; i < N; i++) {
      R[j][i] = i == j ? 1 << 27 : 0;
    }
  }
  int32_t c = cos2t_q27(ia);
  int32_t s = sin2t_q27(ia);
  R[axis1][axis1] = c;
  R[axis1][axis2] = -s;
  R[axis2][axis1] = s;
  R[axis2][axis2] = c;
}
void setR(int axis1, int axis2, float angle) {
  setR(axis1, axis2, (uint32_t)(angle * 11930464.71f)); // (2^32)/360
}

void MxR(void) {
  for (int j = 0; j < N; j++) {
    for (int i = 0; i < N; i++) {
      tmp[j][i] = M[j][i];
    }
  }
  for (int j = 0; j < N; j++) {
    for (int i = 0; i < N; i++) {
      M[j][i] = 0;
      for (int k = 0; k < N; k++) {
        M[j][i] += ___SMMUL(tmp[j][k], R[k][i]) << 5;
      }
    }
  }
}
void rot(int axis1, int axis2, float angle) {
  setR(axis1 % N, axis2 % N, angle);
  MxR();
}
void negLine(int j) {
  for (int i = 0; i < N; i++) {
    M[j][i] = -M[j][i];
  }
}
void negCol(int i) {
  for (int j = 0; j < N; j++) {
    M[j][i] = -M[j][i];
  }
}
void randR(uint32_t seed) {
  uint32_t r = seed * 69069 + 1;
  for (int k = 0; k < 4; k++) {
    for (int j = 0; j < N; j++) {
      for (int i = j + 1; i < N; i++) {
        setR(i, j, r = r * 69069 + 1);
        MxR();
      }
    }
  }
}
void randO(uint32_t seed) {
  uint32_t r = seed * 69069 + 1;
  for (int k = 0; k < 4; k++) {
    for (int j = 0; j < N; j++) {
      for (int i = j + 1; i < N; i++) {
        setR(i, j, r = r * 69069 + 1);
        MxR();
      }
      r = r * 69069 + 1;
      negLine((r >> 24) % N);
    }
  }
}
void transpose(void) {
  for (int j = 0; j < N; j++) {
    for (int i = 0; i < N; i++) {
      tmp[j][i] = M[j][i];
    }
  }
  for (int j = 0; j < N; j++) {
    for (int i = 0; i < N; i++) {
      M[j][i] = tmp[i][j];
    }
  }
}

for (int j = 0; j < N; j++) {
  for (int i = 0; i < N; i++) {
    M[j][i] = (i == j) ? 1 << 27 : 0;
  }
}
{ attr_rotations }

int32_t *out[N];
out[0] = outlet_out0;
out[1] = outlet_out1;
out[2] = outlet_out2;
out[3] = outlet_out3;
out[4] = outlet_out4;
out[5] = outlet_out5;
for (int ni = 0; ni < N; ni++) {
  int32_t *o = out[ni];
  int32_t *m = M[ni];
  {
    int32_t m0 = m[0];
    int32_t m1 = m[1];
    int32_t m2 = m[2];
    for (int i = 0; i < BUFSIZE; i++) {
      o[i] = ___SMMLA(inlet_in2[i], m2,
                      ___SMMLA(inlet_in1[i], m1, ___SMMUL(inlet_in0[i], m0)));
    }
  }
  {
    int32_t m3 = m[3];
    int32_t m4 = m[4];
    int32_t m5 = m[5];
    for (int i = 0; i < BUFSIZE; i++) {
      o[i] =
          ___SMMLA(inlet_in5[i], m5,
                   ___SMMLA(inlet_in4[i], m4, ___SMMLA(inlet_in3[i], m3, o[i])))
          << 5;
    }
  }
}