GRASS GIS 8 Programmer's Manual 8.4.1(2025)-45ca3179ab
Loading...
Searching...
No Matches
shapiro1.c
Go to the documentation of this file.
1#include <stdio.h>
2#include <stdlib.h>
3#include <math.h>
4#include "local_proto.h"
5
6double *Cdhc_shapiro_wilk(double *x, int n)
7{
8 static double y[2];
9 double a[25], s2, *xcopy;
10 double sumb = 0.0, sumx = 0.0, sumx2 = 0.0;
11 int i, k;
12
13 if ((xcopy = (double *)malloc(n * sizeof(double))) == NULL) {
14 fprintf(stderr, "Memory error in Cdhc_shapiro_wilk\n");
15 exit(EXIT_FAILURE);
16 }
17
18 k = n / 2;
19 for (i = 0; i < n; ++i) {
20 xcopy[i] = x[i];
21 sumx += x[i];
22 sumx2 += x[i] * x[i];
23 }
24 s2 = sumx2 - sumx * sumx / n;
25
26 qsort(xcopy, n, sizeof(double), Cdhc_dcmp);
27
28 if (n == 3) {
29 a[0] = (double).7071;
30 }
31 else if (n == 4) {
32 a[0] = (double).6872;
33 a[1] = (double).1677;
34 }
35 else if (n == 5) {
36 a[0] = (double).6646;
37 a[1] = (double).2413;
38 }
39 else if (n == 6) {
40 a[0] = (double).6431;
41 a[1] = (double).2806;
42 a[2] = (double).0875;
43 }
44 else if (n == 7) {
45 a[0] = (double).6233;
46 a[1] = (double).3031;
47 a[2] = (double).1401;
48 }
49 else if (n == 8) {
50 a[0] = (double).6052;
51 a[1] = (double).3164;
52 a[2] = (double).1743;
53 a[3] = (double).0561;
54 }
55 else if (n == 9) {
56 a[0] = (double).5888;
57 a[1] = (double).3244;
58 a[2] = (double).1976;
59 a[3] = (double).0947;
60 }
61 else if (n == 10) {
62 a[0] = (double).5739;
63 a[1] = (double).3291;
64 a[2] = (double).2141;
65 a[3] = (double).1224;
66 a[4] = (double).0399;
67 }
68 else if (n == 11) {
69 a[0] = (double).5601;
70 a[1] = (double).3315;
71 a[2] = (double).226;
72 a[3] = (double).1429;
73 a[4] = (double).0695;
74 }
75 else if (n == 12) {
76 a[0] = (double).5475;
77 a[1] = (double).3325;
78 a[2] = (double).2347;
79 a[3] = (double).1586;
80 a[4] = (double).0922;
81 a[5] = (double).0303;
82 }
83 else if (n == 13) {
84 a[0] = (double).5359;
85 a[1] = (double).3325;
86 a[2] = (double).2412;
87 a[3] = (double).1707;
88 a[4] = (double).1099;
89 a[5] = (double).0539;
90 }
91 else if (n == 14) {
92 a[0] = (double).5251;
93 a[1] = (double).3318;
94 a[2] = (double).246;
95 a[3] = (double).1802;
96 a[4] = (double).124;
97 a[5] = (double).0727;
98 a[6] = (double).024;
99 }
100 else if (n == 15) {
101 a[0] = (double).515;
102 a[1] = (double).3306;
103 a[2] = (double).2495;
104 a[3] = (double).1878;
105 a[4] = (double).1353;
106 a[5] = (double).088;
107 a[6] = (double).0433;
108 }
109 else if (n == 16) {
110 a[0] = (double).5056;
111 a[1] = (double).329;
112 a[2] = (double).2521;
113 a[3] = (double).1939;
114 a[4] = (double).1447;
115 a[5] = (double).1005;
116 a[6] = (double).0593;
117 a[7] = (double).0196;
118 }
119 else if (n == 17) {
120 a[0] = (double).4968;
121 a[1] = (double).3273;
122 a[2] = (double).254;
123 a[3] = (double).1988;
124 a[4] = (double).1524;
125 a[5] = (double).1109;
126 a[6] = (double).0725;
127 a[7] = (double).0359;
128 }
129 else if (n == 18) {
130 a[0] = (double).4886;
131 a[1] = (double).3253;
132 a[2] = (double).2553;
133 a[3] = (double).2027;
134 a[4] = (double).1587;
135 a[5] = (double).1197;
136 a[6] = (double).0837;
137 a[7] = (double).0496;
138 a[8] = (double).0163;
139 }
140 else if (n == 19) {
141 a[0] = (double).4808;
142 a[1] = (double).3232;
143 a[2] = (double).2561;
144 a[3] = (double).2059;
145 a[4] = (double).1641;
146 a[5] = (double).1271;
147 a[6] = (double).0932;
148 a[7] = (double).0612;
149 a[8] = (double).0303;
150 }
151 else if (n == 20) {
152 a[0] = (double).4734;
153 a[1] = (double).3211;
154 a[2] = (double).2565;
155 a[3] = (double).2085;
156 a[4] = (double).1686;
157 a[5] = (double).1334;
158 a[6] = (double).1013;
159 a[7] = (double).0711;
160 a[8] = (double).0422;
161 a[9] = (double).014;
162 }
163 else if (n == 21) {
164 a[0] = (double).4643;
165 a[1] = (double).3185;
166 a[2] = (double).2578;
167 a[3] = (double).2119;
168 a[4] = (double).1736;
169 a[5] = (double).1399;
170 a[6] = (double).1092;
171 a[7] = (double).0804;
172 a[8] = (double).053;
173 a[9] = (double).0263;
174 }
175 else if (n == 22) {
176 a[0] = (double).459;
177 a[1] = (double).3156;
178 a[2] = (double).2571;
179 a[3] = (double).2131;
180 a[4] = (double).1764;
181 a[5] = (double).1443;
182 a[6] = (double).115;
183 a[7] = (double).0878;
184 a[8] = (double).0618;
185 a[9] = (double).0368;
186 a[10] = (double).0122;
187 }
188 else if (n == 23) {
189 a[0] = (double).4542;
190 a[1] = (double).3126;
191 a[2] = (double).2563;
192 a[3] = (double).2139;
193 a[4] = (double).1787;
194 a[5] = (double).148;
195 a[6] = (double).1201;
196 a[7] = (double).0941;
197 a[8] = (double).0696;
198 a[9] = (double).0459;
199 a[10] = (double).0228;
200 }
201 else if (n == 24) {
202 a[0] = (double).4493;
203 a[1] = (double).3098;
204 a[2] = (double).2554;
205 a[3] = (double).2145;
206 a[4] = (double).1807;
207 a[5] = (double).1512;
208 a[6] = (double).1245;
209 a[7] = (double).0997;
210 a[8] = (double).0764;
211 a[9] = (double).0539;
212 a[10] = (double).0321;
213 a[11] = (double).0107;
214 }
215 else if (n == 25) {
216 a[0] = (double).445;
217 a[1] = (double).3069;
218 a[2] = (double).2543;
219 a[3] = (double).2148;
220 a[4] = (double).1822;
221 a[5] = (double).1539;
222 a[6] = (double).1283;
223 a[7] = (double).1046;
224 a[8] = (double).0823;
225 a[9] = (double).061;
226 a[10] = (double).0403;
227 a[11] = (double).02;
228 }
229 else if (n == 26) {
230 a[0] = (double).4407;
231 a[1] = (double).3043;
232 a[2] = (double).2533;
233 a[3] = (double).2151;
234 a[4] = (double).1836;
235 a[5] = (double).1563;
236 a[6] = (double).1316;
237 a[7] = (double).1089;
238 a[8] = (double).0876;
239 a[9] = (double).0672;
240 a[10] = (double).0476;
241 a[11] = (double).0284;
242 a[12] = (double).0094;
243 }
244 else if (n == 27) {
245 a[0] = (double).4366;
246 a[1] = (double).3018;
247 a[2] = (double).2522;
248 a[3] = (double).2152;
249 a[4] = (double).1848;
250 a[5] = (double).1584;
251 a[6] = (double).1346;
252 a[7] = (double).1128;
253 a[8] = (double).0923;
254 a[9] = (double).0728;
255 a[10] = (double).054;
256 a[11] = (double).0358;
257 a[12] = (double).0178;
258 }
259 else if (n == 28) {
260 a[0] = (double).4328;
261 a[1] = (double).2992;
262 a[2] = (double).251;
263 a[3] = (double).2151;
264 a[4] = (double).1857;
265 a[5] = (double).1601;
266 a[6] = (double).1372;
267 a[7] = (double).1162;
268 a[8] = (double).0965;
269 a[9] = (double).0778;
270 a[10] = (double).0598;
271 a[11] = (double).0424;
272 a[12] = (double).0253;
273 a[13] = (double).0084;
274 }
275 else if (n == 29) {
276 a[0] = (double).4291;
277 a[1] = (double).2968;
278 a[2] = (double).2499;
279 a[3] = (double).215;
280 a[4] = (double).1864;
281 a[5] = (double).1616;
282 a[6] = (double).1395;
283 a[7] = (double).1192;
284 a[8] = (double).1002;
285 a[9] = (double).0822;
286 a[10] = (double).065;
287 a[11] = (double).0483;
288 a[12] = (double).032;
289 a[13] = (double).0159;
290 }
291 else if (n == 30) {
292 a[0] = (double).4254;
293 a[1] = (double).2944;
294 a[2] = (double).2487;
295 a[3] = (double).2148;
296 a[4] = (double).187;
297 a[5] = (double).163;
298 a[6] = (double).1415;
299 a[7] = (double).1219;
300 a[8] = (double).1036;
301 a[9] = (double).0862;
302 a[10] = (double).0697;
303 a[11] = (double).0537;
304 a[12] = (double).0381;
305 a[13] = (double).0227;
306 a[14] = (double).0076;
307 }
308 else if (n == 31) {
309 a[0] = (double).422;
310 a[1] = (double).2921;
311 a[2] = (double).2475;
312 a[3] = (double).2145;
313 a[4] = (double).1874;
314 a[5] = (double).1641;
315 a[6] = (double).1433;
316 a[7] = (double).1243;
317 a[8] = (double).1066;
318 a[9] = (double).0899;
319 a[10] = (double).0739;
320 a[11] = (double).0585;
321 a[12] = (double).0435;
322 a[13] = (double).0289;
323 a[14] = (double).0144;
324 }
325 else if (n == 32) {
326 a[0] = (double).4188;
327 a[1] = (double).2898;
328 a[2] = (double).2463;
329 a[3] = (double).2141;
330 a[4] = (double).1878;
331 a[5] = (double).1651;
332 a[6] = (double).1449;
333 a[7] = (double).1265;
334 a[8] = (double).1093;
335 a[9] = (double).0931;
336 a[10] = (double).0777;
337 a[11] = (double).0629;
338 a[12] = (double).0485;
339 a[13] = (double).0344;
340 a[14] = (double).0206;
341 a[15] = (double).0068;
342 }
343 else if (n == 33) {
344 a[0] = (double).4156;
345 a[1] = (double).2876;
346 a[2] = (double).2451;
347 a[3] = (double).2137;
348 a[4] = (double).188;
349 a[5] = (double).166;
350 a[6] = (double).1463;
351 a[7] = (double).1284;
352 a[8] = (double).1118;
353 a[9] = (double).0961;
354 a[10] = (double).0812;
355 a[11] = (double).0669;
356 a[12] = (double).053;
357 a[13] = (double).0395;
358 a[14] = (double).0262;
359 a[15] = (double).0131;
360 }
361 else if (n == 34) {
362 a[0] = (double).4127;
363 a[1] = (double).2854;
364 a[2] = (double).2439;
365 a[3] = (double).2132;
366 a[4] = (double).1882;
367 a[5] = (double).1667;
368 a[6] = (double).1475;
369 a[7] = (double).1301;
370 a[8] = (double).114;
371 a[9] = (double).0988;
372 a[10] = (double).0844;
373 a[11] = (double).0706;
374 a[12] = (double).0572;
375 a[13] = (double).0441;
376 a[14] = (double).0314;
377 a[15] = (double).0187;
378 a[16] = (double).0062;
379 }
380 else if (n == 35) {
381 a[0] = (double).4096;
382 a[1] = (double).2834;
383 a[2] = (double).2427;
384 a[3] = (double).2127;
385 a[4] = (double).1883;
386 a[5] = (double).1673;
387 a[6] = (double).1487;
388 a[7] = (double).1317;
389 a[8] = (double).116;
390 a[9] = (double).1013;
391 a[10] = (double).0873;
392 a[11] = (double).0739;
393 a[12] = (double).061;
394 a[13] = (double).0484;
395 a[14] = (double).0361;
396 a[15] = (double).0239;
397 a[16] = (double).0119;
398 }
399 else if (n == 36) {
400 a[0] = (double).4068;
401 a[1] = (double).2813;
402 a[2] = (double).2415;
403 a[3] = (double).2121;
404 a[4] = (double).1883;
405 a[5] = (double).1678;
406 a[6] = (double).1496;
407 a[7] = (double).1331;
408 a[8] = (double).1179;
409 a[9] = (double).1036;
410 a[10] = (double).09;
411 a[11] = (double).077;
412 a[12] = (double).0645;
413 a[13] = (double).0523;
414 a[14] = (double).0404;
415 a[15] = (double).0287;
416 a[16] = (double).0172;
417 a[17] = (double).0057;
418 }
419 else if (n == 37) {
420 a[0] = (double).404;
421 a[1] = (double).2794;
422 a[2] = (double).2403;
423 a[3] = (double).2116;
424 a[4] = (double).1883;
425 a[5] = (double).1683;
426 a[6] = (double).1505;
427 a[7] = (double).1344;
428 a[8] = (double).1196;
429 a[9] = (double).1056;
430 a[10] = (double).0924;
431 a[11] = (double).0798;
432 a[12] = (double).0677;
433 a[13] = (double).0559;
434 a[14] = (double).0444;
435 a[15] = (double).0331;
436 a[16] = (double).022;
437 a[17] = (double).011;
438 }
439 else if (n == 38) {
440 a[0] = (double).4015;
441 a[1] = (double).2774;
442 a[2] = (double).2391;
443 a[3] = (double).211;
444 a[4] = (double).1881;
445 a[5] = (double).1686;
446 a[6] = (double).1513;
447 a[7] = (double).1356;
448 a[8] = (double).1211;
449 a[9] = (double).1075;
450 a[10] = (double).0947;
451 a[11] = (double).0824;
452 a[12] = (double).0706;
453 a[13] = (double).0592;
454 a[14] = (double).0481;
455 a[15] = (double).0372;
456 a[16] = (double).0264;
457 a[17] = (double).0158;
458 a[18] = (double).0053;
459 }
460 else if (n == 39) {
461 a[0] = (double).3989;
462 a[1] = (double).2755;
463 a[2] = (double).238;
464 a[3] = (double).2104;
465 a[4] = (double).188;
466 a[5] = (double).1689;
467 a[6] = (double).152;
468 a[7] = (double).1366;
469 a[8] = (double).1225;
470 a[9] = (double).1092;
471 a[10] = (double).0967;
472 a[11] = (double).0848;
473 a[12] = (double).0733;
474 a[13] = (double).0622;
475 a[14] = (double).0515;
476 a[15] = (double).0409;
477 a[16] = (double).0305;
478 a[17] = (double).0203;
479 a[18] = (double).0101;
480 }
481 else if (n == 40) {
482 a[0] = (double).3964;
483 a[1] = (double).2737;
484 a[2] = (double).2368;
485 a[3] = (double).2098;
486 a[4] = (double).1878;
487 a[5] = (double).1691;
488 a[6] = (double).1526;
489 a[7] = (double).1376;
490 a[8] = (double).1237;
491 a[9] = (double).1108;
492 a[10] = (double).0986;
493 a[11] = (double).087;
494 a[12] = (double).0759;
495 a[13] = (double).0651;
496 a[14] = (double).0546;
497 a[15] = (double).0444;
498 a[16] = (double).0343;
499 a[17] = (double).0244;
500 a[18] = (double).0146;
501 a[19] = (double).0049;
502 }
503 else if (n == 41) {
504 a[0] = (double).394;
505 a[1] = (double).2719;
506 a[2] = (double).2357;
507 a[3] = (double).2091;
508 a[4] = (double).1876;
509 a[5] = (double).1693;
510 a[6] = (double).1531;
511 a[7] = (double).1384;
512 a[8] = (double).1249;
513 a[9] = (double).1123;
514 a[10] = (double).1004;
515 a[11] = (double).0891;
516 a[12] = (double).0782;
517 a[13] = (double).0677;
518 a[14] = (double).0575;
519 a[15] = (double).0476;
520 a[16] = (double).0379;
521 a[17] = (double).0283;
522 a[18] = (double).0188;
523 a[19] = (double).0094;
524 }
525 else if (n == 42) {
526 a[0] = (double).3917;
527 a[1] = (double).2701;
528 a[2] = (double).2345;
529 a[3] = (double).2085;
530 a[4] = (double).1874;
531 a[5] = (double).1694;
532 a[6] = (double).1535;
533 a[7] = (double).1392;
534 a[8] = (double).1259;
535 a[9] = (double).1136;
536 a[10] = (double).102;
537 a[11] = (double).0909;
538 a[12] = (double).0804;
539 a[13] = (double).0701;
540 a[14] = (double).0602;
541 a[15] = (double).0506;
542 a[16] = (double).0411;
543 a[17] = (double).0318;
544 a[18] = (double).0227;
545 a[19] = (double).0136;
546 a[20] = (double).0045;
547 }
548 else if (n == 43) {
549 a[0] = (double).3894;
550 a[1] = (double).2684;
551 a[2] = (double).2334;
552 a[3] = (double).2078;
553 a[4] = (double).1871;
554 a[5] = (double).1695;
555 a[6] = (double).1539;
556 a[7] = (double).1398;
557 a[8] = (double).1269;
558 a[9] = (double).1149;
559 a[10] = (double).1035;
560 a[11] = (double).0927;
561 a[12] = (double).0824;
562 a[13] = (double).0724;
563 a[14] = (double).0628;
564 a[15] = (double).0534;
565 a[16] = (double).0442;
566 a[17] = (double).0352;
567 a[18] = (double).0263;
568 a[19] = (double).0175;
569 a[20] = (double).0087;
570 }
571 else if (n == 44) {
572 a[0] = (double).3872;
573 a[1] = (double).2667;
574 a[2] = (double).2323;
575 a[3] = (double).2072;
576 a[4] = (double).1868;
577 a[5] = (double).1695;
578 a[6] = (double).1542;
579 a[7] = (double).1405;
580 a[8] = (double).1278;
581 a[9] = (double).116;
582 a[10] = (double).1049;
583 a[11] = (double).0943;
584 a[12] = (double).0842;
585 a[13] = (double).0745;
586 a[14] = (double).0651;
587 a[15] = (double).056;
588 a[16] = (double).0471;
589 a[17] = (double).0383;
590 a[18] = (double).0296;
591 a[19] = (double).0211;
592 a[20] = (double).0126;
593 a[21] = (double).0042;
594 }
595 else if (n == 45) {
596 a[0] = (double).385;
597 a[1] = (double).2651;
598 a[2] = (double).2313;
599 a[3] = (double).2065;
600 a[4] = (double).1865;
601 a[5] = (double).1695;
602 a[6] = (double).1545;
603 a[7] = (double).141;
604 a[8] = (double).1286;
605 a[9] = (double).117;
606 a[10] = (double).1062;
607 a[11] = (double).0959;
608 a[12] = (double).086;
609 a[13] = (double).0765;
610 a[14] = (double).0673;
611 a[15] = (double).0584;
612 a[16] = (double).0497;
613 a[17] = (double).0412;
614 a[18] = (double).0328;
615 a[19] = (double).0245;
616 a[20] = (double).0163;
617 a[21] = (double).0081;
618 }
619 else if (n == 46) {
620 a[0] = (double).383;
621 a[1] = (double).2635;
622 a[2] = (double).2302;
623 a[3] = (double).2058;
624 a[4] = (double).1862;
625 a[5] = (double).1695;
626 a[6] = (double).1548;
627 a[7] = (double).1415;
628 a[8] = (double).1293;
629 a[9] = (double).118;
630 a[10] = (double).1073;
631 a[11] = (double).0972;
632 a[12] = (double).0876;
633 a[13] = (double).0783;
634 a[14] = (double).0694;
635 a[15] = (double).0607;
636 a[16] = (double).0522;
637 a[17] = (double).0439;
638 a[18] = (double).0357;
639 a[19] = (double).0277;
640 a[20] = (double).0197;
641 a[21] = (double).0118;
642 a[22] = (double).0039;
643 }
644 else if (n == 47) {
645 a[0] = (double).3808;
646 a[1] = (double).262;
647 a[2] = (double).2291;
648 a[3] = (double).2052;
649 a[4] = (double).1859;
650 a[5] = (double).1695;
651 a[6] = (double).155;
652 a[7] = (double).142;
653 a[8] = (double).13;
654 a[9] = (double).1189;
655 a[10] = (double).1085;
656 a[11] = (double).0986;
657 a[12] = (double).0892;
658 a[13] = (double).0801;
659 a[14] = (double).0713;
660 a[15] = (double).0628;
661 a[16] = (double).0546;
662 a[17] = (double).0465;
663 a[18] = (double).0385;
664 a[19] = (double).0307;
665 a[20] = (double).0229;
666 a[21] = (double).0153;
667 a[22] = (double).0076;
668 }
669 else if (n == 48) {
670 a[0] = (double).3789;
671 a[1] = (double).2604;
672 a[2] = (double).2281;
673 a[3] = (double).2045;
674 a[4] = (double).1855;
675 a[5] = (double).1693;
676 a[6] = (double).1551;
677 a[7] = (double).1423;
678 a[8] = (double).1306;
679 a[9] = (double).1197;
680 a[10] = (double).1095;
681 a[11] = (double).0998;
682 a[12] = (double).0906;
683 a[13] = (double).0817;
684 a[14] = (double).0731;
685 a[15] = (double).0648;
686 a[16] = (double).0568;
687 a[17] = (double).0489;
688 a[18] = (double).0411;
689 a[19] = (double).0335;
690 a[20] = (double).0259;
691 a[21] = (double).0185;
692 a[22] = (double).0111;
693 a[23] = (double).0037;
694 }
695 else if (n == 49) {
696 a[0] = (double).377;
697 a[1] = (double).2589;
698 a[2] = (double).2271;
699 a[3] = (double).2038;
700 a[4] = (double).1851;
701 a[5] = (double).1692;
702 a[6] = (double).1553;
703 a[7] = (double).1427;
704 a[8] = (double).1312;
705 a[9] = (double).1205;
706 a[10] = (double).1105;
707 a[11] = (double).101;
708 a[12] = (double).0919;
709 a[13] = (double).0832;
710 a[14] = (double).0748;
711 a[15] = (double).0667;
712 a[16] = (double).0588;
713 a[17] = (double).0511;
714 a[18] = (double).0436;
715 a[19] = (double).0361;
716 a[20] = (double).0288;
717 a[21] = (double).0215;
718 a[22] = (double).0143;
719 a[23] = (double).0071;
720 }
721 else if (n == 50) {
722 a[0] = (double).3751;
723 a[1] = (double).2574;
724 a[2] = (double).226;
725 a[3] = (double).2032;
726 a[4] = (double).1847;
727 a[5] = (double).1691;
728 a[6] = (double).1554;
729 a[7] = (double).143;
730 a[8] = (double).1317;
731 a[9] = (double).1212;
732 a[10] = (double).1113;
733 a[11] = (double).102;
734 a[12] = (double).0932;
735 a[13] = (double).0846;
736 a[14] = (double).0764;
737 a[15] = (double).0685;
738 a[16] = (double).0608;
739 a[17] = (double).0532;
740 a[18] = (double).0459;
741 a[19] = (double).0386;
742 a[20] = (double).0314;
743 a[21] = (double).0244;
744 a[22] = (double).0174;
745 a[23] = (double).0104;
746 a[24] = (double).0035;
747 }
748
749 if (n > 50 || n < 3) {
750#ifdef NOISY
751 fprintf(stdout, " THIS IS THE SHAPIRO-WILK TEST FOR SMALL SAMPLES\n");
752 fprintf(stdout, " THE SAMPLE SIZE MUST BE LESS THAN OR EQUAL TO 50\n");
753#endif /* NOISY */
754
755 y[0] = y[1] = 0.0;
756 }
757 else {
758 for (i = 1; i <= k; ++i)
759 sumb += a[i - 1] * (x[n - i + 1] - x[i]);
760
761 y[0] = sumb * sumb / s2;
762 y[1] = s2;
763
764#ifdef NOISY
765 fprintf(stdout, " TEST13 SW(N) =%10.4f\n", y[0]);
766#endif /* NOISY */
767 }
768
769 free(xcopy);
770
771 return y;
772}
#define NULL
Definition ccmath.h:32
int Cdhc_dcmp(const void *i, const void *j)
Definition dcmp.c:1
double * Cdhc_shapiro_wilk(double *x, int n)
Definition shapiro1.c:6
#define x