Subversion Repositories HelenOS

Rev

Rev 230 | Go to most recent revision | Only display areas with differences | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

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