labeller.c

   1 #include <errno.h>
   2
   3 #include <ucw/lib.h>
   4 #include <ucw/gary.h>
   5 #include <ucw/mempool.h>
   6 #include <ucw/eltpool.h>
   7
   8 #include "leo.h"
   9 #include "sym.h"
  10 #include "map.h"
  11 #include "labeller.h"
  12
  13 #define HASH_NODE struct graph_node
  14 #define HASH_PREFIX(x) hash_##x
  15 #define HASH_KEY_ATOMIC id
  16 #define HASH_WANT_FIND
  17 #define HASH_WANT_NEW
  18 #define HASH_WANT_CLEANUP
  19 #include <ucw/hashtable.h>
  20
  21 #include <stdio.h>
  22 #include <stdlib.h>
  23 #include <math.h>
  24 #include <fcntl.h>
  25
  26 #define BLOCK_SIZE 4096
  27
  28 //struct mempool *mpool_requests;
  29
  30 static struct request_point *requests_point;
  31 static struct request_line *requests_line;
  32 static struct request_area *requests_area;
  33
  34 static struct graph_edge **bfs_queue;
  35 static struct longline *longlines; int num_longlines;
  36 static struct buffer_line *buffer_line;
  37 static struct buffer_linelabel *buffer_linelabel;
  38
  39 struct eltpool *ep_individuals;
  40
  41 struct individual **population1;
  42 struct individual **population2;
  43
  44 int dbg_segments = 0;
  45 int dbg_plan = 0;
  46
  47 int page_width_int;
  48 int page_height_int;
  49
  50 int num_edges_dbg;
  51 int num_nodes;
  52 int num_edges = 0;
  53 int dbg_num_hits = 0;
  54
  55 int conf_pop_size = 50;
  56
  57 int conf_penalty_bound = 0;
  58 int conf_stagnation_bound = 0;
  59 int conf_iteration_limit = 4;
  60
  61 int conf_term_cond = TERM_COND_ITERATIONS;
  62
  63 int conf_breed_rbest_perc = 80;
  64 int conf_breed_pop_size_perc = 20;
  65 int conf_breed_perc = 50;                       // Percentage of new pop created by breeding
  66
  67 bool conf_mutate_children = 1;
  68 int conf_mutate_children_prob = 0.3;
  69
  70 int conf_mutate_rbest_perc = 60;
  71 int conf_mutate_pop_size_perc = 20;
  72
  73 int conf_mutate_move_bound = 0.2;
  74 int conf_mutate_regen_bound = 0.1;
  75 int conf_mutate_chvar_bound = 0.1;
  76
  77 int conf_elite_perc = 5;
  78
  79 double conf_max_section_length = 100;
  80 double conf_max_section_overlay = 10;
  81
  82 int old_best = 0; // FIXME: Shall be int max
  83 int iteration = 0;
  84 int pop2_ind;
  85
  86 int conf_part_size = 50;
  87
  88 int move_min = 0;
  89 int move_max = 1;
  90
  91 int num_requests = 0;
  92
  93 void make_graph(void);
  94 void label_graph(void);
  95 void join_edge(struct graph_edge *e, int dir);
  96 void bfs(uns longline);
  97 void make_segments(void);
  98
  99 void make_population(void);
 100 bool shall_terminate(void);
 101 void breed(void);
 102 void mutate(void);
 103 void elite(void);
 104 void rank_population(void);
 105 void plan_individual(struct individual *individual);
 106
 107 void make_bitmap(struct point_variant *v, struct symbol *sym);
 108 void make_bitmap_icon(struct point_variant *v, struct sym_icon *si);
 109 void make_bitmap_point(struct point_variant *v, struct sym_point *sp);
 110 void make_bitmap_label(struct point_variant *v, struct sym_text *text);
 111
 112 void cut_edge(struct graph_edge *e, double dist);
 113 struct request_line *make_new_line(void);
 114 struct request_section *make_new_section(struct request_line *rl);
 115 struct request_segment *make_new_segment(struct request_section *rls, struct symbol *sym);
 116
 117 void dump_bitmaps(struct individual *individual);
 118 void dump_graph(void);
 119 void bfs2(void);
 120 void bfs_edge(struct graph_edge *e, struct graph_node *node, struct graph_node *anode, enum edge_dir dir);
 121 void bfs_wrapper(void);
 122 void oldbfs(void);
 123 void dump_longlines(void);
 124 void dump_linelabel_requests(void);
 125 void dump_individual(struct individual *individual);
 126 void print_label(struct symbol *sym);
 127
 128 double gen_movement(void);
 129 void gen_coords(struct placement *p);
 130 void gen_coords_point(struct placement *p);
 131 void gen_coords_segment(struct placement *p);
 132 void gen_coords_area(struct placement *p);
 133
 134 void make_segments_old(void);
 135
 136 void labeller_cleanup(void);
 137
 138 struct individual **perform_crossover(struct individual *parent1, struct individual *parent2);
 139 void perform_mutation(struct individual *individual);
 140
 141 void hide_segment_labels(struct individual *individual);
 142 void init_placement(struct placement *p, struct individual *individual, struct request *r);
 143 void init_individual(struct individual *i);
 144 struct map_part **get_parts(struct placement *symbol, struct individual *individual);
 145
 146 int randint(int min, int max);
 147
 148 struct placement **get_closure(struct placement *placement, struct individual *parent1, struct individual *parent2);
 149 void copy_symbols(struct placement **closure, struct individual *parent, struct individual *child);
 150 void move_symbol(struct placement *p);
 151 void move_symbol_point(struct placement *p);
 152
 153 struct placement **get_overlapping(struct placement *p);
 154 void filter(struct placement **list, bool *pred);
 155
 156 int flip(int a, int b);
 157 double randdouble(void);
 158
 159 void cleanup(void);
 160
 161 void copy_individual(struct individual *src, struct individual *dest);
 162
 163 int max2(int a, int b);
 164 int min2(int a, int b);
 165 int max4(int a, int b, int c, int d);
 166 int min4(int a, int b, int c, int d);
 167
 168 int max2(int a, int b)
 169 {
 170   return (a > b ? a : b);
 171 }
 172
 173 int min2(int a, int b)
 174 {
 175   return (a < b ? a : b);
 176 }
 177
 178 int max4(int a, int b, int c, int d)
 179 {
 180   return max2(max2(a, b), max2(c, d));
 181 }
 182
 183 int min4(int a, int b, int c, int d)
 184 {
 185   return min2(min2(a, b), min2(c, d));
 186 }
 187
 188 void print_label(struct symbol *sym)
 189 {
 190   switch (sym->type)
 191   {
 192     case SYMBOLIZER_TEXT: ;
 193       struct sym_text *st = (struct sym_text *) sym;
 194       printf("%s\n", osm_val_decode(st->text));
 195     default:
 196       // FIXME
 197       ;
 198   }
 199 }
 200
 201 void labeller_init(void)
 202 {
 203   GARY_INIT(requests_point, 0);
 204   GARY_INIT(requests_line, 0);
 205   GARY_INIT(requests_area, 0);
 206   GARY_INIT(buffer_line, 0);
 207   GARY_INIT(buffer_linelabel, 0);
 208   ep_individuals = ep_new(sizeof(struct individual), 1);
 209
 210   page_width_int = floor(page_width);
 211   page_height_int = floor(page_height);
 212 }
 213
 214 void make_bitmap(struct point_variant *v, struct symbol *sym)
 215 {
 216   switch (sym->type)
 217   {
 218     case SYMBOLIZER_POINT:
 219       make_bitmap_point(v, (struct sym_point *) sym);
 220       break;
 221     case SYMBOLIZER_ICON:
 222       make_bitmap_icon(v, (struct sym_icon *) sym);
 223       break;
 224     case SYMBOLIZER_TEXT:
 225       make_bitmap_label(v, (struct sym_text *) sym);
 226       break;
 227     default:
 228       ASSERT(sym->type != SYMBOLIZER_INVALID);
 229   }
 230 }
 231
 232 void make_bitmap_icon(struct point_variant *v, struct sym_icon *si)
 233 {
 234   v->width = si->sir.icon->width;
 235   v->height = si->sir.icon->height;
 236   v->bitmap = malloc((int) ceil(v->width * v->height * sizeof(bool)));
 237   for (int i=0; i<v->width*v->height; i++) v->bitmap[i] = 1;
 238 }
 239
 240 void make_bitmap_point(struct point_variant *v, struct sym_point *sp)
 241 {
 242   v->width = v->height = sp->size;
 243   v->bitmap = malloc(sp->size*sp->size * sizeof(bool));
 244   // FIXME: Okay, memset would be much nicer here
 245   for (int i=0; i<sp->size*sp->size; i++) v->bitmap[i] = 1;
 246 }
 247
 248 void make_bitmap_label(struct point_variant *v, struct sym_text *text)
 249 {
 250   int x_ld = 0;
 251   int y_ld = 0;
 252   int x_lu = 0;
 253   int y_lu = 0;
 254   int x_rd = 0;
 255   int y_rd = 0;
 256   int x_ru = 0;
 257   int y_ru = 0;
 258
 259   v->width = max4(x_ld, x_lu, x_rd, x_ru) - min4(x_ld, x_lu, x_rd, x_ru);
 260   v->height = max4(y_ld, y_lu, y_rd, y_ru) - min4(y_ld, y_lu, y_rd, y_ru);
 261   //v->bitmap = malloc((int) (ceil(v->width) * ceil(v->height) * sizeof(bool)));
 262
 263   v->width = ceil(text->tw);
 264   v->height = ceil(text->th);
 265   v->bitmap = malloc(v->width * v->height * sizeof(bool));
 266 //  printf("Allocated bitmap of %d bools for %d x %d label\n", v->width * v->height, v->width, v->height);
 267   for (int i=0; i<v->height; i++)
 268     for (int j=0; j<v->width; j++)
 269     {
 270       v->bitmap[i*v->width + j] = 1;
 271 //      printf("Writing at %d\n", i*v->width + j);
 272     }
 273 }
 274
 275 void labeller_add_point(struct symbol *sym, struct osm_object *object, z_index_t zindex)
 276 {
 277 printf("Adding point\n");
 278   if (object->type != OSM_TYPE_NODE)
 279   {
 280     printf("Warning: Point label requested on non-point object\n");
 281     return;
 282   }
 283
 284   struct request_point *r = GARY_PUSH(requests_point);
 285
 286   r->request.type = REQUEST_POINT;
 287   r->request.ind = num_requests++;
 288
 289   r->sym = sym;
 290   r->zindex = zindex;
 291
 292   r->offset_x = 0;
 293   r->offset_y = 0;
 294
 295   r->num_variants = 1;
 296   GARY_INIT(r->variants, 0);
 297
 298   struct point_variant *v = GARY_PUSH(r->variants);
 299
 300   struct osm_node *n = (struct osm_node *) object; // FIXME: Compiler warning
 301   r->x = n->x;
 302   r->y = n->y;
 303   switch (sym->type)
 304   {
 305     case SYMBOLIZER_ICON:
 306       make_bitmap_icon(v, (struct sym_icon *) sym);
 307       r->x = ((struct sym_icon *)sym)->sir.x;
 308       r->y = ((struct sym_icon *)sym)->sir.y;
 309       break;
 310     case SYMBOLIZER_POINT:
 311       make_bitmap_point(v, (struct sym_point *) sym);
 312       break;
 313     case SYMBOLIZER_TEXT: ;
 314       struct sym_text *st = (struct sym_text *) sym;
 315       struct osm_node *n = (struct osm_node *) object;
 316       make_bitmap_label(v, st);
 317     default:
 318       // FIXME
 319       return;
 320   }
 321
 322 //  printf("Inited point to [%.2f; %.2f] on %u\n", r->x, r->y, r->zindex);
 323 }
 324
 325 void labeller_add_line(struct symbol *sym, z_index_t zindex)
 326 {
 327 printf("Adding line on %u\n", zindex);
 328   struct buffer_line *b = GARY_PUSH(buffer_line);
 329   b->line = (struct sym_line *) sym;
 330   b->zindex = zindex;
 331   sym_plan(sym, zindex);
 332 }
 333
 334 void labeller_add_linelabel(struct symbol *sym, struct osm_object *o, z_index_t zindex)
 335 {
 336   if (o->type != OSM_TYPE_WAY)
 337   {
 338     // FIXME
 339     return;
 340   }
 341
 342   printf("[LAB] Labelling way %ju on %u\n", o->id, zindex);
 343   struct buffer_linelabel *ll = GARY_PUSH(buffer_linelabel);
 344   ll->way = (struct osm_way *) o;
 345   ll->label = sym;
 346   ll->zindex = zindex;
 347 }
 348
 349 void labeller_add_arealabel(struct symbol *sym, struct osm_object *o, z_index_t zindex)
 350 {
 351 printf("Adding area on %u\n", zindex);
 352   struct request_area *r = GARY_PUSH(requests_area);
 353
 354   r->request.type = REQUEST_AREA;
 355   r->request.ind = num_requests++;
 356
 357   r->o = (struct osm_multipolygon *) o;
 358   r->zindex = zindex;
 359   r->label = sym;
 360
 361   osm_obj_center(o, &(r->cx), &(r->cy));
 362
 363   GARY_INIT(r->variants, 0);
 364   struct point_variant *v = GARY_PUSH(r->variants);
 365   switch (sym->type)
 366   {
 367     case SYMBOLIZER_ICON:
 368       printf("DEBUG: Icon label\n");
 369       make_bitmap_icon(v, (struct sym_icon *) sym);
 370       break;
 371     case SYMBOLIZER_TEXT:
 372       printf("DEBUG: Text label\n");
 373       make_bitmap_label(v, (struct sym_text *) sym);
 374     default:
 375       // FIXME
 376       ;
 377   }
 378 }
 379
 380 void make_graph(void)
 381 {
 382   hash_init();
 383   struct mempool *mp_edges = mp_new(BLOCK_SIZE);
 384
 385   printf("Extracting nodes, will iterate over %lld ways\n", GARY_SIZE(buffer_line));
 386   for (uns i=0; i<GARY_SIZE(buffer_line); i++)
 387   {
 388     struct osm_way *way = (struct osm_way *) buffer_line[i].line->s.o;
 389     struct graph_node *g_prev = NULL;
 390     struct osm_node *o_prev = NULL;
 391
 392     CLIST_FOR_EACH(struct osm_ref *, ref, way->nodes)
 393     {
 394       // FIXME: Shall osm_object's type be checked here?
 395       struct osm_node *o_node = (struct osm_node *) ref->o;
 396
 397       struct graph_node *g_node = hash_find(ref->o->id);
 398       if (!g_node)
 399       {
 400         g_node = hash_new(ref->o->id);
 401         GARY_INIT(g_node->edges, 0);
 402         g_node->o = o_node;
 403         g_node->id = ref->o->id;
 404         g_node->num = num_nodes++;
 405       }
 406
 407       if (! g_prev)
 408       {
 409         g_prev = g_node;
 410         o_prev = o_node;
 411         continue;
 412       }
 413
 414       struct graph_edge *e = mp_alloc(mp_edges, sizeof(struct graph_edge)); num_edges_dbg++;
 415       e->num = num_edges++;
 416       e->id = buffer_line[i].line->s.o->id;
 417       e->color = buffer_line[i].line->color;
 418       e->length = hypot(abs(o_prev->x - o_node->x), abs(o_prev->y - o_node->y));
 419       e->visited = -1;
 420       e->prev = NULL;
 421       e->next = NULL;
 422       e->n1 = g_prev;
 423       e->n2 = g_node;
 424       e->longline = (uns) -1;
 425       e->line = buffer_line[i].line;
 426       e->dir = DIR_UNSET;
 427       e->label = NULL;
 428
 429       struct graph_edge **edge = GARY_PUSH(g_prev->edges);
 430       *edge = e;
 431       edge = GARY_PUSH(g_node->edges);
 432       *edge = e;
 433
 434       g_prev = g_node;
 435       o_prev = o_node;
 436     }
 437   }
 438
 439   printf("Made graph with %d edges\n", num_edges_dbg);
 440 }
 441
 442 void dump_graph(void)
 443 {
 444   HASH_FOR_ALL(hash, node)
 445   {
 446     printf("* Node: (%d) #%ju [%.2f; %.2f]\n", node->num, node->id, node->o->x, node->o->y);
 447     for (uns i=0; i<GARY_SIZE(node->edges); i++)
 448     {
 449       struct graph_edge *e = node->edges[i];
 450       printf("\t edge (%d) #%ju to ", e->num, e->id);
 451       if (node->edges[i]->n1->id == node->id)
 452         printf("(%d) #%ju [%.2f; %.2f]\n", e->n2->num, e->n2->id, e->n2->o->x, e->n2->o->y);
 453       else if (node->edges[i]->n2->id == node->id)
 454         printf("(%d) #%ju [%.2f; %.2f]\n", e->n1->num, e->n1->id, e->n1->o->x, e->n1->o->y);
 455       else
 456         printf("BEWARE! BEWARE! BEWARE!\n");
 457
 458       printf("\t\t");
 459       if ((node->edges[i]->label)) printf("Labelled\n");
 460       if ((node->edges[i]->label) && (node->edges[i]->label->type == SYMBOLIZER_TEXT)) printf(" labelled %s;", osm_val_decode(((struct sym_text *) node->edges[i]->label)->text));
 461       printf(" colored %d;", node->edges[i]->color);
 462       printf("   length %.2f", node->edges[i]->length);
 463       printf("\n");
 464     }
 465   }
 466   HASH_END_FOR;
 467 }
 468
 469 void label_graph(void)
 470 {
 471 printf("There are %u line labels requested\n", GARY_SIZE(buffer_linelabel));
 472   for (uns i=0; i<GARY_SIZE(buffer_linelabel); i++)
 473   {
 474     if (buffer_linelabel[i].label->type == SYMBOLIZER_TEXT)
 475     printf("Labelling nodes of way %s\n", osm_val_decode(((struct sym_text *) buffer_linelabel[i].label)->text));
 476     CLIST_FOR_EACH(struct osm_ref *, ref, buffer_linelabel[i].way->nodes)
 477     {
 478       printf("Looking for node %ju\n", ref->o->id);
 479       struct graph_node *n = hash_find(ref->o->id);
 480       if (n == NULL)
 481       {
 482         // FIXME: What shall be done?
 483       }
 484       else
 485       {
 486         printf("Searching among %u edges\n", GARY_SIZE(n->edges));
 487         for (uns j=0; j<GARY_SIZE(n->edges); j++)
 488         {
 489           if (n->edges[j]->id == buffer_linelabel[i].way->o.id)
 490           {
 491             printf("Labelling node %ju\n", n->id);
 492             n->edges[j]->label = buffer_linelabel[i].label;
 493             n->edges[j]->zindex = buffer_linelabel[i].zindex;
 494           }
 495         }
 496       }
 497     }
 498   }
 499 }
 500
 501 void bfs_edge(struct graph_edge *e, struct graph_node *node, struct graph_node *anode, enum edge_dir dir)
 502 {
 503 printf("BFS edge called for edge %d (going %d) in direction %d\n", e->num, e->dir, dir);
 504   struct graph_edge *candidate = NULL;
 505
 506   for (uns i=0; i<GARY_SIZE(node->edges); i++)
 507   {
 508     struct graph_edge *other = node->edges[i];
 509 if (e->num == 987) printf("Got label %d\n", e->num);
 510     if ((other->longline != (uns) -1) && (other->longline != e->longline)) continue;
 511 if (e->num == 987) printf("Continuing with edge %d\n", e->num);
 512
 513 printf("Testing %d ?= %d\n", other->visited, e->longline);
 514     if ((uns) other->visited != e->longline) {
 515     printf("Pushing new edge %d / %ju\n", other->num, other->id);
 516     struct graph_edge **e_ptr = GARY_PUSH(bfs_queue);
 517     *e_ptr = other;
 518     other->visited = e->longline;
 519     }
 520
 521     if (((other->n1->id == node->id) && (other->n2->id == anode->id)) ||
 522         ((other->n2->id == node->id) && (other->n1->id == anode->id)))
 523         continue;
 524
 525     if (((other->n1->id == node->id) || (other->n2->id == node->id)) &&
 526         (e->label) && (other->label) &&
 527         (e->label->type == SYMBOLIZER_TEXT) && (other->label->type == SYMBOLIZER_TEXT) &&
 528         (((struct sym_text *) e->label)->text == ((struct sym_text *) other->label)->text))
 529     {
 530       if (! candidate || (other->length > candidate->length))
 531       candidate = other;
 532     }
 533   }
 534
 535   if (candidate)
 536   {
 537 printf("New line in longline %u\n", e->longline);
 538     struct graph_edge *other = candidate;
 539       other->longline = e->longline;
 540       other->dir = dir;
 541       if (((dir == DIR_BWD) && (other->n1->id == node->id)) ||
 542           ((dir == DIR_FWD) && (other->n2->id == node->id)))
 543       {
 544         struct graph_node *swp = other->n2;
 545         other->n2 = other->n1;
 546         other->n1 = swp;
 547       }
 548
 549       switch (dir)
 550       {
 551         case DIR_BWD:
 552           e->prev = other;
 553           other->next = e;
 554           longlines[other->longline].first = other;
 555           break;
 556         case DIR_FWD:
 557           e->next = other;
 558           other->prev = e;
 559           break;
 560         default:
 561           printf("Oops\n");
 562           ASSERT(0);
 563       }
 564   }
 565 }
 566
 567 void bfs(uns longline)
 568 {
 569 printf("BFS called for longline %u\n", longline);
 570 printf("%d longlines are believed to exist, %d exist\n", num_longlines, GARY_SIZE(longlines));
 571   for (uns i=0; i<GARY_SIZE(bfs_queue); i++)
 572   {
 573     struct graph_edge *cur = bfs_queue[i];
 574     printf("Exploring new edge %d; %d remaining\n", cur->num, GARY_SIZE(bfs_queue));
 575     //ASSERT(! cur->visited);
 576
 577     cur->visited = longline;
 578
 579     if (cur->longline == (uns) -1)
 580       continue;
 581
 582     if (cur->dir == DIR_UNSET)
 583     {
 584       cur->dir = DIR_CENTER;
 585       bfs_edge(cur, cur->n1, cur->n2, DIR_BWD);
 586       bfs_edge(cur, cur->n2, cur->n1, DIR_FWD);
 587     }
 588     else
 589     {
 590       switch (cur->dir)
 591       {
 592         case DIR_BWD:
 593           bfs_edge(cur, cur->n1, cur->n2, cur->dir);
 594           break;
 595         case DIR_FWD:
 596           bfs_edge(cur, cur->n2, cur->n1, cur->dir);
 597           break;
 598         default:
 599           // FIXME
 600           ;
 601       }
 602     }
 603   }
 604 }
 605
 606 void bfs_wrapper(void)
 607 {
 608   GARY_INIT(bfs_queue, 0);
 609   GARY_INIT(longlines, 0);
 610
 611   HASH_FOR_ALL(hash, node)
 612   {
 613     for (uns i=0; i<GARY_SIZE(node->edges); i++)
 614     {
 615       if ((node->edges[i]->label) && (node->edges[i]->longline == (uns) -1))
 616       {
 617         GARY_PUSH(longlines);
 618         longlines[num_longlines].first = node->edges[i];
 619         printf("Running new BFS\n");
 620         printf("Creating longline %u\n", num_longlines);
 621         GARY_RESIZE(bfs_queue, 0);
 622         struct graph_edge **e = GARY_PUSH(bfs_queue);
 623         *e = node->edges[i];
 624         node->edges[i]->longline = num_longlines;
 625         bfs(node->edges[i]->longline);
 626         //dump_longlines();
 627         printf("Joined %d edges\n", dbg_num_hits); dbg_num_hits = 0;
 628         printf("Planned %u edges\n", GARY_SIZE(bfs_queue));
 629         num_longlines++;
 630       }
 631     }
 632   }
 633   HASH_END_FOR;
 634
 635   GARY_FREE(bfs_queue);
 636 }
 637
 638 void dump_longlines(void)
 639 {
 640 printf("*** Longlines dump\n");
 641   for (uns i=0; i<GARY_SIZE(longlines); i++)
 642   {
 643 printf("Longline %u:", i);
 644     struct graph_edge *e = longlines[i].first;
 645 if ((e->label) && (e->label->type == SYMBOLIZER_TEXT))
 646   printf(" labelled %s", osm_val_decode(((struct sym_text *) e->label)->text));
 647 printf("\n");
 648
 649     while (e)
 650     {
 651       printf("\t#%ju (%d): [%.2f; %.2f] -- [%.2f; %.2f] (dir %d)\n",
 652              e->id, e->num, e->n1->o->x, e->n1->o->y, e->n2->o->x, e->n2->o->y, e->dir);
 653
 654       e = e->next;
 655     }
 656   }
 657 }
 658
 659 struct request_line *make_new_line(void)
 660 {
 661   struct request_line *rl = GARY_PUSH(requests_line);
 662   rl->request.ind = num_requests++;
 663   rl->request.type = REQUEST_LINE;
 664   GARY_INIT(rl->sections, 0);
 665
 666   return rl;
 667 }
 668
 669 struct request_section *make_new_section(struct request_line *rl)
 670 {
 671   struct request_section *rls = GARY_PUSH(rl->sections);
 672   rls->request.ind = num_requests++;
 673   rls->request.type = REQUEST_SECTION;
 674   rls->num_segments = 0;
 675   GARY_INIT(rls->segments, 0);
 676
 677   return rls;
 678 }
 679
 680 struct request_segment *make_new_segment(struct request_section *rls, struct symbol *sym)
 681 {
 682   struct request_segment *rs = GARY_PUSH(rls->segments);
 683   rls->num_segments++;
 684
 685   rs->request.ind = num_requests++;
 686   rs->request.type = REQUEST_SEGMENT;
 687
 688   struct point_variant *v = malloc(sizeof(struct point_variant));
 689   make_bitmap(v, sym);
 690   rs->variant = v;
 691
 692   return rs;
 693 }
 694
 695 void cut_edge(struct graph_edge *e, double dist)
 696 {
 697   if (dbg_segments)
 698     printf("Cutting [%.2f; %.2f] -- [%.2f; %.2f] to dist %.2f\n", e->n1->o->x, e->n1->o->y, e->n2->o->x, e->n2->o->y, dist);
 699
 700   struct graph_edge *new = malloc(sizeof(struct graph_edge));
 701   *new = *e;
 702   e->next = new;
 703
 704   struct osm_node *n1 = e->n1->o;
 705   struct osm_node *n2 = e->n2->o;
 706
 707   // FIXME
 708   if ((n1->x == n2->x) && (n1->y == n2->y))
 709   {
 710     printf("[%.2f; %.2f] -- [%.2f; %.2f]\n", n1->x, n1->y, n2->x, n2->y);
 711     printf("Won't cut point\n");
 712     return;
 713   }
 714
 715   struct osm_node *n11 = malloc(sizeof(struct osm_node));
 716   struct graph_node *gn = malloc(sizeof(struct graph_node));
 717   gn->o = n11;
 718   double vsize = sqrt(pow(n1->x - n2->x, 2) + pow(n1->y - n2->y, 2));
 719   n11->x = n1->x + (n2->x - n1->x) / vsize * dist;
 720   n11->y = n1->y + (n2->y - n1->y) / vsize * dist;
 721
 722   e->n2 = new->n1 = gn;
 723
 724   e->length = hypot(abs(n1->x - n11->x), abs(n1->y - n11->y));
 725   new->length = hypot(abs(n11->x - n2->x), abs(n11->y - n2->y));
 726 }
 727
 728 void make_segments(void)
 729 {
 730   for (uns i=0; i<GARY_SIZE(longlines); i++)
 731   {
 732     // Skip lines which are not labelled
 733     if (! (longlines[i].first && longlines[i].first->label))
 734       continue;
 735
 736     struct request_line *request = make_new_line();
 737     struct request_section *rls = make_new_section(request);
 738     struct request_segment *rs = NULL;
 739
 740     struct graph_edge *e = longlines[i].first;
 741     double cur_length = 0;
 742
 743     struct sym_text *st = NULL;
 744     if (e->label->type == SYMBOLIZER_TEXT)
 745     {
 746       st = (struct sym_text *) e->label;
 747     }
 748     else
 749     {
 750       printf("Warning: Skipping line\n");
 751       continue;
 752       // FIXME;
 753     }
 754
 755     printf("New longline\n");
 756     while (e)
 757     {
 758       if (cur_length + e->length > conf_max_section_length + conf_max_section_overlay)
 759       {
 760         if (dbg_segments)
 761           printf("Edge too long, length is %.2f; %.2f - %.2f = %.2f\n", e->length, conf_max_section_length, cur_length, conf_max_section_length - cur_length);
 762         cut_edge(e, conf_max_section_length - cur_length);
 763       }
 764
 765       if (cur_length + e->length > conf_max_section_length)
 766       {
 767         if (dbg_segments)
 768           printf("Making new section, new length would be %f, allowed is %.2f / %.2f\n", cur_length + e->length, conf_max_section_length, conf_max_section_overlay);
 769
 770         struct osm_node *n1 = e->n1->o;
 771         struct osm_node *n2 = e->n2->o;
 772         rs = make_new_segment(rls, e->label);
 773         rs->x1 = n1->x;
 774         rs->y1 = n1->y;
 775         rs->x2 = n2->x;
 776         rs->y2 = n2->y;
 777         rs->zindex = e->zindex;
 778
 779         rs->label = malloc(sizeof(struct sym_text));
 780         *((struct sym_text *) rs->label) = *((struct sym_text *) e->label);
 781         rls = make_new_section(request);
 782         cur_length = 0;
 783       }
 784
 785       if (st && (e->length < st->tw))
 786       {
 787         e = e->next;
 788         printf("Warning: Skipping segment\n");
 789         continue;
 790       }
 791
 792       rs = make_new_segment(rls, e->label);
 793       rs->label = malloc(sizeof(struct sym_text));
 794       *((struct sym_text *) rs->label) = *((struct sym_text *) e->label);
 795
 796       rs->x1 = e->n1->o->x;
 797       rs->y1 = e->n1->o->y;
 798       rs->x2 = e->n2->o->x;
 799       rs->y2 = e->n2->o->y;
 800
 801       // FIXME: Set text rotation
 802       rs->angle = atan2(rs->x2 - rs->x1, rs->y2 - rs->y1);
 803       rs->zindex = e->zindex;
 804
 805       cur_length += e->length;
 806       e = e->next;
 807     }
 808
 809     if (request->sections[0].num_segments == 0)
 810     {
 811       // FIXME
 812       printf("WARNING: 0 segment section\n");
 813       GARY_POP(requests_line);
 814       num_requests -= 2;
 815     }
 816   }
 817 }
 818
 819 void dump_linelabel_requests(void)
 820 {
 821   for (uns i=0; i<GARY_SIZE(requests_line); i++)
 822   {
 823     if (requests_line[i].sections[0].num_segments == 0)
 824     {
 825       printf("HEY!\n");
 826       continue;
 827     }
 828     printf("Request for linelabel, %d sections\n", GARY_SIZE(requests_line[i].sections));
 829     print_label(requests_line[i].sections[0].segments[0].label);
 830     for (uns j=0; j<GARY_SIZE(requests_line[i].sections); j++)
 831     {
 832       printf("%d section, %d segments\n", j, GARY_SIZE(requests_line[i].sections[j].segments));
 833       for (uns k=0; k<GARY_SIZE(requests_line[i].sections[j].segments); k++)
 834       {
 835         struct request_segment *rs = &requests_line[i].sections[j].segments[k];
 836         printf("[%.2f; %.2f] -- [%.2f; %.2f]\n", rs->x1, rs->y1, rs->x2, rs->y2);
 837       }
 838     }
 839     printf("\n");
 840   }
 841 }
 842
 843 void dump_bitmaps(struct individual *individual)
 844 {
 845   bool *bitmap = malloc(page_width_int * page_height_int * sizeof(bool));
 846   printf("Bitmap size is %d\n", page_width_int * page_height_int);
 847   for (int i=0; i<page_height_int; i++)
 848     for (int j=0; j<page_width_int; j++)
 849       bitmap[i*page_width_int + j] = 0;
 850
 851   for (uns i=0; i<GARY_SIZE(individual->placements); i++)
 852   {
 853 fprintf(stderr, "%d-th placement\n", i);
 854     struct placement *p = &(individual->placements[i]);
 855     struct point_variant *v = NULL;
 856
 857     switch (p->request->type)
 858     {
 859       case REQUEST_SEGMENT: ;
 860         struct request_segment *rs = (struct request_segment *) p->request;
 861         v = rs->variant;
 862         break;
 863       case REQUEST_POINT: ;
 864         struct request_point *rp = (struct request_point *) p->request;
 865         v = &(rp->variants[p->variant_used]);
 866         break;
 867       case REQUEST_AREA: ;
 868         struct request_area *ra = (struct request_area *) p->request;
 869         printf("Using %d-th of %d variants\n", p->variant_used, GARY_SIZE(ra->variants));
 870         v = &(ra->variants[p->variant_used]);
 871         break;
 872       default:
 873         printf("Testing request type (dump_bitmaps): %d\n", p->request->type);
 874         ASSERT(p->request->type != REQUEST_INVALID);
 875         continue;
 876     }
 877
 878     printf("Got after with %d-th placement of request type %d\n", i, p->request->type);
 879
 880     printf("Rendering %d-th label %d x %d (w x h)\n", i, v->width, v->height);
 881     for (int row = max2(p->y, 0); row < min2(p->y + v->height, page_height_int); row++)
 882     {
 883       for (int col = max2(p->x, 0); col < min2(p->x + v->width, page_width_int); col++)
 884       {
 885         printf("Writing to %d\n", row*page_width_int + col);
 886         bitmap[row * page_width_int + col] = 1;
 887       }
 888     }
 889   }
 890
 891   errno = 0;
 892   FILE *fd_dump = fopen("dump.pbm", "w");
 893   fprintf(fd_dump, "P1\n");
 894   fprintf(fd_dump, "%d %d\n", page_width_int, page_height_int);
 895   for (int i=0; i<page_height_int; i++)
 896   {
 897     for (int j=0; j<page_width_int; j++)
 898     {
 899       fprintf(fd_dump, "%d", bitmap[(int) (i*page_width_int + j)] ? 1 : 0);
 900     }
 901     fprintf(fd_dump, "\n");
 902   }
 903   fclose(fd_dump);
 904 }
 905
 906 void dump_individual(struct individual *individual)
 907 {
 908 printf("*** Dumping INDIVIDUAL ***\n");
 909 printf("(There are %d requests)\n", num_requests);
 910   for (uns i=0; i<GARY_SIZE(individual->placements); i++)
 911   {
 912     struct placement *p = &(individual->placements[i]);
 913
 914     switch (p->request->type)
 915     {
 916       case REQUEST_POINT:
 917         printf("Point at [%.2f; %.2f] on %u\n", p->x, p->y, ((struct request_point *) p->request)->zindex);
 918         break;
 919       case REQUEST_LINE: ;
 920         struct request_line *rl = (struct request_line *) p->request;
 921         printf("Line: ");
 922         print_label(rl->sections[0].segments[0].label);
 923         break;
 924       case REQUEST_SECTION: ;
 925         printf("*");
 926         break;
 927       case REQUEST_SEGMENT: ;
 928         if (p->variant_used >= 0)
 929           printf("Segment placed at [%.2f; %.2f] on %u\n", p->x, p->y, ((struct request_segment *) p->request)->zindex);
 930         else
 931           printf("Segment not placed\n");
 932         break;
 933       case REQUEST_AREA: ;
 934         struct request_area *ra = (struct request_area *) p->request;
 935         printf("Area label ");
 936         print_label(ra->label);
 937         printf(" at [%.2f; %.2f] on %u\n", p->x, p->y, ((struct request_area *) p->request)->zindex);
 938         break;
 939       default:
 940         printf("Testing request type (dump_individual)\n");
 941         ASSERT(p->request->type != 0);
 942     }
 943   }
 944   printf("\nTotal penalty: %d\n", individual->penalty);
 945 }
 946
 947 void plan_individual(struct individual *individual)
 948 {
 949   for (uns i=0; i<GARY_SIZE(individual->placements); i++)
 950   {
 951     struct symbol *s = NULL;
 952     z_index_t zindex = 0;
 953     if (individual->placements[i].variant_used < 0) continue;
 954     switch (individual->placements[i].request->type)
 955     {
 956       case REQUEST_POINT: ;
 957         struct request_point *rp = (struct request_point *) individual->placements[i].request;
 958         s = rp->sym;
 959         zindex = rp->zindex;
 960         break;
 961       case REQUEST_SEGMENT: ;
 962         struct request_segment *rs = (struct request_segment *) individual->placements[i].request;
 963         s = rs->label;
 964         zindex = rs->zindex;
 965         break;
 966       case REQUEST_LINE: ;
 967         break;
 968       case REQUEST_AREA: ;
 969         struct request_area *ra = (struct request_area *) individual->placements[i].request;
 970         s = ra->label;
 971         zindex = ra->zindex;
 972         break;
 973       default:
 974         ASSERT(individual->placements[i].request != REQUEST_INVALID);
 975         continue;
 976     }
 977
 978 if (dbg_plan)
 979   printf("Will plan symbol at [%.2f; %.2f] on %u\n", individual->placements[i].x, individual->placements[i].y, zindex);
 980
 981     if (s) switch (s->type)
 982     {
 983       case SYMBOLIZER_POINT: ;
 984         struct sym_point *sp = (struct sym_point *) s;
 985         sp->x = individual->placements[i].x;
 986         sp->y = individual->placements[i].y;
 987         sym_plan((struct symbol *) sp, zindex);
 988         break;
 989       case SYMBOLIZER_ICON: ;
 990         struct sym_icon *si = (struct sym_icon *) s;
 991         si->sir.x = individual->placements[i].x;
 992         si->sir.y = individual->placements[i].y;
 993         sym_plan((struct symbol *) si, zindex);
 994         break;
 995       case SYMBOLIZER_TEXT: ;
 996         struct sym_text *st = (struct sym_text *) s;
 997         st->x = individual->placements[i].x;
 998         st->y = individual->placements[i].y;
 999         st->next_duplicate = NULL;
1000         if (dbg_plan) printf("Planning text %s at [%.2f; %.2f] on %u, with rotation %.2f\n", osm_val_decode(st->text), st->x, st->y, zindex, st->rotate);
1001         sym_plan((struct symbol *) st, zindex);
1002         break;
1003       default:
1004         ASSERT(s->type != SYMBOLIZER_INVALID);
1005     }
1006   }
1007
1008 }
1009
1010 void labeller_label(void)
1011 {
1012   make_graph();
1013   label_graph();
1014 //dump_graph();
1015   bfs_wrapper();
1016 //dump_longlines();
1017   make_segments();
1018 dump_linelabel_requests();
1019
1020 printf("Having %u point requests, %u line requests and %u area requests\n", GARY_SIZE(requests_point), GARY_SIZE(requests_line), GARY_SIZE(requests_area));
1021
1022   GARY_INIT(population1, conf_pop_size);
1023   GARY_INIT(population2, conf_pop_size);
1024   make_population();
1025
1026   printf("Dealing with %d requests\n", num_requests);
1027
1028 /*
1029   while (! shall_terminate())
1030   {
1031     iteration++;
1032
1033     struct individual **swp = population1;
1034     population1 = population2;
1035     population2 = swp;
1036     pop2_ind = 0;
1037   }
1038 */
1039
1040   dump_individual(population1[0]);
1041 //dump_bitmaps(population1[0]);
1042
1043   plan_individual(population1[0]);
1044
1045   labeller_cleanup();
1046
1047   return;
1048 }
1049
1050 void labeller_cleanup(void)
1051 {
1052 }
1053
1054 void make_population(void)
1055 {
1056   for (int i=0; i<conf_pop_size; i++)
1057   {
1058     printf("Making individual %d\n", i);
1059     struct individual *individual = ep_alloc(ep_individuals); init_individual(individual);
1060     population1[i] = individual;
1061
1062     int p = 0;
1063     for (uns j=0; j<GARY_SIZE(requests_point); j++)
1064     {
1065       init_placement(&(individual->placements[p++]), individual, (struct request *) &requests_point[j]);
1066     }
1067
1068     for (uns j=0; j<GARY_SIZE(requests_line); j++)
1069     {
1070       init_placement(&(individual->placements[p++]), individual, (struct request *) &requests_line[j]);
1071
1072       for (uns k=0; k<GARY_SIZE(requests_line[j].sections); k++)
1073       {
1074         init_placement(&(individual->placements[p++]), individual, (struct request *) &requests_line[j].sections[k]);
1075
1076         for (uns l=0; l<GARY_SIZE(requests_line[j].sections[k].segments); l++)
1077         {
1078           init_placement(&(individual->placements[p++]), individual, (struct request *) &requests_line[j].sections[k].segments[l]);
1079         }
1080       }
1081     }
1082
1083     for (uns j=0; j<GARY_SIZE(requests_area); j++)
1084     {
1085       init_placement(&(individual->placements[p++]), individual, (struct request *) &requests_area[j]);
1086     }
1087
1088     hide_segment_labels(individual);
1089
1090 if (p != num_requests)
1091 {
1092   printf("Say bye\n");
1093   exit(42);
1094 }
1095
1096 printf("Testing p\n");
1097     ASSERT(p == num_requests);
1098   }
1099 }
1100
1101 bool shall_terminate(void)
1102 {
1103   switch (conf_term_cond)
1104   {
1105     case TERM_COND_PENALTY:
1106       return (population1[0]->penalty < conf_penalty_bound);
1107     case TERM_COND_STAGNATION:
1108       return (abs(old_best - population1[0]->penalty) < conf_stagnation_bound);
1109     case TERM_COND_ITERATIONS:
1110       return (iteration >= conf_iteration_limit);
1111     default:
1112       // FIXME: Warn the user that no condition is set
1113       return 1;
1114   }
1115 }
1116
1117 void breed(void)
1118 {
1119   int acc = 0;
1120   int i=0;
1121   printf("%.2f\n", ((double) conf_breed_pop_size_perc/100));
1122   int conf_breed_pop_size = ((double) conf_breed_pop_size_perc/100) * conf_pop_size;
1123   struct individual **breed_buffer;
1124   while (i < conf_breed_pop_size)
1125   {
1126   printf("%d < %d, breeding\n", i, conf_breed_pop_size);
1127     int parent1 = randint(1, conf_breed_pop_size);
1128     int parent2 = randint(1, conf_breed_pop_size);
1129     printf("Will breed %d and %d, chosen of %d best of %d population (intended to be %d)\n", parent1, parent2, conf_breed_pop_size, GARY_SIZE(population1), conf_pop_size);
1130     breed_buffer = perform_crossover(population1[parent1], population1[parent2]);
1131     population2[pop2_ind++] = breed_buffer[0];
1132     population2[pop2_ind++] = breed_buffer[1];
1133     free(breed_buffer);
1134     i++;
1135   }
1136
1137   acc += conf_breed_rbest_perc;
1138
1139   return; // FIXME: DEBUG HACK
1140
1141   int remaining = (1 - acc) * (conf_pop_size * conf_breed_perc);
1142   int step = remaining / conf_pop_size;
1143   for (; i<conf_pop_size; i += 2)
1144   {
1145     printf("Asking for %d and %d of %d\n", i*step, i*(step+1), conf_pop_size);
1146     breed_buffer = perform_crossover(population1[i*step], population1[i*step+1]);
1147     population2[pop2_ind++] = breed_buffer[0];
1148     population2[pop2_ind++] = breed_buffer[1];
1149   }
1150
1151   // FIXME: Could there be one missing individual?
1152 }
1153
1154 struct individual **perform_crossover(struct individual *parent1, struct individual *parent2)
1155 {
1156   struct individual **buffer = malloc(2*sizeof(struct individual));
1157   struct individual *child1 = ep_alloc(ep_individuals); init_individual(child1);
1158   struct individual *child2 = ep_alloc(ep_individuals); init_individual(child2);
1159
1160   printf("Performing crossover\n");
1161
1162   for (uns i=0; i<GARY_SIZE(parent1->placements); i++)
1163   {
1164     printf("%dth placement out of %d\n", i, num_requests);
1165     if (! parent1->placements[i].processed)
1166     {
1167       struct placement **clos_symbols = get_closure(&(parent1->placements[i]), parent1, parent2);
1168       int x = randint(1, 2);
1169
1170       if (x == 1)
1171       {
1172         copy_symbols(clos_symbols, parent1, child1);
1173         copy_symbols(clos_symbols, parent2, child2);
1174       }
1175       else
1176       {
1177         copy_symbols(clos_symbols, parent2, child1);
1178         copy_symbols(clos_symbols, parent1, child2);
1179       }
1180       printf("Symbols copied; %lld\n", GARY_SIZE(clos_symbols));
1181       GARY_FREE(clos_symbols);
1182     }
1183
1184     if (conf_mutate_children)
1185     {
1186       if (randint(1, 1000) < conf_mutate_children_prob * 1000) perform_mutation(child1);
1187       if (randint(1, 1000) < conf_mutate_children_prob * 1000) perform_mutation(child2);
1188     }
1189   }
1190
1191   buffer[0] = child1;
1192   buffer[1] = child2;
1193   return buffer;
1194 }
1195
1196 void mutate(void)
1197 {
1198   int i = 0;
1199   int conf_mutate_pop_size = conf_mutate_pop_size_perc * conf_pop_size;
1200   while (i < conf_mutate_rbest_perc * conf_pop_size)
1201   {
1202     int ind = randint(1, conf_mutate_pop_size);
1203     copy_individual(population2[pop2_ind], population1[ind]);
1204     perform_mutation(population2[pop2_ind]);
1205     pop2_ind++;
1206   }
1207 }
1208
1209 void perform_mutation(struct individual *individual)
1210 {
1211   for (uns i=0; i<GARY_SIZE(individual->placements); i++)
1212   {
1213     int x = randint(1, 1000);
1214     int acc = 0;
1215
1216     if (x <= acc + conf_mutate_move_bound)
1217     {
1218       move_symbol(&(individual->placements[i]));
1219       continue;
1220     }
1221     acc += conf_mutate_move_bound;
1222
1223     if (x <= acc + conf_mutate_regen_bound)
1224     {
1225       gen_coords(&(individual->placements[i]));
1226       continue;
1227     }
1228     acc += conf_mutate_regen_bound;
1229
1230     if (x <= acc + conf_mutate_chvar_bound)
1231     {
1232       if (0) // if num_variants > 1
1233       {
1234         // FIXME: assign new variant
1235       }
1236     }
1237   }
1238 }
1239
1240 void elite(void)
1241 {
1242   for (int i=0; i<conf_elite_perc * conf_pop_size; i++)
1243   {
1244     population2[pop2_ind++] = population1[0];
1245   }
1246 }
1247
1248 void rank_population(void)
1249 {
1250   // FIXME
1251 }
1252
1253 void gen_coords(struct placement *p)
1254 {
1255   switch(p->request->type)
1256   {
1257     case REQUEST_POINT:
1258       gen_coords_point(p);
1259       break;
1260     case REQUEST_AREA:
1261       gen_coords_area(p);
1262       break;
1263     case REQUEST_SEGMENT:
1264       gen_coords_segment(p);
1265       break;
1266     case REQUEST_LINE:
1267       printf("Not yet implemented\n");
1268       break;
1269     default:
1270       printf("Testing request type\n");
1271       ASSERT(p->request->type != REQUEST_INVALID);
1272   }
1273 }
1274
1275 double gen_movement(void)
1276 {
1277   double m = (random() % 1000000) / 10000;
1278   m = pow(m, 1.0/3) * flip(1, -1);
1279   printf("Movement %.2f\n", m);
1280   return m;
1281 }
1282
1283 void gen_coords_point(struct placement *p)
1284 {
1285   p->x = p->x + gen_movement();
1286 }
1287
1288 void gen_coords_segment(struct placement *p)
1289 {
1290   struct request_segment *rs = (struct request_segment *) p->request;
1291   int a = flip(1, 2);
1292   p->x = (a == 1 ? rs->x1 : rs->x2);
1293   p->y = (a == 1 ? rs->y1 : rs->y2);
1294 }
1295
1296 void gen_coords_area(struct placement *p)
1297 {
1298   struct request_area *ra = (struct request_area *) p->request;
1299
1300   p->x = p->x + gen_movement();
1301   p->y = p->y + gen_movement();
1302
1303   printf("Moved label to [%.2f; %.2f] from [%.2f; %.2f]\n", p->x, p->y, ra->cx, ra->cy);
1304 }
1305
1306 struct map_part **get_parts(struct placement *symbol, struct individual *individual)
1307 {
1308   struct map_part **buffer;
1309   GARY_INIT(buffer, 0);
1310   int x_min = symbol->x / conf_part_size;
1311   int x_max = (symbol->x /*+ symbol->bitmap->width*/ + conf_part_size - 1) / conf_part_size;
1312   int y_min = symbol->y / conf_part_size;
1313   int y_max = (symbol->y /*+ symbol->bitmap->height*/ + conf_part_size - 1) / conf_part_size;
1314
1315   for (int x=x_min; x < x_max; x++)
1316     for (int y=y_min; y < y_max; y++)
1317     {
1318       struct map_part *m = GARY_PUSH(buffer);
1319       *m = individual->map[x][y];
1320     }
1321
1322   return buffer;
1323 }
1324
1325 int randint(int min, int max)
1326 {
1327   if (min == max) return min;
1328   int r = random();
1329   //printf("Returning %d + (%d %% (%d - %d)) = %d + %d %% %d = %d + %d = %d\n", min, r, max, min, min, r, max-min, min, r%(max-min), min+(r%(max-min)));
1330   return min + (r % (max - min));
1331   return (r * (max - min));
1332 }
1333
1334 struct placement **get_closure(struct placement *placement, struct individual *parent1, struct individual *parent2 UNUSED)
1335 {
1336   printf("Getting closure\n");
1337   struct placement **closure;
1338   GARY_INIT(closure, 0);
1339   bool *chosen = malloc(GARY_SIZE(parent1->placements) * sizeof(bool));
1340   chosen[placement->request->ind] = 1;
1341
1342   struct placement **p = GARY_PUSH(closure); *p = placement;
1343
1344   uns first = 0;
1345   while (first < GARY_SIZE(closure))
1346   {
1347     printf("Iterating, first is %d\n", first);
1348     struct placement **overlapping = get_overlapping(placement);
1349     filter(overlapping, chosen);
1350     for (uns j=0; j<GARY_SIZE(overlapping); j++)
1351     {
1352       p = GARY_PUSH(closure); *p = overlapping[j];
1353       chosen[overlapping[j]->request->ind] = 1;
1354     }
1355     GARY_FREE(overlapping);
1356     first++;
1357   }
1358
1359   return closure;
1360 }
1361
1362 void copy_symbols(struct placement **closure, struct individual *parent, struct individual *child)
1363 {
1364   //printf("%d\n", child->penalty);
1365   //printf("Closure size: %lld\n", GARY_SIZE(closure));
1366   for (uns i=0; i<GARY_SIZE(closure); i++)
1367   {
1368     int ind = closure[i]->request->ind;
1369     child->placements[ind] = parent->placements[ind];
1370     child->placements[ind].processed = 0;
1371   }
1372 }
1373
1374 void move_symbol(struct placement *p)
1375 {
1376   switch (p->request->type)
1377   {
1378     case REQUEST_POINT:
1379       move_symbol_point(p);
1380     case REQUEST_LINE:
1381     case REQUEST_SEGMENT:
1382     case REQUEST_AREA:
1383       printf("Not yet implemented\n");
1384     default:
1385       ASSERT(p->request->type != REQUEST_INVALID);
1386   }
1387 }
1388
1389 void move_symbol_point(struct placement *p)
1390 {
1391   p->x += (double) (move_min + randdouble()) * flip(1, -1);
1392   p->y += (double) (move_min + randdouble()) * flip(1, -1);
1393 }
1394
1395 void hide_segment_labels(struct individual *individual)
1396 {
1397   // BEWARE: This fully depends on current genetic encoding
1398
1399   int used = -1, num = -1;
1400   for (uns i=0; i<GARY_SIZE(individual->placements); i++)
1401   {
1402     switch (individual->placements[i].request->type)
1403     {
1404       case REQUEST_SECTION:
1405         used = individual->placements[i].variant_used;
1406         num = 0;
1407         break;
1408       case REQUEST_SEGMENT:
1409         if (num == used)
1410           individual->placements[i].variant_used = 0;
1411         else
1412           individual->placements[i].variant_used = -1;
1413         num++;
1414         break;
1415       default:
1416         ;
1417     }
1418   }
1419 }
1420
1421 void init_placement(struct placement *p, struct individual *individual, struct request *r)
1422 {
1423   // FIXME
1424   p->request = r;
1425   p->processed = 0;
1426   p->x = p->y = 0; // To prevent valgrind from complaining
1427   p->variant_used = 0;
1428   p->individual = individual;
1429   switch (r->type)
1430   {
1431     case REQUEST_POINT: ;
1432       struct request_point *rp = (struct request_point *) r;
1433       p->x = rp->x;
1434       p->y = rp->y;
1435       break;
1436     case REQUEST_LINE: ;
1437       break;
1438     case REQUEST_SECTION: ;
1439       struct request_section *rls = (struct request_section *) r;
1440       p->variant_used = randint(0, rls->num_segments);
1441       break;
1442     case REQUEST_SEGMENT: ;
1443       struct request_segment *rs = (struct request_segment *) r;
1444       p->x = rs->x2;
1445       p->y = rs->y2;
1446       break;
1447     case REQUEST_AREA: ;
1448       struct request_area *ra = (struct request_area *) r;
1449       p->x = ra->cx;
1450       p->y = ra->cy;
1451       p->variant_used = 0;
1452       break;
1453     default:
1454       ASSERT(p->request->type != REQUEST_INVALID);
1455       printf("Valid request: %d\n", p->request->type);
1456   }
1457
1458   gen_coords(p);
1459 //  printf("Inited placement to [%.2f; %.2f]\n", p->x, p->y);
1460 }
1461
1462 void init_individual(struct individual *i)
1463 {
1464 //printf("Initing individual\n");
1465   GARY_INIT(i->placements, num_requests);
1466   GARY_INIT(i->map, 0);
1467   i->penalty = 0; // FIXME
1468 }
1469
1470 struct placement **get_overlapping(struct placement *p UNUSED)
1471 {
1472   struct placement **buffer;
1473   GARY_INIT(buffer, 0);
1474   return buffer;
1475 }
1476
1477 void filter(struct placement **list UNUSED, bool *pred UNUSED)
1478 {
1479   // FIXME
1480 }
1481
1482 int flip(int a, int b)
1483 {
1484   return (random() % 2 ? a : b);
1485 }
1486
1487 double randdouble(void)
1488 {
1489   // FIXME: How the hell shall double in range <0, 1> be generated? O:)
1490   return 0.5;
1491 }
1492
1493 void cleanup(void)
1494 {
1495   hash_cleanup();
1496   GARY_FREE(requests_point);
1497   GARY_FREE(requests_line);
1498   GARY_FREE(requests_area);
1499 }
1500
1501 void copy_individual(struct individual *src, struct individual *dest)
1502 {
1503   src->penalty = dest->penalty;
1504   GARY_INIT(dest->placements, GARY_SIZE(src->placements));
1505   for (uns i=0; i<GARY_SIZE(src->placements); i++)
1506   {
1507     dest->placements[i] = src->placements[i];
1508   }
1509 }