You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.

714 lines
18 KiB

Remove exception hooks in favor of real C++ exceptions: * include/cln/exception.h: New file... * include/cln/cln.h: ...included here... * include/cln/rational.h: ...and here. * include/cln/abort.h: Removed. * include/cln/integer.h (ash_exception, exquo_exception): Add exception types. * include/cln/float.h (floating_point_exception, floating_point_nan_exception, floating_point_overflow_exception, floating_point_underflow_exception): Likewise. * include/cln/number_io.h (read_number_exception, read_number_bad_syntax_exception, read_number_junk_exception, read_number_eof_exception): Likewise. * include/cln/object.h (CL_DEFINE_CONVERTER): Replace dynamic assertion with a compile-time assertion. * include/cln/GV.h: Replace nonreturning functions with exceptions. * include/cln/modinteger.h: Likewise. * include/cln/SV.h: Likewise. * include/cln/ring.h: Likewise. * include/cln/string.h: Likewise. * include/cln/univpoly.h: Likewise. * src/base/cl_abort.cc: Removed. * src/base/cl_N_err_d0.cc: Removed. * src/base/cl_d0_exception.cc: New file. * src/base/cl_as_err.cc: Removed. * src/base/cl_as_exception.cc: New file. * src/base/cl_notreached.cc: Removed. * src/base/cl_notreached_exception.cc: New file. * src/base/input/cl_read_err_bad.cc: Removed. * src/base/input/cl_read_bad_syntax_exception.cc: New file. * src/base/input/cl_read_err_junk.cc: Removed. * src/base/input/cl_read_junk_exception.cc: New file. * src/base/input/cl_read_err_eof.cc: Removed. * src/base/input/cl_read_eof_exception.cc: New file. * src/base/cl_N.h (cl_as_error): Removed (see cln/exception.h). * src/base/macros.h (NOTREACHED): Throw. * src/base/cl_malloc.cc (xmalloc): Throw. * src/base/digitseq/cl_2DS_div.cc: Throw. * src/base/digitseq/cl_DS_div.cc: Throw. * src/base/digitseq/cl_DS_mul.cc: Throw. * src/base/digitseq/cl_DS_mul_fftc.h: Throw. * src/base/digitseq/cl_DS_mul_fftcs.h: Throw. * src/base/digitseq/cl_DS_mul_fftm.h: Throw. * src/base/digitseq/cl_DS_mul_fftp.h: Throw. * src/base/digitseq/cl_DS_mul_fftp3.h: Throw. * src/base/digitseq/cl_DS_mul_fftp3m.h: Throw. * src/base/digitseq/cl_DS_mul_fftr.h: Throw. * src/base/digitseq/cl_DS_mul_nuss.h: Throw. * src/base/digitseq/cl_DS_recipsqrt.cc: Throw. * src/base/digitseq/cl_DS_sqrt.cc: Throw. * src/base/hash/cl_hash.h: Throw. * src/base/hash/cl_hash1.h: Throw. * src/base/hash/cl_hash1weak.h: Throw. * src/base/hash/cl_hash2.h: Throw. * src/base/hash/cl_hash2weak.h: Throw. * src/base/hash/cl_hashset.h: Throw. * src/base/hash/cl_hashuniq.h: Throw. * src/base/hash/cl_hashuniqweak.h: Throw. * src/base/proplist/cl_pl_add.cc: Throw. * src/base/ring/cl_no_ring.cc: Throw. * src/base/string/cl_spushstring.h: Throw. * src/base/symbol/cl_symbol.cc: Throw. * src/integer/bitwise/cl_I_ash.h: Removed (see cln/integer.h.) * src/integer/bitwise/cl_I_asherr.cc: Removed. * src/integer/bitwise/cl_I_ash_exception.cc: New file. * src/integer/division/cl_I_exquoerr.cc: Removed. * src/integer/division/cl_I_exquo_exception.cc: New file. * src/integer/cl_I.h: Throw. * src/integer/division/cl_I_exquopos.cc: Throw. * src/integer/bitwise/cl_I_logbitp_I.cc: Throw. * src/integer/bitwise/cl_I_ash.cc: Throw. * src/integer/bitwise/cl_I_ash_I.cc: Throw. * src/integer/division/cl_I_exquo.cc: Throw. * src/integer/gcd/cl_I_gcd_aux2.cc: Throw. * src/integer/conv/cl_I_to_L.cc: Throw. * src/integer/conv/cl_I_to_Q.cc: Throw. * src/integer/conv/cl_I_to_UL.cc: Throw. * src/integer/conv/cl_I_to_UQ.cc: Throw. * src/integer/conv/cl_I_to_digits.cc: Throw. * src/integer/elem/cl_I_div.cc: Throw. * src/integer/algebraic/cl_I_sqrt.cc: Throw. * src/integer/input/cl_I_read.cc: Throw. * src/integer/input/cl_I_read_stream.cc: Throw. * src/integer/misc/cl_I_as.cc: Throw. * src/rational/elem/cl_RA_from_I_I_div.cc: Throw. * src/rational/elem/cl_RA_recip.cc: Throw. * src/rational/input/cl_RA_read.cc: Throw. * src/rational/input/cl_RA_read_stream.cc: Throw. * src/rational/input/cl_RA_readparsed.cc: Throw. * src/rational/misc/cl_RA_as.cc: Throw. * src/float/base/cl_F_err_nan.cc: Removed. * src/float/base/cl_F_nan_exception.cc: New file. * src/float/base/cl_F_err_un.cc: Removed. * src/float/base/cl_F_underflow_exception.cc: New file. * src/float/base/cl_F_err_ov.cc: Removed. * src/float/base/cl_F_overflow_exception.cc: New file. * src/float/cl_F.h (cl_error_floating_point_nan, cl_error_floating_point_overflow, cl_error_floating_point_underflow): Removed (see cln/float.h.) * src/float/sfloat/cl_SF.h: Throw. * src/float/sfloat/elem/cl_SF_div.cc: Throw. * src/float/sfloat/elem/cl_SF_from_RA.cc: Throw. * src/float/sfloat/elem/cl_SF_scale.cc: Throw. * src/float/sfloat/elem/cl_SF_scale_I.cc: Throw. * src/float/sfloat/misc/cl_SF_as.cc: Throw. * src/float/ffloat/cl_FF.h: Throw. * src/float/ffloat/conv/cl_FF_from_float.cc: Throw. * src/float/ffloat/elem/cl_FF_div.cc: Throw. * src/float/ffloat/elem/cl_FF_from_RA.cc: Throw. * src/float/ffloat/elem/cl_FF_scale.cc: Throw. * src/float/ffloat/elem/cl_FF_scale_I.cc: Throw. * src/float/ffloat/misc/cl_FF_as.cc: Throw. * src/float/dfloat/cl_DF.h: Throw. * src/float/dfloat/conv/cl_DF_from_double.cc: Throw. * src/float/dfloat/elem/cl_DF_div.cc: Throw. * src/float/dfloat/elem/cl_DF_from_RA.cc: Throw. * src/float/dfloat/elem/cl_DF_scale.cc: Throw. * src/float/dfloat/elem/cl_DF_scale_I.cc: Throw. * src/float/dfloat/misc/cl_DF_as.cc: Throw. * src/float/lfloat/algebraic/cl_LF_sqrt.cc: Throw. * src/float/lfloat/elem/cl_LF_1plus.cc: Throw. * src/float/lfloat/elem/cl_LF_I_div.cc: Throw. * src/float/lfloat/elem/cl_LF_I_mul.cc: Throw. * src/float/lfloat/elem/cl_LF_div.cc: Throw. * src/float/lfloat/elem/cl_LF_from_I.cc: Throw. * src/float/lfloat/elem/cl_LF_mul.cc: Throw. * src/float/lfloat/elem/cl_LF_scale.cc: Throw. * src/float/lfloat/elem/cl_LF_scale_I.cc: Throw. * src/float/lfloat/elem/cl_LF_square.cc: Throw. * src/float/lfloat/misc/cl_LF_as.cc: Throw. * src/float/lfloat/misc/cl_LF_shorten.cc: Throw. * src/float/lfloat/misc/cl_LF_shortenrel.cc: Throw. * src/float/lfloat/misc/cl_LF_shortenwith.cc: Throw. * src/float/input/cl_F_read.cc: Throw. * src/float/input/cl_F_read_stream.cc: Throw. * src/float/misc/cl_F_as.cc: Throw. * src/float/misc/cl_F_shortenrel.cc: Throw. * src/float/transcendental/cl_LF_coshsinh_aux.cc: Throw. * src/float/transcendental/cl_LF_cossin_aux.cc: Throw. * src/float/transcendental/cl_LF_exp_aux.cc: Throw. * src/float/transcendental/cl_LF_ratseries_a.cc: Throw. * src/float/transcendental/cl_LF_ratseries_ab.cc: Throw. * src/float/transcendental/cl_LF_ratseries_b.cc: Throw. * src/float/transcendental/cl_LF_ratseries_p.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pa.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pab.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pb.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pq.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pqa.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pqab.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pqb.cc: Throw. * src/float/transcendental/cl_LF_ratseries_q.cc: Throw. * src/float/transcendental/cl_LF_ratseries_qa.cc: Throw. * src/float/transcendental/cl_LF_ratseries_qab.cc: Throw. * src/float/transcendental/cl_LF_ratseries_qb.cc: Throw. * src/float/transcendental/cl_LF_ratseries_stream_pq.cc: Throw. * src/float/transcendental/cl_LF_ratseries_stream_pqa.cc: Throw. * src/float/transcendental/cl_LF_ratseries_stream_pqab.cc: Throw. * src/float/transcendental/cl_LF_ratseries_stream_pqb.cc: Throw. * src/float/transcendental/cl_LF_ratsumseries_pqcd_aux.cc: Throw. * src/float/transcendental/cl_LF_ratsumseries_pqd.cc: Throw. * src/float/transcendental/cl_LF_ratsumseries_pqd_aux.cc: Throw. * src/float/transcendental/cl_LF_zeta_int.cc: Throw. * src/real/elem/cl_R_div.cc: Throw. * src/real/format-output/cl_fmt_cardinal.cc: Throw. * src/real/format-output/cl_fmt_newroman.cc: Throw. * src/real/format-output/cl_fmt_oldroman.cc: Throw. * src/real/input/cl_R_read.cc: Throw. * src/real/input/cl_R_read_stream.cc: Throw. * src/real/misc/cl_R_as.cc: Throw. * src/real/random/cl_R_random.cc: Throw. * src/real/transcendental/cl_R_atan2.cc: Throw. * src/real/transcendental/cl_R_log.cc: Throw. * src/complex/input/cl_N_read.cc: Throw. * src/complex/input/cl_N_read_stream.cc: Throw. * src/complex/misc/cl_N_as.cc: Throw. * src/complex/transcendental/cl_C_atanh_aux.cc: Throw. * src/complex/transcendental/cl_C_expt_C.cc: Throw. * src/complex/transcendental/cl_C_log.cc: Throw. * src/complex/transcendental/cl_C_log2.cc: Throw. * src/numtheory/cl_nt_cornacchia1.cc: Throw. * src/numtheory/cl_nt_cornacchia4.cc: Throw. * src/numtheory/cl_nt_isprobprime.cc: Throw. * src/numtheory/cl_nt_jacobi.cc: Throw. * src/numtheory/cl_nt_jacobi_low.cc: Throw. * src/numtheory/cl_nt_sqrtmodp.cc: Throw. * src/modinteger/cl_MI.cc: Throw. * src/modinteger/cl_MI_int.h: Throw. * src/modinteger/cl_MI_montgom.h: Throw. * src/modinteger/cl_MI_pow2.h: Throw. * src/modinteger/cl_MI_rshift.cc: Throw. * src/modinteger/cl_MI_std.h: Throw. * src/polynomial/elem/cl_UP_GF2.h: Throw. * src/polynomial/elem/cl_UP_MI.h: Throw. * src/polynomial/elem/cl_UP_gen.h: Throw. * src/polynomial/elem/cl_UP_named.cc: Throw. * src/polynomial/elem/cl_UP_no_ring.cc (uninitialized_error, uninitialized_ring): Removed (see cln/ring.h.) * src/polynomial/elem/cl_UP_number.h: Throw. * src/polynomial/elem/cl_UP_unnamed.cc: Throw. * src/vector/cl_GV_I.cc: Throw. * src/vector/cl_GV_number.cc: Throw. * tests/timediv2adic-compare.cc: Use default abort(), not cl_abort(). * tests/timeprint-compare.cc: Likewise. * tests/timerecip2adic-compare.cc: Likewise. * doc/cln.tex: Document the exception classes. * examples/contfrac.cc: Use try/catch instead of setjmp/longjmp. * INSTALL: undocument -fno-exceptions. * README: Add exceptions to list of used C++ features.
18 years ago
Extend the exponent range from 32 bits to 64 bits on selected platforms. * include/cln/number.h: Add signatures for operations with long long. * include/cln/complex_class.h: Likewise. * include/cln/real_class.h: Likewise. * include/cln/real.h: Likewise. * include/cln/rational_class.h: Likewise. * include/cln/rational.h: Likewise. * include/cln/integer_class.h: Likewise. * include/cln/integer.h: Likewise. * include/cln/float.h: Likewise. * include/cln/lfloat.h: Likewise. * include/cln/types.h (sintE and uintE): New types for exponents. * include/cln/*float.h: Use the new types for exponents. * include/cln/floatformat.h (float_format_t): Make underlying type compatible with sintE. * doc/cln.tex: Document changed float_exponent return value. * src/float/cl_F.h: Likewise. * src/float/ffloat/misc/cl_FF_exponent.cc: Likewise. * src/float/input/cl_F_read.cc: Likewise. * src/float/lfloat/cl_LF.h: Likewise. * src/float/lfloat/cl_LF_impl.h: Likewise. * src/float/lfloat/algebraic/cl_LF_sqrt.cc: Likewise. * src/float/lfloat/elem/cl_LF_1plus.cc: Likewise. * src/float/lfloat/elem/cl_LF_I_div.cc: Likewise. * src/float/lfloat/elem/cl_LF_I_mul.cc: Likewise. * src/float/lfloat/elem/cl_LF_compare.cc: Likewise. * src/float/lfloat/elem/cl_LF_div.cc: Likewise. * src/float/lfloat/elem/cl_LF_from_I.cc: Likewise. * src/float/lfloat/elem/cl_LF_fround.cc: Likewise. * src/float/lfloat/elem/cl_LF_ftrunc.cc: Likewise. * src/float/lfloat/elem/cl_LF_futrunc.cc: Likewise. * src/float/lfloat/elem/cl_LF_mul.cc: Likewise. * src/float/lfloat/elem/cl_LF_scale.cc: Likewise. * src/float/lfloat/elem/cl_LF_scale_I.cc: Likewise. * src/float/lfloat/elem/cl_LF_square.cc: Likewise. * src/float/lfloat/elem/cl_LF_to_I.cc: Likewise. * src/float/lfloat/misc/cl_LF_decode.cc: Likewise. * src/float/lfloat/misc/cl_LF_exponent.cc: Likewise. * src/float/lfloat/misc/cl_LF_idecode.cc: Likewise. * src/float/lfloat/misc/cl_LF_shortenrel.cc: Likewise. * src/float/lfloat/misc/cl_LF_shortenwith.cc: Likewise. * src/float/misc/cl_F_decode.cc: Likewise. * src/float/misc/cl_F_exponent.cc: Likewise. * src/float/misc/cl_F_shortenrel.cc: Likewise. * src/float/misc/cl_float_format.cc: Likewise. * src/float/output/cl_F_dprint.cc: Likewise. * src/float/sfloat/misc/cl_SF_exponent.cc: Likewise. * src/float/transcendental/cl_F_atanhx.cc: Likewise. * src/float/transcendental/cl_F_atanx.cc: Likewise. * src/float/transcendental/cl_F_cosh.cc: Likewise. * src/float/transcendental/cl_F_expx.cc: Likewise. * src/float/transcendental/cl_F_lnx.cc: Likewise. * src/float/transcendental/cl_F_sinhx.cc: Likewise. * src/float/transcendental/cl_F_sinx.cc: Likewise. * src/float/transcendental/cl_LF_pi.cc: Likewise. * src/integer/cl_I.h: Likewise. * src/complex/algebraic/cl_LF_hypot.cc: Likewise. * src/complex/elem/division/cl_C_LF_recip.cc: Likewise. * src/float/dfloat/misc/cl_DF_exponent.cc: Likewise. * src/integer/conv/cl_I_from_Q2.cc: Added. * src/base/cl_low.h (isqrtC): New function, for 64 bit falls back to... * src/base/low/cl_low_isqrt.cc (isqrt): ...this new implementation. * src/base/cl_macros.h (bitc): Make sure 64 bit is used if required by exponent operations. * examples/pi.cc: Support more than 646456614 decimal digits.
18 years ago
Remove exception hooks in favor of real C++ exceptions: * include/cln/exception.h: New file... * include/cln/cln.h: ...included here... * include/cln/rational.h: ...and here. * include/cln/abort.h: Removed. * include/cln/integer.h (ash_exception, exquo_exception): Add exception types. * include/cln/float.h (floating_point_exception, floating_point_nan_exception, floating_point_overflow_exception, floating_point_underflow_exception): Likewise. * include/cln/number_io.h (read_number_exception, read_number_bad_syntax_exception, read_number_junk_exception, read_number_eof_exception): Likewise. * include/cln/object.h (CL_DEFINE_CONVERTER): Replace dynamic assertion with a compile-time assertion. * include/cln/GV.h: Replace nonreturning functions with exceptions. * include/cln/modinteger.h: Likewise. * include/cln/SV.h: Likewise. * include/cln/ring.h: Likewise. * include/cln/string.h: Likewise. * include/cln/univpoly.h: Likewise. * src/base/cl_abort.cc: Removed. * src/base/cl_N_err_d0.cc: Removed. * src/base/cl_d0_exception.cc: New file. * src/base/cl_as_err.cc: Removed. * src/base/cl_as_exception.cc: New file. * src/base/cl_notreached.cc: Removed. * src/base/cl_notreached_exception.cc: New file. * src/base/input/cl_read_err_bad.cc: Removed. * src/base/input/cl_read_bad_syntax_exception.cc: New file. * src/base/input/cl_read_err_junk.cc: Removed. * src/base/input/cl_read_junk_exception.cc: New file. * src/base/input/cl_read_err_eof.cc: Removed. * src/base/input/cl_read_eof_exception.cc: New file. * src/base/cl_N.h (cl_as_error): Removed (see cln/exception.h). * src/base/macros.h (NOTREACHED): Throw. * src/base/cl_malloc.cc (xmalloc): Throw. * src/base/digitseq/cl_2DS_div.cc: Throw. * src/base/digitseq/cl_DS_div.cc: Throw. * src/base/digitseq/cl_DS_mul.cc: Throw. * src/base/digitseq/cl_DS_mul_fftc.h: Throw. * src/base/digitseq/cl_DS_mul_fftcs.h: Throw. * src/base/digitseq/cl_DS_mul_fftm.h: Throw. * src/base/digitseq/cl_DS_mul_fftp.h: Throw. * src/base/digitseq/cl_DS_mul_fftp3.h: Throw. * src/base/digitseq/cl_DS_mul_fftp3m.h: Throw. * src/base/digitseq/cl_DS_mul_fftr.h: Throw. * src/base/digitseq/cl_DS_mul_nuss.h: Throw. * src/base/digitseq/cl_DS_recipsqrt.cc: Throw. * src/base/digitseq/cl_DS_sqrt.cc: Throw. * src/base/hash/cl_hash.h: Throw. * src/base/hash/cl_hash1.h: Throw. * src/base/hash/cl_hash1weak.h: Throw. * src/base/hash/cl_hash2.h: Throw. * src/base/hash/cl_hash2weak.h: Throw. * src/base/hash/cl_hashset.h: Throw. * src/base/hash/cl_hashuniq.h: Throw. * src/base/hash/cl_hashuniqweak.h: Throw. * src/base/proplist/cl_pl_add.cc: Throw. * src/base/ring/cl_no_ring.cc: Throw. * src/base/string/cl_spushstring.h: Throw. * src/base/symbol/cl_symbol.cc: Throw. * src/integer/bitwise/cl_I_ash.h: Removed (see cln/integer.h.) * src/integer/bitwise/cl_I_asherr.cc: Removed. * src/integer/bitwise/cl_I_ash_exception.cc: New file. * src/integer/division/cl_I_exquoerr.cc: Removed. * src/integer/division/cl_I_exquo_exception.cc: New file. * src/integer/cl_I.h: Throw. * src/integer/division/cl_I_exquopos.cc: Throw. * src/integer/bitwise/cl_I_logbitp_I.cc: Throw. * src/integer/bitwise/cl_I_ash.cc: Throw. * src/integer/bitwise/cl_I_ash_I.cc: Throw. * src/integer/division/cl_I_exquo.cc: Throw. * src/integer/gcd/cl_I_gcd_aux2.cc: Throw. * src/integer/conv/cl_I_to_L.cc: Throw. * src/integer/conv/cl_I_to_Q.cc: Throw. * src/integer/conv/cl_I_to_UL.cc: Throw. * src/integer/conv/cl_I_to_UQ.cc: Throw. * src/integer/conv/cl_I_to_digits.cc: Throw. * src/integer/elem/cl_I_div.cc: Throw. * src/integer/algebraic/cl_I_sqrt.cc: Throw. * src/integer/input/cl_I_read.cc: Throw. * src/integer/input/cl_I_read_stream.cc: Throw. * src/integer/misc/cl_I_as.cc: Throw. * src/rational/elem/cl_RA_from_I_I_div.cc: Throw. * src/rational/elem/cl_RA_recip.cc: Throw. * src/rational/input/cl_RA_read.cc: Throw. * src/rational/input/cl_RA_read_stream.cc: Throw. * src/rational/input/cl_RA_readparsed.cc: Throw. * src/rational/misc/cl_RA_as.cc: Throw. * src/float/base/cl_F_err_nan.cc: Removed. * src/float/base/cl_F_nan_exception.cc: New file. * src/float/base/cl_F_err_un.cc: Removed. * src/float/base/cl_F_underflow_exception.cc: New file. * src/float/base/cl_F_err_ov.cc: Removed. * src/float/base/cl_F_overflow_exception.cc: New file. * src/float/cl_F.h (cl_error_floating_point_nan, cl_error_floating_point_overflow, cl_error_floating_point_underflow): Removed (see cln/float.h.) * src/float/sfloat/cl_SF.h: Throw. * src/float/sfloat/elem/cl_SF_div.cc: Throw. * src/float/sfloat/elem/cl_SF_from_RA.cc: Throw. * src/float/sfloat/elem/cl_SF_scale.cc: Throw. * src/float/sfloat/elem/cl_SF_scale_I.cc: Throw. * src/float/sfloat/misc/cl_SF_as.cc: Throw. * src/float/ffloat/cl_FF.h: Throw. * src/float/ffloat/conv/cl_FF_from_float.cc: Throw. * src/float/ffloat/elem/cl_FF_div.cc: Throw. * src/float/ffloat/elem/cl_FF_from_RA.cc: Throw. * src/float/ffloat/elem/cl_FF_scale.cc: Throw. * src/float/ffloat/elem/cl_FF_scale_I.cc: Throw. * src/float/ffloat/misc/cl_FF_as.cc: Throw. * src/float/dfloat/cl_DF.h: Throw. * src/float/dfloat/conv/cl_DF_from_double.cc: Throw. * src/float/dfloat/elem/cl_DF_div.cc: Throw. * src/float/dfloat/elem/cl_DF_from_RA.cc: Throw. * src/float/dfloat/elem/cl_DF_scale.cc: Throw. * src/float/dfloat/elem/cl_DF_scale_I.cc: Throw. * src/float/dfloat/misc/cl_DF_as.cc: Throw. * src/float/lfloat/algebraic/cl_LF_sqrt.cc: Throw. * src/float/lfloat/elem/cl_LF_1plus.cc: Throw. * src/float/lfloat/elem/cl_LF_I_div.cc: Throw. * src/float/lfloat/elem/cl_LF_I_mul.cc: Throw. * src/float/lfloat/elem/cl_LF_div.cc: Throw. * src/float/lfloat/elem/cl_LF_from_I.cc: Throw. * src/float/lfloat/elem/cl_LF_mul.cc: Throw. * src/float/lfloat/elem/cl_LF_scale.cc: Throw. * src/float/lfloat/elem/cl_LF_scale_I.cc: Throw. * src/float/lfloat/elem/cl_LF_square.cc: Throw. * src/float/lfloat/misc/cl_LF_as.cc: Throw. * src/float/lfloat/misc/cl_LF_shorten.cc: Throw. * src/float/lfloat/misc/cl_LF_shortenrel.cc: Throw. * src/float/lfloat/misc/cl_LF_shortenwith.cc: Throw. * src/float/input/cl_F_read.cc: Throw. * src/float/input/cl_F_read_stream.cc: Throw. * src/float/misc/cl_F_as.cc: Throw. * src/float/misc/cl_F_shortenrel.cc: Throw. * src/float/transcendental/cl_LF_coshsinh_aux.cc: Throw. * src/float/transcendental/cl_LF_cossin_aux.cc: Throw. * src/float/transcendental/cl_LF_exp_aux.cc: Throw. * src/float/transcendental/cl_LF_ratseries_a.cc: Throw. * src/float/transcendental/cl_LF_ratseries_ab.cc: Throw. * src/float/transcendental/cl_LF_ratseries_b.cc: Throw. * src/float/transcendental/cl_LF_ratseries_p.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pa.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pab.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pb.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pq.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pqa.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pqab.cc: Throw. * src/float/transcendental/cl_LF_ratseries_pqb.cc: Throw. * src/float/transcendental/cl_LF_ratseries_q.cc: Throw. * src/float/transcendental/cl_LF_ratseries_qa.cc: Throw. * src/float/transcendental/cl_LF_ratseries_qab.cc: Throw. * src/float/transcendental/cl_LF_ratseries_qb.cc: Throw. * src/float/transcendental/cl_LF_ratseries_stream_pq.cc: Throw. * src/float/transcendental/cl_LF_ratseries_stream_pqa.cc: Throw. * src/float/transcendental/cl_LF_ratseries_stream_pqab.cc: Throw. * src/float/transcendental/cl_LF_ratseries_stream_pqb.cc: Throw. * src/float/transcendental/cl_LF_ratsumseries_pqcd_aux.cc: Throw. * src/float/transcendental/cl_LF_ratsumseries_pqd.cc: Throw. * src/float/transcendental/cl_LF_ratsumseries_pqd_aux.cc: Throw. * src/float/transcendental/cl_LF_zeta_int.cc: Throw. * src/real/elem/cl_R_div.cc: Throw. * src/real/format-output/cl_fmt_cardinal.cc: Throw. * src/real/format-output/cl_fmt_newroman.cc: Throw. * src/real/format-output/cl_fmt_oldroman.cc: Throw. * src/real/input/cl_R_read.cc: Throw. * src/real/input/cl_R_read_stream.cc: Throw. * src/real/misc/cl_R_as.cc: Throw. * src/real/random/cl_R_random.cc: Throw. * src/real/transcendental/cl_R_atan2.cc: Throw. * src/real/transcendental/cl_R_log.cc: Throw. * src/complex/input/cl_N_read.cc: Throw. * src/complex/input/cl_N_read_stream.cc: Throw. * src/complex/misc/cl_N_as.cc: Throw. * src/complex/transcendental/cl_C_atanh_aux.cc: Throw. * src/complex/transcendental/cl_C_expt_C.cc: Throw. * src/complex/transcendental/cl_C_log.cc: Throw. * src/complex/transcendental/cl_C_log2.cc: Throw. * src/numtheory/cl_nt_cornacchia1.cc: Throw. * src/numtheory/cl_nt_cornacchia4.cc: Throw. * src/numtheory/cl_nt_isprobprime.cc: Throw. * src/numtheory/cl_nt_jacobi.cc: Throw. * src/numtheory/cl_nt_jacobi_low.cc: Throw. * src/numtheory/cl_nt_sqrtmodp.cc: Throw. * src/modinteger/cl_MI.cc: Throw. * src/modinteger/cl_MI_int.h: Throw. * src/modinteger/cl_MI_montgom.h: Throw. * src/modinteger/cl_MI_pow2.h: Throw. * src/modinteger/cl_MI_rshift.cc: Throw. * src/modinteger/cl_MI_std.h: Throw. * src/polynomial/elem/cl_UP_GF2.h: Throw. * src/polynomial/elem/cl_UP_MI.h: Throw. * src/polynomial/elem/cl_UP_gen.h: Throw. * src/polynomial/elem/cl_UP_named.cc: Throw. * src/polynomial/elem/cl_UP_no_ring.cc (uninitialized_error, uninitialized_ring): Removed (see cln/ring.h.) * src/polynomial/elem/cl_UP_number.h: Throw. * src/polynomial/elem/cl_UP_unnamed.cc: Throw. * src/vector/cl_GV_I.cc: Throw. * src/vector/cl_GV_number.cc: Throw. * tests/timediv2adic-compare.cc: Use default abort(), not cl_abort(). * tests/timeprint-compare.cc: Likewise. * tests/timerecip2adic-compare.cc: Likewise. * doc/cln.tex: Document the exception classes. * examples/contfrac.cc: Use try/catch instead of setjmp/longjmp. * INSTALL: undocument -fno-exceptions. * README: Add exceptions to list of used C++ features.
18 years ago
Extend the exponent range from 32 bits to 64 bits on selected platforms. * include/cln/number.h: Add signatures for operations with long long. * include/cln/complex_class.h: Likewise. * include/cln/real_class.h: Likewise. * include/cln/real.h: Likewise. * include/cln/rational_class.h: Likewise. * include/cln/rational.h: Likewise. * include/cln/integer_class.h: Likewise. * include/cln/integer.h: Likewise. * include/cln/float.h: Likewise. * include/cln/lfloat.h: Likewise. * include/cln/types.h (sintE and uintE): New types for exponents. * include/cln/*float.h: Use the new types for exponents. * include/cln/floatformat.h (float_format_t): Make underlying type compatible with sintE. * doc/cln.tex: Document changed float_exponent return value. * src/float/cl_F.h: Likewise. * src/float/ffloat/misc/cl_FF_exponent.cc: Likewise. * src/float/input/cl_F_read.cc: Likewise. * src/float/lfloat/cl_LF.h: Likewise. * src/float/lfloat/cl_LF_impl.h: Likewise. * src/float/lfloat/algebraic/cl_LF_sqrt.cc: Likewise. * src/float/lfloat/elem/cl_LF_1plus.cc: Likewise. * src/float/lfloat/elem/cl_LF_I_div.cc: Likewise. * src/float/lfloat/elem/cl_LF_I_mul.cc: Likewise. * src/float/lfloat/elem/cl_LF_compare.cc: Likewise. * src/float/lfloat/elem/cl_LF_div.cc: Likewise. * src/float/lfloat/elem/cl_LF_from_I.cc: Likewise. * src/float/lfloat/elem/cl_LF_fround.cc: Likewise. * src/float/lfloat/elem/cl_LF_ftrunc.cc: Likewise. * src/float/lfloat/elem/cl_LF_futrunc.cc: Likewise. * src/float/lfloat/elem/cl_LF_mul.cc: Likewise. * src/float/lfloat/elem/cl_LF_scale.cc: Likewise. * src/float/lfloat/elem/cl_LF_scale_I.cc: Likewise. * src/float/lfloat/elem/cl_LF_square.cc: Likewise. * src/float/lfloat/elem/cl_LF_to_I.cc: Likewise. * src/float/lfloat/misc/cl_LF_decode.cc: Likewise. * src/float/lfloat/misc/cl_LF_exponent.cc: Likewise. * src/float/lfloat/misc/cl_LF_idecode.cc: Likewise. * src/float/lfloat/misc/cl_LF_shortenrel.cc: Likewise. * src/float/lfloat/misc/cl_LF_shortenwith.cc: Likewise. * src/float/misc/cl_F_decode.cc: Likewise. * src/float/misc/cl_F_exponent.cc: Likewise. * src/float/misc/cl_F_shortenrel.cc: Likewise. * src/float/misc/cl_float_format.cc: Likewise. * src/float/output/cl_F_dprint.cc: Likewise. * src/float/sfloat/misc/cl_SF_exponent.cc: Likewise. * src/float/transcendental/cl_F_atanhx.cc: Likewise. * src/float/transcendental/cl_F_atanx.cc: Likewise. * src/float/transcendental/cl_F_cosh.cc: Likewise. * src/float/transcendental/cl_F_expx.cc: Likewise. * src/float/transcendental/cl_F_lnx.cc: Likewise. * src/float/transcendental/cl_F_sinhx.cc: Likewise. * src/float/transcendental/cl_F_sinx.cc: Likewise. * src/float/transcendental/cl_LF_pi.cc: Likewise. * src/integer/cl_I.h: Likewise. * src/complex/algebraic/cl_LF_hypot.cc: Likewise. * src/complex/elem/division/cl_C_LF_recip.cc: Likewise. * src/float/dfloat/misc/cl_DF_exponent.cc: Likewise. * src/integer/conv/cl_I_from_Q2.cc: Added. * src/base/cl_low.h (isqrtC): New function, for 64 bit falls back to... * src/base/low/cl_low_isqrt.cc (isqrt): ...this new implementation. * src/base/cl_macros.h (bitc): Make sure 64 bit is used if required by exponent operations. * examples/pi.cc: Support more than 646456614 decimal digits.
18 years ago
19 years ago
24 years ago
24 years ago
24 years ago
24 years ago
24 years ago
25 years ago
  1. 2007-mm-dd, version 1.2.0
  2. =========================
  3. Implementation changes
  4. ----------------------
  5. * Added support for huge numbers...
  6. * CLN now uses C++ exceptions for error handling.
  7. * Removed cl_boolean. Use built-in C++ bool instead.
  8. * Dropped the cln-config script. Please use pkg-config instead.
  9. 2006-08-08, version 1.1.13
  10. ==========================
  11. * Compilation fixes for 64-bit brokenness introduced in last release.
  12. 2006-08-06, version 1.1.12
  13. ==========================
  14. Implementation changes
  15. ----------------------
  16. * Fix rare assertion when printing quite large floats.
  17. Other changes
  18. -------------
  19. * Compilation fixes for several platforms: *BSD, Intel Mac, and MinGW.
  20. 2005-11-23, version 1.1.11
  21. ==========================
  22. Algorithmic changes
  23. -------------------
  24. * Considerably improved performance of number input.
  25. 2005-10-22, version 1.1.10
  26. ==========================
  27. Implementation changes
  28. ----------------------
  29. * Removed the vptr of modular integer and univariate polynomial
  30. classes in order to fix some crashes in that sector of CLN. Code
  31. using this needs to be recompiled.
  32. * Many more operator signatures make the integration of CLN types with
  33. float/double easier.
  34. Other changes
  35. -------------
  36. * Several compilation fixes.
  37. * Made it possible to cross-compile CLN.
  38. 2004-11-03, version 1.1.9
  39. =========================
  40. Algorithmic changes
  41. -------------------
  42. * Input of numbers in bases 2, 4, 8, 16 and 32 is now done in linear bit
  43. complexity as opposed to O(N^2). Useful for all kinds of persistency.
  44. Implementation changes
  45. ----------------------
  46. * Fixed several bugs in the integer input and output routines that could be
  47. blamed for random crashes in the following cases: output in base 32 for
  48. quite large numbers, input in base 2 for fixnums and input in base 3 for
  49. fixnums on a 64 bit target.
  50. * Fixed crash when radix specifiers were used in input streams.
  51. * Speed up on x86_64 and ia64 by adding some inline assembly.
  52. Other changes
  53. -------------
  54. * Fixes for compilation on MacOS X and little endian Mips.
  55. 2004-07-01, version 1.1.8
  56. =========================
  57. Implementation changes
  58. ----------------------
  59. * Fix bug in log(a,b) when b is 1/n.
  60. * Fix crash in shared library on AMD64.
  61. Other changes
  62. -------------
  63. * Fixes for some language conformance issues.
  64. 2004-05-02, version 1.1.7
  65. =========================
  66. Implementation changes
  67. ----------------------
  68. * Fixed two corner case bugs in cl_LF multiplication and sqrt.
  69. * Workaround GCC 3.[123] compiler bug.
  70. * Parallel build should work reliably now.
  71. 2004-01-01, version 1.1.6
  72. =========================
  73. Functional changes
  74. ------------------
  75. * New function
  76. `sintL ldegree(const cl_UP& x)'.
  77. (Suggested by Munagala Ramanath.)
  78. Implementation changes
  79. ----------------------
  80. * Fixed a bug in polynomial subtraction.
  81. (Reported by Munagala Ramanath.)
  82. Other changes
  83. -------------
  84. * Provide a package control file cln.pc.
  85. * Updated infrastructure to autoconf-2.57.
  86. * Fixed compilation issues with GCC-3.4 prereleases and on Mac OSX.
  87. 2002-05-28, version 1.1.5
  88. =========================
  89. Other changes
  90. -------------
  91. * Fix compilation issues with GCC-3.1.
  92. 2002-01-04, version 1.1.4
  93. =========================
  94. Other changes
  95. -------------
  96. * Fix compilation issues with GMP-4.0 and on IBM s390.
  97. * Updated infrastructure to autoconf-2.52.
  98. 2001-11-05, version 1.1.3
  99. =========================
  100. Implementation changes
  101. ----------------------
  102. * Fix a computation error of sin and sinh for small arguments.
  103. (Reported by Christopher Kennedy.)
  104. * Module ordering works with gcc-3.0 -fuse-cxa-atexit now.
  105. * -DNO_ASM works properly on Sparc.
  106. 2001-07-25, version 1.1.2
  107. Implementation changes
  108. ----------------------
  109. * Minor cleanups and portability issues.
  110. (No need to upgrade if you have been happy so far.)
  111. 2001-05-31, version 1.1.1
  112. =========================
  113. Implementation changes
  114. ----------------------
  115. * Minor cleanups for GCC 3.0 and libstdc++-v3.
  116. * Fixes several little bugs and portability issues.
  117. 2000-12-14, version 1.1.0
  118. =========================
  119. Functional changes
  120. ------------------
  121. * ISO/IEC 14882 fashion adjustments:
  122. Put everything into namespace cln. All fundamental data types still
  123. carry their old names. Other non-macro identifiers are now written as
  124. cln::foo instead of cl_foo, except where the cl_ comes from a data type
  125. (as in cl_I_to_int()). Headers are installed into a separate
  126. directory, #include <cln/foo.h> instead of <cl_foo.h>. Applications
  127. must be manually ported to the new style. We apologize for the
  128. inconvenience.
  129. Implementation changes
  130. ----------------------
  131. * Removed internal copy of mpn subdirectory tree from libgmp-2.0.2.
  132. Configuring with --with-gmp now relies on an installed libgmp library
  133. version 3 or above. We recommend libgmp-3.1 or newer.
  134. * Adjusted several break-even points to make better use of a faster libgmp and
  135. better match present-day CPUs.
  136. * Fix several errors for architectures with 64-bit wide words.
  137. (Others than Alpha, which worked already.)
  138. * Fix compilation errors with current GCC-3.0 snapshots.
  139. Other changes
  140. -------------
  141. * Added package tools: script cln-config and autoconf macro AC_PATH_CLN
  142. (in file cln.m4).
  143. 2000-01-13, version 1.0.3
  144. =========================
  145. Functional changes
  146. ------------------
  147. * New function
  148. `cl_I doublefactorial (uintL n)'.
  149. Implementation changes
  150. ----------------------
  151. * Fix several little configuration errors.
  152. * Fix some compilation errors with gcc-2.96 prereleases.
  153. 1999-09-07, version 1.0.2
  154. =========================
  155. Functional changes
  156. ------------------
  157. * New functions
  158. `cl_I numerator (const cl_RA&)'
  159. `cl_I denominator (const cl_RA&)'.
  160. Suggested by Richard Kreckel and Sylvain Pion.
  161. * New function `cl_equal_hashcode' for the classes cl_N, cl_R, cl_RA, cl_I,
  162. cl_F, cl_SF, cl_FF, cl_DF, cl_LF.
  163. Suggested by Richard Kreckel.
  164. Implementation changes
  165. ----------------------
  166. * Fix an endless loop when either of the functions
  167. `cl_boolean rootp (const cl_RA&, uintL, cl_RA*)',
  168. `cl_boolean rootp (const cl_RA&, const cl_I&, cl_RA*)'
  169. was called with an integer argument. Thanks to Richard Kreckel.
  170. * Fix a bug in the addition and subtraction of rational numbers which could
  171. cause unnormalized rational numbers like 3/1 to be created.
  172. 1999-06-09, version 1.0.1
  173. =========================
  174. Algorithmic changes
  175. -------------------
  176. * Speed up the functions `logand', `lognand', `logandc2'.
  177. Implementation changes
  178. ----------------------
  179. * Fix compilation errors with gcc-2.7.2, egcs-1.1.2 and gcc-2.95.
  180. * Fix compilation errors on HPPA, MIPS, some versions of DEC Alpha, OpenBSD,
  181. and SPARC with LiDIA.
  182. * Fix compilation errors on Win32. Thanks to Mumit Khan.
  183. 1999-01-12, version 1.0
  184. =======================
  185. Functional changes
  186. ------------------
  187. * New include files, containing I/O functions:
  188. <cl_integer_io.h>
  189. <cl_rational_io.h>
  190. <cl_float_io.h>
  191. <cl_sfloat_io.h> <cl_ffloat_io.h> <cl_dfloat_io.h> <cl_lfloat_io.h>
  192. <cl_real_io.h>
  193. <cl_complex_io.h>
  194. The file <cl_io.h> now contains only I/O functions for characters and
  195. C integers.
  196. * To access the base ring of a univariate polynomial ring, now use
  197. `R->basering()' instead of `R->basering'.
  198. * Implement `plusp', `max', `min' for the classes cl_F, cl_SF, cl_FF, cl_DF,
  199. cl_LF, cl_RA, cl_I.
  200. * Implement `abs' for the class cl_RA.
  201. * Replaced `read_number' with specialized functions `read_complex',
  202. `read_real', `read_float', `read_rational', `read_integer'.
  203. * Replaced the functions
  204. `void fprint (cl_ostream stream, int x)'
  205. `void fprint (cl_ostream stream, unsigned int x)'
  206. `void fprint (cl_ostream stream, long x)'
  207. `void fprint (cl_ostream stream, unsigned long x)'
  208. with
  209. `void fprintdecimal (cl_ostream stream, int x)'
  210. `void fprintdecimal (cl_ostream stream, unsigned int x)'
  211. `void fprintdecimal (cl_ostream stream, long x)'
  212. `void fprintdecimal (cl_ostream stream, unsigned long x)'
  213. Algorithmic changes
  214. -------------------
  215. * The function `cl_I xgcd (const cl_I& a, const cl_I& b, cl_I* u, cl_I* v)'
  216. now normalizes the resulting u and v to be of minimum absolute value.
  217. (Suggested by the LiDIA group.)
  218. * The conversion from string to cl_F, cl_R, cl_N now chooses the format of
  219. the resulting float depending on the number of significant digits present
  220. in the input string.
  221. * Speed up the functions
  222. `cl_R operator* (const cl_RA&, const cl_LF&)'
  223. `cl_R operator* (const cl_LF&, const cl_RA&)'
  224. `cl_R operator/ (const cl_RA&, const cl_LF&)'
  225. `cl_LF operator/ (const cl_LF&, const cl_RA&)'
  226. Implementation changes
  227. ----------------------
  228. * Fix for `cl_I ash (const cl_I& x, const cl_I& n)' when |n| is between
  229. 2^29 and 2^37.
  230. (Reported by Thomas Papanikolaou.)
  231. * Fixed a memory leak in the long-float truncate function.
  232. (Reported by Keith Briggs.)
  233. * Fixed a memory leak in the modular integer exponentiation.
  234. * Reduced the stack space consumption. The default stack size (typically
  235. 8 MB on Unix) should now be sufficient in order to run all CLN programs.
  236. * Portability fixes for m68k platforms.
  237. (Reported and fixed by Roman Hodek.)
  238. Other changes
  239. -------------
  240. * Restructured the sources. The subpackages
  241. - base
  242. - base + integer
  243. - base + integer + rational
  244. - base + integer + rational + float
  245. - base + integer + rational + float + real
  246. - base + integer + rational + float + real + complex
  247. - base + integer + modinteger
  248. are now self-contained.
  249. * The configuration script can be called with option `--enable-shared',
  250. to build CLN as a shared library.
  251. * Updated the LiDIA interface.
  252. 1998-05-07, version 0.98
  253. ========================
  254. Implementation changes
  255. ----------------------
  256. * Portability fixes for 64-bit platforms like DEC Alpha.
  257. (Reported by John Cremona and Thomas Papanikolaou.)
  258. 1998-02-27
  259. ==========
  260. Other changes
  261. -------------
  262. * Portability fixes for new ISO/ANSI C++ compilers.
  263. * Updated the LiDIA interface for LiDIA version 1.3 and fixed a couple of
  264. bugs.
  265. 1997-09-06
  266. ==========
  267. Implementation changes
  268. ----------------------
  269. * Portability fix for i686 platforms.
  270. (Reported by Michael Pfeifer.)
  271. 1997-09-01
  272. ==========
  273. Functional changes
  274. ------------------
  275. * New include files:
  276. <cl_ring.h>, declaring general rings,
  277. <cl_GV.h>, <cl_SV.h>, defining vectors,
  278. <cl_univpoly.h>, defining univariate polynomial rings.
  279. Algorithmic changes
  280. -------------------
  281. * Speed up the internal computation of ln(2) and ln(10) by about 20%.
  282. * Speed up the computation of exp (for 1000 digits or more) by about 10%.
  283. Implementation changes
  284. ----------------------
  285. * Portability fix for 64-bit platforms like DEC Alpha:
  286. Fixed a bug: endless loop in `sqrt' for large precision.
  287. (Reported by Blake Jones.)
  288. * Portability fix for RS/6000 platforms.
  289. 1997-02-21
  290. ==========
  291. Algorithmic changes
  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. Implementation changes
  299. ----------------------
  300. * Portability fix for SunOS 4 and other BSD platforms.
  301. 1997-01-31
  302. ==========
  303. Algorithmic changes
  304. -------------------
  305. * Fixed a bug: endless recursion in `cl_boolean sqrtp (const cl_RA&, cl_RA*)'.
  306. * Fixed a bug: Only the first 85% of the digits returned by `cl_eulerconst'
  307. for moderate precision (< 11370 digits) were correct.
  308. (Reported by Thomas Papanikolaou.)
  309. Implementation changes
  310. ----------------------
  311. * Fix static initialization order bug. Requires g++ 2.7.0 or newer
  312. (g++ 2.8.0 or newer on HPPA).
  313. * New method `debug_print()', for pretty printing of CLN objects, intended
  314. to be called from the debugger.
  315. 1997-01-07
  316. ==========
  317. Functional changes
  318. ------------------
  319. * New functions `float cl_float_approx (const cl_R&)',
  320. `double cl_double_approx (const cl_R&)'
  321. for converting a real number to a C `float' or `double'.
  322. Algorithmic changes
  323. -------------------
  324. * Speed up `cl_zeta' using Cohen-Villegas-Zagier convergence acceleration.
  325. Implementation changes
  326. ----------------------
  327. * New function `cl_print', for pretty printing of CLN objects, intended to be
  328. called from the debugger.
  329. * Portability fixes for 64-bit platforms like DEC Alpha.
  330. 1996-11-29
  331. ==========
  332. Functional changes
  333. ------------------
  334. * New include files:
  335. <cl_string.h>, defining strings,
  336. <cl_symbol.h>, defining symbols,
  337. <cl_proplist.h>, defining property lists,
  338. <cl_condition.h>, defining conditions/exceptions.
  339. * New functions `cl_F cl_catalanconst ()',
  340. `cl_F cl_catalanconst (const cl_F&)',
  341. `cl_F cl_catalanconst (cl_float_format_t)'
  342. which return Catalan's constant.
  343. (Suggested by Thomas Papanikolaou.)
  344. * New functions `cl_F cl_zeta (int s)',
  345. `cl_F cl_zeta (int s, const cl_F&)',
  346. `cl_F cl_zeta (int s, cl_float_format_t)'
  347. which return the Riemann zeta function at an integral point s>1.
  348. * New functions `cl_F cl_exp1 ()',
  349. `cl_F cl_exp1 (const cl_F&)',
  350. `cl_F cl_exp1 (cl_float_format_t)'
  351. which return e = exp(1).
  352. * New function `cl_I binomial (uintL n, uintL k)' which returns the
  353. binomial coefficient (n choose k).
  354. * New functions `int cl_I_to_int (const cl_I&)',
  355. `unsigned int cl_I_to_uint (const cl_I&)',
  356. `long cl_I_to_long (const cl_I&)',
  357. `unsigned long cl_I_to_ulong (const cl_I&)'
  358. for converting an integer to a C `int' or `long'.
  359. * New functions `float cl_float_approx (const cl_I&)',
  360. `float cl_float_approx (const cl_RA&)',
  361. `double cl_double_approx (const cl_I&)',
  362. `double cl_double_approx (const cl_RA&)'
  363. for converting a rational number to a C `float' or `double'.
  364. Implementation changes
  365. ----------------------
  366. * Moved the sources to subdirectories.
  367. (Suggested by J�rg Arndt.)
  368. Other changes
  369. -------------
  370. * New benchmark for better comparing LiDIA, Pari and CLN.
  371. * Added an interface to LiDIA, allows using CLN as a replacement of libI.
  372. (Suggested by Thomas Papanikolaou.)
  373. * Added an ILOG Talk interface, for interactive use of the library.
  374. 1996-10-13
  375. ==========
  376. Functional changes
  377. ------------------
  378. * New include file <cl_modinteger.h>, defining modular integer rings.
  379. * New functions `cl_F cl_eulerconst ()', `cl_F cl_eulerconst (const cl_F&)',
  380. `cl_F cl_eulerconst (cl_float_format_t)' which return Euler's constant.
  381. Algorithmic changes
  382. -------------------
  383. * Speed up square root of large numbers by use of Newton's algorithm.
  384. * Speed up multiplication and division of large long-floats by small integers.
  385. * Speed up the computation of pi, e, exp, log, sin, cos, tan, sinh, cosh, tanh,
  386. asin, acos, atan, asinh, acosh, atanh.
  387. All of them are now computed with asymptotically fast algorithms, of
  388. bit complexity O(log(N)^2*M(N)) = O(N^(1+epsilon)) for N digits.
  389. * Fixed several bugs in the transcendental functions routines. In particular,
  390. the `log' function went into an endless loop.
  391. * Fixed a bug: The cl_I -> cl_DF and cl_I -> cl_LF conversions didn't work
  392. correctly for negative integers.
  393. 1996-09-08
  394. ==========
  395. Functional changes
  396. ------------------
  397. * New include file <cl_numtheory.h>, containing number theoretic functions,
  398. for now only the Jacobi symbol.
  399. Algorithmic changes
  400. -------------------
  401. * Speed up squaring of large numbers by 30%.
  402. * Speed up division of large numbers by use of Newton's algorithm. The
  403. improvement is noticeable for numbers with at least about 1000 words =
  404. 32000 bits = 10000 decimal digits.
  405. * Speed up the binary-to-decimal conversion of large numbers by use of
  406. a divide-and-conquer method. The improvement is noticeable for numbers
  407. of at least 1250 words = 40000 bits = 12000 decimal digits.
  408. * Fixed a bug: The functions `cl_boolean minusp (const type&)' and
  409. `cl_boolean zerop (const type&)' were missing for type=cl_I and type=cl_RA.
  410. * Fixed a bug: The function `cl_boolean logtest (const cl_I&, const cl_I&)'
  411. didn't work correctly if both arguments were bignums.
  412. * Fixed a bug: The function `cl_I mod (const cl_I&, const cl_I&)', when called
  413. 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.
  416. Implementation changes
  417. ----------------------
  418. * Speed up CLN by using the low-level addition/subtraction, shift and
  419. multiplication loops from GNU gmp-2.0.2.
  420. * Fixed a severe reference counting bug in the constructors
  421. `type::type (const long)' and `type::type (const unsigned long)'
  422. for type=cl_number, type=cl_N, type=cl_RA, type=cl_I that most often
  423. caused a core dump.
  424. * CLN's word sequences can be chosen to be stored big-endian in memory
  425. (like it was up to now) or little-endian (needed for interfacing to GMP).
  426. 1996-05-20
  427. ==========
  428. Functional changes
  429. ------------------
  430. * New include file <cl_timing.h>, contains timing facilities.
  431. Algorithmic changes
  432. -------------------
  433. * Speed up the multiplication of large numbers by use of an FFT based
  434. algorithm (Sch�nhage-Strassen). The improvement is noticeable when both
  435. factors have at least 1250 words = 40000 bits = 12000 decimal digits.
  436. * Speed up the functions `cl_I gcd (const cl_I&, const cl_I&)' and
  437. `cl_I xgcd (const cl_I&, const cl_I&, cl_I*, cl_I*)' by a factor
  438. of 1.5. Simple trick.
  439. * Speed up the function `cl_boolean sqrtp (const cl_I&, cl_I*)' using a
  440. trick from Henri Cohen.
  441. * Fixed an endless loop in the function `uint32 gcd (uint32, uint32)'
  442. which occurred when one of the arguments is zero.
  443. * Fixed an endless loop in the function `cl_I factorial (uintL)'.
  444. Implementation changes
  445. ----------------------
  446. * <cl_io.h> now uses iostreams by default instead of stdio FILE pointers.
  447. (Reported by Tito Flagella.)
  448. * Fixed linking error when compiling without -O.
  449. (Reported by Tito Flagella.)
  450. * Fixed syntax error in <cl_input.h>.
  451. (Reported by Tito Flagella.)
  452. * Fixed syntax errors in src/cl_DS.h and src/cl_FF_plus.h.
  453. (Reported by Marcus Daniels.)
  454. * More portability fixes.
  455. * Configuration now supports the --srcdir option.
  456. (Reported by Marcus Daniels.)
  457. 1996-01-03
  458. ==========
  459. * Compiles well on Sparc, using gcc.
  460. * Mail to beta testers.
  461. 1995-12-06
  462. ==========
  463. * Compiles well on Linux, using gcc.