Compare commits

..

2 Commits

Author SHA1 Message Date
Landon Curt Noll
97e9429000 Release calc version 2.11.0t4 2017-05-21 15:38:27 -07:00
Landon Curt Noll
1ce630ac19 Release calc version 2.11.0t3 2017-05-21 15:38:26 -07:00
34 changed files with 479 additions and 343 deletions

143
BUGS
View File

@@ -1,8 +1,8 @@
If you notice something wrong, strange or broken, try rereading:
README.FIRST
README
BUGS (in particular the bottom problems or mis-features section)
HOWTO.INSTALL
BUGS (this file)
If that does not help, cd to the calc source directory and try:
@@ -18,13 +18,13 @@ If it does not, then something is really broken!
If you made and modifications to calc beyond the simple Makefile
configuration, try backing them out and see if things get better.
Check to see if the version of calc you are using is current. Calc
distributions may be obtained from the official calc repository:
To be sure that your version of calc is up to date, check out:
ftp://ftp.uu.net/pub/calc
http://reality.sgi.com/chongo/calc/calc-download.html
If you are an alpha or beta tester, you may have a special pre-released
version that is more advanced than what is in the ftp archive.
The calc web site is located at:
http://reality.sgi.com/chongo/calc
=-=
@@ -46,6 +46,8 @@ When you send your report, please include the following information:
* the type of compiler you were using
* any compiler warnings or erros that you saw
* cd to the calc source directory, and type:
make debug > debug.out 2>&1 (sh, ksh, bash users)
@@ -57,22 +59,8 @@ Stack traces from core dumps are useful to send as well.
=-=
The official calc repository is located in:
ftp://ftp.uu.net/pub/calc
If you don't have ftp access to that site, or if your version is more
recent than what has been released to the ftp archive, you may, as a
last resort, send EMail to:
chongo@toad.com
Indicate the version you have and that you would like a more up to date version.
=-=
Send any comments, suggestions and most importantly, fixes (in the form
of a context diff patch) to:
Send any comments, compiler warning messages, suggestions and most
importantly, fixes (in the form of a context diff patch) to:
calc-tester@postofc.corp.sgi.com
@@ -118,103 +106,22 @@ Known problems or mis-features:
8 getcommands(toplevel = 536864560) ["codegen.c":122, 0x1200284dc]
9 main() ["calc.c":279, 0x12001f61c]
However On 25-Sep-99 9:18, Ernest Bowen wrote:
Michel van der List reports that he was using the Dec compiler.
Last night I downloaded 2.11.0t1 and without changing Makefile
or anything else did a make and make check. Absolutely no
reported errors or warnings and the check "passed all tests".
The machine I'm using is a Dec 20 alpha running Digitial
UNIX V4.0B (Rev 564).
However On 25-Sep-99 9:18, Ernest Bowen wrote:
Any ideas on how to reconcile these different results?
Last night I downloaded 2.11.0t1 and without changing Makefile
or anything else did a make and make check. Absolutely no
reported errors or warnings and the check "passed all tests".
The machine I'm using is a Dec 20 alpha running Digitial
UNIX V4.0B (Rev 564).
* On Sun SPARCstation running Solaris 2.5.1 compiling with
cc version 4.0 with -Xc option 'make check' fails on the hash
functions. However when compiled with gcc everything seems to
work just fine.
Ernest Bowen was using gcc.
From: "Dr.D.J.Picton" <dave@aps5.ph.bham.ac.uk>
Date: Fri, 8 May 1998 10:54:43 +0100
We definitely want to hunt down this Dec's cc compiler issue. It
could very well be that gcc is letting something slide where as
Dec's cc compiler is compiling it faithfully into core dumping code.
The fault could very well be the calc code and not the compiler.
Platform: Sun SPARCstation
Operating system: Solaris 2.5.1
Compiler: cc version 4.0 with -Xc option.
The code compiled OK but the sha() and sha1() functions give errors with
make check. These functions clearly are broken. For example,
sha(sha()) gives
0xea19f96cdd568ad1084a17acc6103d92c0ef7708 not
0xf96cea198ad1dd5617ac084a3d92c6107708c0ef
Interestingly, the consistency checks work OK e.g.
>B = blk() = {"a", "b", "c"};
>sha(B) == sha("abc");
1
**** Non-true result (0): 7108: sha(sha()) ==
0xf96cea198ad1dd5617ac084a3d92c6107708c0ef
**** Non-true result (0): 7109:
sha(sha("a"))==0x37f297772fae4cb1ba39b6cf9cf0381180bd62f2
**** Non-true result (0): 7110:
sha(sha("ab"))==0x488373d362684af3d3f7a6a408b59dfe85419e09
**** Non-true result (0): 7111:
sha(sha("abc"))==0x0164b8a914cd2a5e74c4f7ff082c4d97f1edf880
**** Non-true result (0): 7112:
sha(sha("abcd"))==0x082c73b06f71185d840fb4b28eb3abade67714bc
**** Non-true result (0): 7113:
sha(sha("abcde"))==0xd624e34951bb800f0acae773001df8cffe781ba8
**** Non-true result (0): 7114:
sha(sha("abcdef"))==0x2a589f7750598dc0ea0a608719e04327f609279a
**** Non-true result (0): 7115:
sha(sha("abcdefg"))==0x5bdf01f9298e9d19d3f8d15520fd74eed600b497
**** Non-true result (0): 7116:
sha(sha("abcdefgh"))==0x734ba8b31975d0dbae4d6e249f4e8da270796c94
**** Non-true result (0): 7117: sha(sha(1)) ==
0x864c8d09e828c7c31d62693736a5a9302c282777
**** Non-true result (0): 7118: sha(sha(2)) ==
0x2c0b59c512cb20fad6bb0883b69c9f5a46545808
**** Non-true result (0): 7119:
sha(sha(22/7)==0x7ddb7f8a7e9d70757f157744fddea7a6c6a6bcc6
**** Non-true result (0): 7120:
sha(sha(isqrt(2e1000)==0x6db8d9cf0b018b8f9cbbf5aa1edb8066d19e1bb0
7121: sha("x", "y", "z") == sha("xyz")
**** Non-true result (0): 7122: sha(sha("this is", 7^19-8, ..., "hash")) ==
0x21e4...
7123: z = sha(list(1,2,3), "curds and whey", 2^21701-1, pi());
**** Non-true result (0): 7124: sha(z) ==
0x36dcca3e51865c30a2cf738023cda446f1368340
**** Non-true result (0): 7207: sha1(sha1()) ==
0xda39a3ee5e6b4b0d3255bfef95601890afd80709
7208: sha1("x", "y", "z") == sha1("xyz")
**** Non-true result (0): 7209: sha1(sha1("this is",7^19-8,"a
composit",3i+4.5,"hash")) == ...
7210: z = sha1(list(1,2,3), "curds and whey", 2^21701-1, pi());
**** Non-true result (0): 7211: sha1(z) ==
0xc19e7317675dbf71e293b4c41e117169e9da5b6f
**** Non-true result (0): 7217:
sha1(sha1("a"))==0x86f7e437faa5a7fce15d1ddcb9eaeaea377667b8
**** Non-true result (0): 7218:
sha1(sha1("ab"))==0xda23614e02469a0d7c7bd1bdab5c9c474b1904dc
**** Non-true result (0): 7219:
sha1(sha1("abc"))==0xa9993e364706816aba3e25717850c26c9cd0d89d
**** Non-true result (0): 7220:
sha1(sha1("abcd"))==0x81fe8bfe87576c3ecb22426f8e57847382917acf
**** Non-true result (0): 7221:
sha1(sha1("abcde"))==0x03de6c570bfe24bfc328ccd7ca46b76eadaf4334
**** Non-true result (0): 7222:
sha1(sha1("abcdef"))==0x1f8ac10f23c5b5bc1167bda84b833e5c057a77d2
**** Non-true result (0): 7223:
sha1(sha1("abcdefg"))==0x2fb5e13419fc89246865e7a324f476ec624e8740
**** Non-true result (0): 7224:
sha1(sha1("abcdefgh"))==0x425af12a0743502b322e93a015bcf868e324d56a
**** Non-true result (0): 7225:
sha1(sha1(1))==0x53dd4e1734ad47d45e41c23e4ce42d7f1f98c1e8
**** Non-true result (0): 7226:
sha1(sha1(22/7))==0xf8e2510f85f7b9bf088b321188e9f70620f44246
**** Non-true result (0): 7227:
sha1(sha1(isqrt(2e1000)))==0x6852a1365c51050c3d039e3c5d9cf29c12283ef4
Any more information on calc core dumps compiled with Dec's cc compiler
would be most appreciated!

69
CHANGES
View File

@@ -1,15 +1,13 @@
Following is the change from calc version 2.11.0t1 to date:
Removed the makefile symbol MAIN. Now forcing all funcions to correctly
be reclared main. To passify some old broken compilers, a return 0;
Removed the makefile symbol MAIN. Now forcing all functions to correctly
be declared main. To satisfy some old broken compilers, a return 0;
(instead of an exit(0);) is used at the end of main().
A few of files that were added to calc used 4 character indentation
whereas most of calc uses 8 character indentation. These imported
sources have been changed to conform better with the calc style.
Fixed misc compile warnings.
Added the program calc_errno.c and the Makefile symbol ERRNO_DECL.
If ERRNO_DECL is empty, calc_errno.c will try various ways to
declare errno, sys_errlist and sys_nerr. On success or when
@@ -21,6 +19,69 @@ Following is the change from calc version 2.11.0t1 to date:
The func.c file now includes the constructed header file calc_errno.h
to ensure that errno, sys_errlist and sys_nerr are declared correctly.
Changed check.awk to be more 'old awk' friendly.
Made some of the source a little more ++ friendly. We are NOT
porting calc to C++! We will NOT support C++ compilation of calc.
Calc will written ANSI C. We just compiled with a suggestion from
Love-Jensen, John <jlove-jensen@globalmt.com> to make calc's version
of C a little more to C++ compilers. Avoiding symbols such as new
or try for example.
Thanks to the efforts of Ernest Bowen <ernie@neumann.une.edu.au> and
Dr.D.J.Picton <dave@aps5.ph.bham.ac.uk>, a nasty endian-ness bug
in the sha and sha1 hash functions that showed up on machines such
as the Sparc was fixed.
Updated README.FIRST and BUGS to reflect new URLs and addresses.
Added a HOWTO.INSTALL file.
Reordered cc Makefile variable sets in the main Makefile.
Fixed a bug and applied a fix that was reported by Ernest Bowen
<ernie@turing.une.edu.au>. Added regression tests 1103 to 1112.
Fixed a bug in version.c related to MINOR_PATCHs in both the
empty and non-empty MINOR_PATCH cases.
Fixed malloc and bad storage issues reported by Michel van der List
<vanderlistmj@sbphrd.com>.
Fixed some problems related to path processing while opening files.
Under extreme cases, an excessively long filename or CALCPATH value
could create problems. Placed guards in opensearchfile() function
in input.c to catch these cases.
Fixed cases were malloc failures were silently ignored in input.c.
Eliminated the PATHSIZE limit and the PATHSIZE symbol.
Added MAX_CALCRC to limit the length of the $CALCRC environment
variable to 1024 chars.
Fixed the magic number relating to the initial number of constants
declared by initconstants(). It is now related to the length
of the initnumbs[] NUMBER array.
Added a 'Dec Alpha / Compaq Tru64 cc (non-gnu) compiler set'
section to the main Makefile.
Fixed a string handling bug discovered by Dr.D.J.Picton
<dave@aps5.ph.bham.ac.uk> in the custom demo code.
Added FORCE_STDC symbol. When defined it will force __STDC__ like
conditions. Thus for compilers with as the Solaris cc compiler
that are ANSI-like but still define __STDC__ as 0, one can use
-DFORCE_STDC and make use of ANSI-like features.
Removed the CCSHS symbol from the Makefile. The shs.c and shs1.c
files are now compiled with full ${CFLAGS}.
The custom.c file is now compiled with full ${CFLAGS}.
Fixed misc compile warnings and notices.
Following is the change from calc version 2.10.3t5.38 to 2.11.0t0:

53
HOWTO.INSTALL Normal file
View File

@@ -0,0 +1,53 @@
Installing calc in 4 easy steps:
1) Look at the makefile, and adjust it to suit your needs.
Here are some Makefile hints:
Select a compiler set by commenting in the appropriate
set of cc options. As shipped the Makefile assumes
a gcc-like environment such as Linux.
If a more appropriate cc set if found below, comment
out the Linux set and comment in that set or edit
the gcc set or the common cc set as needed.
You may or may not need RANLIB when building libraries.
As shipped the Makefile assumes RANLIB is needed.
Comment the in/out the RANLIB value if ranlib does
not work or does not exist.
You may want to change the default pager used by calc.
As shipped the Makefile assumes 'more'. On your system
you may find 'less' to be a better pager.
The CALCBINDINGS is matter of personal taste. As shipped
the Makefile assumes a default quasi-emacs-like command
line editor. Changing CALCBINDINGS= altbind will cause ^D
to end calc in a fashion similar to that of the bc(1) command.
Set TOPDIR to be the place under which help files, calc,
include files and calc libs are to be installed. As shipped
the Makefile assumes a TOPDIR of /usr/local/lib.
Set BINDIR to the place where calc is installed. As shipped
the Makefile assumes a BINDIR /usr/local/bin.
Adjust other Makefile variables as needed.
2) build calc:
make all
==> We are interested in any compiler warnings (and errors) that
you may find. See the BUGS file if you find any.
3) test calc:
make check
==> If you run into problems, follow the instructions in the BUGS file
4) install calc:
make install

154
Makefile
View File

@@ -314,7 +314,8 @@ CATEXT= 1
# The cat page is not built or installed
#
# If in doubt and you don't want to fool with man pages, set MANDIR
# and CATDIR to empty and ignore the lines below.
# and CATDIR to empty and ignore the NROFF, NROFF_ARG and MANMAKE
# lines below.
#
NROFF= nroff
#NROFF=
@@ -516,7 +517,6 @@ ALLOW_CUSTOM= -DCUSTOM
# ICFLAGS are given to ${CC} for intermediate progs
#
# CCMAIN are flags for ${CC} when files with main() instead of CFLAGS
# CCSHS are flags given to ${CC} for compiling shs.c & shs1.c instead of CFLAGS
#
# LCFLAGS are CC-style flags for ${LINT}
# LDFLAGS are flags given to ${CC} for linking .o files
@@ -526,6 +526,53 @@ ALLOW_CUSTOM= -DCUSTOM
#
###
#
# Linux set
#
# Tested on Red Hat 6.0 Linux but should run on almost any Linux release.
#
# for better performance, set the following above:
# DEBUG= -O2
#
CCWARN= -Wall -Wno-implicit -Wno-comment
CCOPT= ${DEBUG} ${NO_SHARED}
CCMISC=
#
CFLAGS= ${CCWARN} ${CCOPT} ${CCMISC}
ICFLAGS= ${CCWARN} ${CCMISC}
#
CCMAIN= ${ICFLAGS}
#
LCFLAGS=
LDFLAGS= ${NO_SHARED} ${LD_NO_SHARED}
ILDFLAGS=
#
CC= ${PURIFY} gcc
#
###
#
# gcc set (some call it gcc2, some call it gcc)
#
# for better performance, set the following above:
# DEBUG= -O2
#
#CCWARN= -Wall -Wno-implicit -Wno-comment
#CCOPT= ${DEBUG} ${NO_SHARED}
#CCMISC=
#
#CFLAGS= ${CCWARN} ${CCOPT} ${CCMISC}
#ICFLAGS= ${CCWARN} ${CCMISC}
#
#CCMAIN= ${ICFLAGS}
#
#LCFLAGS=
#LDFLAGS= ${NO_SHARED} ${LD_NO_SHARED}
#ILDFLAGS=
#
#CC= ${PURIFY} gcc
#CC= ${PURIFY} gcc2
#
###
#
# common cc set
#
#CCWARN=
@@ -536,38 +583,6 @@ ALLOW_CUSTOM= -DCUSTOM
#ICFLAGS= ${CCWARN} ${CCMISC}
#
#CCMAIN= ${ICFLAGS}
#CCSHS= ${CFLAGS}
#
#LCFLAGS=
#LDFLAGS= ${NO_SHARED} ${LD_NO_SHARED}
#ILDFLAGS=
#
#CC= ${PURIFY} cc
#
###
#
# SGI IRIX5.3 (or earlier) -o32 C Compiler
#
# You must set above:
# RANLIB=:
# LONGLONG_BITS= 0
#
# for better performance, set the following above:
# DEBUG= -O2 -g3
#
# If you have the directory /usr/lib/nonshared, then set the following above:
# NO_SHARED= -non_shared
# LD_NO_SHARED= -Wl,-rdata_shared
#
#CCWARN= -fullwarn -woff 835
#CCOPT= ${DEBUG} ${NO_SHARED}
#CCMISC= -use_readonly_const
#
#CFLAGS= ${CCWARN} ${CCOPT} ${CCMISC}
#ICFLAGS= ${CCWARN} ${CCMISC}
#
#CCMAIN= ${ICFLAGS}
#CCSHS= ${CFLAGS}
#
#LCFLAGS=
#LDFLAGS= ${NO_SHARED} ${LD_NO_SHARED}
@@ -599,7 +614,6 @@ ALLOW_CUSTOM= -DCUSTOM
#ICFLAGS= ${CCWARN} ${CCMISC}
#
#CCMAIN= ${ICFLAGS}
#CCSHS=
#
#LCFLAGS=
#LDFLAGS= ${NO_SHARED} ${LD_NO_SHARED}
@@ -625,7 +639,6 @@ ALLOW_CUSTOM= -DCUSTOM
#ICFLAGS= ${CCWARN} ${CCMISC}
#
#CCMAIN= ${ICFLAGS}
#CCSHS= ${CFLAGS}
#
#LCFLAGS=
#LDFLAGS= ${NO_SHARED} ${LD_NO_SHARED}
@@ -645,7 +658,6 @@ ALLOW_CUSTOM= -DCUSTOM
#ICFLAGS= ${CCWARN} ${CCMISC}
#
#CCMAIN= ${ICFLAGS}
#CCSHS= ${CFLAGS}
#
#LCFLAGS=
#LDFLAGS= ${NO_SHARED} ${LD_NO_SHARED}
@@ -655,20 +667,22 @@ ALLOW_CUSTOM= -DCUSTOM
#
###
#
# Solaris 2.x Sun cc compiler
# Solaris Sun cc compiler set
#
# for better performance, set the following above:
# DEBUG= -O
#
# We need -DFORCE_STDC to make use of ANSI-C like features and
# to avoid the use of -Xc (which as a lose performance wise).
#
#CCWARN=
#CCOPT= ${DEBUG} ${NO_SHARED}
#CCMISC=-Xc
#CCMISC= -DFORCE_STDC
#
#CFLAGS= ${CCWARN} ${CCOPT} ${CCMISC}
#ICFLAGS= ${CCWARN} ${CCMISC}
#
#CCMAIN= ${ICFLAGS}
#CCSHS= ${CFLAGS}
#
#LCFLAGS=
#LDFLAGS= ${NO_SHARED} ${LD_NO_SHARED}
@@ -678,30 +692,7 @@ ALLOW_CUSTOM= -DCUSTOM
#
###
#
# Red Hat Linux 6.0 and BSDI's BSD/OS 4.0 (or later) set
#
# for better performance, set the following above:
# DEBUG= -O2
#
CCWARN= -Wall -Wno-implicit -Wno-comment
CCOPT= ${DEBUG} ${NO_SHARED}
CCMISC=
#
CFLAGS= ${CCWARN} ${CCOPT} ${CCMISC}
ICFLAGS= ${CCWARN} ${CCMISC}
#
CCMAIN= ${ICFLAGS}
CCSHS= ${CFLAGS}
#
LCFLAGS=
LDFLAGS= ${NO_SHARED} ${LD_NO_SHARED}
ILDFLAGS=
#
CC= ${PURIFY} gcc
#
###
#
# gcc set (some call it gcc2, some call it gcc)
# BSDI's BSD/OS 4.0 (or later) set
#
# for better performance, set the following above:
# DEBUG= -O2
@@ -714,14 +705,34 @@ CC= ${PURIFY} gcc
#ICFLAGS= ${CCWARN} ${CCMISC}
#
#CCMAIN= ${ICFLAGS}
#CCSHS= ${CFLAGS}
#
#LCFLAGS=
#LDFLAGS= ${NO_SHARED} ${LD_NO_SHARED}
#ILDFLAGS=
#
#CC= ${PURIFY} gcc
#CC= ${PURIFY} gcc2
#
###
#
# Dec Alpha / Compaq Tru64 cc (non-gnu) compiler set
#
# for better performance, set the following above:
# DEBUG= -O2
#
#CCWARN=
#CCOPT= ${DEBUG} ${NO_SHARED}
#CCMISC=
#
#CFLAGS= ${CCWARN} ${CCOPT} ${CCMISC}
#ICFLAGS= ${CCWARN} ${CCMISC}
#
#CCMAIN= ${ICFLAGS}
#
#LCFLAGS=
#LDFLAGS= ${NO_SHARED} ${LD_NO_SHARED}
#ILDFLAGS=
#
#CC= ${PURIFY} cc
##############################################################################
#-=-=-=-=-=-=-=-=- Be careful if you change something below -=-=-=-=-=-=-=-=-#
@@ -920,7 +931,7 @@ C_SRC= ${LIBSRC} ${CALCSRC} ${UTIL_C_SRC}
# These files are found (but not built) in the distribution
#
DISTLIST= ${C_SRC} ${H_SRC} ${MAKE_FILE} BUGS CHANGES LIBRARY README \
calc.man lint.sed README.FIRST ${UTIL_MISC_SRC}
calc.man lint.sed README.FIRST HOWTO.INSTALL ${UTIL_MISC_SRC}
# complete list of .o files
#
@@ -973,17 +984,11 @@ calc.o: calc.c ${MAKE_FILE}
${CC} ${CCMAIN} ${CCOPT} ${ALLOW_CUSTOM} -c calc.c
custom.o: custom.c ${MAKE_FILE}
${CC} ${CCOPT} ${ALLOW_CUSTOM} -c custom.c
${CC} ${CFLAGS} ${ALLOW_CUSTOM} -c custom.c
hist.o: hist.c ${MAKE_FILE}
${CC} ${CFLAGS} ${TERMCONTROL} -c hist.c
shs.o: shs.c ${MAKE_FILE}
${CC} ${CCSHS} -c shs.c
shs1.o: shs1.c ${MAKE_FILE}
${CC} ${CCSHS} -c shs1.c
func.o: func.c ${MAKE_FILE}
${CC} ${CFLAGS} ${ALLOW_CUSTOM} -c func.c
@@ -2328,7 +2333,6 @@ env:
@echo "CFLAGS=${CFLAGS}"; echo ""
@echo "ICFLAGS=${ICFLAGS}"; echo ""
@echo "CCMAIN=${CCMAIN}"; echo ""
@echo "CCSHS=${CCSHS}"; echo ""
@echo "LCFLAGS=${LCFLAGS}"; echo ""
@echo "LDFLAGS=${LDFLAGS}"; echo ""
@echo "ILDFLAGS=${ILDFLAGS}"; echo ""

5
README
View File

@@ -4,6 +4,11 @@
#
# Arbitrary precision calculator.
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
NOTE: This is an old historic README. We recommend that you
read README.FIRST and HOWTO.INSTALL for more info.
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
I am allowing this calculator to be freely distributed for your enjoyment.
Like all multi-precision programs, you should not depend absolutely on
its results, since bugs in such programs can be insidious and only rarely

View File

@@ -1,12 +1,17 @@
Dear alpha tester,
Dear calc user,
Thanks for taking the time to try out this alpha version of calc! We are
interested in any/all feedback that you may have on this version. In
particular we would like to hear about:
See the HOWTO.INSTALL file for information on how to build and install calc.
To be sure that your version of calc is up to date, check out:
http://reality.sgi.com/chongo/calc/calc-download.html
We are interested in any/all feedback on recent versions of calc.
In particular we would like to hear about:
* compiler warnings
* compile problems
* regression test problems (try: make check)
* compiler warnings
* special compile flags/options that you needed
* Makefile problems
* help file problems
@@ -15,38 +20,35 @@ particular we would like to hear about:
We would like to offer a clean compile across a wide verity of platforms,
so if you can test on several, so much the better!
Calc distributions may be obtained from:
ftp://ftp.uu.net/pub/calc
If you don't have ftp access to that site, or if you do not find a more
recent version (you may have a special pre-released version that is
more advanced than what is in the ftp archive) send EMail to:
chongo@toad.com
Indicate the version you have and that you would like a more up
to date version.
If you run into problems, see the BUGS file.
=-=
Misc items TODO before Beta release:
See the file:
* improve the coverage in the 'SEE ALSO' help file lists
help/todo
* where reasonable, be sure that regress.cal tests builtin functions
or run:
* add the Blum-Blum-Shub random() generator code
calc 'help todo'
* add code to allow of the reading, writing and processing of binary data
for a wish/todo list. Code contributions are welcome.
* add shs, shs-1 and md5 hashing functions. Use align32.h.
=-=
* add mod h*2^n+/-1 function for integers
To join the calc-tester mailing list. Send a request to:
* be sure that CHANGES is up to date,
look over the help/todo file and update as needed,
revisit issues in the BUGS file and
change this file :-)
calc-tester-request@postofc.corp.sgi.com
* clean the source code and document it better
Your message body (not the subject) should consist of:
subscribe calc-tester address
end
name your_full_name
where ``address'' is your EMail address and ``your_full_name'' is
your full name.
The calc web site is located at:
http://reality.sgi.com/chongo/calc

View File

@@ -17,7 +17,7 @@
#ifdef HAVE_MALLOC_H
# include <malloc.h>
#else
# if defined(__STDC__) && __STDC__ != 0
#if defined(FORCE_STDC) || (defined(__STDC__) && __STDC__ != 0) || defined(__cplusplus)
extern void *malloc();
extern void *realloc();
extern void free();
@@ -36,7 +36,7 @@
# if defined(HAVE_NEWSTR)
extern void *memcpy();
extern void *memset();
# if defined(__STDC__) && __STDC__ != 0
#if defined(FORCE_STDC) || (defined(__STDC__) && __STDC__ != 0) || defined(__cplusplus)
extern size_t strlen();
# else
extern long strlen();
@@ -65,7 +65,7 @@ extern int strcmp();
#if !defined(HAVE_MEMMOVE)
# undef CALC_SIZE_T
# if defined(__STDC__) && __STDC__ != 0
#if defined(FORCE_STDC) || (defined(__STDC__) && __STDC__ != 0) || defined(__cplusplus)
# define CALC_SIZE_T size_t
# else
# define CALC_SIZE_T long

48
block.c
View File

@@ -71,7 +71,7 @@ static void blkchk(BLOCK*);
BLOCK *
blkalloc(int len, int chunk)
{
BLOCK *new; /* new block allocated */
BLOCK *nblk; /* new block allocated */
/*
* firewall
@@ -84,8 +84,8 @@ blkalloc(int len, int chunk)
/*
* allocate BLOCK
*/
new = (BLOCK *)malloc(sizeof(BLOCK));
if (new == NULL) {
nblk = (BLOCK *)malloc(sizeof(BLOCK));
if (nblk == NULL) {
math_error("cannot allocate block");
/*NOTREACHED*/
}
@@ -93,23 +93,23 @@ blkalloc(int len, int chunk)
/*
* initialize BLOCK
*/
new->blkchunk = chunk;
new->maxsize = ((len+chunk)/chunk)*chunk;
new->data = (USB8*)malloc(new->maxsize);
if (new->data == NULL) {
nblk->blkchunk = chunk;
nblk->maxsize = ((len+chunk)/chunk)*chunk;
nblk->data = (USB8*)malloc(nblk->maxsize);
if (nblk->data == NULL) {
math_error("cannot allocate block data storage");
/*NOTREACHED*/
}
memset(new->data, 0, new->maxsize);
new->datalen = len;
memset(nblk->data, 0, nblk->maxsize);
nblk->datalen = len;
/*
* return BLOCK
*/
if (conf->calc_debug > 0) {
blkchk(new);
blkchk(nblk);
}
return new;
return nblk;
}
@@ -227,7 +227,7 @@ blkchk(BLOCK *blk)
BLOCK *
blkrealloc(BLOCK *blk, int newlen, int newchunk)
{
USB8 *new; /* realloced storage */
USB8 *nblk; /* realloced storage */
int newmax; /* new maximum stoage size */
/*
@@ -258,20 +258,20 @@ blkrealloc(BLOCK *blk, int newlen, int newchunk)
if (newmax != blk->maxsize) {
/* reallocate new storage */
new = (USB8*)realloc(blk->data, newmax);
if (new == NULL) {
nblk = (USB8*)realloc(blk->data, newmax);
if (nblk == NULL) {
math_error("cannot reallocate block storage");
/*NOTREACHED*/
}
/* clear any new storage */
if (newmax > blk->maxsize) {
memset(new + blk->maxsize, 0, (newmax - blk->maxsize));
memset(nblk+blk->maxsize, 0, (newmax-blk->maxsize));
}
blk->maxsize = newmax;
/* restore the data pointers */
blk->data = new;
blk->data = nblk;
}
/*
@@ -391,13 +391,13 @@ blktrunc(BLOCK *blk)
BLOCK *
blk_copy(BLOCK *blk)
{
BLOCK *new; /* copy of blk */
BLOCK *nblk; /* copy of blk */
/*
* malloc new block
*/
new = (BLOCK *)malloc(sizeof(BLOCK));
if (new == NULL) {
nblk = (BLOCK *)malloc(sizeof(BLOCK));
if (nblk == NULL) {
math_error("blk_copy: cannot malloc BLOCK");
/*NOTREACHED*/
}
@@ -405,18 +405,18 @@ blk_copy(BLOCK *blk)
/*
* duplicate most of the block
*/
*new = *blk;
*nblk = *blk;
/*
* duplicate block data
*/
new->data = (USB8 *)malloc(blk->maxsize);
if (new->data == NULL) {
nblk->data = (USB8 *)malloc(blk->maxsize);
if (nblk->data == NULL) {
math_error("blk_copy: cannot duplicate block data");
/*NOTREACHED*/
}
memcpy(new->data, blk->data, blk->maxsize);
return new;
memcpy(nblk->data, blk->data, blk->maxsize);
return nblk;
}

5
calc.c
View File

@@ -212,7 +212,7 @@ main(int argc, char **argv)
/*
* if allowed or needed, print version and setup bindings
*/
if (*str == '\0') {
if (str >= cmdbuf+MAXCMD || *str == '\0') {
/*
* check for pipe mode and/or non-tty stdin
*/
@@ -349,7 +349,8 @@ math_error(char *fmt, ...)
if (post_init) {
longjmp(jmpbuf, 1);
} else {
fprintf(stderr, "no jmpbuf jumpback point - ABORTING!!!\n");
fprintf(stderr, "It is too early provide a command line prompt "
"so we must simply exit. Sorry!\n");
/*
* don't call libcalc_call_me_last() -- we might loop
* and besides ... this is an unusual internal error case

2
calc.h
View File

@@ -26,7 +26,7 @@
#define DEFAULTCALCHELP "help" /* help file that -h prints */
#define DEFAULTSHELL "sh" /* default shell to use */
#define CALCEXT ".cal" /* extension for files read in */
#define PATHSIZE 1024 /* maximum length of path name */
#define MAX_CALCRC 1024 /* maximum length of $CALCRC */
#define HOMECHAR '~' /* char which indicates home directory */
#define DOTCHAR '.' /* char which indicates current directory */
#define PATHCHAR '/' /* char which separates path components */

View File

@@ -71,5 +71,9 @@ END {
if (error > 0 && havebuf0) {
print buf0;
}
exit (error > 0);
if (error > 0) {
exit(1);
} else {
exit(0);
}
}

View File

@@ -81,8 +81,13 @@ static long getinitlist(void);
void
getcommands(BOOL toplevel)
{
char name[PATHSIZE+1]; /* program name */
char name[MAXCMD+1+1]; /* program name */
/* firewall */
name[0] = '\0';
name[MAXCMD+1] = '\0';
/* getcommands */
if (!toplevel)
enterfilescope();
for (;;) {

14
const.c
View File

@@ -27,10 +27,11 @@ initconstants(void)
math_error("Unable to allocate constant table");
/*NOTREACHED*/
}
for (i = 0; i < 8; i++)
for (i = 0; i < INITCONSTCOUNT; i++)
consttable[i] = initnumbs[i];
constcount = 8;
constavail = CONSTALLOCSIZE - 8;
consttable[INITCONSTCOUNT] = NULL; /* firewall */
constcount = INITCONSTCOUNT;
constavail = CONSTALLOCSIZE - INITCONSTCOUNT;
}
@@ -187,10 +188,9 @@ freeconstant(unsigned long index)
void
trimconstants(void)
{
NUMBER **qp;
qp = &consttable[constcount];
while (constcount > 0 && (*--qp)->links == 0) {
while (constcount > 0 &&
(consttable[constcount-1] == NULL ||
consttable[constcount-1]->links == 0)) {
constcount--;
constavail++;
}

View File

@@ -99,6 +99,7 @@ static struct infoname sys_info[] = {
{"LONGLONG_BITS", "length of a long long, or 0", NULL, (FULL)LONGLONG_BITS},
{"LONG_BITS", "bit length of a long", NULL, (FULL)LONG_BITS},
{"MAP_POPCNT", "number of odd primes in pr_map", NULL, (FULL)MAP_POPCNT},
{"MAX_CALCRC", "maximum allowed length of $CALCRC", NULL, (FULL)MAX_CALCRC},
{"MAXCMD", "max length of command invocation", NULL, (FULL)MAXCMD},
{"MAXDIM", "max number of dimensions in matrices", NULL, (FULL)MAXDIM},
{"MAXERROR", "max length of error message string", NULL, (FULL)MAXERROR},
@@ -128,8 +129,6 @@ static struct infoname sys_info[] = {
{"NXT_MAP_PRIME", "smallest odd prime not in pr_map", NULL, (FULL)NXT_MAP_PRIME},
{"NXT_PFACT_VAL", "next prime for higher pfact values", NULL, (FULL)NXT_PFACT_VAL},
{"OFF_T_BITS", "file offset size in bits", NULL, (FULL)OFF_T_BITS},
{"PATHSIZE", "max length of path name", NULL, (FULL)PATHSIZE},
{"PATHSIZE", "max length of path name", NULL, (FULL)PATHSIZE},
{"PIX_32B", "max pix() value", NULL, (FULL)PIX_32B},
{"POW_ALG2", "default size for using REDC for powers", NULL, (FULL)POW_ALG2},
{"REDC_ALG2", "default size using alternative REDC alg", NULL, (FULL)REDC_ALG2},
@@ -238,8 +237,9 @@ c_sysinfo(char *name, int count, VALUE **vals)
} else if (vals[0]->v_type == V_STR) {
/* convert vals[0] to upper case string */
buf = (char *)malloc(strlen((char *)vals[0]->v_str)+1);
for (q = (char *)vals[0]->v_str, r = buf; *q; ++q, ++r) {
buf = (char *)malloc(strlen((char *)vals[0]->v_str->s_str)+1);
for (q = (char *)vals[0]->v_str->s_str, r = buf; *q; ++q, ++r)
{
if (isascii(*q) && islower(*q)) {
*r = *q - 'a' + 'A';
} else {
@@ -262,7 +262,7 @@ c_sysinfo(char *name, int count, VALUE **vals)
/* return value as string */
result.v_type = V_STR;
result.v_subtype = V_NOSUBTYPE;
result.v_str = (STRING *)p->str;
result.v_str = makestring(p->str);
}
/* return found infotype as value */

10
hash.c
View File

@@ -147,13 +147,13 @@ hash_free(HASH *state)
HASH *
hash_copy(HASH *state)
{
HASH *new; /* copy of state */
HASH *hnew; /* copy of state */
/*
* malloc new state
*/
new = (HASH *)malloc(sizeof(HASH));
if (new == NULL) {
hnew = (HASH *)malloc(sizeof(HASH));
if (hnew == NULL) {
math_error("hash_init: cannot malloc HASH");
/*NOTREACHED*/
}
@@ -161,8 +161,8 @@ hash_copy(HASH *state)
/*
* duplicate state
*/
memcpy((void *)new, (void *)state, sizeof(HASH));
return new;
memcpy((void *)hnew, (void *)state, sizeof(HASH));
return hnew;
}

View File

@@ -54,7 +54,7 @@ char buf[BUFSIZ];
void
try(char *fmt, ...)
try_this(char *fmt, ...)
{
va_list ap;
@@ -79,7 +79,7 @@ main(void)
/*
* test variable args and vsprintf/sprintf
*/
try("@%d:%s:%d@", 1, "hi", 2);
try_this("@%d:%s:%d@", 1, "hi", 2);
if (strcmp(buf, "@1:hi:2@") != 0) {
#if !defined(DONT_HAVE_VSPRINTF)
/* <stdarg.h> with vsprintf() didn't work */
@@ -88,7 +88,7 @@ main(void)
#endif
exit(1);
}
try("%s %d%s%d%d %s",
try_this("%s %d%s%d%d %s",
"Landon Noll 1st proved that", 2, "^", 23209, -1, "was prime");
if (strcmp(buf,
"Landon Noll 1st proved that 2^23209-1 was prime") != 0) {

View File

@@ -52,7 +52,7 @@ char buf[BUFSIZ];
#include <varargs.h>
void
try(char *fmt, ...)
try_this(char *fmt, ...)
{
va_list ap;
@@ -68,7 +68,7 @@ try(char *fmt, ...)
#else
void
try(char *a, int b, char *c, int d)
try_this(char *a, int b, char *c, int d)
{
return;
}
@@ -87,7 +87,7 @@ main(void)
/*
* test variable args and vsprintf/sprintf
*/
try("@%d:%s:%d@", 1, "hi", 2);
try_this("@%d:%s:%d@", 1, "hi", 2);
if (strcmp(buf, "@1:hi:2@") != 0) {
#if !defined(DONT_HAVE_VSPRINTF)
/* <varargs.h> with vsprintf() didn't work */
@@ -96,7 +96,7 @@ main(void)
#endif
exit(1);
}
try("%s %d%s%d%d %s",
try_this("%s %d%s%d%d %s",
"Landon Noll 1st proved that", 2, "^", 23209, -1, "was prime");
if (strcmp(buf,
"Landon Noll 1st proved that 2^23209-1 was prime") != 0) {

View File

@@ -235,3 +235,9 @@ Needed enhancements
One should make available a the fundimental math operations
on ZVALUE, NUMBER and perhaps COMPLEX (without all of the
other stuff) in a separate library.
* improve the coverage in the 'SEE ALSO' help file lists
* where reasonable, be sure that regress.cal tests builtin functions
* clean the source code and document it better

1
hist.h
View File

@@ -30,7 +30,6 @@
#define DOTCHAR '.' /* char which indicates current directory */
#define PATHCHAR '/' /* char which separates path components */
#define LISTCHAR ':' /* char which separates paths in a list */
#define PATHSIZE 1024 /* maximum length of path name */
/*

78
input.c
View File

@@ -91,20 +91,51 @@ opensearchfile(char *name, char *pathlist, char *extension, int rd_once)
{
int i;
char *cp;
char path[PATHSIZE+1]; /* name being searched for */
char *path; /* name being searched for */
struct stat statbuf; /* stat of the path */
unsigned long namelen; /* length of name */
unsigned long extlen; /* length of the extension if non-NULL or 0 */
unsigned long pathlen; /* length of the pathlist if non-NULL or 0 */
/*
* allocate storage for the longest name being searched for
*
* We will allocate more than we could ever want/need.
* The longest we could ever need would be:
*
* pathlist (as a single long string)
* /
* name
* .
* extenstion
* \0
* guard byte
*/
namelen = strlen(name);
if (extension != NULL) {
extlen = strlen(extension);
} else {
extlen = 0;
}
if (pathlist != NULL) {
pathlen = strlen(pathlist);
} else {
pathlen = 0;
}
path = malloc(pathlen+1 + 1 + namelen+1 + extlen+1 + 1 + 1);
if (path == NULL) {
math_error("Cannot allocate filename path buffer");
/*NOTREACHED*/
}
/*
* Don't try the extension if the filename already contains it.
*/
if (extension) {
unsigned long namelen = strlen(name);
unsigned long extlen = strlen(extension);
if (namelen >= extlen &&
strcmp(&name[namelen-extlen], extension) == 0)
extension = NULL;
if (extension != NULL && namelen >= extlen &&
strcmp(&name[namelen-extlen], extension) == 0) {
extension = NULL;
}
/*
* If the name is absolute, or if there is no path list, then
* make one which just searches for the name straight. Then
@@ -127,21 +158,26 @@ opensearchfile(char *name, char *pathlist, char *extension, int rd_once)
*cp++ = PATHCHAR;
strcpy(cp, name);
i = openfile(path);
if ((i < 0) && (extension != NULL && extension[0] != '\0')) {
if ((i < 0) &&
(extension != NULL && extension[0] != '\0')) {
strcat(path, extension);
i = openfile(path);
}
} while ((i < 0) && *pathlist);
/* examine what our search produced */
if (i < 0)
if (i < 0) {
free(path);
return i;
}
if (cip->i_fp == NULL) {
/* cannot find a file to open */
free(path);
return -3;
}
if (fstat(fileno(cip->i_fp), &statbuf) < 0) {
/* unable to fstat the open file */
free(path);
return -4;
}
@@ -149,6 +185,7 @@ opensearchfile(char *name, char *pathlist, char *extension, int rd_once)
if (rd_once == TRUE && isinoderead(&statbuf) >= 0) {
/* file is in readset and reopen is false */
closeinput();
free(path);
return 1;
}
@@ -156,10 +193,12 @@ opensearchfile(char *name, char *pathlist, char *extension, int rd_once)
if (addreadset(name, path, &statbuf) < 0) {
/* cannot add to readset */
closeinput();
free(path);
return -1;
}
/* file was added to/updated in readset */
free(path);
return 0;
}
@@ -189,7 +228,7 @@ homeexpand(char *name)
char *home2; /* fullpath of the home directory */
char *fullpath; /* the malloced expanded path */
char *after; /* after the ~user or ~ */
char username[PATHSIZE+1]; /* extratced username */
char *username; /* extratced username */
/* firewall */
if (name[0] != HOMECHAR)
@@ -217,11 +256,15 @@ homeexpand(char *name)
}
/* just malloc the home directory and return it */
fullpath = (char *)malloc(strlen(ent->pw_dir)+1);
if (fullpath == NULL) {
return NULL;
}
strcpy(fullpath, ent->pw_dir);
return fullpath;
}
if (after-name > PATHSIZE+1) {
/* username is too big */
username = (char *) malloc(after-name + 1 + 1);
if (username == NULL) {
/* failed to malloc username */
return NULL;
}
strncpy(username, name+1, after-name-1);
@@ -229,6 +272,7 @@ homeexpand(char *name)
/* get that user's home directory */
ent = (struct passwd *)getpwnam(username);
free(username);
if (ent == NULL) {
/* unknown user */
return NULL;
@@ -241,6 +285,9 @@ homeexpand(char *name)
* build the fullpath given the home directory
*/
fullpath = (char *)malloc(strlen(home2)+strlen(after)+1);
if (fullpath == NULL) {
return NULL;
}
sprintf(fullpath, "%s%s", home2, after);
return fullpath;
}
@@ -333,6 +380,9 @@ openfile(char *name)
cip->i_fp = fp;
cip->i_line = 1;
cip->i_name = (char *)malloc(strlen(name) + 1);
if (cip->i_name == NULL) {
return -1;
}
strcpy(cip->i_name, name);
return 0;
}
@@ -658,7 +708,7 @@ reread(void)
void
runrcfiles(void)
{
char path[PATHSIZE+1]; /* name being searched for */
char path[MAX_CALCRC+1+1]; /* name being searched for */
char *cp;
char *newcp;
char *p;

View File

@@ -702,6 +702,8 @@ define test_functions()
{
local a, b;
local pi;
local h, n, r, m, v;
local n2, m2, v2;
print '700: Beginning test_functions';
@@ -1163,7 +1165,31 @@ define test_functions()
vrfy(popcnt(pi(1e-20),0) == 69, '1101: popcnt(pi(1e-20),0) == 69');
vrfy(popcnt(17^777) == 1593, '1102: popcnt(17^777) == 1593');
print '1103: Ending test_functions';
/*
* more hnrmod testing
*/
vrfy(hnrmod(21<<100+5,3,100,1) == (21<<100+5)%(3<<100+1),
'1103: hnrmod(21<<100+5,3,100,1) == (21<<100+5)%(3<<100+1)');
vrfy(hnrmod(21<<500+7,3,500,1) == (21<<500+7)%(3<<500+1),
'1104: hnrmod(21<<500+7,3,500,1) == (21<<500+7)%(3<<500+1)');
vrfy(hnrmod(-767256,84,1,0) == (-767256)%(84<<1+0),
'1105: hnrmod(-767256,84,1,0) == (-767256)%(84<<1+0)');
vrfy(hnrmod(-831150,75,1,0) == (-831150)%(75<<1+0),
'1106: hnrmod(-831150,75,1,0) == (-831150)%(75<<1+0)');
vrfy(hnrmod(-767256,84,1,1) == (-767256)%(84<<1+1),
'1107: hnrmod(-767256,84,1,1) == (-767256)%(84<<1+1)');
vrfy(hnrmod(-831150,75,1,1) == (-831150)%(75<<1+1),
'1108: hnrmod(-831150,75,1,1) == (-831150)%(75<<1+1)');
vrfy(hnrmod(-767256,84,1,-1) == (-767256)%(84<<1-1),
'1109: hnrmod(-767256,84,1,-1) == (-767256)%(84<<1-1)');
vrfy(hnrmod(-831150,75,1,-1) == (-831150)%(75<<1-1),
'1110: hnrmod(-831150,75,1,-1) == (-831150)%(75<<1-1)');
vrfy(hnrmod(21<<100+5,3,100,0) == (21<<100+5)%(3<<100+0),
'1111: hnrmod(21<<100+5,3,100,0) == (21<<100+5)%(3<<100+0)');
vrfy(hnrmod(21<<500+7,3,500,-1) == (21<<500+7)%(3<<500-1),
'1112: hnrmod(21<<500+7,3,500,-1) == (21<<500+7)%(3<<500-1)');
print '1113: Ending test_functions';
}
print '017: parsed test_functions()';

View File

@@ -65,7 +65,7 @@ BOOL inputwait; /* TRUE if in a terminal input wait */
jmp_buf jmpbuf; /* for errors */
int start_done = FALSE; /* TRUE => start up processing finished */
char *program = "calc"; /* our name */
char cmdbuf[MAXCMD+1]; /* command line expression */
char cmdbuf[MAXCMD+1+1]; /* command line expression */
/*
@@ -281,6 +281,11 @@ initenv(void)
if (calcrc == NULL) {
calcrc = DEFAULTCALCRC;
}
if (strlen(calcrc) > MAX_CALCRC) {
math_error("The $CALCRC variable is longer than %d chars",
MAX_CALCRC);
/*NOTREACHED*/
}
/* determine the $CALCBINDINGS value */
calcbindings = (no_env ? NULL : getenv(CALCBINDINGS));

View File

@@ -206,7 +206,7 @@ main(int argc, char **argv)
'/', "* signed 64 bits *", '/');
putchar('\n');
printf("%c%s%c\n", '/',"* how to form 64 bit constants *",'/');
#if defined(__STDC__) && __STDC__ != 0
#if defined(FORCE_STDC) || (defined(__STDC__) && __STDC__ != 0) || defined(__cplusplus)
printf("#define U(x) x ## ULL\n");
printf("#define L(x) x ## LL\n");
#else
@@ -322,7 +322,7 @@ main(int argc, char **argv)
'/', "* signed 64 bits *", '/');
putchar('\n');
printf("%c%s%c\n", '/',"* how to form 64 bit constants *",'/');
#if defined(__STDC__) && __STDC__ != 0
#if defined(FORCE_STDC) || (defined(__STDC__) && __STDC__ != 0) || defined(__cplusplus)
printf("#define U(x) x ## UL\n");
printf("#define L(x) x ## L\n");
#else
@@ -341,7 +341,7 @@ main(int argc, char **argv)
'/', "* signed 64 bits *", '/');
putchar('\n');
printf("%c%s%c\n", '/',"* how to form 64 bit constants *",'/');
#if defined(__STDC__) && __STDC__ != 0
#if defined(FORCE_STDC) || (defined(__STDC__) && __STDC__ != 0) || defined(__cplusplus)
printf("#define U(x) x ## ULL\n");
printf("#define L(x) x ## LL\n");
#else

14
qmath.c
View File

@@ -21,8 +21,6 @@ NUMBER _qnegone_ = { { _oneval_, 1, 1 }, { _oneval_, 1, 0 }, 1, NULL };
NUMBER _qonehalf_ = { { _oneval_, 1, 0 }, { _twoval_, 1, 0 }, 1, NULL };
NUMBER _qonesqbase_ = { { _oneval_, 1, 0 }, { _sqbaseval_, 2, 0 }, 1, NULL };
#define INITCONSTCOUNT 8
NUMBER * initnumbs[INITCONSTCOUNT] = {&_qzero_, &_qone_, &_qtwo_, &_qthree_,
&_qfour_, &_qten_, &_qnegone_, &_qonehalf_};
@@ -1365,8 +1363,8 @@ qcmpi(NUMBER *q, long n)
#define NNALLOC 1000
static NUMBER *freeNum;
static NUMBER **firstNums;
static NUMBER *freeNum = NULL;
static NUMBER **firstNums = NULL;
static long blockcount = 0;
@@ -1389,8 +1387,12 @@ qalloc(void)
temp->links = 0;
}
blockcount++;
newfn = (NUMBER **)
realloc(firstNums, blockcount * sizeof(NUMBER *));
if (firstNums == NULL) {
newfn = (NUMBER **) malloc(blockcount * sizeof(NUMBER *));
} else {
newfn = (NUMBER **)
realloc(firstNums, blockcount * sizeof(NUMBER *));
}
if (newfn == NULL) {
math_error("Cannot allocate new number block");
/*NOTREACHED*/

View File

@@ -13,6 +13,8 @@
#include "zmath.h"
#define INITCONSTCOUNT 8 /* number of initnumbs[] pre-defined constants */
/*
* Rational arithmetic definitions.

2
shs.c
View File

@@ -598,7 +598,7 @@ shs_final_state(HASH *state)
/*
* load ZVALUE
*/
#if BASEB == 16
#if BASEB == 16 && CALC_BYTE_ORDER == LITTLE_ENDIAN
for (i=0; i < ret.len; i+=2) {
ret.v[ret.len-i-1] = ((HALF*)dig->digest)[i+1];
ret.v[ret.len-i-2] = ((HALF*)dig->digest)[i];

2
shs1.c
View File

@@ -574,7 +574,7 @@ shs1_final_state(HASH *state)
/*
* load ZVALUE
*/
#if BASEB == 16
#if BASEB == 16 && CALC_BYTE_ORDER == LITTLE_ENDIAN
for (i=0; i < ret.len; i+=2) {
ret.v[ret.len-i-1] = ((HALF*)dig->digest)[i+1];
ret.v[ret.len-i-2] = ((HALF*)dig->digest)[i];

View File

@@ -960,8 +960,8 @@ stringrsearch(STRING *s1, STRING *s2, long start, long end, ZVALUE *index)
#define STRALLOC 100
static STRING *freeStr;
static STRING **firstStrs;
static STRING *freeStr = NULL;
static STRING **firstStrs = NULL;
static long blockcount = 0;
@@ -984,8 +984,12 @@ stralloc(void)
temp->s_links = 0;
}
blockcount++;
newfn = (STRING **)
realloc(firstStrs, blockcount * sizeof(STRING *));
if (firstStrs == NULL) {
newfn = (STRING **) malloc( blockcount * sizeof(STRING *));
} else {
newfn = (STRING **)
realloc(firstStrs, blockcount * sizeof(STRING *));
}
if (newfn == NULL) {
math_error("Cannot allocate new string block");
/*NOTREACHED*/

View File

@@ -12,7 +12,7 @@
#define MAJOR_VER 2 /* major version */
#define MINOR_VER 11 /* minor version */
#define MAJOR_PATCH 0 /* patch level or 0 if no patch */
#define MINOR_PATCH "2" /* test number or empty string if no patch */
#define MINOR_PATCH "4" /* test number or empty string if no patch */
/*
* calc version constants
@@ -35,7 +35,7 @@ static char *stored_version = NULL; /* version formed if != NULL */
* This function returns a malloced version string. This version
* string does not contain the title, just:
*
* x.y.ztsomething
* x.y.ztw
* x.y.z
* x.y
*/
@@ -56,17 +56,17 @@ version(void)
*/
if (sizeof(MINOR_PATCH) > 1) {
sprintf(verbuf,
"%d.%d.%d%s", calc_major_ver, calc_minor_ver,
"%d.%d.%dt%s", calc_major_ver, calc_minor_ver,
calc_major_patch, calc_minor_patch);
} else if (MAJOR_PATCH > 0) {
sprintf(verbuf,
"%d.%d.%s", calc_major_ver, calc_minor_ver, calc_minor_patch);
"%d.%d.%d", calc_major_ver, calc_minor_ver, calc_major_patch);
} else {
sprintf(verbuf, "%d.%d", calc_major_ver, calc_minor_ver);
}
/*
* same the versions string into a newly malloced buffer
* save the versions string into a newly malloced buffer
*/
stored_version = (char *)malloc(strlen(verbuf)+1);
if (stored_version == NULL) {

36
zfunc.c
View File

@@ -1683,7 +1683,7 @@ done: if (s == 0) {
void
zroot(ZVALUE z1, ZVALUE z2, ZVALUE *dest)
{
ZVALUE try, quo, old, temp, temp2;
ZVALUE ztry, quo, old, temp, temp2;
ZVALUE k1; /* holds k - 1 */
int sign;
long i;
@@ -1733,12 +1733,12 @@ zroot(ZVALUE z1, ZVALUE z2, ZVALUE *dest)
* is unimportant.
*/
highbit = (highbit + k - 1) / k;
try.len = (highbit / BASEB) + 1;
try.v = alloc(try.len);
zclearval(try);
try.v[try.len-1] = ((HALF)1 << (highbit % BASEB));
try.sign = 0;
old.v = alloc(try.len);
ztry.len = (highbit / BASEB) + 1;
ztry.v = alloc(ztry.len);
zclearval(ztry);
ztry.v[ztry.len-1] = ((HALF)1 << (highbit % BASEB));
ztry.sign = 0;
old.v = alloc(ztry.len);
old.len = 1;
zclearval(old);
old.sign = 0;
@@ -1746,10 +1746,10 @@ zroot(ZVALUE z1, ZVALUE z2, ZVALUE *dest)
* Main divide and average loop
*/
for (;;) {
zpowi(try, k1, &temp);
zpowi(ztry, k1, &temp);
zquo(z1, temp, &quo, 0);
zfree(temp);
i = zrel(try, quo);
i = zrel(ztry, quo);
if (i <= 0) {
/*
* Current try is less than or equal to the root since it is
@@ -1758,24 +1758,24 @@ zroot(ZVALUE z1, ZVALUE z2, ZVALUE *dest)
* we are done since no improvement occurred.
* If not, save the improved value and loop some more.
*/
if ((i == 0) || (zcmp(old, try) == 0)) {
if ((i == 0) || (zcmp(old, ztry) == 0)) {
zfree(quo);
zfree(old);
try.sign = (BOOL)((HALF)sign);
zquicktrim(try);
*dest = try;
ztry.sign = (BOOL)((HALF)sign);
zquicktrim(ztry);
*dest = ztry;
return;
}
old.len = try.len;
zcopyval(try, old);
old.len = ztry.len;
zcopyval(ztry, old);
}
/* average current try and quotent for the new try */
zmul(try, k1, &temp);
zfree(try);
zmul(ztry, k1, &temp);
zfree(ztry);
zadd(quo, temp, &temp2);
zfree(temp);
zfree(quo);
zquo(temp2, z2, &try, 0);
zquo(temp2, z2, &ztry, 0);
zfree(temp2);
}
}

View File

@@ -1117,7 +1117,7 @@ zmodi(ZVALUE z, long n)
val = 0;
while (len--)
val = ((val << BASEB) + ((FULL) *h1--)) % n;
if (z.sign)
if (val && z.sign)
val = n - val;
return (long)val;
}

View File

@@ -577,7 +577,7 @@ typedef struct {
*/
#if FULL_BITS == 64
# if defined(__STDC__) && __STDC__ != 0
# if defined(FORCE_STDC) || (defined(__STDC__) && __STDC__ != 0) || defined(__cplusplus)
# define HVAL(a,b) (HALF)(0x ## a ## b)
# else
# define HVAL(a,b) (HALF)(0x/**/a/**/b)
@@ -585,7 +585,7 @@ typedef struct {
#elif 2*FULL_BITS == 64
# if defined(__STDC__) && __STDC__ != 0
# if defined(FORCE_STDC) || (defined(__STDC__) && __STDC__ != 0) || defined(__cplusplus)
# define HVAL(a,b) (HALF)0x##b, (HALF)0x##a
# else
/* NOTE: Due to a SunOS cc bug, don't put spaces in the HVAL call! */

2
zmod.c
View File

@@ -2337,7 +2337,7 @@ zhnrmod(ZVALUE v, ZVALUE zh, ZVALUE zn, ZVALUE zr, ZVALUE *res)
/* ret = ((ret < 0) ? ret+modlus : ((ret == modulus) ? 0 : ret)); */
if (ret.sign) {
zadd(ret, modulus, &t);
zfree(t);
zfree(ret);
ret = t;
} else if (zrelval == 0) {
zfree(ret);

View File

@@ -142,7 +142,7 @@
#if FULL_BITS == SBITS
# define SLEN 1 /* a 64 bit slot can be held in a FULL */
# if defined(__STDC__) && __STDC__ != 0
# if defined(FORCE_STDC) || (defined(__STDC__) && __STDC__ != 0) || defined(__cplusplus)
# define SVAL(a,b) (FULL)U(0x ## a ## b)
# define SHVAL(a,b,c,d) (HALF)0x ## c ## d, (HALF)0x ## a ## b
# else
@@ -166,7 +166,7 @@
#elif 2*FULL_BITS == SBITS
# define SLEN 2 /* a 64 bit slot needs 2 FULLs */
# if defined(__STDC__) && __STDC__ != 0
# if defined(FORCE_STDC) || (defined(__STDC__) && __STDC__ != 0) || defined(__cplusplus)
# define SVAL(a,b) (FULL)0x##b, (FULL)0x##a
# define SHVAL(a,b,c,d) (HALF)0x##d, (HALF)0x##c, \
(HALF)0x##b, (HALF)0x##a