28 #define min(a, b)   (a) < (b) ? a : b  argument
44 #define SWAPINIT(a, es) swaptype = ((char *)a - (char *)0) % sizeof(long) || \  argument
48 swapfunc(char *a, char *b, int n, int swaptype)  in swapfunc()  argument
51         swapcode(long, a, b, n)  in swapfunc()
53         swapcode(char, a, b, n)  in swapfunc()
56 #define swap(a, b)                  \  argument
58         long t = *(long *)(a);          \
59         *(long *)(a) = *(long *)(b);        \
62         swapfunc(a, b, es, swaptype)
64 #define vecswap(a, b, n)    if ((n) > 0) swapfunc(a, b, n, swaptype)  argument
67 med3(char *a, char *b, char *c, int (*cmp)(const void *, const void *))  in med3()  argument
69     return cmp(a, b) < 0 ?  in med3()
70            (cmp(b, c) < 0 ? b : (cmp(a, c) < 0 ? c : a ))  in med3()
71               :(cmp(b, c) > 0 ? b : (cmp(a, c) < 0 ? a : c ));  in med3()
79     char *a = (char*)aa;  in qsort_local()  local
81 loop:   SWAPINIT(a, es);  in qsort_local()
84         for (pm = (char *)a + es; pm < (char *) a + n * es; pm += es)  in qsort_local()
85             for (pl = pm; pl > (char *) a && cmp(pl - es, pl) > 0;  in qsort_local()
90     pm = (char *)a + (n / 2) * es;  in qsort_local()
92         pl = (char *)a;  in qsort_local()
93         pn = (char *)a + (n - 1) * es;  in qsort_local()
102     swap(a, pm);  in qsort_local()
103     pa = pb = (char *)a + es;  in qsort_local()
105     pc = pd = (char *)a + (n - 1) * es;  in qsort_local()
107         while (pb <= pc && (r = cmp(pb, a)) <= 0) {  in qsort_local()
115         while (pb <= pc && (r = cmp(pc, a)) >= 0) {  in qsort_local()
131         for (pm = (char *) a + es; pm < (char *) a + n * es; pm += es)  in qsort_local()
132             for (pl = pm; pl > (char *) a && cmp(pl - es, pl) > 0;  in qsort_local()
138     pn = (char *)a + n * es;  in qsort_local()
139     r = min(pa - (char *)a, pb - pa);  in qsort_local()
140     vecswap(a, pb - r, r);  in qsort_local()
144         qsort_local(a, r / es, es, cmp);  in qsort_local()
147         a = pn - r;  in qsort_local()
225                 float a = A[i * n + k];  in doMatrixMultiplication()  local
226                 sum[0] += a * B[k * n + j];  in doMatrixMultiplication()
227                 sum[1] += a * B[k * n + j + 1];  in doMatrixMultiplication()
228                 sum[2] += a * B[k * n + j + 2];  in doMatrixMultiplication()
229                 sum[3] += a * B[k * n + j + 3];  in doMatrixMultiplication()
230                 sum[4] += a * B[k * n + j + 4];  in doMatrixMultiplication()
231                 sum[5] += a * B[k * n + j + 5];  in doMatrixMultiplication()
232                 sum[6] += a * B[k * n + j + 6];  in doMatrixMultiplication()
233                 sum[7] += a * B[k * n + j + 7];  in doMatrixMultiplication()