NEWS

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