/arch/mips/math-emu/sp_frexp.c

https://bitbucket.org/evzijst/gittest · C · 53 lines · 22 code · 4 blank · 27 comment · 1 complexity · 517cf95b22e0cd6c19dcb2cb60e06048 MD5 · raw file 

    

  1. /* IEEE754 floating point arithmetic
    2. 

  2.  * single precision
    3. 

  3.  */
    4. 

  4. /*
    5. 

  5.  * MIPS floating point support
    6. 

  6.  * Copyright (C) 1994-2000 Algorithmics Ltd.
    7. 

  7.  * http://www.algor.co.uk
    8. 

  8.  *
    9. 

  9.  * ########################################################################
    10. 

  10.  *
    11. 

  11.  *  This program is free software; you can distribute it and/or modify it
    12. 

  12.  *  under the terms of the GNU General Public License (Version 2) as
    13. 

  13.  *  published by the Free Software Foundation.
    14. 

  14.  *
    15. 

  15.  *  This program is distributed in the hope it will be useful, but WITHOUT
    16. 

  16.  *  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
    17. 

  17.  *  FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
    18. 

  18.  *  for more details.
    19. 

  19.  *
    20. 

  20.  *  You should have received a copy of the GNU General Public License along
    21. 

  21.  *  with this program; if not, write to the Free Software Foundation, Inc.,
    22. 

  22.  *  59 Temple Place - Suite 330, Boston MA 02111-1307, USA.
    23. 

  23.  *
    24. 

  24.  * ########################################################################
    25. 

  25.  */
    26. 

  26. 

  27. 

  28. #include "ieee754sp.h"
    29. 

  29. 

  30. /* close to ieeep754sp_logb
    31. 

  31. */
    32. 

  32. ieee754sp ieee754sp_frexp(ieee754sp x, int *eptr)
    33. 

  33. {
    34. 

  34. 	COMPXSP;
    35. 

  35. 	CLEARCX;
    36. 

  36. 	EXPLODEXSP;
    37. 

  37. 

  38. 	switch (xc) {
    39. 

  39. 	case IEEE754_CLASS_SNAN:
    40. 

  40. 	case IEEE754_CLASS_QNAN:
    41. 

  41. 	case IEEE754_CLASS_INF:
    42. 

  42. 	case IEEE754_CLASS_ZERO:
    43. 

  43. 		*eptr = 0;
    44. 

  44. 		return x;
    45. 

  45. 	case IEEE754_CLASS_DNORM:
    46. 

  46. 		SPDNORMX;
    47. 

  47. 		break;
    48. 

  48. 	case IEEE754_CLASS_NORM:
    49. 

  49. 		break;
    50. 

  50. 	}
    51. 

  51. 	*eptr = xe + 1;
    52. 

  52. 	return buildsp(xs, -1 + SP_EBIAS, xm & ~SP_HIDDEN_BIT);
    53. 

  53. }
    54.