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
|
#include <stdio.h>
#include <math.h>
extern float MAXNUMF, MAXLOGF, MINLOGF;
int
main()
{
float exp1, minnum, x, y, z, e;
exp1 = expf(1.0F);
minnum = powif(2.0F,-149);
x = exp1;
y = MINLOGF + logf(0.501);
/*y = MINLOGF - 0.405;*/
z = powf(x,y);
e = (z - minnum) / minnum;
printf("%.16e %.16e\n", z, e);
x = exp1;
y = MAXLOGF;
z = powf(x,y);
e = (z - MAXNUMF) / MAXNUMF;
printf("%.16e %.16e\n", z, e);
x = MAXNUMF;
y = 1.0F/MAXLOGF;
z = powf(x,y);
e = (z - exp1) / exp1;
printf("%.16e %.16e\n", z, e);
x = exp1;
y = MINLOGF;
z = powf(x,y);
e = (z - minnum) / minnum;
printf("%.16e %.16e\n", z, e);
exit(0);
}
|
