* Fixed bug in example

[ginac.git] / check / times.cpp

diff --git a/check/times.cpp b/check/times.cpp
index ef038f3089599d7d84a9ba6ac8df217ee2d17f85..ce5730b9f5d7bf48707d849fb3ab6f9ab324a652 100644 (file)
--- a/check/times.cpp
+++ b/check/times.cpp
@@ -3,7 +3,7 @@
  *  Main program that calls the individual timings. */
 
 /*
- *  GiNaC Copyright (C) 1999-2001 Johannes Gutenberg University Mainz, Germany
+ *  GiNaC Copyright (C) 1999-2003 Johannes Gutenberg University Mainz, Germany
  *
  *  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
@@ -23,8 +23,35 @@
 #include <stdexcept>
 #include "times.h"
 
+/** Generate a random amount of symbols and destroy them again immediatly.
+ *  This operation effectively makes the serial numbers of all subsequent
+ *  symbols unpredictable.  If the serials are unpredictable, then so are
+ *  their hash values.  If the hash values are unpredictable, then so are
+ *  the canonical orderings.  If the canonical orderings are unpredictable,
+ *  all subsequent times are subject to some variation.  This variation,
+ *  however is natural and desireable for two reasons: First, we cannot know
+ *  how many symbols have been generated before in real world computations.
+ *  Second, the following timings are subject to some semi-random variation
+ *  anyways because short timings need to be repeated until enough time has
+ *  gone by for the measurement to be reliable.  During this process the serial
+ *  numbers will be shifted anyways in a semi-random way.  It is better not
+ *  to lull the user in a false sense of reproducibility and instead confront
+ *  her with the normal variation to be expected.
+ */
+void randomify_symbol_serials()
+{
+       srand((unsigned)time(NULL));
+       const int m = rand() % 666;
+       for (int s=0; s<m; ++s ) {
+               symbol* tmp = new symbol;
+               delete tmp;
+       }
+}
+
 int main()
 {
+       randomify_symbol_serials();
+
        unsigned result = 0;
        
 #define TIME(which) \
@@ -39,6 +66,7 @@ try { \
        TIME(gammaseries)
        TIME(vandermonde)
        TIME(toeplitz)
+       TIME(hashmap)
        TIME(lw_A)
        TIME(lw_B)
        TIME(lw_C)
@@ -57,6 +85,7 @@ try { \
        TIME(lw_Q)
        TIME(lw_Qprime)
        TIME(antipode)
+       TIME(fateman_expand)
        
        if (result) {
                cout << "Error: something went wrong. ";