Import of the latest gdtoa sources, which include fixes for minor

problems relating to NaNs and rounding.

1 files changed, 10 insertions, 5 deletions

diff --git a/contrib/gdtoa/makefile b/contrib/gdtoa/makefile
index 3387a65..209fefc 100644
--- a/contrib/gdtoa/makefile
+++ b/contrib/gdtoa/makefile
@@ -30,13 +30,18 @@ CFLAGS = -g -DINFNAN_CHECK
 .c.o:
 	$(CC) -c $(CFLAGS) $*.c
 
-all: arith.h gdtoa.a
+all: arith.h gd_qnan.h gdtoa.a
 
 arith.h: arithchk.c
 	$(CC) $(CFLAGS) arithchk.c || $(CC) -DNO_LONG_LONG $(CFLAGS) arithchk.c
 	./a.out >arith.h
 	rm -f a.out arithchk.o
 
+gd_qnan.h: arith.h qnan.c
+	$(CC) $(CFLAGS) qnan.c
+	./a.out >gd_qnan.h
+	rm -f a.out qnan.o
+
 gdtoa.a: dmisc.c dtoa.c g_Qfmt.c g__fmt.c g_ddfmt.c g_dfmt.c g_ffmt.c\
 	 g_xLfmt.c g_xfmt.c gdtoa.c gethex.c gmisc.c hd_init.c hexnan.c\
 	 misc.c smisc.c strtoIQ.c strtoId.c strtoIdd.c strtoIf.c strtoIg.c\
@@ -51,10 +56,10 @@ gdtoa.a: dmisc.c dtoa.c g_Qfmt.c g__fmt.c g_ddfmt.c g_dfmt.c g_ffmt.c\
 
 xs0 = README arithchk.c dmisc.c dtoa.c g_Qfmt.c g__fmt.c g_ddfmt.c g_dfmt.c\
 	 g_ffmt.c g_xLfmt.c g_xfmt.c gdtoa.c gdtoa.h gdtoaimp.h gethex.c\
-	 gmisc.c hd_init.c hexnan.c makefile misc.c smisc.c strtoIQ.c\
+	 gmisc.c hd_init.c hexnan.c makefile misc.c qnan.c smisc.c strtoIQ.c\
 	 strtoId.c strtoIdd.c strtoIf.c strtoIg.c strtoIx.c strtoIxL.c\
-	 strtod.c strtodI.c strtodg.c strtof.c strtopQ.c strtopd.c strtopdd.c\
-	 strtopf.c strtopx.c strtopxL.c strtorQ.c strtord.c strtordd.c\
+	 strtod.c strtodI.c strtodg.c strtodnrp.c strtof.c strtopQ.c strtopd.c\
+	 strtopdd.c strtopf.c strtopx.c strtopxL.c strtorQ.c strtord.c strtordd.c\
 	 strtorf.c strtorx.c strtorxL.c sum.c ulp.c
 
 # "make xsum.out" to check for transmission errors; source for xsum is
@@ -66,4 +71,4 @@ xsum.out: xsum0.out $(xs0)
 	cmp xsum0.out xsum1.out && mv xsum1.out xsum.out || diff xsum[01].out
 
 clean:
-	rm -f arith.h *.[ao] xsum.out xsum1.out
+	rm -f arith.h gd_qnan.h *.[ao] xsum.out xsum1.out