![[BACK]](/cvsweb/icons/back.gif){kind=icon}
Return to d45.c CVS log ![[TXT]](/cvsweb/icons/text.gif){kind=icon}
![[DIR]](/cvsweb/icons/dir.gif){kind=icon}
Up to [Research Unix] / researchv10no / cmd / lcc / ph|
Return to d45.c CVS log | Up to [Research Unix] / researchv10no / cmd / lcc / ph |
1.1 root 1: /* The Plum Hall Validation Suite for C
2: * Unpublished copyright (c) 1986-1991, Chiron Systems Inc and Plum Hall Inc.
3: * VERSION: 4
4: * DATE: 1993-01-01
5: * The "ANSI" mode of this suite corresponds to official ANSI C, X3.159-1989.
6: * As per your license agreement, your distribution is not to be moved or copied outside the Designated Site
7: * without specific permission from Plum Hall Inc.
8: */
9:
10: #define LIB_TEST 1
11: #include "defs.h"
12: #if !ANSI && !V7
13: #define SKIP45 /* prior to V7, too little commonality in math to be tested */
14: #endif
15:
16: #include "flags.h"
17: #ifndef SKIP45
18: /*
19: * 4.5 - Mathematics
20: */
21: #include <math.h>
22: #if !ANSI
23: #include <errno.h>
24: extern int errno;
25: #ifndef HUGE_VAL
26: #define HUGE_VAL HUGE /* the most common other name */
27: #endif /* of HUGE_VAL */
28: #else /* if ANSI */
29: #if ANSI8703
30: #include <errno.h>
31: #else /* ANSI, but not 8703 */
32: #include <stddef.h>
33: #endif /* ANSI8703 */
34: #endif /* !ANSI */
35: static void d4_5_2();
36: static void d4_5_3();
37: static void d4_5_4();
38: static void d4_5_5();
39: static void d4_5_6();
40: /* 4.5 - Mathmatics <math.h>
41: * This version is not intended to verify the "correctness" of the
42: * results of the math library, but rather to verify the existence
43: * of the functions and some domain/range properties.
44: */
45: void d4_5()
46: {
47: Filename = "d45.c";
48: checkthat(__LINE__, EDOM != ERANGE);
49: checkthat(__LINE__, EDOM != 0);
50: checkthat(__LINE__, ERANGE != 0);
51: checkthat(__LINE__, HUGE_VAL > 0);
52: d4_5_2();
53: d4_5_3();
54: d4_5_4();
55: d4_5_5();
56: d4_5_6();
57: }
58:
59: /*
60: * 4.5.2 - Trignometric functions
61: */
62: static void d4_5_2()
63: {
64: #ifndef SKIP452
65: double pi = acos(-1.);
66: double dummy;
67:
68: /* 4.5.2.1 acos
69: * has domain on [-1,1]
70: */
71: dequals(__LINE__, acos(1.), 0.0);
72: dequals(__LINE__, acos(0.), pi/2);
73: dequals(__LINE__, acos(-1.), pi);
74: dequals(__LINE__, acos(cos(pi/4)), pi/4);
75:
76: /* 4.5.2.2 asin
77: * has domain on [-1,1]
78: */
79: dequals(__LINE__, asin(1.), pi/2);
80: dequals(__LINE__, asin(0.), 0.0);
81: dequals(__LINE__, asin(-1.), -pi/2);
82: dequals(__LINE__, asin(cos(pi/4)), pi/4);
83:
84: /* 4.5.2.3 atan
85: * has domain (-HUGE_VAL, HUGE_VAL)
86: */
87: dequals(__LINE__, atan(0.0), 0.0);
88: dequals(__LINE__, atan(1.0), pi/4);
89: dequals(__LINE__, atan(tan(pi/7)), pi/7);
90:
91:
92: /* 4.5.2.4 atan2
93: * figures out the correct quadrant
94: */
95: dequals(__LINE__, atan2(sin(pi/4), cos(pi/4)), pi/4);
96: dequals(__LINE__, atan2(sin(3*pi/4), cos(3*pi/4)), 3*pi/4);
97: dequals(__LINE__, atan2(sin(-pi/4), cos(-pi/4)), -pi/4);
98: dequals(__LINE__, atan2(sin(-3*pi/4), cos(-3*pi/4)), -3*pi/4);
99: errno = 0;
100: atan2(0.0, 0.0);
101: iequals( - __LINE__, errno, EDOM);
102:
103: /* 4.5.2.5 cos
104: * works in radians
105: */
106: dequals(__LINE__, cos(0.0), 1.0);
107: dequals(__LINE__, cos(pi/2), 0.0);
108: dequals(__LINE__, cos(pi), -1.0);
109: dequals(__LINE__, cos(3*pi/2), 0.0);
110: dequals(__LINE__, cos(pi/4), cos(-pi/4));
111: dequals(__LINE__, cos(pi/4), sin(pi/4));
112:
113: /* 4.5.2.6 sin
114: * works in radians
115: */
116: dequals(__LINE__, sin(0.0), 0.0);
117: dequals(__LINE__, sin(pi/2), 1.0);
118: dequals(__LINE__, sin(pi), 0.0);
119: dequals(__LINE__, sin(3*pi/2), -1.0);
120: dequals(__LINE__, sin(pi/2), -sin(-pi/2));
121:
122: /* 4.5.2.7 tan
123: * also works on radians
124: */
125: dequals(__LINE__, tan(0.0), 0.0);
126: dequals(__LINE__, tan(pi/4), 1.0);
127: dequals(__LINE__, tan(-pi/4), -1.0);
128: #endif /* SKIP452 */
129: }
130:
131: /*
132: * 4.5.3 - Hyperbolic functions
133: */
134: static void d4_5_3()
135: {
136: #ifndef SKIP453
137: double pi = acos(-1.);
138:
139: /* 4.5.3.1 cosh
140: * range error returns HUGE_VAL and sets errno
141: */
142: errno = 0;
143: dequals(__LINE__, cosh(HUGE_VAL), HUGE_VAL);
144: iequals( - __LINE__, errno, ERANGE);
145: dequals(__LINE__, cosh(1.234), cosh(-1.234));
146: errno = 0;
147:
148: /* 4.5.3.2 sinh
149: */
150: dequals(__LINE__, sinh(HUGE_VAL), HUGE_VAL);
151: iequals( - __LINE__, errno, ERANGE);
152: dequals(__LINE__, sinh(1.234), -sinh(-1.234));
153:
154: /* 4.5.3.3 tanh
155: */
156: dequals(__LINE__, tanh(1.234), sinh(1.234)/cosh(1.234));
157: dequals(__LINE__, tanh(0.0), 0.0);
158: dequals(__LINE__, tanh(1.234), -tanh(-1.234));
159: #endif /* SKIP453 */
160: }
161:
162: /*
163: * 4.5.4 - Exponential and logrithmic functions
164: */
165: static void d4_5_4()
166: {
167: #ifndef SKIP454
168: int i;
169: double dummy;
170:
171: /* 4.5.4.1 exp
172: */
173: dequals(__LINE__, exp(HUGE_VAL), HUGE_VAL);
174: iequals( - __LINE__, errno, ERANGE);
175: dequals(__LINE__, exp(0.0), 1.0);
176: dequals(__LINE__, exp(-1.1), 1.0/exp(1.1));
177:
178: /* 4.5.4.2 frexp
179: */
180: dequals(__LINE__, frexp(1.234, &i), .617);
181: iequals(__LINE__, i, 1);
182: dequals(__LINE__, frexp(12.34, &i), .77125);
183: iequals(__LINE__, i, 4);
184: dequals(__LINE__, frexp(0.0, &i), 0.0);
185: iequals(__LINE__, i, 0);
186:
187: /* 4.5.4.3 ldexp
188: */
189: dequals(__LINE__, ldexp(3.0, 2), 12.0);
190: errno = 0;
191: dequals(__LINE__, ldexp(1e37, 32767), HUGE_VAL);
192: iequals( - __LINE__, errno, ERANGE);
193:
194: /* 4.5.4.4 log
195: */
196: errno = 0;
197: dummy = log(-1.0);
198: iequals(__LINE__, errno, EDOM);
199: errno = 0;
200: iequals( - __LINE__, errno, ERANGE);
201: dequals(__LINE__, exp(log(1.234)), 1.234);
202:
203: /* 4.5.4.5 log10
204: */
205: errno = 0;
206: dummy = log10(-1.0);
207: iequals(__LINE__, errno, EDOM);
208: errno = 0;
209: iequals( - __LINE__, errno, ERANGE);
210: dequals(__LINE__, log10(1.0), 0.0);
211: dequals(__LINE__, log10(1.e10), 10.0);
212:
213: /* 4.5.4.6 modf
214: */
215: dequals(__LINE__, modf(12.345, &dummy), .345);
216: dequals(__LINE__, dummy, 12.);
217: dequals(__LINE__, modf(0.0, &dummy), 0.0);
218: dequals(__LINE__, dummy, 0.0);
219: #endif /* SKIP454 */
220: }
221:
222: /*
223: * 4.5.5 - Power functions
224: */
225: static void d4_5_5()
226: {
227: #ifndef SKIP455
228: double dummy;
229:
230: /* 4.5.5.1 pow
231: */
232: errno = 0;
233: dummy = pow(0.0, 0.0);
234: iequals( - __LINE__, errno, EDOM);
235: errno = 0;
236: dummy = pow(0.0, -1.0);
237: iequals( - __LINE__, errno, EDOM);
238: errno = 0;
239: dummy = pow(-4.3, 1.1);
240: iequals(__LINE__, errno, EDOM);
241: dequals(__LINE__, log10(pow(10.0, 1.234)), 1.234);
242: dequals(__LINE__, pow(-3.0, 3.0), -27.0);
243:
244: /* range errors */
245: errno = 0;
246: dequals(__LINE__, pow(10.0, HUGE_VAL), HUGE_VAL);
247: iequals( - __LINE__, errno, ERANGE);
248: errno = 0;
249: /* underflow */
250: dequals(__LINE__, pow(10.0, -HUGE_VAL), 0.0);
251: iequals( - __LINE__, errno, ERANGE);
252:
253: /* 4.5.5.2 sqrt
254: */
255: dequals(__LINE__, pow(sqrt(1.234), 2.0), 1.234);
256: #endif /* SKIP455 */
257: }
258:
259: /*
260: * 4.5.6 - Nearest integer, absolute value, and remainder functions
261: */
262: static void d4_5_6()
263: {
264: #ifndef SKIP456
265: /* 4.5.6.1 ceil
266: */
267: dequals(__LINE__, ceil(0.0), 0.0);
268: dequals(__LINE__, ceil(12.345), 13.0);
269: dequals(__LINE__, ceil(-12.345), -12.0);
270:
271: /* 4.5.6.2 fabs
272: */
273: dequals(__LINE__, fabs(12.345), 12.345);
274: dequals(__LINE__, fabs(-12.345), 12.345);
275:
276: /* 4.5.6.3 floor
277: */
278: dequals(__LINE__, floor(0.0), 0.0);
279: dequals(__LINE__, floor(12.345), 12.0);
280: dequals(__LINE__, floor(-12.345), -13.0);
281: #if ANSI
282: /* 4.5.6.4 fmod
283: */
284: dequals(__LINE__, fmod(13.1, 5.0), 3.1);
285: dequals(__LINE__, fmod(-13.1, 5.0), -3.1);
286: dequals(__LINE__, fmod(13.1, -5.0), 3.1);
287: dequals(__LINE__, fmod(-13.1, -5.0), -3.1);
288: #if ANSI8612
289: dequals(- __LINE__, fmod(1., 0.), 0.); /* DEC86 */
290: checkthat(__LINE__, fmod(1E38, 1E-38) <= 1E-38); /* DEC86 */
291: #endif
292: #endif
293: #endif /* SKIP456 */
294: }
295:
296: #else /* if SKIP45 */
297: void d4_5() { pr_skip("d4_5: SKIPPED ENTIRELY\n"); }
298: #endif
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.