-#define _ISOC99_SOURCE
-#define _POSIX_C_SOURCE 200112L
+/**
+ * libsortnetwork - src/sn-evolution.c
+ * Copyright (C) 2008-2010 Florian octo Forster
+ *
+ * This program is free software; you can redistribute it and/or modify it
+ * under the terms of the GNU General Public License as published by the
+ * Free Software Foundation; only version 2 of the License is applicable.
+ *
+ * This program is distributed in the hope that it will be useful, but
+ * WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+ * General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License along
+ * with this program; if not, write to the Free Software Foundation, Inc.,
+ * 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
+ *
+ * Authors:
+ * Florian octo Forster <ff at octo.it>
+ **/
+
+#ifndef _ISOC99_SOURCE
+# define _ISOC99_SOURCE
+#endif
+#ifndef _POSIX_C_SOURCE
+# define _POSIX_C_SOURCE 200809L
+#endif
+#ifndef _XOPEN_SOURCE
+# define _XOPEN_SOURCE 700
+#endif
#include <stdlib.h>
#include <stdio.h>
+#include <stdint.h>
#include <string.h>
#include <sys/types.h>
#include <unistd.h>
#include <signal.h>
#include <assert.h>
+#include <limits.h>
+
+#include <pthread.h>
+
+#include <population.h>
#include "sn_network.h"
+#include "sn_random.h"
-struct population_entry_s
-{
- sn_network_t *network;
- int rating;
-};
-typedef struct population_entry_s population_entry_t;
+#if !defined(__GNUC__) || !__GNUC__
+# define __attribute__(x) /**/
+#endif
+
+static uint64_t iteration_counter = 0;
+static int inputs_num = 16;
+static int inputs_num_is_power_of_two = 1;
+
+static char *initial_input_file = NULL;
+static char *best_output_file = NULL;
+
+static int stats_interval = 0;
-static int iterations_num = 1000000;
static int max_population_size = 128;
-static int inputs_num = 16;
+static population_t *population;
+
+static int evolution_threads_num = 4;
-static population_entry_t *population = NULL;
-int population_size = 0;
+static int do_loop = 0;
-static void sigint_handler (int signal)
+static void sigint_handler (int signal __attribute__((unused)))
{
- iterations_num = 0;
+ do_loop++;
} /* void sigint_handler */
-static int init_random (void)
+static void exit_usage (const char *name)
{
- int fd;
- unsigned int r;
+ printf ("Usage: %s [options]\n"
+ "\n"
+ "Valid options are:\n"
+ " -i <file> Initial input file (REQUIRED)\n"
+ " -o <file> Write the current best solution to <file>\n"
+ " -p <num> Size of the population (default: 128)\n"
+ " -P <peer> Send individuals to <peer> (may be repeated)\n"
+ " -t <num> Number of threads (default: 4)\n"
+ "\n",
+ name);
+ exit (1);
+} /* void exit_usage */
- fd = open ("/dev/random", O_RDONLY);
- if (fd < 0)
+int read_options (int argc, char **argv)
+{
+ int option;
+
+ while ((option = getopt (argc, argv, "i:o:p:P:s:t:h")) != -1)
{
- perror ("open");
- return (-1);
- }
+ switch (option)
+ {
+ case 'i':
+ {
+ if (initial_input_file != NULL)
+ free (initial_input_file);
+ initial_input_file = strdup (optarg);
+ break;
+ }
- read (fd, (void *) &r, sizeof (r));
- close (fd);
+ case 'o':
+ {
+ if (best_output_file != NULL)
+ free (best_output_file);
+ best_output_file = strdup (optarg);
+ break;
+ }
- srand (r);
+ case 'p':
+ {
+ int tmp = atoi (optarg);
+ if (tmp > 0)
+ {
+ max_population_size = tmp;
+ population_set_size (population, (size_t) max_population_size);
+ }
+ break;
+ }
- return (0);
-} /* int init_random */
+ case 'P':
+ {
+ int status;
+
+ status = population_add_peer (population, optarg, /* port = */ NULL);
+ if (status != 0)
+ {
+ fprintf (stderr, "population_add_peer failed with status %i.\n",
+ status);
+ }
+ break;
+ }
-static int bounded_random (int upper_bound)
-{
- double r = ((double) rand ()) / ((double) RAND_MAX);
- return ((int) (r * upper_bound));
-}
+ case 's':
+ {
+ int tmp = atoi (optarg);
+ if (tmp > 0)
+ stats_interval = tmp;
+ break;
+ }
-static void exit_usage (const char *name)
-{
- printf ("%s <file0>\n", name);
- exit (1);
-} /* void exit_usage */
+ case 't':
+ {
+ int tmp = atoi (optarg);
+ if (tmp >= 1)
+ evolution_threads_num = tmp;
+ break;
+ }
+
+ case 'h':
+ default:
+ exit_usage (argv[0]);
+ }
+ }
+
+ return (0);
+} /* int read_options */
static int rate_network (const sn_network_t *n)
{
return (rate);
} /* int rate_network */
-static int population_print_stats (int iterations)
+static int mutate_network (sn_network_t *n)
{
- int best = -1;
- int total = 0;
- int i;
-
- for (i = 0; i < population_size; i++)
+ sn_network_t *n_copy;
+ int stage_index;
+ sn_stage_t *s;
+ int comparator_index;
+ int status;
+
+ n_copy = sn_network_clone (n);
+ if (n_copy == NULL)
{
- if ((best == -1) || (best > population[i].rating))
- best = population[i].rating;
- total += population[i].rating;
+ fprintf (stderr, "mutate_network: sn_network_clone failed.\n");
+ return (-1);
}
- printf ("Iterations: %6i; Best: %i; Average: %.2f;\n",
- iterations, best, ((double) total) / ((double) population_size));
-
- return (0);
-} /* int population_print_stats */
+ stage_index = sn_bounded_random (0, SN_NETWORK_STAGE_NUM (n_copy) - 1);
+ s = SN_NETWORK_STAGE_GET (n_copy, stage_index);
-static int insert_into_population (sn_network_t *n)
-{
- int rating;
- int i;
+ comparator_index = sn_bounded_random (0, SN_STAGE_COMP_NUM (s) - 1);
+ sn_stage_comparator_remove (s, comparator_index);
- rating = rate_network (n);
+ status = sn_network_brute_force_check (n_copy);
+
+ sn_network_destroy (n_copy);
- if (population_size < max_population_size)
- {
- population[population_size].network = n;
- population[population_size].rating = rating;
- population_size++;
- return (0);
- }
-
- for (i = 0; i < population_size; i++)
- if (population[i].rating > rating)
- break;
+ if (status < 0)
+ return (-1);
+ else if (status > 0) /* Mutated network does not sort anymore. */
+ return (1);
- if (i < population_size)
- {
- sn_network_destroy (population[i].network);
- population[i].network = n;
- population[i].rating = rating;
- }
- else
- {
- sn_network_destroy (n);
- }
+ /* We saved one comparator \o/ Let's do the same change on the original
+ * network. */
+ s = SN_NETWORK_STAGE_GET (n, stage_index);
+ sn_stage_comparator_remove (s, comparator_index);
return (0);
-} /* int insert_into_population */
+} /* int mutate_network */
static int create_offspring (void)
{
- int p0;
- int p1;
+ sn_network_t *p0;
+ sn_network_t *p1;
sn_network_t *n;
- p0 = bounded_random (population_size);
- p1 = bounded_random (population_size);
+ p0 = population_get_random (population);
+ p1 = population_get_random (population);
+
+ assert (p0 != NULL);
+ assert (p1 != NULL);
+
+ /* combine the two parents */
+ n = sn_network_combine (p0, p1);
- n = sn_network_combine (population[p0].network, population[p1].network);
+ sn_network_destroy (p0);
+ sn_network_destroy (p1);
+ /* randomly chose an input and do a min- or max-cut. */
while (SN_NETWORK_INPUT_NUM (n) > inputs_num)
{
int pos;
enum sn_network_cut_dir_e dir;
- pos = bounded_random (SN_NETWORK_INPUT_NUM (n));
- dir = (bounded_random (2) == 0) ? DIR_MIN : DIR_MAX;
+ pos = sn_bounded_random (0, SN_NETWORK_INPUT_NUM (n) - 1);
+ dir = (sn_bounded_random (0, 1) == 0) ? DIR_MIN : DIR_MAX;
assert ((pos >= 0) && (pos < SN_NETWORK_INPUT_NUM (n)));
sn_network_cut_at (n, pos, dir);
}
+ /* compress the network to get a compact representation */
sn_network_compress (n);
assert (SN_NETWORK_INPUT_NUM (n) == inputs_num);
- insert_into_population (n);
+ if ((SN_NETWORK_INPUT_NUM (n) <= 16) && (sn_bounded_random (0, 100) <= 1))
+ mutate_network (n);
+
+ population_insert (population, n);
+
+ sn_network_destroy (n);
return (0);
} /* int create_offspring */
-static int start_evolution (void)
+static void *evolution_thread (void *arg __attribute__((unused)))
+{
+ while (do_loop == 0)
+ {
+ create_offspring ();
+ /* XXX: Not synchronized! */
+ iteration_counter++;
+ }
+
+ return ((void *) 0);
+} /* int start_evolution */
+
+static int evolution_start (int threads_num)
{
+ pthread_t threads[threads_num]; /* C99 ftw! */
int i;
- for (i = 0; i < iterations_num; i++)
+ for (i = 0; i < threads_num; i++)
{
- if ((i % 1000) == 0)
- population_print_stats (i);
+ int status;
- create_offspring ();
+ status = pthread_create (&threads[i], /* attr = */ NULL,
+ evolution_thread, /* arg = */ NULL);
+ if (status != 0)
+ {
+ fprintf (stderr, "evolution_start: pthread_create[%i] failed "
+ "with status %i.\n",
+ i, status);
+ threads[i] = 0;
+ }
+ }
+
+ while (do_loop == 0)
+ {
+ int status;
+
+ status = sleep (1);
+ if (status == 0)
+ {
+ sn_network_t *n;
+ int stages_num;
+ int comparators_num;
+ int rating;
+ int iter;
+
+ iter = iteration_counter;
+
+ n = population_get_fittest (population);
+
+ rating = rate_network (n);
+
+ stages_num = SN_NETWORK_STAGE_NUM (n);
+ comparators_num = 0;
+ for (i = 0; i < stages_num; i++)
+ {
+ sn_stage_t *s;
+
+ s = SN_NETWORK_STAGE_GET (n, i);
+ comparators_num += SN_STAGE_COMP_NUM (s);
+ }
+
+ sn_network_destroy (n);
+
+ printf ("Best after approximately %i iterations: "
+ "%i comparators in %i stages. Rating: %i.\n",
+ iter, comparators_num, stages_num, rating);
+ }
+ }
+
+ for (i = 0; i < threads_num; i++)
+ {
+ if (threads[i] == 0)
+ continue;
+ pthread_join (threads[i], /* value_ptr = */ NULL);
}
return (0);
-} /* int start_evolution */
+} /* int evolution_start */
int main (int argc, char **argv)
{
struct sigaction sigint_action;
+ struct sigaction sigterm_action;
- if (argc != 2)
- exit_usage (argv[0]);
+ population = population_create ((pi_rate_f) rate_network,
+ (pi_copy_f) sn_network_clone,
+ (pi_free_f) sn_network_destroy);
+ if (population == NULL)
+ {
+ fprintf (stderr, "population_create failed.\n");
+ return (1);
+ }
+
+ population_set_serialization (population,
+ (pi_serialize_f) sn_network_serialize,
+ (pi_unserialize_f) sn_network_unserialize);
- init_random ();
+ read_options (argc, argv);
+ if (initial_input_file == NULL)
+ exit_usage (argv[0]);
memset (&sigint_action, '\0', sizeof (sigint_action));
sigint_action.sa_handler = sigint_handler;
sigaction (SIGINT, &sigint_action, NULL);
- population = (population_entry_t *) malloc (max_population_size
- * sizeof (population_entry_t));
- if (population == NULL)
- {
- printf ("Malloc failed.\n");
- return (-1);
- }
- memset (population, '\0', max_population_size
- * sizeof (population_entry_t));
+ memset (&sigterm_action, '\0', sizeof (sigterm_action));
+ sigterm_action.sa_handler = sigint_handler;
+ sigaction (SIGTERM, &sigterm_action, NULL);
+
+ population_start_listen_thread (population, NULL, NULL);
{
- sn_network_t *n = sn_network_read_file (argv[1]);
+ sn_network_t *n;
+ int tmp;
+
+ n = sn_network_read_file (initial_input_file);
if (n == NULL)
{
printf ("n == NULL\n");
return (1);
}
- population[0].network = n;
- population[0].rating = rate_network (n);
- population_size++;
- }
-
- start_evolution ();
- {
- int i;
- int best_rate = -1;
- int best_index = -1;
+ inputs_num = SN_NETWORK_INPUT_NUM(n);
- for (i = 0; i < population_size; i++)
+ /* Determine if `inputs_num' is a power of two. If so, more merge
+ * algorithms can be used. */
+ tmp = inputs_num;
+ inputs_num_is_power_of_two = 1;
+ while (tmp > 0)
{
- if ((best_rate == -1) || (best_rate > population[i].rating))
+ if ((tmp % 2) != 0)
{
- best_rate = population[i].rating;
- best_index = i;
+ if (tmp == 1)
+ inputs_num_is_power_of_two = 1;
+ else
+ inputs_num_is_power_of_two = 0;
+ break;
}
+ tmp = tmp >> 1;
}
- sn_network_show (population[best_index].network);
+ population_insert (population, n);
+ sn_network_destroy (n);
+ }
+
+ printf ("Current configuration:\n"
+ " Initial network: %s\n"
+ " Number of inputs: %3i\n"
+ " Population size: %3i\n"
+ "=======================\n",
+ initial_input_file, inputs_num, max_population_size);
+
+ evolution_start (evolution_threads_num);
+
+ printf ("Exiting after %llu iterations.\n",
+ (unsigned long long) iteration_counter);
+
+ {
+ sn_network_t *n;
+
+ n = population_get_fittest (population);
+ if (n != NULL)
+ {
+ if (best_output_file != NULL)
+ sn_network_write_file (n, best_output_file);
+ sn_network_show (n);
+ sn_network_destroy (n);
+ }
}
+ population_destroy (population);
+
return (0);
} /* int main */