#ifdef FLOAT_ASM_MCS51
-// float __fsmul (float a, float b) reentrant
-static void dummy(void) _naked
+// float __fsmul (float a, float b) __reentrant
+static void dummy(void) __naked
{
- _asm
+ __asm
.globl ___fsmul
___fsmul:
// extract the two inputs, placing them into:
- // sign exponent mantiassa
- // ---- -------- ---------
+ // sign exponent mantissa
+ // ---- -------- --------
// a: sign_a exp_a r4/r3/r2
// b: sign_b exp_b r7/r6/r5
cpl sign_a
00003$:
- // add the exponents
+ // check if either input is infinity
+ mov a, exp_b
+ cjne a, #0xFF, 00004$
+ ljmp fs_return_inf
+00004$:
mov a, exp_a
+ cjne a, #0xFF, 00005$
+ ljmp fs_return_inf
+00005$:
+
+ // add the exponents
add a, exp_b
+ // if carry then no underflow
+ jc 00006$
add a, #130
- mov exp_a, a
+ jc 00007$
+ ljmp fs_return_zero
+
+00006$:
+ add a, #131
+ dec a
+ jnc 00007$
+ ljmp fs_return_inf
+00007$:
+ mov exp_a, a
// now we need to multipy r4/r3/r2 * r7/r6/r5
// ------------------------------------------
// r4 * r7 << 32
//
// This adds quite a bit of code, but it is a LOT faster
- // that three calls to __mululong...
+ // than three calls to __mululong...
// output goes into r4/r3/r2/r1/r0/xx
00010$:
ljmp fs_round_and_return
- _endasm;
+ __endasm;
}
#else
-
/*
** libgcc support for software floating point.
** Copyright (C) 1991 by Pipeline Associates, Inc. All rights reserved.
/* (c)2000/2001: hacked a little by johan.knol@iduna.nl for sdcc */
-
union float_long
{
float f;
volatile unsigned long result;
volatile int exp;
char sign;
-
+
fl1.f = a1;
fl2.f = a2;
result &= ~HIDDEN;
/* pack up and go home */
- fl1.l = PACK (sign ? SIGNBIT : 0 , (unsigned long)exp, result);
+ if (exp >= 0x100)
+ fl1.l = (sign ? SIGNBIT : 0) | __INFINITY;
+ else if (exp < 0)
+ fl1.l = 0;
+ else
+ fl1.l = PACK (sign ? SIGNBIT : 0 , exp, result);
return (fl1.f);
}
#endif
-