Login

Subversion Repositories NedoOS

Rev

Blame | Last modification | View Log | Download | RSS feed

# extended
This is an 80-bit float format with a 1-bit sign, 15-bit exponent and 64-bit
significand.

# Format
`extended` floats are stored in little endian. The "most significant bit" is
sign, the next 15 bits are exponent, and the next 64 bits encode the 64-bit
significand (note that the top bit of the significand **is** stored explicitly).
The exponent is stored with a bias of +16384 (so an exponent of 0 is stored as
the literal value 16384). Like the so-called "single" format this library uses,
special values +0 ,-0, +inf, -inf, and NaN are stored with an exponent of -16384
(literal stored as a 0), and the sign and top two bits of the significand
determine which special value it is.

```
m is significand
e is exponent
s is sign
- is any value (0 or 1, doesn't matter)

        seeeeeee eeeeeeee mmmmmmmm mmmmmmmm mmmmmmmm mmmmmmmm mmmmmmmm mmmmmmmm mmmmmmmm mmmmmmmm
+0      00000000 00000000 00------ -------- -------- -------- -------- -------- -------- --------
-0      10000000 00000000 00------ -------- -------- -------- -------- -------- -------- --------
+inf    00000000 00000000 1------- -------- -------- -------- -------- -------- -------- --------
+inf    10000000 00000000 1------- -------- -------- -------- -------- -------- -------- --------
NaN     -0000000 00000000 01------ -------- -------- -------- -------- -------- -------- --------
 1      01000000 00000000 00000000 00000000 00000000 00000000 00000000 00000000 00000000 00000000
 2      01000000 00000001 00000000 00000000 00000000 00000000 00000000 00000000 00000000 00000000
 1      11000000 00000000 00000000 00000000 00000000 00000000 00000000 00000000 00000000 00000000
pi      01000000 00000001 11001001 00001111 11011010 10100010 00100001 01101000 11000010 00110101

```

## Input/Output
```
  x: HL points to the first operand (if any)
  y: DE points to the second operand (if any)
  t: IX points to the third operand (if any)
  z: BC points to the output
```
There is one notable exception:

* `xcmp` takes float inputs, but returns output in the zero flag and carry flag

## Routines
These float routines are for the most part prefixed with `x` to distinguish
them from the other float routines in this library. Most of them are in their
own file, but some files contain several routines. These cases are noted in the
chart below.

```
xabs       |x|
xacos      acos(x)
xacosh     acosh(x)
xadd       x+y
xamean     (x+y)/2
xasin      asin(x)
xasinh     asinh(x)
xatan      atan(x)
xatanh     atanh(x)
xbg        1/BG(x,y)  (reciprocal Borchardt-Gauss mean)
xcis       {cos(x), sin(x)}
xcmp       compare x to y
xcos       cos(x)
xcosh      cosh(x)
xdiv       x/y
xdiv2      x/2
xexp       e^x
xfma       x*y+t
xgeomean   sqrt(x*y)
xinv       1/x
xlg        log2(x)
xln        log(x)     (a.k.a. ln(x))
xlog       log_y(x)
xlog10     log10(x)
xmod1      x mod 1
xmul       x*y
xmul2      x*2  (found in the extended/mul directory)
xmul3      x*3  (found in the extended/mul directory)
xmul5      x*5  (found in the extended/mul directory)
xmul7      x*7  (found in the extended/mul directory)
xmul10     x*10 (found in the extended/mul directory)
xmul11     x*11 (found in the extended/mul directory)
xmul13     x*13 (found in the extended/mul directory)
xmul15     x*15 (found in the extended/mul directory)
xmul17     x*17 (found in the extended/mul directory)
xmul31     x*31 (found in the extended/mul directory)
xneg       -x
xpow       x^y
xpow2      2^x
xpow10     10^x
xrand      random number on [0,1), uniform distribution
xrsub      -x+y
xsin       sin(x)
xsinh      sinh(x)
xsqrt      sqrt(x)
xsub       x-y
xtan       tan(x)
xtanh      tanh(x)
strtox     "string" ==> extended float
xtostr     string(x)
xtoTI      TIFloat(x)
TItox      TIFloat ==> extended
```

## Example Usage
Calculate `10^(pi+e)`:

```
  ; Add pi+e
  ld hl, const_pi
  ld de, const_e
  ld bc, xOP1
  call xadd

  ; BC points to the output, so want to load BC into HL, then call pow10
  ld h, b
  ld l, c
  call xpow10
  ; result is at BC (which is the same as HL and xOP1 here)
```