File: pow-nan.c

package info (click to toggle)
avr-libc 1%3A1.6.2.cvs20080610-2
  • links: PTS
  • area: main
  • in suites: lenny
  • size: 14,848 kB
  • ctags: 55,619
  • sloc: ansic: 92,267; asm: 6,692; sh: 4,131; makefile: 2,481; python: 976; pascal: 426; perl: 116
file content (152 lines) | stat: -rw-r--r-- 4,450 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
/* Test of pow() function. Invalid args.
   $Id: pow-nan.c,v 1.1 2007/02/05 21:35:58 dmix Exp $
 */
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include "progmem.h"

union lofl_u {
    long lo;
    float fl;
};

volatile union lofl_u v = { 1 };

PROGMEM const struct {		/* Table of test cases.	*/
    union lofl_u x, y;		/* args.	*/
} t[] = {

    /* pow (NaN, finite) --> NaN	*/
    { { 0x7f800001 }, { 0x3f800000 } },
    { { 0x7f800100 }, { 0x3f800000 } },
    { { 0x7f810000 }, { 0x3f800000 } },
    { { 0x7fc00000 }, { 0x3f800000 } },
    { { 0x7fffffff }, { 0x3f800000 } },

    { { 0xff800001 }, { 0x3f800000 } },
    { { 0xff800100 }, { 0x3f800000 } },
    { { 0xff810000 }, { 0x3f800000 } },
    { { 0xffc00000 }, { 0x3f800000 } },
    { { 0xffffffff }, { 0x3f800000 } },

/*  { { 0x7fc00000 }, { 0x00000000 } },	    Glibc: pow(NaN,0) --> 1	*/
    { { 0x7fc00000 }, { 0x00000001 } },
    { { 0x7fc00000 }, { 0x00800000 } },
    { { 0x7fc00000 }, { 0x7f7fffff } },

/*  { { 0x7fc00000 }, { 0x80000000 } },	    Glibc: pow(NaN,-0) --> 1	*/
    { { 0x7fc00000 }, { 0x80000001 } },
    { { 0x7fc00000 }, { 0x80800000 } },
    { { 0x7fc00000 }, { 0xff7fffff } },

    /* pow (finite, NaN) --> NaN	*/
    { { 0x3f800001 }, { 0x7f800001 } },
    { { 0x3f800001 }, { 0x7f800100 } },
    { { 0x3f800001 }, { 0x7f810000 } },
    { { 0x3f800001 }, { 0x7fc00000 } },
    { { 0x3f800001 }, { 0x7fffffff } },

    { { 0x3f800001 }, { 0xff800001 } },
    { { 0x3f800001 }, { 0xff800100 } },
    { { 0x3f800001 }, { 0xff810000 } },
    { { 0x3f800001 }, { 0xffc00000 } },
    { { 0x3f800001 }, { 0xffffffff } },
    
    { { 0x00000001 }, { 0x7fc00000 } },
    { { 0x00800000 }, { 0x7fc00000 } },
    { { 0x7f7fffff }, { 0x7fc00000 } },

    { { 0x80000001 }, { 0x7fc00000 } },
    { { 0x80800000 }, { 0x7fc00000 } },
    { { 0xff7fffff }, { 0x7fc00000 } },
    
/*  { { 0x3f800000 }, { 0x7fc00000 } },	    Glibc: pow(1,NaN) --> 1	*/
    
    /* pow (-1, NaN) --> NaN	*/
    { { 0xbf800000 }, { 0x7fc00000 } },

    /* pow (0, NaN) --> NaN	*/
    { { 0x00000000 }, { 0x7fc00000 } },
    { { 0x80000000 }, { 0x7fc00000 } },

    /* pow (NaN, NaN) --> NaN	*/
    { { 0x7fc00000 }, { 0x7fc00000 } },
    
    /* pow (negative, nonintegral)  -->  NaN	*/
    { { 0x80000001 }, { .fl= 1.1 } },
    { { 0x807fffff }, { .fl= 1.1 } },
    { { 0x80800000 }, { .fl= 1.1 } },
    { { 0xbf800000 }, { .fl= 1.1 } },
    { { 0xff7fffff }, { .fl= 1.1 } },

#if 0
    /* fmod (finite, 0) --> NaN	*/
    { { 0x00000001 }, { 0x00000000 } },
    { { 0x007fffff }, { 0x00000000 } },
    { { 0x00800000 }, { 0x00000000 } },
    { { 0x7f7fffff }, { 0x00000000 } },
    { { 0x3f800000 }, { 0x80000000 } },
    { { 0xbf800000 }, { 0x00000000 } },
    { { 0xbf800000 }, { 0x80000000 } },
    
    /* fmod (Inf, 0) --> NaN	*/
    { { 0x7f800000 }, { 0x00000000 } },
    { { 0x7f800000 }, { 0x80000000 } },
    { { 0xff800000 }, { 0x00000000 } },
    { { 0xff800000 }, { 0x80000000 } },
    
    /* fmod (Inf, Inf) --> NaN	*/
    { { 0x7f800000 }, { 0x7f800000 } },
    { { 0x7f800000 }, { 0xff800000 } },
    { { 0xff800000 }, { 0x7f800000 } },
    { { 0xff800000 }, { 0xff800000 } },

    /* fmod (NaN, 1) --> NaN	*/
    { { 0x7f800001 }, { 0x3f800000 } },
    { { 0x7f800100 }, { 0x3f800000 } },
    { { 0x7f810000 }, { 0x3f800000 } },
    { { 0x7fc00000 }, { 0x3f800000 } },
    { { 0x7fffffff }, { 0x3f800000 } },
    { { 0xff800001 }, { 0x3f800000 } },
    { { 0xff800100 }, { 0x3f800000 } },
    { { 0xff810000 }, { 0x3f800000 } },
    { { 0xffc00000 }, { 0x3f800000 } },
    { { 0xffffffff }, { 0x3f800000 } },

    /* fmod (1, NaN) --> NaN	*/
    { { 0x3f800000 }, { 0x7f800001 } },
    { { 0x3f800000 }, { 0x7f800100 } },
    { { 0x3f800000 }, { 0x7f810000 } },
    { { 0x3f800000 }, { 0x7fc00000 } },
    { { 0x3f800000 }, { 0x7fffffff } },
    { { 0x3f800000 }, { 0xff800001 } },
    { { 0x3f800000 }, { 0xff800100 } },
    { { 0x3f800000 }, { 0xff810000 } },
    { { 0x3f800000 }, { 0xffc00000 } },
    { { 0x3f800000 }, { 0xffffffff } },
#endif
};

void x_exit (int index)
{
#ifndef	__AVR__
    fprintf (stderr, "t[%d]:  %#lx\n", index - 1, v.lo);
#endif
    exit (index ? index : -1);
}

int main ()
{
    union lofl_u x, y;
    int i;
    
    for (i = 0; i < (int) (sizeof(t) / sizeof(t[0])); i++) {
	x.lo = pgm_read_dword (& t[i].x);
	y.lo = pgm_read_dword (& t[i].y);
	v.fl = pow (x.fl, y.fl);
	if (!isnan(v.fl))
	    x_exit (i+1);
    }
    return 0;
}