1 2000-05-xx, version 1.1.0
2 =========================
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.
10 * Adjusted several break-even points to make better use of a faster libgmp and
11 better match present-day CPUs.
14 2000-01-13, version 1.0.3
15 =========================
21 `cl_I doublefactorial (uintL n)'.
23 Implementation changes
24 ----------------------
26 * Fix several little configuration errors.
28 * Fix some compilation errors with gcc-2.96 prereleases.
31 1999-09-07, version 1.0.2
32 =========================
38 `cl_I numerator (const cl_RA&)'
39 `cl_I denominator (const cl_RA&)'.
40 Suggested by Richard Kreckel and Sylvain Pion.
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.
46 Implementation changes
47 ----------------------
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.
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.
58 1999-06-09, version 1.0.1
59 =========================
64 * Speed up the functions `logand', `lognand', `logandc2'.
66 Implementation changes
67 ----------------------
69 * Fix compilation errors with gcc-2.7.2, egcs-1.1.2 and gcc-2.95.
71 * Fix compilation errors on HPPA, MIPS, some versions of DEC Alpha, OpenBSD,
74 * Fix compilation errors on Win32. Thanks to Mumit Khan.
77 1999-01-12, version 1.0
78 =======================
83 * New include files, containing I/O functions:
87 <cl_sfloat_io.h> <cl_ffloat_io.h> <cl_dfloat_io.h> <cl_lfloat_io.h>
90 The file <cl_io.h> now contains only I/O functions for characters and
93 * To access the base ring of a univariate polynomial ring, now use
94 `R->basering()' instead of `R->basering'.
96 * Implement `plusp', `max', `min' for the classes cl_F, cl_SF, cl_FF, cl_DF,
99 * Implement `abs' for the class cl_RA.
101 * Replaced `read_number' with specialized functions `read_complex',
102 `read_real', `read_float', `read_rational', `read_integer'.
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)'
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)'
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.)
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
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&)'
132 Implementation changes
133 ----------------------
135 * Fix for `cl_I ash (const cl_I& x, const cl_I& n)' when |n| is between
137 (Reported by Thomas Papanikolaou.)
139 * Fixed a memory leak in the long-float truncate function.
140 (Reported by Keith Briggs.)
142 * Fixed a memory leak in the modular integer exponentiation.
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.
147 * Portability fixes for m68k platforms.
148 (Reported and fixed by Roman Hodek.)
153 * Restructured the sources. The subpackages
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.
163 * The configuration script can be called with option `--enable-shared',
164 to build CLN as a shared library.
166 * Updated the LiDIA interface.
169 1998-05-07, version 0.98
170 ========================
172 Implementation changes
173 ----------------------
175 * Portability fixes for 64-bit platforms like DEC Alpha.
176 (Reported by John Cremona and Thomas Papanikolaou.)
185 * Portability fixes for new ISO/ANSI C++ compilers.
187 * Updated the LiDIA interface for LiDIA version 1.3 and fixed a couple of
194 Implementation changes
195 ----------------------
197 * Portability fix for i686 platforms.
198 (Reported by Michael Pfeifer.)
208 <cl_ring.h>, declaring general rings,
209 <cl_GV.h>, <cl_SV.h>, defining vectors,
210 <cl_univpoly.h>, defining univariate polynomial rings.
215 * Speed up the internal computation of ln(2) and ln(10) by about 20%.
217 * Speed up the computation of exp (for 1000 digits or more) by about 10%.
219 Implementation changes
220 ----------------------
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.)
226 * Portability fix for RS/6000 platforms.
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%.
241 Implementation changes
242 ----------------------
244 * Portability fix for SunOS 4 and other BSD platforms.
253 * Fixed a bug: endless recursion in `cl_boolean sqrtp (const cl_RA&, cl_RA*)'.
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.)
259 Implementation changes
260 ----------------------
262 * Fix static initialization order bug. Requires g++ 2.7.0 or newer
263 (g++ 2.8.0 or newer on HPPA).
265 * New method `debug_print()', for pretty printing of CLN objects, intended
266 to be called from the debugger.
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'.
282 * Speed up `cl_zeta' using Cohen-Villegas-Zagier convergence acceleration.
284 Implementation changes
285 ----------------------
287 * New function `cl_print', for pretty printing of CLN objects, intended to be
288 called from the debugger.
290 * Portability fixes for 64-bit platforms like DEC Alpha.
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.
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.)
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.
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).
321 * New function `cl_I binomial (uintL n, uintL k)' which returns the
322 binomial coefficient (n choose k).
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'.
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'.
336 Implementation changes
337 ----------------------
339 * Moved the sources to subdirectories.
340 (Suggested by Jörg Arndt.)
345 * New benchmark for better comparing LiDIA, Pari and CLN.
347 * Added an interface to LiDIA, allows using CLN as a replacement of libI.
348 (Suggested by Thomas Papanikolaou.)
350 * Added an ILOG Talk interface, for interactive use of the library.
359 * New include file <cl_modinteger.h>, defining modular integer rings.
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.
367 * Speed up square root of large numbers by use of Newton's algorithm.
369 * Speed up multiplication and division of large long-floats by small integers.
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.
376 * Fixed several bugs in the transcendental functions routines. In particular,
377 the `log' function went into an endless loop.
379 * Fixed a bug: The cl_I -> cl_DF and cl_I -> cl_LF conversions didn't work
380 correctly for negative integers.
389 * New include file <cl_numtheory.h>, containing number theoretic functions,
390 for now only the Jacobi symbol.
395 * Speed up squaring of large numbers by 30%.
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.
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.
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.
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.
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.
414 * Fixed a bug: The function `uintL power2p (const cl_I&)' returned wrong
415 values for arguments x = 2^n, n >= 29.
417 Implementation changes
418 ----------------------
420 * Speed up CLN by using the low-level addition/subtraction, shift and
421 multiplication loops from GNU gmp-2.0.2.
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
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).
438 * New include file <cl_timing.h>, contains timing facilities.
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.
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.
451 * Speed up the function `cl_boolean sqrtp (const cl_I&, cl_I*)' using a
452 trick from Henri Cohen.
454 * Fixed an endless loop in the function `uint32 gcd (uint32, uint32)'
455 which occurred when one of the arguments is zero.
457 * Fixed an endless loop in the function `cl_I factorial (uintL)'.
459 Implementation changes
460 ----------------------
462 * <cl_io.h> now uses iostreams by default instead of stdio FILE pointers.
463 (Reported by Tito Flagella.)
465 * Fixed linking error when compiling without -O.
466 (Reported by Tito Flagella.)
468 * Fixed syntax error in <cl_input.h>.
469 (Reported by Tito Flagella.)
471 * Fixed syntax errors in src/cl_DS.h and src/cl_FF_plus.h.
472 (Reported by Marcus Daniels.)
474 * More portability fixes.
476 * Configuration now supports the --srcdir option.
477 (Reported by Marcus Daniels.)
483 * Compiles well on Sparc, using gcc.
485 * Mail to beta testers.
491 * Compiles well on Linux, using gcc.