]> www.ginac.de Git - cln.git/blobdiff - benchmarks/timebench1.cc
* src/integer/conv/cl_I_from_digits.cc: Made input of all numbers in
[cln.git] / benchmarks / timebench1.cc
index e2bad84a6d5ac0f5d7f2b540da1f7fbe44a529ba..a626bf3f9da091d0df5c0cc10d09cad79a0bed5f 100644 (file)
@@ -1,17 +1,17 @@
 // Benchmarks from the LiDIA home page
 
-#include <cl_number.h>
-#include <cl_io.h>
-#include <cl_integer.h>
-#include <cl_float.h>
-#include <cl_float_io.h>
-#include <cl_real.h>
-#include <cl_real_io.h>
-#include <cl_complex.h>
-#include <cl_complex_io.h>
-#include <stdlib.h>
-#include <string.h>
-#include <cl_timing.h>
+#include <cln/number.h>
+#include <cln/io.h>
+#include <cln/integer.h>
+#include <cln/float.h>
+#include <cln/float_io.h>
+#include <cln/real.h>
+#include <cln/real_io.h>
+#include <cln/complex.h>
+#include <cln/complex_io.h>
+#include <cstdlib>
+#include <cstring>
+#include <cln/timing.h>
 
 // Timings on Linux i486 33 MHz, 1000 decimal places = 104 32-bit words.
 // Function              LiDIA       Pari       CLISP         CLN
@@ -61,42 +61,42 @@ int main (int argc, char * argv[])
        if (argc < 1)
                exit(1);
 
-       fprint(cl_stderr, "Number of repetitions: ");
-       fprintdecimal(cl_stderr, repetitions);
-       fprint(cl_stderr, "\n");
+       fprint(std::cerr, "Number of repetitions: ");
+       fprintdecimal(std::cerr, repetitions);
+       fprint(std::cerr, "\n");
 
-       cl_float_format_t prec = cl_float_format(1000);
+       float_format_t prec = float_format(1000);
 
-       fprint(cl_stderr, "pi\n");
+       fprint(std::cerr, "pi\n");
        { cl_F p;
-         { CL_TIMING; p = cl_pi(prec); }
+         { CL_TIMING; p = pi(prec); }
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
-             { cl_F p = cl_pi(prec); }
+             { cl_F p = pi(prec); }
          }
          cout << p << endl << endl;
        }
 
-       fprint(cl_stderr, "gamma\n");
+       fprint(std::cerr, "gamma\n");
        { cl_F p;
-         { CL_TIMING; p = cl_eulerconst(prec); }
+         { CL_TIMING; p = eulerconst(prec); }
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
-             { cl_F p = cl_eulerconst(prec); }
+             { cl_F p = eulerconst(prec); }
          }
          cout << p << endl << endl;
        }
 
-       fprint(cl_stderr, "e\n");
-       { cl_F p = cl_exp1(prec);
+       fprint(std::cerr, "e\n");
+       { cl_F p = exp1(prec);
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
-             { cl_F p = cl_exp1(prec); }
+             { cl_F p = exp1(prec); }
          }
          cout << p << endl << endl;
        }
 
-       fprint(cl_stderr, "sqrt(3)\n");
+       fprint(std::cerr, "sqrt(3)\n");
        { cl_R p = sqrt(cl_float(3,prec));
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
@@ -105,7 +105,7 @@ int main (int argc, char * argv[])
          cout << p << endl << endl;
        }
 
-       fprint(cl_stderr, "exp(log(2))\n");
+       fprint(std::cerr, "exp(log(2))\n");
        { cl_N p = exp(log(cl_float(2,prec)));
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
@@ -114,7 +114,7 @@ int main (int argc, char * argv[])
          cout << p << endl << endl;
        }
 
-       fprint(cl_stderr, "log(exp(2))\n");
+       fprint(std::cerr, "log(exp(2))\n");
        { cl_N p = log(exp(cl_float(2,prec)));
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
@@ -123,25 +123,25 @@ int main (int argc, char * argv[])
          cout << p << endl << endl;
        }
 
-       fprint(cl_stderr, "sin(pi/3)\n");
-       { cl_R p = sin(cl_pi(prec)/3);
+       fprint(std::cerr, "sin(pi/3)\n");
+       { cl_R p = sin(pi(prec)/3);
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
-             { cl_R p = sin(cl_pi(prec)/3); }
+             { cl_R p = sin(pi(prec)/3); }
          }
          cout << p << endl << endl;
        }
 
-       fprint(cl_stderr, "cos(pi/3)\n");
-       { cl_R p = cos(cl_pi(prec)/3);
+       fprint(std::cerr, "cos(pi/3)\n");
+       { cl_R p = cos(pi(prec)/3);
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
-             { cl_R p = cos(cl_pi(prec)/3); }
+             { cl_R p = cos(pi(prec)/3); }
          }
          cout << p << endl << endl;
        }
 
-       fprint(cl_stderr, "arcsin(sqrt(3)/2)\n");
+       fprint(std::cerr, "arcsin(sqrt(3)/2)\n");
        { cl_N p = asin(sqrt(cl_float(3,prec))/2);
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
@@ -150,7 +150,7 @@ int main (int argc, char * argv[])
          cout << p << endl << endl;
        }
 
-       fprint(cl_stderr, "arccos(sqrt(3)/2)\n");
+       fprint(std::cerr, "arccos(sqrt(3)/2)\n");
        { cl_N p = acos(sqrt(cl_float(3,prec))/2);
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
@@ -159,7 +159,7 @@ int main (int argc, char * argv[])
          cout << p << endl << endl;
        }
 
-       fprint(cl_stderr, "sinh(log(2))\n");
+       fprint(std::cerr, "sinh(log(2))\n");
        { cl_N p = sinh(log(cl_float(2,prec)));
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
@@ -168,7 +168,7 @@ int main (int argc, char * argv[])
          cout << p << endl << endl;
        }
 
-       fprint(cl_stderr, "cosh(log(2))\n");
+       fprint(std::cerr, "cosh(log(2))\n");
        { cl_N p = cosh(log(cl_float(2,prec)));
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
@@ -177,20 +177,20 @@ int main (int argc, char * argv[])
          cout << p << endl << endl;
        }
 
-       fprint(cl_stderr, "arsinh(pi)\n");
-       { cl_N p = asinh(cl_pi(prec));
+       fprint(std::cerr, "arsinh(pi)\n");
+       { cl_N p = asinh(pi(prec));
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
-             { cl_N p = asinh(cl_pi(prec)); }
+             { cl_N p = asinh(pi(prec)); }
          }
          cout << p << endl << endl;
        }
 
-       fprint(cl_stderr, "arcosh(pi)\n");
-       { cl_N p = acosh(cl_pi(prec));
+       fprint(std::cerr, "arcosh(pi)\n");
+       { cl_N p = acosh(pi(prec));
          { CL_TIMING;
            for (int rep = repetitions; rep > 0; rep--)
-             { cl_N p = acosh(cl_pi(prec)); }
+             { cl_N p = acosh(pi(prec)); }
          }
          cout << p << endl << endl;
        }