Subscribe to particular warning flags, not to -Wextra.

-Wextra now entails -Wcast-function-type, which is incompatible with
various Unix APIs like sigaction.

Fix reader NaN notation: -nan.0/+nan.0, not just nan.0.

Some trivial tests for read/write invariance.

One xfail: The `nan.0' notation reads as NaN, not as a symbol, and
`+nan.0' is not recognized.

More atan2 edge cases.

Fix edge cases of flonum divide and atan2.

Test edge cases of flonum-divide primitive and of atan.

Let the system math library elicit IEEE 754 exceptions.

Don't trap in our ucode primitive wrappers if the floating-point
exception isn't trapped.

Test the flonum-sqrt ucode primitive explicitly.

The compiler open-codes it in flo:sqrt on amd64, which hardware gives
the correct behaviour, but the ucode does not right now.

Follow IEEE 754 semantics: disable floating-point traps by default.

Add tests for floating-point exceptions.

Mind my i's and j\s.

I broke random-bytevector for large inputs.

Rewrite random number generator.

New one has 32-byte state s, produces output x by splitting 64-byte
ChaCha20_s(0) into 32-byte halves s' and x and replacing the state by
s'.

I added two alternate implementations of flo:random-unit, one which
samples real numbers uniformly from [0,1] and rounds them to
floating-point numbers, and one which samples real numbers uniformly
from (2^{-emin - p - 1}, 1 - ulp(1)/4) and rounds them to
floating-point numbers.  The latter is wrong for various reasons but
it is closer to what we historically provided, so it is what we use
for now.

I removed the fallback in case get-entropy (/dev/urandom) fails,
which means this won't work on Windows until someone teaches the
microcode to call CryptGenRandom there, and won't work in a chroot
unless someone teaches it to use getentropy(2) or getrandom(2) or
whatever.

If this causes any problems, feel free to back out this commit --
aside from refusing to fall back to getting `entropy' from the clock
and having a different export format, this is intended to be a
drop-in replacement for the old random.scm (hence the recent tests),
so if I made a mistake just back it out, let me know what went wrong,
and I'll add more tests before re-merging it.

Test that make-random-source is deterministic.

Randomize the dispatch tag cache; don't copy the system state.

Use get-entropy primitive to initialize random states.

Remember to close /dev/urandom when done.

Check for fenv.h and ieeefp.h functions only if the header is there.

Makes it easier to disable these and test our own version by just
setting ac_cv_header_fenv_h=no ac_cv_header_ieeefp_h=no.

Add -Werror after all autoconf checks.

If we do it before, the autoconf checks -- which are full of compiler
warnings -- all fail.

Need "osentropy.h" for OS_get_entropy.

Use "config.h" to get uint8_t.  No more need of <stddef.h> for size_t.

Fix prototypes for x87_read/write_environment.

Use a fresh random state for randomized tests.

More floating-point parameters.

New primitive (get-entropy bv) randomly fills a 32-byte vector.

Turn off -Werror here; breaks configure.

Should turn it back on somewhere else.

Missed a run-sub-test.

New procedure (bytevector-zero-explicit! <bv> [<start> [<end>]]).

Intended to zero the memory backing a bytevector even if the values
will never be used again and an aggressively optimizing compiler can
prove that.

Doesn't actually work (GC can move stuff without zeroing it) but it
may help to have it in order to tag where it _would_ be needed if we
could make it work in the future.

Fix indentation.

Make some of my stupid old randomized tests a little less slow.

Would be better to rewrite these so they're not as stupid.

Halfway fix some broken tests.

- Use with-test-properties, not run-sub-test.

- Disable run-sub-test in the test environment; it doesn't make sense
  for tests to call it.

- Add a rudimentary expect-failure.  (Feel free to spruce this up
  with a message detailing the nature of the expectation, a note in
  the final report, &c.)

- Partly fix test-char-set.scm, which was silently nonfunctional
  because of run-sub-test.  Leave an xfail for the part I don't know
  how to fix.  (Chris -- if you could fix this, that would be great!)

Depend on x11-const directly, not via a stamp and submake.

Going via the submake means parallel make can get confused and cause
two processes to try to build the same thing at the same time and
stomp on each other's toes.

Make the C compiler warnings that we subscribe to fatal.

Confirm the primitive respects the bounds too.

Simplify construction of chacha core test vector tests.

Test fenceposts in chacha core primitives.

Unfix the fix for an off-by-zero.

arg_ulong_index_integer(argno, n) takes n as _exclusive_ upper bound.

In this case, n - 64 is a valid starting index.  So the exclusive
upper bound is n - 64 + 1 = n - 63.

Fix more formatting.

Eliminate unused variable warnings.

Fix a bug in prchacha: compiler warned of unused variable.

Also there was an off-by-one thinko, and incorrect PRIMITIVE_HEADER arguments.
While I was at it I tweaked the formatting a little.

Fix sign of (flo:ulp -infinity).  Simplify.

flo:ulp is supposed to return a positive magnitude, in particular the
absolute distance from x to the next larger floating-point number in
magnitude.

Tidy up ieee754.scm, add some tests, and export some of it.

Allow sequences to contain zero elements.

Don't use upper-case symbols in loader.

The system might be in a non-case-fold mode.

Don't print internals of promises.

Fix printing of promises.

Merge branch 'release-10'

Make macos dmg file read only.

Fix problem with building and installing x11 plugin from native dist.

Comment what log-logistic and logit-exp are for.

Define (flo:ulp x) to be distance to next float in direction of x.

Thus, (flo:ulp 1.) = (flo:ulp -1.) = flo:ulp-of-one.  Bottoms out at
infinity: (flo:ulp x) = x if x is infinite.  This definition applies
to zero, too, giving the smallest subnormal.

Add Mode: C line on request from mejja.

(Not keen to otherwise reformat this code from an upstream, though.)

Clarify some comments.

Reorganize logsumexp for clarity.  No functional change.

Mask underflow exception in logsumexp because it doesn't matter.

Add (logsumexp (list 1 2 3)) = (log (+ (exp 1) (exp 2) (exp 3))).

Don't assume that it's possible to run MIT_SCHEME_EXE during configuration.

This won't work with a native distribution, since there may not be a working
scheme on the host.  Instead use other ways to do the necessary work.

Bump revision to 10.1.1 to fix stupid bug.

Reduce default heap size for 32-bit systems.

During build was getting sporadic "unable to mmap heap" errors.

Fix typo in previous change.

Eliminate flo:epsilon.

The name is confusing.  It means two different things: the distance
from 1 to the next positive floating-point number, and the maximum
relative error of rounding to the nearest floating-point number.
Provide the separate names flo:ulp-of-one and flo:error-bound for
these so that you're not even tempted to be confused by the name
`epsilon'.

(But continue to say `eps' in the proofs in the comments where it's
obvious that we're talking about relative error bounds and a short
name for that is essential for legibility.)

Tidy up error analysis of logit and logistic.

Add (logit1/2+ p) = (logit (+ 1/2 p)) and (logistic-1/2 x) = (-
(logistic x) 1/2), for reasons like log1p and expm1.  Add some
trivial tests of the edge cases they cover where logit and logistic
are ill-conditioned.

Merge branch 'release-10'

No need to separately configure plugins.

Change "module" to "plugin".

Eliminate compiler warning.

Use create-dmg script to automate building of macOS image.

Fix a couple of compiler warnings from 32-bit system.

Merge branch 'release-10' of git.sv.gnu.org:/srv/git/mit-scheme into release-10

Fix thinko.

Eliminate unavailable (and unused) load option.

Forgot to define this.  Ummm, I guess cref accidentally made it work?

Expose flo:radix = 2 too, just for good measure, and for assertions.

Export primitive flonum-denormalize as flo:ldexp and flo:scalbn.

These are the familiar -- and when FLT_RADIX = 2, equivalent -- names
for the operation in libm.

`Denormalize' is not right: neither does it turn a normal number into
a subnormal, nor does it put anything into a canonical interval like
[1/2, 1) or [1, 2).  I personally favour `shift' -- it is, in Z, the
same as the bitwise shift operation.  But unless anyone else wants to
join me against this windmill I'll just stay with the familiar if
ugly libm names.

Disable errno checking in libm functions.

If there's a floating-point exception, we have a mechanism for
dealing with it -- and continuing without a trap if the user so
desires.

Leave in the restriction on flo:log1p and flo:expm1 because those are
implementation limitations in the i386 compiler; the full functions
are implemented in the runtime.

log1pexp(x) is never exactly x even if that's the best approximation.

Tabify.

Use named constants rather than magic numbers.

Define and export some more floating-point parameters.

Fix typo.

Update make-doc-files to work with plugins.

test-driver is autocrud too.

Change way that plugins are registered so that lib and doc are separate.

This is needed when building just the documentation, as for the release.

Simple-minded criterion to compute x^-y by 1/x^y was too simple.

For any |x| > 1, we can't do this without comparing the size of
log(x) and y well enough to discern what might yield subnormal, so
might as well just compute the general case then.

Omit unnecessary case for logistic.  Add a test for this case.

The relative error is small in naive evaluation of 1/(1 + e^{-x})
even if x is very negative.

Fix potential shell portability issues.

Extend expt to generate subnormals with exact integer exponents.

Expose the minimum and maximum finite normal floating-point exponents.

Add some .gitignore files so I can use git status.

New primitives for ChaCha stream cipher core.

Should maybe provide a version that produces many blocks of output,
and perhaps xor's them into an existing message, but this will do for
my needs now.  Plan to use this for a PRNG.

Fix unparser->print method fallout.

Avoid using uninitialized stack garbage.

clock() returns (clock_t)-1 on error, not necessarily negative.

If clock_t is unsigned, n < 0 is always false.  Compare for equality
with (clock_t)-1.

Add plugin docs to doc files.

More version numbers.

Remove no-longer-needed args to configure.

More version numbers.

Update some version numbers.

No need to separately configure plugins.

Change "module" to "plugin".