git.gag.com Git - debian/cc1111/blob - device/lib/_ulong2fs.c

   1 /* Floating point library in optimized assembly for 8051
   2  * Copyright (c) 2004, Paul Stoffregen, paul@pjrc.com
   3  *
   4  * This program is free software; you can redistribute it and/or
   5  * modify it under the terms of the GNU Library General Public License
   6  * as published by the Free Software Foundation; either version 2
   7  * of the License, or (at your option) any later version.
   8  *
   9  * This library is distributed in the hope that it will be useful,
  10  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  11  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12  * GNU General Public License for more details.
  13  *
  14  * You should have received a copy of the GNU General Public License
  15  * along with this program; if not, write to the Free Software
  16  * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA  02111-1307, USA.
  17  */
  18
  19
  20 #define SDCC_FLOAT_LIB
  21 #include <float.h>
  22
  23
  24 #ifdef FLOAT_ASM_MCS51
  25
  26 // float long __ulong2fs (float x)
  27 static void dummy(void) __naked
  28 {
  29         __asm
  30         .globl  ___ulong2fs
  31 ___ulong2fs:
  32         mov     r4, a
  33         mov     r3, b
  34         mov     r2, dph
  35         mov     r1, dpl
  36         mov     a, #158
  37         .globl  ulong2fs_doit
  38 ulong2fs_doit:
  39         clr     sign_a
  40 long2fs_doit:
  41         mov     exp_a, a
  42         lcall   fs_normalize_a
  43         ljmp    fs_round_and_return
  44         __endasm;
  45 }
  46
  47 #else
  48
  49 /*
  50 ** libgcc support for software floating point.
  51 ** Copyright (C) 1991 by Pipeline Associates, Inc.  All rights reserved.
  52 ** Permission is granted to do *anything* you want with this file,
  53 ** commercial or otherwise, provided this message remains intact.  So there!
  54 ** I would appreciate receiving any updates/patches/changes that anyone
  55 ** makes, and am willing to be the repository for said changes (am I
  56 ** making a big mistake?).
  57 **
  58 ** Pat Wood
  59 ** Pipeline Associates, Inc.
  60 ** pipeline!phw@motown.com or
  61 ** sun!pipeline!phw or
  62 ** uunet!motown!pipeline!phw
  63 */
  64
  65 /* (c)2000/2001: hacked a little by johan.knol@iduna.nl for sdcc */
  66
  67 union float_long
  68   {
  69     float f;
  70     long l;
  71   };
  72
  73 float __ulong2fs (unsigned long a )
  74 {
  75   int exp = 24 + EXCESS;
  76   volatile union float_long fl;
  77
  78   if (!a)
  79     {
  80       return 0.0;
  81     }
  82
  83   while (a & NORM)
  84     {
  85       // we lose accuracy here
  86       a >>= 1;
  87       exp++;
  88     }
  89
  90   while (a < HIDDEN)
  91     {
  92       a <<= 1;
  93       exp--;
  94     }
  95
  96 #if 1
  97   if ((a&0x7fffff)==0x7fffff) {
  98     a=0;
  99     exp++;
 100   }
 101 #endif
 102
 103   a &= ~HIDDEN ;
 104   /* pack up and go home */
 105   fl.l = PACK(0,(unsigned long)exp, a);
 106
 107   return (fl.f);
 108 }
 109
 110 #endif