/**
- * collectd - src/sn-evolution.c
- * Copyright (C) 2008 Florian octo Forster
+ * 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
* 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
*
* Authors:
- * Florian octo Forster <octo at verplant.org>
+ * Florian octo Forster <ff at octo.it>
**/
#ifndef _ISOC99_SOURCE
# define _ISOC99_SOURCE
#endif
#ifndef _POSIX_C_SOURCE
-# define _POSIX_C_SOURCE 200112L
+# define _POSIX_C_SOURCE 200809L
+#endif
+#ifndef _XOPEN_SOURCE
+# define _XOPEN_SOURCE 700
#endif
#include <stdlib.h>
#include "sn_network.h"
#include "sn_random.h"
-/* Yes, this is ugly, but the GNU libc doesn't export it with the above flags.
- * */
-char *strdup (const char *s);
+#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 max_population_size = 128;
static population_t *population;
+static int evolution_threads_num = 4;
+
static int do_loop = 0;
-static void sigint_handler (int signal)
+static void sigint_handler (int signal __attribute__((unused)))
{
do_loop++;
} /* void sigint_handler */
" -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);
{
int option;
- while ((option = getopt (argc, argv, "i:o:p:P:s:h")) != -1)
+ while ((option = getopt (argc, argv, "i:o:p:P:s:t:h")) != -1)
{
switch (option)
{
{
int tmp = atoi (optarg);
if (tmp > 0)
+ {
max_population_size = tmp;
+ population_set_size (population, (size_t) max_population_size);
+ }
+ break;
+ }
+
+ 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;
}
break;
}
+ case 't':
+ {
+ int tmp = atoi (optarg);
+ if (tmp >= 1)
+ evolution_threads_num = tmp;
+ break;
+ }
+
case 'h':
default:
exit_usage (argv[0]);
assert (SN_NETWORK_INPUT_NUM (n) == inputs_num);
- if (sn_bounded_random (0, 100) <= 1)
+ if ((SN_NETWORK_INPUT_NUM (n) <= 16) && (sn_bounded_random (0, 100) <= 1))
mutate_network (n);
population_insert (population, n);
return (0);
} /* int create_offspring */
-static void *evolution_thread (void *arg)
+static void *evolution_thread (void *arg __attribute__((unused)))
{
while (do_loop == 0)
{
if (status == 0)
{
sn_network_t *n;
+ int stages_num;
+ int comparators_num;
int rating;
+ int iter;
- i = iteration_counter;
+ 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 ("After approximately %i iterations: "
- "Currently best rating: %i\n",
- i, rating);
+ printf ("Best after approximately %i iterations: "
+ "%i comparators in %i stages. Rating: %i.\n",
+ iter, comparators_num, stages_num, rating);
}
}
struct sigaction sigint_action;
struct sigaction sigterm_action;
+ 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);
+
read_options (argc, argv);
if (initial_input_file == NULL)
exit_usage (argv[0]);
sigterm_action.sa_handler = sigint_handler;
sigaction (SIGTERM, &sigterm_action, NULL);
- 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_start_listen_thread (population, NULL, NULL);
{
sn_network_t *n;
+ int tmp;
n = sn_network_read_file (initial_input_file);
if (n == NULL)
inputs_num = SN_NETWORK_INPUT_NUM(n);
+ /* 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 ((tmp % 2) != 0)
+ {
+ if (tmp == 1)
+ inputs_num_is_power_of_two = 1;
+ else
+ inputs_num_is_power_of_two = 0;
+ break;
+ }
+ tmp = tmp >> 1;
+ }
+
population_insert (population, n);
sn_network_destroy (n);
}
"=======================\n",
initial_input_file, inputs_num, max_population_size);
- evolution_start (3);
+ evolution_start (evolution_threads_num);
printf ("Exiting after %llu iterations.\n",
(unsigned long long) iteration_counter);
{
if (best_output_file != NULL)
sn_network_write_file (n, best_output_file);
- else
- sn_network_show (n);
+ sn_network_show (n);
sn_network_destroy (n);
}
}