
1997-02-21
==========

Algorithmic changes
-------------------

* Speed up the computation of the elementary transcendental functions:
  exp, sinh, cosh, asinh, atanh (for 100 digits) by about 30%,
  log (for 100 digits) by about 25%,
  sin, cos (for 1000 or 10000 digits) by about 20%,
  asin, acos, atan (for 1000 digits) by about 25%.

Implementation changes
----------------------

* Portability fix for SunOS 4 and other BSD platforms.


1997-01-31
==========

Algorithmic changes
-------------------

* Fixed a bug: endless recursion in `cl_boolean sqrtp (const cl_RA&, cl_RA*)'.

* Fixed a bug: Only the first 85% of the digits returned by `cl_eulerconst'
  for moderate precision (< 11370 digits) were correct.
  (Reported by Thomas Papanikolaou.)

Implementation changes
----------------------

* Fix static initialization order bug. Requires g++.

* New method `debug_print()', for pretty printing of CLN objects, intended
  to be called from the debugger.


1997-01-07
==========

Functional changes
------------------

* New functions `float cl_float_approx (const cl_R&)',
                `double cl_double_approx (const cl_R&)'
  for converting a real number to a C `float' or `double'.

Algorithmic changes
-------------------

* Speed up `cl_zeta' using Cohen-Villegas-Zagier convergence acceleration.

Implementation changes
----------------------

* New function `cl_print', for pretty printing of CLN objects, intended to be
  called from the debugger.

* Portability fixes for 64-bit platforms like DEC Alpha.


1996-11-29
==========

Functional changes
------------------

* New include files:
  <cl_string.h>, defining strings,
  <cl_symbol.h>, defining symbols,
  <cl_proplist.h>, defining property lists,
  <cl_condition.h>, defining conditions/exceptions.

* New functions `cl_F cl_catalanconst ()',
                `cl_F cl_catalanconst (const cl_F&)',
                `cl_F cl_catalanconst (cl_float_format_t)'
  which return Catalan's constant.
  (Suggested by Thomas Papanikolaou.)

* New functions `cl_F cl_zeta (int s)',
                `cl_F cl_zeta (int s, const cl_F&)',
                `cl_F cl_zeta (int s, cl_float_format_t)'
  which return the Riemann zeta function at an integral point s>1.

* New functions `cl_F cl_exp1 ()',
                `cl_F cl_exp1 (const cl_F&)',
                `cl_F cl_exp1 (cl_float_format_t)'
  which return e = exp(1).

* New function `cl_I binomial (uintL n, uintL k)' which returns the
  binomial coefficient (n choose k).

* New functions `int cl_I_to_int (const cl_I&)',
       `unsigned int cl_I_to_uint (const cl_I&)',
                `long cl_I_to_long (const cl_I&)',
       `unsigned long cl_I_to_ulong (const cl_I&)'
  for converting an integer to a C `int' or `long'.

* New functions `float cl_float_approx (const cl_I&)',
                `float cl_float_approx (const cl_RA&)',
                `double cl_double_approx (const cl_I&)',
                `double cl_double_approx (const cl_RA&)'
  for converting a rational number to a C `float' or `double'.

Implementation changes
----------------------

* Moved the sources to subdirectories.
  (Suggested by Jörg Arndt.)

Other changes
-------------

* New benchmark for better comparing LiDIA, Pari and CLN.

* Added an interface to LiDIA, allows using CLN as a replacement of libI.
  (Suggested by Thomas Papanikolaou.)

* Added an ILOG Talk interface, for interactive use of the library.


1996-10-13
==========

Functional changes
------------------

* New include file <cl_modinteger.h>, defining modular integer rings.

* New functions `cl_F cl_eulerconst ()', `cl_F cl_eulerconst (const cl_F&)',
  `cl_F cl_eulerconst (cl_float_format_t)' which return Euler's constant.

Algorithmic changes
-------------------

* Speed up square root of large numbers by use of Newton's algorithm.

* Speed up multiplication and division of large long-floats by small integers.

* Speed up the computation of pi, e, exp, log, sin, cos, tan, sinh, cosh, tanh,
  asin, acos, atan, asinh, acosh, atanh.
  All of them are now computed with asymptotically fast algorithms, of
  bit complexity O(log(N)^2*M(N)) = O(N^(1+epsilon)) for N digits.

* Fixed several bugs in the transcendental functions routines. In particular,
  the `log' function went into an endless loop.

* Fixed a bug: The cl_I -> cl_DF and cl_I -> cl_LF conversions didn't work
  correctly for negative integers.


1996-09-08
==========

Functional changes
------------------

* New include file <cl_numtheory.h>, containing number theoretic functions,
  for now only the Jacobi symbol.

Algorithmic changes
-------------------

* Speed up squaring of large numbers by 30%.

* Speed up division of large numbers by use of Newton's algorithm. The
  improvement is noticeable for numbers with at least about 1000 words =
  32000 bits = 10000 decimal digits.

* Speed up the binary-to-decimal conversion of large numbers by use of
  a divide-and-conquer method. The improvement is noticeable for numbers
  of at least 1250 words = 40000 bits = 12000 decimal digits.

* Fixed a bug: The functions `cl_boolean minusp (const type&)' and
  `cl_boolean zerop (const type&)' were missing for type=cl_I and type=cl_RA.

* Fixed a bug: The function `cl_boolean logtest (const cl_I&, const cl_I&)'
  didn't work correctly if both arguments were bignums.

* Fixed a bug: The function `cl_I mod (const cl_I&, const cl_I&)', when called
  as `mod(-n*m,m)' (n>0), returned m instead of 0.

* Fixed a bug: The function `uintL power2p (const cl_I&)' returned wrong
  values for arguments x = 2^n, n >= 29.

Implementation changes
----------------------

* Speed up CLN by using the low-level addition/subtraction, shift and
  multiplication loops from GNU gmp-2.0.2.

* Fixed a severe reference counting bug in the constructors
  `type::type (const long)' and `type::type (const unsigned long)'
  for type=cl_number, type=cl_N, type=cl_RA, type=cl_I that most often
  caused a core dump.

* CLN's word sequences can be chosen to be stored big-endian in memory
  (like it was up to now) or little-endian (needed for interfacing to GMP).


1996-05-20
==========

Functional changes
------------------

* New include file <cl_timing.h>, contains timing facilities.

Algorithmic changes
-------------------

* Speed up the multiplication of large numbers by use of an FFT based
  algorithm (Schönhage-Strassen). The improvement is noticeable when both
  factors have at least 1250 words = 40000 bits = 12000 decimal digits.

* Speed up the functions `cl_I gcd (const cl_I&, const cl_I&)' and
  `cl_I xgcd (const cl_I&, const cl_I&, cl_I*, cl_I*)' by a factor
  of 1.5. Simple trick.

* Speed up the function `cl_boolean sqrtp (const cl_I&, cl_I*)' using a
  trick from Henri Cohen.

* Fixed an endless loop in the function `uint32 gcd (uint32, uint32)'
  which occurred when one of the arguments is zero.

* Fixed an endless loop in the function `cl_I factorial (uintL)'.

Implementation changes
----------------------

* <cl_io.h> now uses iostreams by default instead of stdio FILE pointers.
  (Reported by Tito Flagella.)

* Fixed linking error when compiling without -O.
  (Reported by Tito Flagella.)

* Fixed syntax error in <cl_input.h>.
  (Reported by Tito Flagella.)

* Fixed syntax errors in src/cl_DS.h and src/cl_FF_plus.h.
  (Reported by Marcus Daniels.)

* More portability fixes.

* Configuration now supports the --srcdir option.
  (Reported by Marcus Daniels.)


1996-01-03
==========

* Compiles well on Sparc, using gcc.

* Mail to beta testers.


1995-12-06
==========

* Compiles well on Linux, using gcc.

