NEWS

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