* as/link/lkar.h: sgetl and sputl are independent of endianness
[fw/sdcc] / device / lib / _fslt.c
index 1661a648659d98a4ae10b6bbd771140489bbcfb0..854d43da296281deb2560a180cd6723eaa828d86 100644 (file)
@@ -1,3 +1,75 @@
+/* Floating point library in optimized assembly for 8051
+ * Copyright (c) 2004, Paul Stoffregen, paul@pjrc.com
+ *
+ * This program is free software; you can redistribute it and/or
+ * modify it under the terms of the GNU Library General Public License
+ * as published by the Free Software Foundation; either version 2
+ * of the License, or (at your option) any later version.
+ *
+ * This library is distributed in the hope that it will be useful,
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
+ * GNU General Public License for more details.
+ *
+ * You should have received a copy of the GNU General Public License
+ * along with this program; if not, write to the Free Software
+ * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
+ */
+
+
+#define SDCC_FLOAT_LIB
+#include <float.h>
+
+
+#ifdef FLOAT_ASM_MCS51
+
+// char __fslt (float a, float b)
+static void dummy(void) __naked
+{
+       __asm
+       .globl  ___fslt
+___fslt:
+       mov     r7, a
+       mov     r0, sp
+       dec     r0
+       dec     r0
+       lcall   fs_check_negative_zeros
+       setb    sign_a
+       rlc     a
+       mov     a, @r0
+       jc      a_negative
+a_positive:
+       jnb     acc.7, ab_positive
+       // a is positive and b is negative, so a > b
+       mov     dpl, #0
+       ret
+a_negative:
+       jb      acc.7, ab_negative
+       // a is negative and b is positive, so a < b
+       mov     dpl, #1
+       ret
+ab_positive:
+       clr     sign_a
+ab_negative:
+       lcall   fs_compare_uint32
+       mov     a, r1
+       jnz     ab_different
+       // a and b are equal
+       mov     dpl, a
+       ret
+ab_different:
+       jb      sign_a, skip_invert
+       cpl     c
+skip_invert:
+       clr     a
+       rlc     a
+       mov     dpl, a
+       ret
+       __endasm;
+}
+
+#else
+
 /*
 ** libgcc support for software floating point.
 ** Copyright (C) 1991 by Pipeline Associates, Inc.  All rights reserved.
@@ -16,8 +88,6 @@
 
 /* (c)2000/2001: hacked a little by johan.knol@iduna.nl for sdcc */
 
-#include <float.h>
-
 union float_long
   {
     float f;
@@ -32,15 +102,15 @@ char __fslt (float a1, float a2)
   fl1.f = a1;
   fl2.f = a2;
 
-#if 0
-  if (fl1.l<0 && fl2.l<0)
-    {
-      fl1.l ^= SIGNBIT;
-      fl2.l ^= SIGNBIT;
-    }
-#endif
+  if (fl1.l<0 && fl2.l<0) {
+    if (fl2.l < fl1.l)
+      return (1);
+    return (0);
+  }
 
   if (fl1.l < fl2.l)
     return (1);
   return (0);
 }
+
+#endif