crmf 0.1.1 → 0.1.2

data/ext/crlibm-1.0beta5/Makefile.am

@@ -0,0 +1,134 @@
+## Process this file with automake to produce Makefile.in
+
+lib_LIBRARIES = libcrlibm.a
+
+# shared libraries need libtool. We'll do it later.
+#lib_LTLIBRARIES = libcrlibm.la
+
+
+# Note that at the moment, asking for intervals disables double-extended
+if USE_HARDWARE_DE
+libcrlibm_a_SOURCES =  \
+	crlibm.h crlibm_private.h crlibm_private.c double-extended.h\
+	triple-double.h triple-double.c\
+	exp-td.h exp-td.c exp-td-standalone.c expm1-standalone.c  \
+	expm1.h expm1.c \
+	log-de.c log-de.h \
+	log1p.c \
+	log10-td.h log10-td.c \
+	log2-td.h  log2-td.c \
+	rem_pio2_accurate.h rem_pio2_accurate.c \
+	trigo_fast.c trigo_fast.h trigo_accurate.c trigo_accurate.h \
+	trigpi.c trigpi.h \
+	asincos.h asincos.c \
+	pow.h pow.c\
+	atan_fast.c atan_fast.h atan_accurate.h atan_accurate.c \
+	csh_fast.h csh_fast.c
+else 
+libcrlibm_a_SOURCES =  \
+	crlibm.h crlibm_private.h crlibm_private.c \
+	triple-double.h triple-double.c\
+	exp-td.h exp-td.c exp-td-standalone.c expm1-standalone.c  \
+	expm1.h expm1.c \
+	log.c log.h \
+	log1p.c \
+	rem_pio2_accurate.h rem_pio2_accurate.c \
+	trigo_fast.c trigo_fast.h trigo_accurate.c trigo_accurate.h \
+	trigpi.c trigpi.h \
+	asincos.h asincos.c \
+	pow.h pow.c\
+	atan_fast.c atan_fast.h atan_accurate.h atan_accurate.c \
+	csh_fast.h csh_fast.c 
+endif
+
+
+# add all the scs_lib useful functions
+libcrlibm_a_LIBADD = scs_lib/scs_private.o scs_lib/addition_scs.o  \
+	scs_lib/division_scs.o scs_lib/print_scs.o\
+	scs_lib/double2scs.o scs_lib/zero_scs.o\
+	scs_lib/multiplication_scs.o scs_lib/scs2double.o
+
+include_HEADERS = crlibm.h
+
+SUBDIRS = scs_lib .
+
+# Add your compiler here. Default is assumed to be gcc.
+if COMPILER_ICC
+# mp means maintain precision. Should be useless with std=c99, this is an ICC bug
+AM_CFLAGS =  -std=c99 -mp -Qoption,cpp,--extended_float_types -IPF_fp_speculationsafe
+else # default assumed to be gcc
+if USE_SSE2
+AM_CFLAGS = -std=c99 -Wall -msse2 -mfpmath=sse 
+else
+AM_CFLAGS = -std=c99 -Wall -O3
+endif
+endif
+
+tripledoubleprocedures.pdf: docs/latex/tripledoubleprocedures.tex
+	cd docs/latex; \
+	pdflatex tripledoubleprocedures; bibtex tripledoubleprocedures; \
+	pdflatex tripledoubleprocedures; pdflatex tripledoubleprocedures;\
+	mv tripledoubleprocedures.pdf ..
+
+crlibm.pdf: docs/latex/crlibm.tex
+	cd docs/latex; pdflatex crlibm; bibtex crlibm; pdflatex crlibm; pdflatex crlibm;\
+	mv crlibm.pdf .. 
+
+doc: crlibm.pdf tripledoubleprocedures.pdf
+
+EXTRA_DIST = VERSION TODO\
+	maple/common-procedures.mpl \
+	maple/double-extended.mpl maple/triple-double.mpl\
+	maple/exp-td.mpl maple/csh.mpl \
+	maple/log-td.mpl  maple/log-de.mpl \
+	maple/log2-td.mpl maple/log10-td.mpl \
+	maple/trigo.mpl maple/atan.mpl maple/atan-pentium.mpl \
+	maple/asin-td.mpl maple/acos-td.mpl \
+	maple/sqrt.mpl \
+	maple/expm1.mpl \
+	maple/gal.mpl \
+	maple/pow.mpl \
+	maple/atan-itanium.mpl \
+	gappa/trigos/SinCosCase3.gappa gappa/trigos/TanCase2.gappa \
+	gappa/trigos/TanCase2.gappa \
+	gappa/trigos/SinCosCase3.gappa gappa/trigos/TanCase2.gappa \
+	gappa/sqrt12.gappa gappa/sqrt13.gappa gappa/sqrt-proof.sh gappa/sqrt.sed \
+	gappa/log-de/log-de-acc-index0-E0.gappa    gappa/log-de/log-de-index0-E0.gappa \
+	gappa/log-de/log-de-acc-index0-E1N.gappa   gappa/log-de/log-de-index0-E1N.gappa \
+	gappa/log-de/log-de-acc-index1N-E0.gappa   gappa/log-de/log-de-index1N-E0N.gappa \
+	gappa/log-de/log-de-acc-index1N-E1N.gappa  gappa/log-de/README \
+	gappa/log-td/log-td.gappa             gappa/log-td/log-td-E0.gappa \
+	gappa/log-td/log-td-E0-logir0.gappa   gappa/log-td/log-td-accurate.gappa \
+	gappa/log-td/log-td-accurate-E0.gappa gappa/log-td/log-td-accurate-E0-logir0.gappa \
+	gappa/exp-td/exp-td-accurate1.gappa gappa/exp-td/exp-td-accurate2.gappa \
+	gappa/exp-td/exp-td-accurate3.gappa gappa/exp-td/exp-td-accurate4.gappa \
+	gappa/exp-td/exp-td-proof.sh \
+	gappa/trigpi/cospi-accurate.gappa  gappa/trigpi/cospi-quick.gappa \
+	gappa/trigpi/sinpi-accurate.gappa  gappa/trigpi/sinpi-quick.gappa \
+	gappa/asin/asinAccuR0.gappa gappa/asin/asinAccuR1.gappa  \
+	gappa/asin/asinAccuR2.gappa gappa/asin/asinAccuR3.gappa \
+	gappa/asin/asinAccuR4.gappa gappa/asin/asinAccuR5.gappa \
+	gappa/asin/asinAccuR6.gappa gappa/asin/asinAccuR7.gappa \
+	gappa/asin/asinAccuR8.gappa gappa/asin/asinAccuR9.gappa \
+	gappa/asin/asinQuickR0.gappa gappa/asin/asinQuickR1.gappa \
+	gappa/asin/asinQuickR2.gappa gappa/asin/asinQuickR3.gappa \
+	gappa/asin/asinQuickR4.gappa gappa/asin/asinQuickR5.gappa \
+	gappa/asin/asinQuickR6.gappa gappa/asin/asinQuickR7.gappa \
+	gappa/asin/asinQuickR8.gappa gappa/asin/asinQuickR9.gappa \
+	docs/latex/crlibm.tex \
+	docs/latex/0_intro.tex docs/latex/0_getting-started.tex\
+	docs/latex/1_common.tex docs/latex/sqrt.tex \
+	docs/latex/exp.tex docs/latex/expm1.tex \
+	docs/latex/log.tex docs/latex/log1p.tex docs/latex/log2.tex docs/latex/log10.tex\
+	docs/latex/asin.tex docs/latex/acos.tex docs/latex/atan.tex \
+	docs/latex/csh.tex docs/latex/trigo.tex docs/latex/trigpi.tex \
+	docs/latex/pow.tex \
+	docs/latex/fig_scs docs/latex/fig_exp\
+	docs/latex/elem-fun.bib\
+	docs/latex/tripledoubleprocedures.tex \
+	scs_lib/tests/tbx_timing.h \
+	tests/testperfs.sh\
+	triple-double.h double-extended.h \
+	log-de.c log-de.h log-td.c log-td.h\
+	atan-pentium.c atan-itanium.c 
+ 

data/ext/crlibm-1.0beta5/README

@@ -0,0 +1,31 @@
+cr-libm, a portable, efficient, correctly rounded mathematical library.
+
+To compile: 
+./configure; make;
+
+(for options see ./configure --help)
+
+To run the selftest:
+make check;
+
+
+
+Beta release note: 
+
+This version contains one state-of-the-art implementation. The log
+comes in two versions, log-td (portable) and log-de (optimised for
+processors with double-extended hardware). The ./configure should
+enable the one that is best for your system, except on Itanium: see
+below. 
+
+There is no known bug, but several shortcomings are listed in
+the TODO file.
+
+
+Note for Itanium systems: Out-of-the box it should compile and
+selftest on Linux under both gcc and icc (feedback on HPUX is
+welcome). However, optimised versions of log and atan are
+disabled. They currently require icc 8.1 or later. See the beginning
+of log-de.c and atan-itanium.c for instructions on how to enable them.
+
+

data/ext/crlibm-1.0beta5/README.DEV

@@ -0,0 +1,23 @@
+
+To compile crlibm after a cvs checkout, type
+
+./prepare
+make
+
+or, if you want to perform accuracy tests against the MPFR library,
+
+./prepare --enable-mpfr --with-gmp=/usr/local
+make
+
+(if necessary replace /usr/local/ with the actual location of the MPFR
+source)
+
+The prepare script requires to be invoked only once. After that, use:
+./configure 
+
+The prepare script ends with a call to ./configure, and the options to
+prepare are the options of configure. For a list of these options, type
+ 
+./configure --help
+
+

data/ext/crlibm-1.0beta5/README.md

@@ -0,0 +1,5 @@
+
+Retreived from <https://github.com/taschini/crlibm>, commit eb30637.
+
+Removed `docs/`, `gappa/`, `mapple/` and `tests/`.
+

data/ext/crlibm-1.0beta5/TODO

@@ -0,0 +1,66 @@
+TODO before 1.0 final release
+
+
+Serious bugs
+------------
+No known serious bug left so far. Exception raising should be improved.
+
+New functions
+-------------
+Document the cases on which power is guaranteed correctly rounded values.
+Possibly write 2^x and 10^x (in harmonious combination with power)
+
+Proof stuff
+-----------
+Update chapter asin and write chapter acos. 
+   Add the scripts generating the polynomials.
+Write Gappa proofs for log2-td, log10-td, expm1, log1p
+Check all the existing .gappa against last version of Gappa
+ + set the paths right (where there are paths). 
+Write proper proofs of the DE integer rounding tests, or remove DE altogether.
+
+
+Cleaning up stuff
+-----------------
+Check .mpl with underscored variables (some versions of Maple don't like them)
+Purge obsolete .mpl
+Ensure compilation on icc9 and SunStudio11
+Solution: check recent icc first. Then add FMADE macros everywhere FMA intended
+add make timings
+
+Autotest stuff
+--------------
+Add worst cases for all which are currently written TODO
+A lot of the  Lefevre worst cases are still missing
+You never have enough test cases. In particular asin and acos
+More cases around the various boundaries would be welcome.
+
+Optim stuff
+-----------
+Remove remaining SCS: in atan, in trigo.
+Replace the current test for non-IEEE subnormal (on x86) in the code
+    with a test taking into account that x86-64 is OK
+Minor improvements to atanpi as mentionned in the doc. 
+Fuse atan-pentium and atan-itanium into atan-de, and add in atanpi. 
+Or maybe remove all the DE stuff altogether for cleaning up.
+
+===========================================================================
+
+TODO some day for a later release
+
+New functions
+-------------
+Write  inverse hyperbolic functions
+
+Optimisation stuff
+------------------
+Write TD versions of: trigo
+Write DE versions of: most functions
+
+
+Use technique in paper by Karp and Markstein for sqrt12 and sqrt13
+
+
+Interval stuff
+--------------
+Write interval functions for each CR function

data/ext/crlibm-1.0beta5/VERSION

@@ -0,0 +1 @@
+1.0beta5