NEWS

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