NEWS

   1 2002-05-28, version 1.1.5
   2 =========================
   3
   4 Other changes
   5 -------------
   6
   7 * Fix compilation issues with GCC-3.1.
   8
   9
  10 2002-01-04, version 1.1.4
  11 =========================
  12
  13 Other changes
  14 -------------
  15
  16 * Fix compilation issues with GMP-4.0 and on IBM s390.
  17
  18 * Updated infrastructure to autoconf-2.52.
  19
  20
  21 2001-11-05, version 1.1.3
  22 =========================
  23
  24 Implementation changes
  25 ----------------------
  26
  27 * Fix a computation error of sin and sinh for small arguments.
  28   (Reported by Christopher Kennedy.)
  29
  30 * Module ordering works with gcc-3.0 -fuse-cxa-atexit now.
  31
  32 * -DNO_ASM works properly on Sparc.
  33
  34
  35 2001-07-25, version 1.1.2
  36
  37 Implementation changes
  38 ----------------------
  39
  40 * Minor cleanups and portability issues.
  41   (No need to upgrade if you have been happy so far.)
  42
  43
  44 2001-05-31, version 1.1.1
  45 =========================
  46
  47 Implementation changes
  48 ----------------------
  49
  50 * Minor cleanups for GCC 3.0 and libstdc++-v3.
  51
  52 * Fixes several little bugs and portability issues.
  53
  54
  55 2000-12-14, version 1.1.0
  56 =========================
  57
  58 Functional changes
  59 ------------------
  60
  61 * ISO/IEC 14882 fashion adjustments:
  62   Put everything into namespace cln.  All fundamental data types still
  63   carry their old names.  Other non-macro identifiers are now written as
  64   cln::foo instead of cl_foo, except where the cl_ comes from a data type
  65   (as in cl_I_to_int()).  Headers are installed into a separate
  66   directory, #include <cln/foo.h> instead of <cl_foo.h>.  Applications
  67   must be manually ported to the new style.  We apologize for the
  68   inconvenience.
  69
  70 Implementation changes
  71 ----------------------
  72
  73 * Removed internal copy of mpn subdirectory tree from libgmp-2.0.2.
  74   Configuring with --with-gmp now relies on an installed libgmp library
  75   version 3 or above.  We recommend libgmp-3.1 or newer.
  76
  77 * Adjusted several break-even points to make better use of a faster libgmp and
  78   better match present-day CPUs.
  79
  80 * Fix several errors for architectures with 64-bit wide words.
  81   (Others than Alpha, which worked already.)
  82
  83 * Fix compilation errors with current GCC-3.0 snapshots.
  84
  85 Other changes
  86 -------------
  87
  88 * Added package tools: script cln-config and autoconf macro AC_PATH_CLN
  89   (in file cln.m4).
  90
  91
  92 2000-01-13, version 1.0.3
  93 =========================
  94
  95 Functional changes
  96 ------------------
  97
  98 * New function
  99   `cl_I doublefactorial (uintL n)'.
 100
 101 Implementation changes
 102 ----------------------
 103
 104 * Fix several little configuration errors.
 105
 106 * Fix some compilation errors with gcc-2.96 prereleases.
 107
 108
 109 1999-09-07, version 1.0.2
 110 =========================
 111
 112 Functional changes
 113 ------------------
 114
 115 * New functions
 116   `cl_I numerator (const cl_RA&)'
 117   `cl_I denominator (const cl_RA&)'.
 118   Suggested by Richard Kreckel and Sylvain Pion.
 119
 120 * New function `cl_equal_hashcode' for the classes cl_N, cl_R, cl_RA, cl_I,
 121   cl_F, cl_SF, cl_FF, cl_DF, cl_LF.
 122   Suggested by Richard Kreckel.
 123
 124 Implementation changes
 125 ----------------------
 126
 127 * Fix an endless loop when either of the functions
 128   `cl_boolean rootp (const cl_RA&, uintL, cl_RA*)',
 129   `cl_boolean rootp (const cl_RA&, const cl_I&, cl_RA*)'
 130   was called with an integer argument. Thanks to Richard Kreckel.
 131
 132 * Fix a bug in the addition and subtraction of rational numbers which could
 133   cause unnormalized rational numbers like 3/1 to be created.
 134
 135
 136 1999-06-09, version 1.0.1
 137 =========================
 138
 139 Algorithmic changes
 140 -------------------
 141
 142 * Speed up the functions `logand', `lognand', `logandc2'.
 143
 144 Implementation changes
 145 ----------------------
 146
 147 * Fix compilation errors with gcc-2.7.2, egcs-1.1.2 and gcc-2.95.
 148
 149 * Fix compilation errors on HPPA, MIPS, some versions of DEC Alpha, OpenBSD,
 150   and SPARC with LiDIA.
 151
 152 * Fix compilation errors on Win32. Thanks to Mumit Khan.
 153
 154
 155 1999-01-12, version 1.0
 156 =======================
 157
 158 Functional changes
 159 ------------------
 160
 161 * New include files, containing I/O functions:
 162     <cl_integer_io.h>
 163     <cl_rational_io.h>
 164     <cl_float_io.h>
 165     <cl_sfloat_io.h> <cl_ffloat_io.h> <cl_dfloat_io.h> <cl_lfloat_io.h>
 166     <cl_real_io.h>
 167     <cl_complex_io.h>
 168   The file <cl_io.h> now contains only I/O functions for characters and
 169   C integers.
 170
 171 * To access the base ring of a univariate polynomial ring, now use
 172   `R->basering()' instead of `R->basering'.
 173
 174 * Implement `plusp', `max', `min' for the classes cl_F, cl_SF, cl_FF, cl_DF,
 175   cl_LF, cl_RA, cl_I.
 176
 177 * Implement `abs' for the class cl_RA.
 178
 179 * Replaced `read_number' with specialized functions `read_complex',
 180   `read_real', `read_float', `read_rational', `read_integer'.
 181
 182 * Replaced the functions
 183   `void fprint (cl_ostream stream, int x)'
 184   `void fprint (cl_ostream stream, unsigned int x)'
 185   `void fprint (cl_ostream stream, long x)'
 186   `void fprint (cl_ostream stream, unsigned long x)'
 187   with
 188   `void fprintdecimal (cl_ostream stream, int x)'
 189   `void fprintdecimal (cl_ostream stream, unsigned int x)'
 190   `void fprintdecimal (cl_ostream stream, long x)'
 191   `void fprintdecimal (cl_ostream stream, unsigned long x)'
 192
 193 Algorithmic changes
 194 -------------------
 195
 196 * The function `cl_I xgcd (const cl_I& a, const cl_I& b, cl_I* u, cl_I* v)'
 197   now normalizes the resulting u and v to be of minimum absolute value.
 198   (Suggested by the LiDIA group.)
 199
 200 * The conversion from string to cl_F, cl_R, cl_N now chooses the format of
 201   the resulting float depending on the number of significant digits present
 202   in the input string.
 203
 204 * Speed up the functions
 205   `cl_R operator* (const cl_RA&, const cl_LF&)'
 206   `cl_R operator* (const cl_LF&, const cl_RA&)'
 207   `cl_R operator/ (const cl_RA&, const cl_LF&)'
 208   `cl_LF operator/ (const cl_LF&, const cl_RA&)'
 209
 210 Implementation changes
 211 ----------------------
 212
 213 * Fix for `cl_I ash (const cl_I& x, const cl_I& n)' when |n| is between
 214   2^29 and 2^37.
 215   (Reported by Thomas Papanikolaou.)
 216
 217 * Fixed a memory leak in the long-float truncate function.
 218   (Reported by Keith Briggs.)
 219
 220 * Fixed a memory leak in the modular integer exponentiation.
 221
 222 * Reduced the stack space consumption. The default stack size (typically
 223   8 MB on Unix) should now be sufficient in order to run all CLN programs.
 224
 225 * Portability fixes for m68k platforms.
 226   (Reported and fixed by Roman Hodek.)
 227
 228 Other changes
 229 -------------
 230
 231 * Restructured the sources. The subpackages
 232     - base
 233     - base + integer
 234     - base + integer + rational
 235     - base + integer + rational + float
 236     - base + integer + rational + float + real
 237     - base + integer + rational + float + real + complex
 238     - base + integer + modinteger
 239   are now self-contained.
 240
 241 * The configuration script can be called with option `--enable-shared',
 242   to build CLN as a shared library.
 243
 244 * Updated the LiDIA interface.
 245
 246
 247 1998-05-07, version 0.98
 248 ========================
 249
 250 Implementation changes
 251 ----------------------
 252
 253 * Portability fixes for 64-bit platforms like DEC Alpha.
 254   (Reported by John Cremona and Thomas Papanikolaou.)
 255
 256
 257 1998-02-27
 258 ==========
 259
 260 Other changes
 261 -------------
 262
 263 * Portability fixes for new ISO/ANSI C++ compilers.
 264
 265 * Updated the LiDIA interface for LiDIA version 1.3 and fixed a couple of
 266   bugs.
 267
 268
 269 1997-09-06
 270 ==========
 271
 272 Implementation changes
 273 ----------------------
 274
 275 * Portability fix for i686 platforms.
 276   (Reported by Michael Pfeifer.)
 277
 278
 279 1997-09-01
 280 ==========
 281
 282 Functional changes
 283 ------------------
 284
 285 * New include files:
 286   <cl_ring.h>, declaring general rings,
 287   <cl_GV.h>, <cl_SV.h>, defining vectors,
 288   <cl_univpoly.h>, defining univariate polynomial rings.
 289
 290 Algorithmic changes
 291 -------------------
 292
 293 * Speed up the internal computation of ln(2) and ln(10) by about 20%.
 294
 295 * Speed up the computation of exp (for 1000 digits or more) by about 10%.
 296
 297 Implementation changes
 298 ----------------------
 299
 300 * Portability fix for 64-bit platforms like DEC Alpha:
 301   Fixed a bug: endless loop in `sqrt' for large precision.
 302   (Reported by Blake Jones.)
 303
 304 * Portability fix for RS/6000 platforms.
 305
 306
 307 1997-02-21
 308 ==========
 309
 310 Algorithmic changes
 311 -------------------
 312
 313 * Speed up the computation of the elementary transcendental functions:
 314   exp, sinh, cosh, asinh, atanh (for 100 digits) by about 30%,
 315   log (for 100 digits) by about 25%,
 316   sin, cos (for 1000 or 10000 digits) by about 20%,
 317   asin, acos, atan (for 1000 digits) by about 25%.
 318
 319 Implementation changes
 320 ----------------------
 321
 322 * Portability fix for SunOS 4 and other BSD platforms.
 323
 324
 325 1997-01-31
 326 ==========
 327
 328 Algorithmic changes
 329 -------------------
 330
 331 * Fixed a bug: endless recursion in `cl_boolean sqrtp (const cl_RA&, cl_RA*)'.
 332
 333 * Fixed a bug: Only the first 85% of the digits returned by `cl_eulerconst'
 334   for moderate precision (< 11370 digits) were correct.
 335   (Reported by Thomas Papanikolaou.)
 336
 337 Implementation changes
 338 ----------------------
 339
 340 * Fix static initialization order bug. Requires g++ 2.7.0 or newer
 341   (g++ 2.8.0 or newer on HPPA).
 342
 343 * New method `debug_print()', for pretty printing of CLN objects, intended
 344   to be called from the debugger.
 345
 346
 347 1997-01-07
 348 ==========
 349
 350 Functional changes
 351 ------------------
 352
 353 * New functions `float cl_float_approx (const cl_R&)',
 354                 `double cl_double_approx (const cl_R&)'
 355   for converting a real number to a C `float' or `double'.
 356
 357 Algorithmic changes
 358 -------------------
 359
 360 * Speed up `cl_zeta' using Cohen-Villegas-Zagier convergence acceleration.
 361
 362 Implementation changes
 363 ----------------------
 364
 365 * New function `cl_print', for pretty printing of CLN objects, intended to be
 366   called from the debugger.
 367
 368 * Portability fixes for 64-bit platforms like DEC Alpha.
 369
 370
 371 1996-11-29
 372 ==========
 373
 374 Functional changes
 375 ------------------
 376
 377 * New include files:
 378   <cl_string.h>, defining strings,
 379   <cl_symbol.h>, defining symbols,
 380   <cl_proplist.h>, defining property lists,
 381   <cl_condition.h>, defining conditions/exceptions.
 382
 383 * New functions `cl_F cl_catalanconst ()',
 384                 `cl_F cl_catalanconst (const cl_F&)',
 385                 `cl_F cl_catalanconst (cl_float_format_t)'
 386   which return Catalan's constant.
 387   (Suggested by Thomas Papanikolaou.)
 388
 389 * New functions `cl_F cl_zeta (int s)',
 390                 `cl_F cl_zeta (int s, const cl_F&)',
 391                 `cl_F cl_zeta (int s, cl_float_format_t)'
 392   which return the Riemann zeta function at an integral point s>1.
 393
 394 * New functions `cl_F cl_exp1 ()',
 395                 `cl_F cl_exp1 (const cl_F&)',
 396                 `cl_F cl_exp1 (cl_float_format_t)'
 397   which return e = exp(1).
 398
 399 * New function `cl_I binomial (uintL n, uintL k)' which returns the
 400   binomial coefficient (n choose k).
 401
 402 * New functions `int cl_I_to_int (const cl_I&)',
 403        `unsigned int cl_I_to_uint (const cl_I&)',
 404                 `long cl_I_to_long (const cl_I&)',
 405        `unsigned long cl_I_to_ulong (const cl_I&)'
 406   for converting an integer to a C `int' or `long'.
 407
 408 * New functions `float cl_float_approx (const cl_I&)',
 409                 `float cl_float_approx (const cl_RA&)',
 410                 `double cl_double_approx (const cl_I&)',
 411                 `double cl_double_approx (const cl_RA&)'
 412   for converting a rational number to a C `float' or `double'.
 413
 414 Implementation changes
 415 ----------------------
 416
 417 * Moved the sources to subdirectories.
 418   (Suggested by Jörg Arndt.)
 419
 420 Other changes
 421 -------------
 422
 423 * New benchmark for better comparing LiDIA, Pari and CLN.
 424
 425 * Added an interface to LiDIA, allows using CLN as a replacement of libI.
 426   (Suggested by Thomas Papanikolaou.)
 427
 428 * Added an ILOG Talk interface, for interactive use of the library.
 429
 430
 431 1996-10-13
 432 ==========
 433
 434 Functional changes
 435 ------------------
 436
 437 * New include file <cl_modinteger.h>, defining modular integer rings.
 438
 439 * New functions `cl_F cl_eulerconst ()', `cl_F cl_eulerconst (const cl_F&)',
 440   `cl_F cl_eulerconst (cl_float_format_t)' which return Euler's constant.
 441
 442 Algorithmic changes
 443 -------------------
 444
 445 * Speed up square root of large numbers by use of Newton's algorithm.
 446
 447 * Speed up multiplication and division of large long-floats by small integers.
 448
 449 * Speed up the computation of pi, e, exp, log, sin, cos, tan, sinh, cosh, tanh,
 450   asin, acos, atan, asinh, acosh, atanh.
 451   All of them are now computed with asymptotically fast algorithms, of
 452   bit complexity O(log(N)^2*M(N)) = O(N^(1+epsilon)) for N digits.
 453
 454 * Fixed several bugs in the transcendental functions routines. In particular,
 455   the `log' function went into an endless loop.
 456
 457 * Fixed a bug: The cl_I -> cl_DF and cl_I -> cl_LF conversions didn't work
 458   correctly for negative integers.
 459
 460
 461 1996-09-08
 462 ==========
 463
 464 Functional changes
 465 ------------------
 466
 467 * New include file <cl_numtheory.h>, containing number theoretic functions,
 468   for now only the Jacobi symbol.
 469
 470 Algorithmic changes
 471 -------------------
 472
 473 * Speed up squaring of large numbers by 30%.
 474
 475 * Speed up division of large numbers by use of Newton's algorithm. The
 476   improvement is noticeable for numbers with at least about 1000 words =
 477   32000 bits = 10000 decimal digits.
 478
 479 * Speed up the binary-to-decimal conversion of large numbers by use of
 480   a divide-and-conquer method. The improvement is noticeable for numbers
 481   of at least 1250 words = 40000 bits = 12000 decimal digits.
 482
 483 * Fixed a bug: The functions `cl_boolean minusp (const type&)' and
 484   `cl_boolean zerop (const type&)' were missing for type=cl_I and type=cl_RA.
 485
 486 * Fixed a bug: The function `cl_boolean logtest (const cl_I&, const cl_I&)'
 487   didn't work correctly if both arguments were bignums.
 488
 489 * Fixed a bug: The function `cl_I mod (const cl_I&, const cl_I&)', when called
 490   as `mod(-n*m,m)' (n>0), returned m instead of 0.
 491
 492 * Fixed a bug: The function `uintL power2p (const cl_I&)' returned wrong
 493   values for arguments x = 2^n, n >= 29.
 494
 495 Implementation changes
 496 ----------------------
 497
 498 * Speed up CLN by using the low-level addition/subtraction, shift and
 499   multiplication loops from GNU gmp-2.0.2.
 500
 501 * Fixed a severe reference counting bug in the constructors
 502   `type::type (const long)' and `type::type (const unsigned long)'
 503   for type=cl_number, type=cl_N, type=cl_RA, type=cl_I that most often
 504   caused a core dump.
 505
 506 * CLN's word sequences can be chosen to be stored big-endian in memory
 507   (like it was up to now) or little-endian (needed for interfacing to GMP).
 508
 509
 510 1996-05-20
 511 ==========
 512
 513 Functional changes
 514 ------------------
 515
 516 * New include file <cl_timing.h>, contains timing facilities.
 517
 518 Algorithmic changes
 519 -------------------
 520
 521 * Speed up the multiplication of large numbers by use of an FFT based
 522   algorithm (Schönhage-Strassen). The improvement is noticeable when both
 523   factors have at least 1250 words = 40000 bits = 12000 decimal digits.
 524
 525 * Speed up the functions `cl_I gcd (const cl_I&, const cl_I&)' and
 526   `cl_I xgcd (const cl_I&, const cl_I&, cl_I*, cl_I*)' by a factor
 527   of 1.5. Simple trick.
 528
 529 * Speed up the function `cl_boolean sqrtp (const cl_I&, cl_I*)' using a
 530   trick from Henri Cohen.
 531
 532 * Fixed an endless loop in the function `uint32 gcd (uint32, uint32)'
 533   which occurred when one of the arguments is zero.
 534
 535 * Fixed an endless loop in the function `cl_I factorial (uintL)'.
 536
 537 Implementation changes
 538 ----------------------
 539
 540 * <cl_io.h> now uses iostreams by default instead of stdio FILE pointers.
 541   (Reported by Tito Flagella.)
 542
 543 * Fixed linking error when compiling without -O.
 544   (Reported by Tito Flagella.)
 545
 546 * Fixed syntax error in <cl_input.h>.
 547   (Reported by Tito Flagella.)
 548
 549 * Fixed syntax errors in src/cl_DS.h and src/cl_FF_plus.h.
 550   (Reported by Marcus Daniels.)
 551
 552 * More portability fixes.
 553
 554 * Configuration now supports the --srcdir option.
 555   (Reported by Marcus Daniels.)
 556
 557
 558 1996-01-03
 559 ==========
 560
 561 * Compiles well on Sparc, using gcc.
 562
 563 * Mail to beta testers.
 564
 565
 566 1995-12-06
 567 ==========
 568
 569 * Compiles well on Linux, using gcc.
 570