oa 2010
[aversive.git] / projects / microb2010 / tests / oa / main.c
diff --git a/projects/microb2010/tests/oa/main.c b/projects/microb2010/tests/oa/main.c
new file mode 100644 (file)
index 0000000..5d05ae9
--- /dev/null
@@ -0,0 +1,467 @@
+#include <stdio.h>
+#include <stdlib.h>
+#include <stdint.h>
+#include <math.h>
+
+/* XXX TODO
+static
+const
+change x,y -> i,j to avoid confusion with coords
+could be optimized in mem space: it is not needed to store the x,y coord,
+   we can process it from idx. however it will be less optimized for speed
+
+*/
+
+#define OFFSET_CORN_X 150
+#define OFFSET_CORN_Y 222
+#define STEP_CORN_X 225
+#define STEP_CORN_Y 250
+
+#define CORN_NB 18
+
+#define WAYPOINTS_NBX 13
+#define WAYPOINTS_NBY 8
+
+/* enum is better */
+#define TYPE_WAYPOINT 0
+#define TYPE_DANGEROUS 1
+#define TYPE_WHITE_CORN 2
+#define TYPE_BLACK_CORN 3
+#define TYPE_OBSTACLE 4
+
+/* XXX enum possible ? else just rename */
+#define START      0
+#define UP         1
+#define UP_RIGHT   2
+#define DOWN_RIGHT 3
+#define DOWN       4
+#define DOWN_LEFT  5
+#define UP_LEFT    6
+#define END        7
+
+struct point {
+       int32_t x;
+       int32_t y;
+};
+
+struct djpoint {
+       struct point pos;
+       uint16_t weight;
+       struct djpoint *parent;
+       
+       uint8_t type:3;
+       uint8_t parent_pos:3;
+       uint8_t updated:1;
+       uint8_t todo:1;
+};
+
+uint8_t corn_table[CORN_NB];
+static struct djpoint djpoints[WAYPOINTS_NBX][WAYPOINTS_NBY];
+
+/* table to find the symetric idx */
+uint8_t corn_sym[] = {
+       15, 16, 17, 12, 13, 14, 10, 11, 8, 9, 6, 7, 3, 4, 5, 0, 1, 2 
+};
+
+uint8_t corn_side_confs[9][2] = {
+       { 1, 4 },
+       { 0, 4 },
+       { 2, 4 },
+       { 2, 3 },
+       { 0, 3 },
+       { 1, 3 },
+       { 1, 6 },
+       { 0, 6 },
+       { 2, 6 },
+};
+uint8_t corn_center_confs[4][2] = {
+       { 5, 8 },
+       { 7, 8 },
+       { 5, 9 },
+       { 7, 8 },
+};
+
+
+/* return index from neigh pointer */
+#define PT2IDX(neigh) ( ((void *)(neigh)-(void *)(&djpoints)) / sizeof(*neigh) )
+
+void dump(void)
+{
+       int8_t i, j;
+       struct djpoint *pt;
+
+       printf("         ");
+       for (i=0; i<WAYPOINTS_NBX; i++) {
+               printf(" %2d ", i);
+       }
+       printf("\n");
+                       
+       for (j=WAYPOINTS_NBY*2-1; j>=0; j--) {
+               printf("%3d   ", j/2);
+
+               if ((j&1) == 0)
+                       printf("    ");
+
+               for (i=0; i<WAYPOINTS_NBX; i++) {
+                       pt = &djpoints[i][j/2];
+                       
+                       if (((i+j) & 1) == 0)
+                               continue;
+
+                       if (pt->type == TYPE_OBSTACLE)
+                               printf("     X  ");
+                       else if (pt->type == TYPE_DANGEROUS)
+                               printf("     D  ");
+                       else if (pt->type == TYPE_WHITE_CORN)
+                               printf("     W  ");
+                       else if (pt->type == TYPE_BLACK_CORN)
+                               printf("     B  ");
+                       else if (pt->type == TYPE_WAYPOINT)
+                               printf(" %5d  ", pt->weight);
+                       else
+                               printf("     ?  "); 
+               }
+               printf("\n");
+       }
+}
+
+static inline uint8_t opposite_position(uint8_t pos)
+{
+       pos += 3;
+       if (pos > UP_LEFT)
+               pos -= 6;
+       return pos;
+}
+
+/* return coord of the entry in the table from the pointer */
+static void djpoint2ij(struct djpoint *pt, int8_t *x, int8_t *y)
+{
+       int8_t idx = PT2IDX(pt);
+       *x = idx / WAYPOINTS_NBY;
+       *y = idx % WAYPOINTS_NBY;
+}
+
+/* get the neighbour of the point at specified position */
+static struct djpoint *get_neigh(struct djpoint *pt,
+                                uint8_t position)
+{
+       int8_t i,j;
+       struct djpoint *neigh;
+
+       djpoint2ij(pt, &i, &j);
+
+       switch (position) {
+       case UP:
+               j++;
+               break;
+       case UP_RIGHT:
+               if (!(i & 1)) j++;
+               i++;
+               break;
+       case DOWN_RIGHT:
+               if (i & 1) j--;
+               i++;
+               break;
+       case DOWN:
+               j--;
+               break;
+       case DOWN_LEFT:
+               if (i & 1) j--;
+               i--;
+               break;
+       case UP_LEFT:
+               if (!(i & 1)) j++;
+               i--;
+               break;
+       default:
+               return NULL;
+       }
+       if (i < 0 || j < 0 || i >= WAYPOINTS_NBX || j >= WAYPOINTS_NBY)
+               return NULL;
+
+       neigh = &djpoints[i][j];
+
+       if (neigh->type != TYPE_WAYPOINT)
+               return NULL;
+
+       return neigh;
+}
+
+static struct djpoint *get_next_neigh(struct djpoint *pt, uint8_t *position)
+{
+       struct djpoint *neigh = NULL;
+       uint8_t pos = *position + 1;
+
+       for (pos = *position + 1; pos < END; pos++) {
+               neigh = get_neigh(pt, pos);
+               if (neigh != NULL)
+                       break;
+       }
+
+       *position = pos;
+       return neigh;
+}
+
+/* browse all points */
+#define POINT_FOREACH(cur)                                             \
+       for (cur = &djpoints[0][0];                                     \
+            cur < &djpoints[WAYPOINTS_NBX][WAYPOINTS_NBY];             \
+            cur ++)
+
+/* XXX comment */
+#define NEIGH_FOREACH(neigh, pos, point)                       \
+       for (pos = START, neigh = get_next_neigh(point, &pos);  \
+            neigh;                                             \
+            neigh = get_next_neigh(point, &pos))
+
+int dijkstra_init(void)
+{
+       return 0;
+}
+
+static uint16_t dist(struct djpoint *p1, struct djpoint *p2)
+{
+       double vx, vy;
+       vx = p2->pos.x - p1->pos.x;
+       vy = p2->pos.y - p1->pos.y;
+       return sqrt(vx * vx + vy * vy);
+}
+
+void dijkstra_process_neighs(struct djpoint *pt)
+{
+       struct djpoint *neigh;
+       uint8_t pos, parent_pos;
+       uint16_t weight;
+       int8_t i,j;
+
+       djpoint2ij(pt, &i, &j);
+       printf("at %d %d:\n", i, j);
+
+       NEIGH_FOREACH(neigh, pos, pt) {
+               weight = pt->weight + dist(pt, neigh);
+               parent_pos = opposite_position(pos);
+
+               /* bonus if we keep the same direction */
+               if (parent_pos == pt->parent_pos ||
+                   pt->parent_pos == END)
+                       weight = weight - 1;
+
+               printf("  pos=%d: ppos=%d opos=%d nw=%d w=%d\n", pos,
+                      pt->parent_pos, parent_pos,
+                      neigh->weight, weight);
+               if (neigh->weight == 0 || weight < neigh->weight) {
+                       djpoint2ij(neigh, &i, &j);
+                       //printf("    %d,%d updated\n", i, j);
+                       neigh->weight = weight;
+                       neigh->parent_pos = parent_pos;
+                       neigh->updated = 1;
+               }
+       }
+}
+
+int dijkstra(struct djpoint *start)
+{
+       struct djpoint *cur;
+       uint8_t todolist = 1;
+
+       start->todo = 1;
+
+       while (todolist) {
+               printf("\n");
+               dump();
+               /* process all neighbours of todo list */
+               POINT_FOREACH(cur) {
+                       if (!cur->todo)
+                               continue;
+                       dijkstra_process_neighs(cur);
+                       cur->todo = 0;
+               }
+
+               /* convert updated list in todo list */
+               todolist = 0;
+               POINT_FOREACH(cur) {
+                       if (!cur->updated)
+                               continue;
+                       todolist = 1;
+                       cur->todo = 1;
+                       cur->updated = 0;
+               }               
+       }
+       return 0; /* XXX */
+}
+
+int8_t coord_to_corn_idx(int8_t i, int8_t j)
+{
+       if (i == 0 && j == 2) return 0;
+       if (i == 0 && j == 4) return 1;
+       if (i == 0 && j == 6) return 2;
+       if (i == 2 && j == 3) return 3;
+       if (i == 2 && j == 5) return 4;
+       if (i == 2 && j == 7) return 5;
+       if (i == 4 && j == 4) return 6;
+       if (i == 4 && j == 6) return 7;
+       if (i == 6 && j == 5) return 8;
+       if (i == 6 && j == 7) return 9;
+       if (i == 8 && j == 4) return 10;
+       if (i == 8 && j == 6) return 11;
+       if (i == 10 && j == 3) return 12;
+       if (i == 10 && j == 5) return 13;
+       if (i == 10 && j == 7) return 14;
+       if (i == 12 && j == 2) return 15;
+       if (i == 12 && j == 4) return 16;
+       if (i == 12 && j == 6) return 17;
+       return -1;
+}
+
+int8_t corn_get_sym(int8_t i)
+{
+       if (i >= CORN_NB)
+               return -1;
+       return corn_sym[i];
+}
+
+void init_corn_table(uint8_t conf_side, uint8_t conf_center)
+{
+       int8_t sym, i;
+
+       printf("confs = %d, %d\n", conf_side, conf_center);
+       for (i=0; i<CORN_NB; i++) {
+               if (i == corn_side_confs[conf_side][0] ||
+                   i == corn_side_confs[conf_side][1]) {
+                       corn_table[i] = TYPE_BLACK_CORN;
+                       continue;
+               }
+               sym = corn_get_sym(i);
+               if (sym == corn_side_confs[conf_side][0] ||
+                   sym == corn_side_confs[conf_side][1]) {
+                       corn_table[i] = TYPE_BLACK_CORN;
+                       continue;
+               }
+               if (i == corn_center_confs[conf_center][0] ||
+                   i == corn_center_confs[conf_center][1]) {
+                       corn_table[i] = TYPE_BLACK_CORN;
+                       continue;
+               }
+               sym = corn_get_sym(i);
+               if (sym == corn_center_confs[conf_center][0] ||
+                   sym == corn_center_confs[conf_center][1]) {
+                       corn_table[i] = TYPE_BLACK_CORN;
+                       continue;
+               }
+               corn_table[i] = TYPE_WHITE_CORN;
+       }
+}
+
+/* init waypoints position */
+void init_waypoints(void)
+{
+       int8_t i, j;
+       int32_t x, y;
+       struct djpoint *pt;
+
+       x = OFFSET_CORN_X;
+       for (i=0; i<WAYPOINTS_NBX; i++) {
+
+               if (i & 1)
+                       y = OFFSET_CORN_Y - STEP_CORN_Y/2;
+               else
+                       y = OFFSET_CORN_Y;                      
+
+               for (j=0; j<WAYPOINTS_NBY; j++) {
+                       pt = &djpoints[i][j];
+
+                       pt->pos.x = x;
+                       pt->pos.y = y;
+
+                       pt->type = TYPE_WAYPOINT;
+                       pt->parent_pos = END;
+                       pt->updated = 0;
+                       pt->todo = 0;
+
+                       y += STEP_CORN_Y;
+               }
+
+               x += STEP_CORN_X;
+       }
+}
+
+/* update the type and weight of waypoints, before starting
+ * dijkstra */
+void update_waypoints(void)
+{
+       int8_t i, j, c;
+       struct djpoint *pt;
+
+       for (i=0; i<WAYPOINTS_NBX; i++) {
+
+               for (j=0; j<WAYPOINTS_NBY; j++) {
+                       pt = &djpoints[i][j];
+
+                       pt->weight = 0;
+
+                       /* corn */
+                       c = coord_to_corn_idx(i, j);
+                       if (c >= 0) {
+                               pt->type = corn_table[c];
+                               continue;
+                       }
+                       /* too close of border */
+                       if ((i & 1) == 1 && j == (WAYPOINTS_NBY-1)) {
+                               pt->type = TYPE_OBSTACLE;
+                               continue;
+                       }
+                       /* hill */
+                       if (i >= 2 && i < (WAYPOINTS_NBX-2) && j < 2) {
+                               pt->type = TYPE_OBSTACLE;
+                               continue;
+                       }
+                       /* dangerous points */
+                       if (i == 0 || i == (WAYPOINTS_NBX-1)) {
+                               pt->type = TYPE_DANGEROUS;
+                               continue;
+                       }
+                       if ( (i&1) == 0 && j == (WAYPOINTS_NBY-1)) {
+                               pt->type = TYPE_DANGEROUS;
+                               continue;
+                       }
+                       pt->type = TYPE_WAYPOINT;
+               }
+       }
+}
+
+int get_path(struct djpoint *start, struct djpoint *end)
+{
+       struct djpoint *cur;
+       uint8_t prev_direction = 0;
+
+       for (cur=start;
+            cur != NULL && cur->parent_pos != END && cur != end;
+            cur = get_neigh(cur, cur->parent_pos)) {
+               if (prev_direction != cur->parent_pos)
+                       printf("%d %d (%d)\n", cur->pos.x, cur->pos.y, cur->parent_pos);
+               prev_direction = cur->parent_pos;
+       }
+       printf("%d %d\n", end->pos.x, end->pos.y);
+
+       return 0; /* XXX */
+}
+
+int main(void)
+{
+       struct djpoint *start;
+       struct djpoint *end;
+       
+       start = &djpoints[1][1];
+       end = &djpoints[12][1];
+
+       init_corn_table(0, 0);
+       init_waypoints();
+       update_waypoints();
+
+       dijkstra(end);
+       dump();
+
+       get_path(start, end);
+
+       return 0;
+}