//#include <cln/modinteger.h>
//#include <cln/numtheory.h>
#include <cln/timing.h>
+using namespace cln;
+#include <iostream>
+using namespace std;
#define DUMP(expr) \
- fprint(stdout, #expr" = "); fprint(stdout, expr); fprint(stdout, "\n");
+ fprint(cout, #expr" = "); fprint(cout, expr); fprint(cout, "\n");
int main (int argc, char* argv[])
{
#if 0
cl_F archimedes = pi((float_format_t)10000);
- fprint(stdout, archimedes);
- fprint(stdout, "\n");
+ fprint(cout, archimedes);
+ fprint(cout, "\n");
#endif
#if 0
cl_FF x1 = "-0.2173f0";
cl_FF x2 = "5.5084f9";
cl_FF y = "-1.19698f9";
- fprint(stdout, "x1 = "); print_float_binary(stdout,x1); fprint(stdout, " = "); fprint(stdout,x1); fprint(stdout, "\n");
- fprint(stdout, "x2 = "); print_float_binary(stdout,x2); fprint(stdout, " = "); fprint(stdout,x2); fprint(stdout, "\n");
- fprint(stdout, "y = "); print_float_binary(stdout,y); fprint(stdout, " = "); fprint(stdout,y); fprint(stdout, "\n");
+ fprint(cout, "x1 = "); print_float_binary(cout,x1); fprint(cout, " = "); fprint(cout,x1); fprint(cout, "\n");
+ fprint(cout, "x2 = "); print_float_binary(cout,x2); fprint(cout, " = "); fprint(cout,x2); fprint(cout, "\n");
+ fprint(cout, "y = "); print_float_binary(cout,y); fprint(cout, " = "); fprint(cout,y); fprint(cout, "\n");
cl_FF x = x1*x2;
- fprint(stdout, "x1*x2 = "); print_float_binary(stdout,x); fprint(stdout, " = "); fprint(stdout,x); fprint(stdout, "\n");
+ fprint(cout, "x1*x2 = "); print_float_binary(cout,x); fprint(cout, " = "); fprint(cout,x); fprint(cout, "\n");
#endif
#if 0
cl_I y = ++x;
x *= 2;
x++;
- fprint(stdout, "x = "); fprint(stdout, x); fprint(stdout, "\n");
- fprint(stdout, "y = "); fprint(stdout, y); fprint(stdout, "\n");
+ fprint(cout, "x = "); fprint(cout, x); fprint(cout, "\n");
+ fprint(cout, "y = "); fprint(cout, y); fprint(cout, "\n");
#endif
#if 0
cl_I a = "77371252437321868671713407";
cl_I w;
- cl_boolean squarep = isqrt(a,&w);
+ bool squarep = isqrt(a,&w);
DUMP(squarep);
DUMP(w);
DUMP(expt_pos(w,2) <= a);
#if 0
float_format_t f = float_format(atoi(argv[1]));
extern cl_LF zeta3 (uintC len);
- uintC len = (uintL)f/intDsize+1;
+ uintC len = (uintC)f/intDsize+1;
{ CL_TIMING; cout << zeta(2,f) << endl; }
{ CL_TIMING; cout << expt(pi(f),2)/6 << endl; }
{ CL_TIMING; cout << zeta(3,f) << endl; }