Subversion Repositories HelenOS-historic

Rev

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

Rev 647 Rev 652
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<softfloat.h>
 29 
#include<softfloat.h>
30 
#include<sftypes.h>
 30 
#include<sftypes.h>
31 
#include<arithmetic.h>
 31 
#include<arithmetic.h>
32 
#include<conversion.h>
 32 
#include<conversion.h>
33 
#include<comparison.h>
 33 
#include<comparison.h>
34 
#include<other.h>
 34 
#include<other.h>
35 
 
 35 
 
36 
/* Arithmetic functions */
 36 
/* Arithmetic functions */
37 
 
 37 
 
38 
float __addsf3(float a, float b)
 38 
float __addsf3(float a, float b)
39 
{
 39 
{
40 
    float32 fa, fb;
 40 
    float32 fa, fb;
41 
    fa.f=a;
 41 
    fa.f=a;
42 
    fb.f=b;
 42 
    fb.f=b;
43 
    if (fa.parts.sign!=fb.parts.sign) {
 43 
    if (fa.parts.sign!=fb.parts.sign) {
44 
        if (fa.parts.sign) {
 44 
        if (fa.parts.sign) {
45 
            fa.parts.sign=0;
 45 
            fa.parts.sign=0;
46 
            return subFloat32(fb,fa).f;
 46 
            return subFloat32(fb,fa).f;
47 
        };
 47 
        };
48 
        fb.parts.sign=0;
 48 
        fb.parts.sign=0;
49 
        return subFloat32(fa,fb).f;
 49 
        return subFloat32(fa,fb).f;
50 
    }
 50 
    }
51 
    return addFloat32(fa,fb).f;
 51 
    return addFloat32(fa,fb).f;
52 
};
 52 
};
53 
 
 53 
 
54 
float __subsf3(float a, float b)
 54 
float __subsf3(float a, float b)
55 
{
 55 
{
56 
    float32 fa, fb;
 56 
    float32 fa, fb;
57 
    fa.f=a;
 57 
    fa.f=a;
58 
    fb.f=b;
 58 
    fb.f=b;
59 
    if (fa.parts.sign!=fb.parts.sign) {
 59 
    if (fa.parts.sign!=fb.parts.sign) {
60 
        fb.parts.sign!=fb.parts.sign;
 60 
        fb.parts.sign!=fb.parts.sign;
61 
        return addFloat32(fa,fb).f;
 61 
        return addFloat32(fa,fb).f;
62 
    }
 62 
    }
63 
    return subFloat32(fa,fb).f;
 63 
    return subFloat32(fa,fb).f;
64 
};
 64 
};
65 
 
 65 
 
66 
float __negsf2(float a)
 66 
float __negsf2(float a)
67 
{
 67 
{
68 
    float32 fa;
 68 
    float32 fa;
69 
    fa.f=a;
 69 
    fa.f=a;
70 
    fa.parts.sign=!fa.parts.sign;
 70 
    fa.parts.sign=!fa.parts.sign;
71 
    return fa.f;
 71 
    return fa.f;
72 
};
 72 
};
73 
 
 73 
 
74 
double __negdf2(double a)
 74 
double __negdf2(double a)
75 
{
 75 
{
76 
    float64 fa;
 76 
    float64 fa;
77 
    fa.d=a;
 77 
    fa.d=a;
78 
    fa.parts.sign=!fa.parts.sign;
 78 
    fa.parts.sign=!fa.parts.sign;
79 
    return fa.d;
 79 
    return fa.d;
80 
};
 80 
};
81 
 
 81 
 
82 
/* Conversion functions */
 82 
/* Conversion functions */
83 
 
 83 
 
84 
/* Comparison functions */
 84 
/* Comparison functions */
85 
 
 85 
 
- 
 
 86 
/* a<b .. -1
- 
 
 87 
 * a=b ..  0
- 
 
 88 
 * a>b ..  1
- 
 
 89 
 * */
- 
 
 90 
 
- 
 
 91 
int __cmpsf2(double a, double b)
- 
 
 92 
{
- 
 
 93 
    float32 fa,fb;
- 
 
 94 
    fa.f=a;
- 
 
 95 
    fb.f=b;
- 
 
 96 
    if ((isFloat32NaN(fa))||(isFloat32NaN(fb))) {
- 
 
 97 
        return 1; /* no special constant for unordered - maybe signaled? */
- 
 
 98 
    };
- 
 
 99 
 
- 
 
 100 
   
- 
 
 101 
    if (isFloat32eq(fa,fb)) {
- 
 
 102 
        return 0;
- 
 
 103 
    };
- 
 
 104 
   
- 
 
 105 
    if (isFloat32lt(fa,fb)) {
- 
 
 106 
        return -1;
- 
 
 107 
        };
- 
 
 108 
    return 1;
- 
 
 109 
}
- 
 
 110 
 
- 
 
 111 
int __unordsf2(float a, float b)
- 
 
 112 
{
- 
 
 113 
    float32 fa,fb;
- 
 
 114 
    fa.f=a;
- 
 
 115 
    fb.f=b;
- 
 
 116 
    return ((isFloat32NaN(fa))||(isFloat32NaN(fb)));
- 
 
 117 
};
- 
 
 118 
 
- 
 
 119 
/**
- 
 
 120 
 * @return zero, if neither argument is a NaN and are equal
- 
 
 121 
 * */
- 
 
 122 
int __eqsf2(float a, float b)
- 
 
 123 
{
- 
 
 124 
    float32 fa,fb;
- 
 
 125 
    fa.f=a;
- 
 
 126 
    fb.f=b;
- 
 
 127 
    if ((isFloat32NaN(fa))||(isFloat32NaN(fb))) {
- 
 
 128 
        /* TODO: sigNaNs*/
- 
 
 129 
        return 1;
- 
 
 130 
        };
- 
 
 131 
    return isFloat32eq(fa,fb)-1;
- 
 
 132 
};
- 
 
 133 
 
- 
 
 134 
/* strange behavior, but it was in gcc documentation */
- 
 
 135 
int __nesf2(float a, float b)
- 
 
 136 
{
- 
 
 137 
    return __eqsf2(a,b);
- 
 
 138 
};
86 
/* Other functions */
 139 
/* Other functions */
87 
 
 140 
 
88 
 
 141