t_chaos.cpp

Go to the documentation of this file.

/*****************************************************************************\
*                                                                             *
*  Name   : test_chaos                                                        *
*  Author : Chris Koeritz                                                     *
*                                                                             *
*******************************************************************************
* Copyright (c) 1992-$now By Author.  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; either version 2 of   *
* the License or (at your option) any later version.  This is online at:      *
*     http://www.fsf.org/copyleft/gpl.html                                    *
* Please send any updates to: fred@gruntose.com                               *
\*****************************************************************************/

//#define DEBUG_CHAOS

#include <basis/function.h>
#include <basis/chaos.h>
#include <basis/guards.h>
#include <basis/istring.h>
#include <opsystem/application_shell.h>
#include <loggers/console_logger.h>
#include <data_struct/static_memory_gremlin.h>

#define MAX_RANDOM_NUMBER 500
#define MAX_TEST_CYCLES 1083414
#define AVG_EXPECTED (double(MAX_TEST_CYCLES) / double(MAX_RANDOM_NUMBER))
#define VARIATION_ALLOWED (AVG_EXPECTED * 0.05)

#define LOG(s) EMERGENCY_LOG(program_wide_logger(), s)

class test_chaos : public application_shell
{
public:
  test_chaos() : application_shell(class_name()) {}
  IMPLEMENT_CLASS_NAME("test_chaos");
  virtual int execute();
};

int test_chaos::execute()
{
#ifdef DEBUG_CHAOS
  LOG(istring(istring::SPRINTF, "average expected=%f, "
      "variation allowed=%f", AVG_EXPECTED, VARIATION_ALLOWED));
#endif
  int results[MAX_RANDOM_NUMBER];
  for (int k = 0; k < MAX_RANDOM_NUMBER; k++) results[k] = 0;
  chaos randomizer;

  for (int i = 0; i < MAX_TEST_CYCLES; i++) {
    // first test if exclusivity is ensured...
    int res = randomizer.exclusive(0, MAX_RANDOM_NUMBER-1);
    if ( (res <= 0) || (res >= MAX_RANDOM_NUMBER-1) )
      deadly_error(__WHERE__.s(), "exclusive test", "out of bounds");

    // then test for our statistics.
    int base = randomizer.inclusive(-1000, 1000);
      // pick a base for the number below.
    res = randomizer.inclusive(base, base + MAX_RANDOM_NUMBER - 1);
    if ( (res < base) || (res > base + MAX_RANDOM_NUMBER - 1) )
      deadly_error(__WHERE__.s(), "exclusive test", "out of bounds");
    results[res - base]++;
  }
#ifdef DEBUG_CHAOS
  LOG("Anomalies:");
#endif
  int printed_any = false;
  for (int j = 0; j < MAX_RANDOM_NUMBER; j++) {
    if (absolute_value(results[j] - AVG_EXPECTED) > VARIATION_ALLOWED) {
      printed_any++;
#ifdef DEBUG_CHAOS
      LOG(istring(istring::SPRINTF, "%d: difference=%f",
          j, double(results[j] - AVG_EXPECTED)));
#endif
    }
  }
#ifdef DEBUG_CHAOS
  if (!printed_any) LOG("None");
#endif

  if (printed_any > MAX_RANDOM_NUMBER / 20)
    deadly_error(__WHERE__.s(), "probability anomalies",
        "more than the allowed number exist");
  guards::alert_message("chaos:: works for those functions tested.");
  return 0;
}

HOOPLE_MAIN(test_chaos, )

---