mirror of https://github.com/ksherlock/mpw.git
381 lines
5.5 KiB
C
381 lines
5.5 KiB
C
#include <Math.h>
|
|
#include <SANE.h>
|
|
|
|
#include <stdio.h>
|
|
|
|
//extern pascal extended NAN(short x);
|
|
|
|
void dump_decimal(const decimal *d)
|
|
{
|
|
fprintf(stdout, "%d : %d : %.*s\n",
|
|
d->sgn, d->exp, d->sig.length, d->sig.text);
|
|
}
|
|
|
|
|
|
void test_inf(void)
|
|
{
|
|
|
|
decimal d;
|
|
decform df;
|
|
char buffer[80+1];
|
|
|
|
printf("\nInfinity\n");
|
|
|
|
df.style = 1;
|
|
df.digits = 10;
|
|
num2dec(&df,-inf(), &d);
|
|
dump_decimal(&d);
|
|
//
|
|
dec2str(&df, &d, buffer);
|
|
fprintf(stdout, "%s\n", buffer);
|
|
|
|
df.style = 0;
|
|
df.digits = 10;
|
|
num2dec(&df, -inf(), &d);
|
|
dump_decimal(&d);
|
|
//
|
|
dec2str(&df, &d, buffer);
|
|
fprintf(stdout, "%s\n", buffer);
|
|
|
|
df.style = 1;
|
|
df.digits = 0;
|
|
num2dec(&df, inf(), &d);
|
|
dump_decimal(&d);
|
|
//
|
|
dec2str(&df, &d, buffer);
|
|
fprintf(stdout, "%s\n", buffer);
|
|
|
|
df.style = 0;
|
|
df.digits = 0;
|
|
num2dec(&df, inf(), &d);
|
|
dump_decimal(&d);
|
|
//
|
|
dec2str(&df, &d, buffer);
|
|
fprintf(stdout, "%s\n", buffer);
|
|
|
|
}
|
|
|
|
void test_nan(void)
|
|
{
|
|
|
|
|
|
// num2dec generates N[16-digit]
|
|
// str2dec generates NAN(001) -- NAN(255)
|
|
|
|
decimal d;
|
|
decform df;
|
|
char buffer[80+1];
|
|
|
|
printf("\nNaN\n");
|
|
|
|
df.style = 1;
|
|
df.digits = 0;
|
|
num2dec(&df,nan(1), &d);
|
|
dump_decimal(&d);
|
|
// 7 : N4001000000000000
|
|
dec2str(&df, &d, buffer);
|
|
fprintf(stdout, "%s\n", buffer);
|
|
|
|
df.style = 0;
|
|
df.digits = 0;
|
|
num2dec(&df, nan(2), &d);
|
|
dump_decimal(&d);
|
|
//
|
|
dec2str(&df, &d, buffer);
|
|
fprintf(stdout, "%s\n", buffer);
|
|
|
|
|
|
df.style = 1;
|
|
df.digits = 10;
|
|
num2dec(&df,-nan(3), &d);
|
|
dump_decimal(&d);
|
|
//
|
|
dec2str(&df, &d, buffer);
|
|
fprintf(stdout, "%s\n", buffer);
|
|
|
|
df.style = 0;
|
|
df.digits = 10;
|
|
num2dec(&df, nan(4), &d);
|
|
dump_decimal(&d);
|
|
//
|
|
dec2str(&df, &d, buffer);
|
|
fprintf(stdout, "%s\n", buffer);
|
|
|
|
df.style = 0;
|
|
df.digits = 10;
|
|
num2dec(&df, nan(255), &d);
|
|
dump_decimal(&d);
|
|
// 7 : N40FF000000000000
|
|
|
|
dec2str(&df, &d, buffer);
|
|
fprintf(stdout, "%s\n", buffer);
|
|
}
|
|
|
|
void test_fxc2dec(void)
|
|
{
|
|
decimal d;
|
|
decform df;
|
|
|
|
df.style = 0;
|
|
df.digits = 100;
|
|
num2dec(&df, 1.125, &d);
|
|
dump_decimal(&d);
|
|
// 0 : -18 : 1125000000000000000
|
|
|
|
df.style = 0;
|
|
df.digits = 10;
|
|
num2dec(&df, 1.125, &d);
|
|
dump_decimal(&d);
|
|
// 0 : -9 : 1125000000
|
|
|
|
df.style = 0;
|
|
df.digits = 2;
|
|
num2dec(&df, 1.125, &d);
|
|
dump_decimal(&d);
|
|
// 0 : -1 : 11
|
|
|
|
|
|
|
|
df.style = 1;
|
|
df.digits = 100;
|
|
num2dec(&df, 1.125, &d);
|
|
dump_decimal(&d);
|
|
// 0 : -18 : 1125000000000000000
|
|
|
|
df.style = 1;
|
|
df.digits = 10;
|
|
num2dec(&df, 1.125, &d);
|
|
dump_decimal(&d);
|
|
// 0 : -10 : 11250000000
|
|
|
|
|
|
df.style = 1;
|
|
df.digits = 2;
|
|
num2dec(&df, 1.125, &d);
|
|
dump_decimal(&d);
|
|
//0 : -2 : 112
|
|
|
|
//
|
|
|
|
df.style = 0;
|
|
df.digits = 100;
|
|
num2dec(&df, 1e25, &d);
|
|
dump_decimal(&d);
|
|
// 0 : 7 : 1000000000000000000
|
|
|
|
df.style = 0;
|
|
df.digits = 10;
|
|
num2dec(&df, 1e25, &d);
|
|
dump_decimal(&d);
|
|
// 0 : 16 : 1000000000
|
|
|
|
df.style = 0;
|
|
df.digits = 2;
|
|
num2dec(&df, 1e25, &d);
|
|
dump_decimal(&d);
|
|
// 0 : 24 : 10
|
|
|
|
|
|
|
|
df.style = 1;
|
|
df.digits = 100;
|
|
num2dec(&df, 1e25, &d);
|
|
dump_decimal(&d);
|
|
// 0 : 7 : 1000000000000000000
|
|
|
|
df.style = 1;
|
|
df.digits = 10;
|
|
num2dec(&df, 1e25, &d);
|
|
dump_decimal(&d);
|
|
// 0 : 7 : 1000000000000000000
|
|
|
|
|
|
df.style = 1;
|
|
df.digits = 2;
|
|
num2dec(&df, 1e25, &d);
|
|
dump_decimal(&d);
|
|
//0 : 7 : 1000000000000000000
|
|
|
|
|
|
|
|
df.style = 1;
|
|
df.digits = 10;
|
|
num2dec(&df, 0, &d);
|
|
dump_decimal(&d);
|
|
//
|
|
|
|
df.style = 0;
|
|
df.digits = 10;
|
|
num2dec(&df, 0, &d);
|
|
dump_decimal(&d);
|
|
//
|
|
|
|
|
|
df.style = 1;
|
|
df.digits = 0;
|
|
num2dec(&df, -0.0, &d);
|
|
dump_decimal(&d);
|
|
//
|
|
|
|
df.style = 0;
|
|
df.digits = 0;
|
|
num2dec(&df, -0.0, &d);
|
|
dump_decimal(&d);
|
|
//
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
df.style = 0;
|
|
df.digits = 10;
|
|
num2dec(&df, 0.125, &d);
|
|
dump_decimal(&d);
|
|
//
|
|
|
|
df.style = 1;
|
|
df.digits = 10;
|
|
num2dec(&df, 0.125, &d);
|
|
dump_decimal(&d);
|
|
//
|
|
|
|
df.style = 0;
|
|
df.digits = 4;
|
|
num2dec(&df, 0.00000125, &d);
|
|
dump_decimal(&d);
|
|
// s/b -9 1250
|
|
|
|
df.style = 1;
|
|
df.digits = 4;
|
|
num2dec(&df, 0.00000125, &d);
|
|
dump_decimal(&d);
|
|
// s/b -4 0
|
|
|
|
df.style = 1;
|
|
df.digits = 6;
|
|
num2dec(&df, 0.0625, &d);
|
|
dump_decimal(&d);
|
|
// s/b -6 62500
|
|
|
|
df.style = 1;
|
|
df.digits = 6;
|
|
num2dec(&df, 1.125, &d);
|
|
dump_decimal(&d);
|
|
// s/b -6 1125000
|
|
|
|
df.style = 1;
|
|
df.digits = 6;
|
|
num2dec(&df, 1.0625, &d);
|
|
dump_decimal(&d);
|
|
// s/b -6 1062500
|
|
}
|
|
|
|
pascal void fp68k_3(void *, void *, unsigned short) = 0xA9EB;
|
|
|
|
void test_fx2l(void) {
|
|
|
|
long double x;
|
|
long int l;
|
|
int i;
|
|
|
|
static long double data[] = {
|
|
1.25,
|
|
1.5,
|
|
1.75,
|
|
2.25,
|
|
2.5,
|
|
2.75,
|
|
-1.25,
|
|
-1.5,
|
|
-1.75,
|
|
-2.25,
|
|
-2.5,
|
|
-2.75,
|
|
};
|
|
|
|
|
|
|
|
|
|
x = inf(); // 1.0 / 0.0;
|
|
fp68k_3(&x, &l, 0x2810);
|
|
printf("fx2l(inf) = %lx\n", l);
|
|
|
|
|
|
x = -inf(); // -1.0 / 0.0;
|
|
fp68k_3(&x, &l, 0x2810);
|
|
printf("fx2l(-inf) = %lx\n", l);
|
|
|
|
x = nan(1); // -1.0 / 0.0;
|
|
fp68k_3(&x, &l, 0x2810);
|
|
printf("fx2l(nan) = %lx\n", l);
|
|
|
|
|
|
x = 1e21;
|
|
fp68k_3(&x, &l, 0x2810);
|
|
printf("fx2l(1e21) = %lx\n", l);
|
|
|
|
x = -1e21;
|
|
fp68k_3(&x, &l, 0x2810);
|
|
printf("fx2l(-1e21) = %lx\n", l);
|
|
|
|
|
|
|
|
setround(UPWARD);
|
|
for (i = 0; i < 12; ++i) {
|
|
x = data[i];
|
|
fp68k_3(&x, &l, 0x2810);
|
|
printf("fx2l(%f) = %ld\n", x, l);
|
|
}
|
|
setround(DOWNWARD);
|
|
for (i = 0; i < 12; ++i) {
|
|
x = data[i];
|
|
fp68k_3(&x, &l, 0x2810);
|
|
printf("fx2l(%f) = %ld\n", x, l);
|
|
}
|
|
|
|
setround(TONEAREST);
|
|
for (i = 0; i < 12; ++i) {
|
|
x = data[i];
|
|
fp68k_3(&x, &l, 0x2810);
|
|
printf("fx2l(%f) = %ld\n", x, l);
|
|
}
|
|
setround(TOWARDZERO);
|
|
for (i = 0; i < 12; ++i) {
|
|
x = data[i];
|
|
fp68k_3(&x, &l, 0x2810);
|
|
printf("fx2l(%f) = %ld\n", x, l);
|
|
}
|
|
|
|
|
|
|
|
|
|
}
|
|
|
|
int main(int argc, char **argv)
|
|
{
|
|
//extended x;
|
|
//decimal d;
|
|
|
|
//short index;
|
|
//short valid;
|
|
|
|
(void)argc;
|
|
(void)argv;
|
|
|
|
//index = 0;
|
|
//valid = 0;
|
|
//str2dec("1.125", &index, &d, &valid);
|
|
//x = dec2num(&d);
|
|
|
|
|
|
test_inf();
|
|
test_nan();
|
|
|
|
test_fxc2dec();
|
|
test_fx2l();
|
|
|
|
return 0;
|
|
}
|