| 1 | //
 | 
|---|
| 2 | // Cforall Version 1.0.0 Copyright (C) 2016 University of Waterloo
 | 
|---|
| 3 | //
 | 
|---|
| 4 | // The contents of this file are covered under the licence agreement in the
 | 
|---|
| 5 | // file "LICENCE" distributed with Cforall.
 | 
|---|
| 6 | //
 | 
|---|
| 7 | // bits/algorithms.h -- Builtins for exception handling.
 | 
|---|
| 8 | //
 | 
|---|
| 9 | // Author           : Thierry Delisle
 | 
|---|
| 10 | // Created On       : Mon Oct 30 13:37:34 2017
 | 
|---|
| 11 | // Last Modified By : --
 | 
|---|
| 12 | // Last Modified On : --
 | 
|---|
| 13 | // Update Count     : 0
 | 
|---|
| 14 | //
 | 
|---|
| 15 | 
 | 
|---|
| 16 | #pragma once
 | 
|---|
| 17 | 
 | 
|---|
| 18 | #ifdef SAFE_SORT
 | 
|---|
| 19 | forall( otype T | {  int ?<?( T, T ); int ?>?( T, T ); } ) static inline void __libcfa_small_sort2( T * arr );
 | 
|---|
| 20 | forall( otype T | {  int ?<?( T, T ); int ?>?( T, T ); } ) static inline void __libcfa_small_sort3( T * arr );
 | 
|---|
| 21 | forall( otype T | {  int ?<?( T, T ); int ?>?( T, T ); } ) static inline void __libcfa_small_sort4( T * arr );
 | 
|---|
| 22 | forall( otype T | {  int ?<?( T, T ); int ?>?( T, T ); } ) static inline void __libcfa_small_sort5( T * arr );
 | 
|---|
| 23 | forall( otype T | {  int ?<?( T, T ); int ?>?( T, T ); } ) static inline void __libcfa_small_sort6( T * arr );
 | 
|---|
| 24 | forall( otype T | {  int ?<?( T, T ); int ?>?( T, T ); } ) static inline void __libcfa_small_sortN( T * arr, size_t dim );
 | 
|---|
| 25 | 
 | 
|---|
| 26 | forall( otype T | {  int ?<?( T, T ); int ?>?( T, T ); } )
 | 
|---|
| 27 | static inline void __libcfa_small_sort( T * arr, size_t dim ) {
 | 
|---|
| 28 |         switch( dim ) {
 | 
|---|
| 29 |                 case 1 : return;
 | 
|---|
| 30 |                 case 2 : __libcfa_small_sort2( arr ); return;
 | 
|---|
| 31 |                 case 3 : __libcfa_small_sort3( arr ); return;
 | 
|---|
| 32 |                 case 4 : __libcfa_small_sort4( arr ); return;
 | 
|---|
| 33 |                 case 5 : __libcfa_small_sort5( arr ); return;
 | 
|---|
| 34 |                 case 6 : __libcfa_small_sort6( arr ); return;
 | 
|---|
| 35 |                 default: __libcfa_small_sortN( arr, dim ); return;
 | 
|---|
| 36 |         }
 | 
|---|
| 37 | }
 | 
|---|
| 38 | 
 | 
|---|
| 39 | #define min(x, y) (y > x ? x : y)
 | 
|---|
| 40 | #define max(x, y) (y > x ? y : x)
 | 
|---|
| 41 | #define SWAP(x,y) { T a = min(arr[x], arr[y]); T b = max(arr[x], arr[y]); arr[x] = a; arr[y] = b;}
 | 
|---|
| 42 | 
 | 
|---|
| 43 | forall( otype T | {  int ?<?( T, T ); int ?>?( T, T ); } )
 | 
|---|
| 44 | static inline void __libcfa_small_sort2( T * arr ) {
 | 
|---|
| 45 |         SWAP(0, 1);
 | 
|---|
| 46 | }
 | 
|---|
| 47 | 
 | 
|---|
| 48 | forall( otype T | {  int ?<?( T, T ); int ?>?( T, T ); } )
 | 
|---|
| 49 | static inline void __libcfa_small_sort3( T * arr ) {
 | 
|---|
| 50 |         SWAP(1, 2);
 | 
|---|
| 51 |         SWAP(0, 2);
 | 
|---|
| 52 |         SWAP(0, 1);
 | 
|---|
| 53 | }
 | 
|---|
| 54 | 
 | 
|---|
| 55 | forall( otype T | {  int ?<?( T, T ); int ?>?( T, T ); } )
 | 
|---|
| 56 | static inline void __libcfa_small_sort4( T * arr ) {
 | 
|---|
| 57 |         SWAP(0, 1);
 | 
|---|
| 58 |         SWAP(2, 3);
 | 
|---|
| 59 |         SWAP(0, 2);
 | 
|---|
| 60 |         SWAP(1, 3);
 | 
|---|
| 61 |         SWAP(1, 2);
 | 
|---|
| 62 | }
 | 
|---|
| 63 | 
 | 
|---|
| 64 | forall( otype T | {  int ?<?( T, T ); int ?>?( T, T ); } )
 | 
|---|
| 65 | static inline void __libcfa_small_sort5( T * arr ) {
 | 
|---|
| 66 |         SWAP(0, 1);
 | 
|---|
| 67 |         SWAP(3, 4);
 | 
|---|
| 68 |         SWAP(2, 4);
 | 
|---|
| 69 |         SWAP(2, 3);
 | 
|---|
| 70 |         SWAP(0, 3);
 | 
|---|
| 71 |         SWAP(0, 2);
 | 
|---|
| 72 |         SWAP(1, 4);
 | 
|---|
| 73 |         SWAP(1, 3);
 | 
|---|
| 74 |         SWAP(1, 2);
 | 
|---|
| 75 | }
 | 
|---|
| 76 | 
 | 
|---|
| 77 | forall( otype T | {  int ?<?( T, T ); int ?>?( T, T ); } )
 | 
|---|
| 78 | static inline void __libcfa_small_sort6( T * arr ) {
 | 
|---|
| 79 |         SWAP(1, 2);
 | 
|---|
| 80 |         SWAP(4, 5);
 | 
|---|
| 81 |         SWAP(0, 2);
 | 
|---|
| 82 |         SWAP(3, 5);
 | 
|---|
| 83 |         SWAP(0, 1);
 | 
|---|
| 84 |         SWAP(3, 4);
 | 
|---|
| 85 |         SWAP(1, 4);
 | 
|---|
| 86 |         SWAP(0, 3);
 | 
|---|
| 87 |         SWAP(2, 5);
 | 
|---|
| 88 |         SWAP(1, 3);
 | 
|---|
| 89 |         SWAP(2, 4);
 | 
|---|
| 90 |         SWAP(2, 3);
 | 
|---|
| 91 | }
 | 
|---|
| 92 | 
 | 
|---|
| 93 | forall( otype T | {  int ?<?( T, T ); int ?>?( T, T ); } )
 | 
|---|
| 94 | static inline void __libcfa_small_sortN( T * arr, size_t dim ) {
 | 
|---|
| 95 |         int i, j;
 | 
|---|
| 96 |         for (i = 1; i < dim; i++) {
 | 
|---|
| 97 |                 T tmp = arr[i];
 | 
|---|
| 98 |                 for (j = i; j >= 1 && tmp < arr[j-1]; j--) {
 | 
|---|
| 99 |                         arr[j] = arr[j-1];
 | 
|---|
| 100 |                 }
 | 
|---|
| 101 |                 arr[j] = tmp;
 | 
|---|
| 102 |         }
 | 
|---|
| 103 | }
 | 
|---|
| 104 | 
 | 
|---|
| 105 | #else
 | 
|---|
| 106 | 
 | 
|---|
| 107 | static inline void __libcfa_small_sort2( void* * arr );
 | 
|---|
| 108 | static inline void __libcfa_small_sort3( void* * arr );
 | 
|---|
| 109 | static inline void __libcfa_small_sort4( void* * arr );
 | 
|---|
| 110 | static inline void __libcfa_small_sort5( void* * arr );
 | 
|---|
| 111 | static inline void __libcfa_small_sort6( void* * arr );
 | 
|---|
| 112 | static inline void __libcfa_small_sortN( void* * arr, size_t dim );
 | 
|---|
| 113 | 
 | 
|---|
| 114 | forall( dtype T )
 | 
|---|
| 115 | static inline void __libcfa_small_sort( T* * arr, size_t dim ) {
 | 
|---|
| 116 |         switch( dim ) {
 | 
|---|
| 117 |                 case 1 : return;
 | 
|---|
| 118 |                 case 2 : __libcfa_small_sort2( (void **) arr ); return;
 | 
|---|
| 119 |                 case 3 : __libcfa_small_sort3( (void **) arr ); return;
 | 
|---|
| 120 |                 case 4 : __libcfa_small_sort4( (void **) arr ); return;
 | 
|---|
| 121 |                 case 5 : __libcfa_small_sort5( (void **) arr ); return;
 | 
|---|
| 122 |                 case 6 : __libcfa_small_sort6( (void **) arr ); return;
 | 
|---|
| 123 |                 default: __libcfa_small_sortN( (void **) arr, dim ); return;
 | 
|---|
| 124 |         }
 | 
|---|
| 125 | }
 | 
|---|
| 126 | 
 | 
|---|
| 127 | #define min(x, y) (y > x ? x : y)
 | 
|---|
| 128 | #define max(x, y) (y > x ? y : x)
 | 
|---|
| 129 | #define SWAP(x,y) { void* a = min(arr[x], arr[y]); void* b = max(arr[x], arr[y]); arr[x] = a; arr[y] = b;}
 | 
|---|
| 130 | 
 | 
|---|
| 131 | static inline void __libcfa_small_sort2( void* * arr ) {
 | 
|---|
| 132 |         SWAP(0, 1);
 | 
|---|
| 133 | }
 | 
|---|
| 134 | 
 | 
|---|
| 135 | static inline void __libcfa_small_sort3( void* * arr ) {
 | 
|---|
| 136 |         SWAP(1, 2);
 | 
|---|
| 137 |         SWAP(0, 2);
 | 
|---|
| 138 |         SWAP(0, 1);
 | 
|---|
| 139 | }
 | 
|---|
| 140 | 
 | 
|---|
| 141 | static inline void __libcfa_small_sort4( void* * arr ) {
 | 
|---|
| 142 |         SWAP(0, 1);
 | 
|---|
| 143 |         SWAP(2, 3);
 | 
|---|
| 144 |         SWAP(0, 2);
 | 
|---|
| 145 |         SWAP(1, 3);
 | 
|---|
| 146 |         SWAP(1, 2);
 | 
|---|
| 147 | }
 | 
|---|
| 148 | 
 | 
|---|
| 149 | static inline void __libcfa_small_sort5( void* * arr ) {
 | 
|---|
| 150 |         SWAP(0, 1);
 | 
|---|
| 151 |         SWAP(3, 4);
 | 
|---|
| 152 |         SWAP(2, 4);
 | 
|---|
| 153 |         SWAP(2, 3);
 | 
|---|
| 154 |         SWAP(0, 3);
 | 
|---|
| 155 |         SWAP(0, 2);
 | 
|---|
| 156 |         SWAP(1, 4);
 | 
|---|
| 157 |         SWAP(1, 3);
 | 
|---|
| 158 |         SWAP(1, 2);
 | 
|---|
| 159 | }
 | 
|---|
| 160 | 
 | 
|---|
| 161 | static inline void __libcfa_small_sort6( void* * arr ) {
 | 
|---|
| 162 |         SWAP(1, 2);
 | 
|---|
| 163 |         SWAP(4, 5);
 | 
|---|
| 164 |         SWAP(0, 2);
 | 
|---|
| 165 |         SWAP(3, 5);
 | 
|---|
| 166 |         SWAP(0, 1);
 | 
|---|
| 167 |         SWAP(3, 4);
 | 
|---|
| 168 |         SWAP(1, 4);
 | 
|---|
| 169 |         SWAP(0, 3);
 | 
|---|
| 170 |         SWAP(2, 5);
 | 
|---|
| 171 |         SWAP(1, 3);
 | 
|---|
| 172 |         SWAP(2, 4);
 | 
|---|
| 173 |         SWAP(2, 3);
 | 
|---|
| 174 | }
 | 
|---|
| 175 | 
 | 
|---|
| 176 | static inline void __libcfa_small_sortN( void* * arr, size_t dim ) {
 | 
|---|
| 177 |         int i, j;
 | 
|---|
| 178 |         for (i = 1; i < dim; i++) {
 | 
|---|
| 179 |                 void* tmp = arr[i];
 | 
|---|
| 180 |                 for (j = i; j >= 1 && tmp < arr[j-1]; j--) {
 | 
|---|
| 181 |                         arr[j] = arr[j-1];
 | 
|---|
| 182 |                 }
 | 
|---|
| 183 |                 arr[j] = tmp;
 | 
|---|
| 184 |         }
 | 
|---|
| 185 | }
 | 
|---|
| 186 | 
 | 
|---|
| 187 | #endif
 | 
|---|
| 188 | 
 | 
|---|
| 189 | #undef SWAP
 | 
|---|
| 190 | #undef min
 | 
|---|
| 191 | #undef max
 | 
|---|