lib/dihedral.c


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40


#include "dihedral.h"

dihedral_t *dihedral_compose(q_t q, q_t g1i, const dihedral_t *g2) {
  // we only need to consider the action of reflections
  dihedral_t *g3 = (dihedral_t *)malloc(1 * sizeof(dihedral_t));

  g3->r = !g2->r;
  g3->i = (g1i + q - g2->i) % q;

  return g3;
}

q_t dihedral_act(q_t q, q_t gi, q_t s) {
  // we only need to consider the action of reflections

  return (gi + q - s) % q;
}

q_t dihedral_inverse_act(q_t q, const dihedral_t *g, q_t s) {
  if (g->r) {
    return (q - ((q + s - g->i) % q)) % q;
  } else {
    return (q + s - g->i) % q;
  }
}

q_t *dihedral_gen_transformations(q_t q) {
  q_t *transformations = (q_t *)malloc(q * q * sizeof(q_t));

  for (q_t i = 0; i < q; i++) {
    for (q_t j = 0; j < q; j++) {
      transformations[q * i + j] = dihedral_act(q, i, j);
    }
  }

  return transformations;
}