Subversion Repositories HelenOS

Rev

Rev 1962 | Blame | Compare with Previous | Last modification | View Log | Download | RSS feed

  1. /*
  2.  * Copyright (C) 2005 Josef Cejka
  3.  * All rights reserved.
  4.  *
  5.  * Redistribution and use in source and binary forms, with or without
  6.  * modification, are permitted provided that the following conditions
  7.  * are met:
  8.  *
  9.  * - Redistributions of source code must retain the above copyright
  10.  *   notice, this list of conditions and the following disclaimer.
  11.  * - Redistributions in binary form must reproduce the above copyright
  12.  *   notice, this list of conditions and the following disclaimer in the
  13.  *   documentation and/or other materials provided with the distribution.
  14.  * - The name of the author may not be used to endorse or promote products
  15.  *   derived from this software without specific prior written permission.
  16.  *
  17.  * THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY EXPRESS OR
  18.  * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES
  19.  * OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED.
  20.  * IN NO EVENT SHALL THE AUTHOR BE LIABLE FOR ANY DIRECT, INDIRECT,
  21.  * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT
  22.  * NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
  23.  * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
  24.  * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
  25.  * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
  26.  * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
  27.  */
  28.  
  29. #include <arch/types.h>
  30.  
  31. #define FMATH_EXPONENT_BIAS 1023
  32.  
  33. typedef unsigned char fmath_ld_descr_t[8];
  34.  
  35. typedef union { double bf; unsigned char ldd[8]; }  fmath_ld_union_t;
  36.  
  37. /**returns exponent in binary encoding*/
  38. signed short fmath_get_binary_exponent(double num);
  39.  
  40. /**returns exponent in decimal encoding*/
  41. double fmath_get_decimal_exponent(double num);
  42.  
  43. /**returns mantisa in binary encoding */
  44. __u64 fmath_get_binary_mantisa(double num) ;
  45.  
  46. /** Function for extract integer part from double
  47. * @param num input value
  48. * @param intp integer part of num
  49. * @return non-integer part
  50. */
  51. double fmath_fint(double num, double *intp);
  52.  
  53. /** count base^exponent from positive exponent
  54. * @param base
  55. * @param exponent - Must be > 0.0
  56. * @return base^exponent or 0.0 (if exponent <=0.0)
  57. */
  58. double fmath_dpow(double base, double exponent) ;
  59.  
  60. /** return 1, if num is NaN */
  61. int fmath_is_nan(double num);
  62.  
  63. /** return 1, if fmath is a infinity */
  64. int fmath_is_infinity(double num);
  65.