696 lines
21 KiB
C
696 lines
21 KiB
C
/*
|
|
*
|
|
* $Id: blur.c,v 1.6 2002-04-17 02:28:33 bdenney Exp $
|
|
*
|
|
*/
|
|
|
|
#include <stdio.h>
|
|
#include <sys/time.h>
|
|
#include <unistd.h>
|
|
#include <assert.h>
|
|
|
|
#define MAX 128
|
|
int array[MAX][MAX];
|
|
int array2[MAX][MAX];
|
|
#define BLUR_SIZE 3
|
|
#define BLUR_WINDOW_HALF 1
|
|
|
|
#define DEFAULT_TIMES 1000
|
|
|
|
void blur_simple()
|
|
{
|
|
int sum;
|
|
int x,y,x2,y2;
|
|
for (x=1; x<MAX-1; x++)
|
|
for (y=1; y<MAX-1; y++)
|
|
{
|
|
sum = 0;
|
|
for (x2=x-BLUR_WINDOW_HALF; x2<=x+BLUR_WINDOW_HALF; x2++)
|
|
for (y2=y-BLUR_WINDOW_HALF; y2<=y+BLUR_WINDOW_HALF; y2++)
|
|
sum += array[x2][y2];
|
|
array2[x][y] = sum;
|
|
}
|
|
}
|
|
|
|
#ifdef BLUR_UNROLL_INNER
|
|
void blur_unroll_inner()
|
|
{
|
|
int sum;
|
|
int x,y,x2,y2;
|
|
for (x=1; x<MAX-1; x++)
|
|
for (y=1; y<MAX-1; y++)
|
|
{
|
|
sum = array[x-1][y-1] + array[x-1][y] + array[x-1][y+1]
|
|
+ array[x][y-1] + array[x][y] + array[x][y+1]
|
|
+ array[x+1][y-1] + array[x+1][y] + array[x+1][y+1];
|
|
array2[x][y] = sum;
|
|
}
|
|
}
|
|
#endif
|
|
|
|
#ifdef BLUR_USE_FUNCTION_CALL
|
|
void blur_func (int *sum, int *data)
|
|
{
|
|
*sum += *data;
|
|
}
|
|
|
|
void blur_funcall()
|
|
{
|
|
int sum;
|
|
int x,y,x2,y2;
|
|
for (x=1; x<MAX-1; x++)
|
|
for (y=1; y<MAX-1; y++)
|
|
{
|
|
sum = 0;
|
|
for (x2=x-BLUR_WINDOW_HALF; x2<=x+BLUR_WINDOW_HALF; x2++)
|
|
for (y2=y-BLUR_WINDOW_HALF; y2<=y+BLUR_WINDOW_HALF; y2++)
|
|
blur_func (&sum, &array[x2][y2]);
|
|
array2[x][y] = sum;
|
|
}
|
|
}
|
|
#endif //ifdef BLUR_USE_FUNCTION_CALL
|
|
|
|
#ifdef BLUR_USE_SWITCH
|
|
/// With BLUR_USE_SWITCH, don't expect the result to match. I had to make
|
|
//all the cases in the switch do something slightly different, so that they
|
|
//weren't all merged together by a compiler optimization.\n");
|
|
void blur_switch()
|
|
{
|
|
int sum;
|
|
int x,y,x2,y2;
|
|
for (x=1; x<MAX-1; x++)
|
|
for (y=1; y<MAX-1; y++)
|
|
{
|
|
sum = 0;
|
|
for (x2=x-BLUR_WINDOW_HALF; x2<=x+BLUR_WINDOW_HALF; x2++)
|
|
for (y2=y-BLUR_WINDOW_HALF; y2<=y+BLUR_WINDOW_HALF; y2++) {
|
|
switch (y2)
|
|
{
|
|
// all cases are slightly different
|
|
case 0:
|
|
sum += array[x][y2]; break;
|
|
case 1:
|
|
sum += array[x][y2]; break;
|
|
case 2:
|
|
sum += array[x][y2]; break;
|
|
case 3:
|
|
sum += array[y][y2]; break;
|
|
case 4:
|
|
sum += array[y2][y2]; break;
|
|
case 5:
|
|
sum += array[y2][y2]; break;
|
|
case 6:
|
|
sum += array[y2][y2]; break;
|
|
case 7:
|
|
sum += array[x2][y2]; break;
|
|
case 8:
|
|
sum += array[x][x2]; break;
|
|
case 9:
|
|
sum += array[x][x2]; break;
|
|
case 10:
|
|
sum += array[y][x2]; break;
|
|
case 11:
|
|
sum += array[x][y]; break;
|
|
case 12:
|
|
sum += array[x][y]; break;
|
|
case 13:
|
|
sum += array[y2][x]; break;
|
|
case 14:
|
|
sum += array[y2][x]; break;
|
|
case 15:
|
|
sum += array[x2][x]; break;
|
|
case 16:
|
|
sum += array[y][x]; break;
|
|
case 17:
|
|
sum += array[y][x]; break;
|
|
case 18:
|
|
sum += array[y][x]; break;
|
|
case 19:
|
|
sum += array[x][x]; break;
|
|
case 20:
|
|
sum += array[x][x]; break;
|
|
case 21:
|
|
sum += array[x2][x2]; break;
|
|
case 22:
|
|
sum += array[x2][x2]; break;
|
|
case 23:
|
|
sum += array[x2][x2]; break;
|
|
case 24:
|
|
sum += array[x2][x2]; break;
|
|
case 25:
|
|
sum += array[y2][x2]; break;
|
|
case 26:
|
|
sum += array[x2][y]; break;
|
|
case 27:
|
|
sum += array[y2][y]; break;
|
|
case 28:
|
|
sum += array[y2][y]; break;
|
|
case 29:
|
|
sum += array[y2][y]; break;
|
|
case 30:
|
|
sum += array[y][y]; break;
|
|
case 31:
|
|
sum += array[y][y]; break;
|
|
default:
|
|
sum += array[x2][y2]; break;
|
|
}
|
|
}
|
|
array2[x][y] = sum;
|
|
}
|
|
}
|
|
|
|
#endif // ifdef BLUR_USE_SWITCH
|
|
|
|
#ifdef BLUR_USE_SWITCH_CALL
|
|
void blur_func00 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func01 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func02 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func03 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func04 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func05 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func06 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func07 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func08 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func09 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func10 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func11 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func12 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func13 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func14 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func15 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func16 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func17 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func18 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func19 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func20 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func21 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func22 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func23 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func24 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func25 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func26 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func27 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func28 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func29 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func30 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func31 (int *sum, int *data) { *sum += *data; }
|
|
|
|
void blur_switch_call()
|
|
{
|
|
int sum;
|
|
int x,y,x2,y2;
|
|
for (x=1; x<MAX-1; x++)
|
|
for (y=1; y<MAX-1; y++)
|
|
{
|
|
sum = 0;
|
|
for (x2=x-BLUR_WINDOW_HALF; x2<=x+BLUR_WINDOW_HALF; x2++)
|
|
for (y2=y-BLUR_WINDOW_HALF; y2<=y+BLUR_WINDOW_HALF; y2++) {
|
|
switch (y2)
|
|
{
|
|
case 0: blur_func00 (&sum, &array[x2][y2]); break;
|
|
case 1: blur_func01 (&sum, &array[x2][y2]); break;
|
|
case 2: blur_func02 (&sum, &array[x2][y2]); break;
|
|
case 3: blur_func03 (&sum, &array[x2][y2]); break;
|
|
case 4: blur_func04 (&sum, &array[x2][y2]); break;
|
|
case 5: blur_func05 (&sum, &array[x2][y2]); break;
|
|
case 6: blur_func06 (&sum, &array[x2][y2]); break;
|
|
case 7: blur_func07 (&sum, &array[x2][y2]); break;
|
|
case 8: blur_func08 (&sum, &array[x2][y2]); break;
|
|
case 9: blur_func09 (&sum, &array[x2][y2]); break;
|
|
case 10: blur_func10 (&sum, &array[x2][y2]); break;
|
|
case 11: blur_func11 (&sum, &array[x2][y2]); break;
|
|
case 12: blur_func12 (&sum, &array[x2][y2]); break;
|
|
case 13: blur_func13 (&sum, &array[x2][y2]); break;
|
|
case 14: blur_func14 (&sum, &array[x2][y2]); break;
|
|
case 15: blur_func15 (&sum, &array[x2][y2]); break;
|
|
case 16: blur_func16 (&sum, &array[x2][y2]); break;
|
|
case 17: blur_func17 (&sum, &array[x2][y2]); break;
|
|
case 18: blur_func18 (&sum, &array[x2][y2]); break;
|
|
case 19: blur_func19 (&sum, &array[x2][y2]); break;
|
|
case 20: blur_func20 (&sum, &array[x2][y2]); break;
|
|
case 21: blur_func21 (&sum, &array[x2][y2]); break;
|
|
case 22: blur_func22 (&sum, &array[x2][y2]); break;
|
|
case 23: blur_func23 (&sum, &array[x2][y2]); break;
|
|
case 24: blur_func24 (&sum, &array[x2][y2]); break;
|
|
case 25: blur_func25 (&sum, &array[x2][y2]); break;
|
|
case 26: blur_func26 (&sum, &array[x2][y2]); break;
|
|
case 27: blur_func27 (&sum, &array[x2][y2]); break;
|
|
case 28: blur_func28 (&sum, &array[x2][y2]); break;
|
|
case 29: blur_func29 (&sum, &array[x2][y2]); break;
|
|
case 30: blur_func30 (&sum, &array[x2][y2]); break;
|
|
case 31: blur_func31 (&sum, &array[x2][y2]); break;
|
|
default:
|
|
sum += array[x2][y2];
|
|
}
|
|
}
|
|
array2[x][y] = sum;
|
|
}
|
|
}
|
|
#endif //ifdef BLUR_USE_SWITCH_CALL
|
|
|
|
#ifdef BLUR_FNPTR_SWITCH
|
|
void blur_func00 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func01 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func02 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func03 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func04 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func05 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func06 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func07 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func08 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func09 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func10 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func11 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func12 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func13 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func14 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func15 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func16 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func17 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func18 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func19 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func20 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func21 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func22 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func23 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func24 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func25 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func26 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func27 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func28 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func29 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func30 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func31 (int *sum, int *data) { *sum += *data; }
|
|
|
|
void bluf_fnptr_switch()
|
|
{
|
|
int sum;
|
|
int x,y,x2,y2;
|
|
void (*ptr)(int *sum, int *data) = NULL;
|
|
for (x=1; x<MAX-1; x++)
|
|
for (y=1; y<MAX-1; y++)
|
|
{
|
|
sum = 0;
|
|
for (x2=x-BLUR_WINDOW_HALF; x2<=x+BLUR_WINDOW_HALF; x2++)
|
|
for (y2=y-BLUR_WINDOW_HALF; y2<=y+BLUR_WINDOW_HALF; y2++) {
|
|
switch (y2) {
|
|
case 0: ptr = blur_func00; break;
|
|
case 1: ptr = blur_func01; break;
|
|
case 2: ptr = blur_func02; break;
|
|
case 3: ptr = blur_func03; break;
|
|
case 4: ptr = blur_func04; break;
|
|
case 5: ptr = blur_func05; break;
|
|
case 6: ptr = blur_func06; break;
|
|
case 7: ptr = blur_func07; break;
|
|
case 8: ptr = blur_func08; break;
|
|
case 9: ptr = blur_func09; break;
|
|
case 10: ptr = blur_func10; break;
|
|
case 11: ptr = blur_func11; break;
|
|
case 12: ptr = blur_func12; break;
|
|
case 13: ptr = blur_func13; break;
|
|
case 14: ptr = blur_func14; break;
|
|
case 15: ptr = blur_func15; break;
|
|
case 16: ptr = blur_func16; break;
|
|
case 17: ptr = blur_func17; break;
|
|
case 18: ptr = blur_func18; break;
|
|
case 19: ptr = blur_func19; break;
|
|
case 20: ptr = blur_func20; break;
|
|
case 21: ptr = blur_func21; break;
|
|
case 22: ptr = blur_func22; break;
|
|
case 23: ptr = blur_func23; break;
|
|
case 24: ptr = blur_func24; break;
|
|
case 25: ptr = blur_func25; break;
|
|
case 26: ptr = blur_func26; break;
|
|
case 27: ptr = blur_func27; break;
|
|
case 28: ptr = blur_func28; break;
|
|
case 29: ptr = blur_func29; break;
|
|
case 30: ptr = blur_func30; break;
|
|
case 31: ptr = blur_func31; break;
|
|
case 32: ptr = blur_func00; break;
|
|
case 33: ptr = blur_func01; break;
|
|
case 34: ptr = blur_func02; break;
|
|
case 35: ptr = blur_func03; break;
|
|
case 36: ptr = blur_func04; break;
|
|
case 37: ptr = blur_func05; break;
|
|
case 38: ptr = blur_func06; break;
|
|
case 39: ptr = blur_func07; break;
|
|
case 40: ptr = blur_func08; break;
|
|
case 41: ptr = blur_func09; break;
|
|
case 42: ptr = blur_func10; break;
|
|
case 43: ptr = blur_func11; break;
|
|
case 44: ptr = blur_func12; break;
|
|
case 45: ptr = blur_func13; break;
|
|
case 46: ptr = blur_func14; break;
|
|
case 47: ptr = blur_func15; break;
|
|
case 48: ptr = blur_func16; break;
|
|
case 49: ptr = blur_func17; break;
|
|
case 50: ptr = blur_func18; break;
|
|
case 51: ptr = blur_func19; break;
|
|
case 52: ptr = blur_func20; break;
|
|
case 53: ptr = blur_func21; break;
|
|
case 54: ptr = blur_func22; break;
|
|
case 55: ptr = blur_func23; break;
|
|
case 56: ptr = blur_func24; break;
|
|
case 57: ptr = blur_func25; break;
|
|
case 58: ptr = blur_func26; break;
|
|
case 59: ptr = blur_func27; break;
|
|
case 60: ptr = blur_func28; break;
|
|
case 61: ptr = blur_func29; break;
|
|
case 62: ptr = blur_func30; break;
|
|
case 63: ptr = blur_func31; break;
|
|
case 64: ptr = blur_func00; break;
|
|
case 65: ptr = blur_func01; break;
|
|
case 66: ptr = blur_func02; break;
|
|
case 67: ptr = blur_func03; break;
|
|
case 68: ptr = blur_func04; break;
|
|
case 69: ptr = blur_func05; break;
|
|
case 70: ptr = blur_func06; break;
|
|
case 71: ptr = blur_func07; break;
|
|
case 72: ptr = blur_func08; break;
|
|
case 73: ptr = blur_func09; break;
|
|
case 74: ptr = blur_func10; break;
|
|
case 75: ptr = blur_func11; break;
|
|
case 76: ptr = blur_func12; break;
|
|
case 77: ptr = blur_func13; break;
|
|
case 78: ptr = blur_func14; break;
|
|
case 79: ptr = blur_func15; break;
|
|
case 80: ptr = blur_func16; break;
|
|
case 81: ptr = blur_func17; break;
|
|
case 82: ptr = blur_func18; break;
|
|
case 83: ptr = blur_func19; break;
|
|
case 84: ptr = blur_func20; break;
|
|
case 85: ptr = blur_func21; break;
|
|
case 86: ptr = blur_func22; break;
|
|
case 87: ptr = blur_func23; break;
|
|
case 88: ptr = blur_func24; break;
|
|
case 89: ptr = blur_func25; break;
|
|
case 90: ptr = blur_func26; break;
|
|
case 91: ptr = blur_func27; break;
|
|
case 92: ptr = blur_func28; break;
|
|
case 93: ptr = blur_func29; break;
|
|
case 94: ptr = blur_func30; break;
|
|
case 95: ptr = blur_func31; break;
|
|
case 96: ptr = blur_func00; break;
|
|
case 97: ptr = blur_func01; break;
|
|
case 98: ptr = blur_func02; break;
|
|
case 99: ptr = blur_func03; break;
|
|
case 100: ptr = blur_func04; break;
|
|
case 101: ptr = blur_func05; break;
|
|
case 102: ptr = blur_func06; break;
|
|
case 103: ptr = blur_func07; break;
|
|
case 104: ptr = blur_func08; break;
|
|
case 105: ptr = blur_func09; break;
|
|
case 106: ptr = blur_func10; break;
|
|
case 107: ptr = blur_func11; break;
|
|
case 108: ptr = blur_func12; break;
|
|
case 109: ptr = blur_func13; break;
|
|
case 110: ptr = blur_func14; break;
|
|
case 111: ptr = blur_func15; break;
|
|
case 112: ptr = blur_func16; break;
|
|
case 113: ptr = blur_func17; break;
|
|
case 114: ptr = blur_func18; break;
|
|
case 115: ptr = blur_func19; break;
|
|
case 116: ptr = blur_func20; break;
|
|
case 117: ptr = blur_func21; break;
|
|
case 118: ptr = blur_func22; break;
|
|
case 119: ptr = blur_func23; break;
|
|
case 120: ptr = blur_func24; break;
|
|
case 121: ptr = blur_func25; break;
|
|
case 122: ptr = blur_func26; break;
|
|
case 123: ptr = blur_func27; break;
|
|
case 124: ptr = blur_func28; break;
|
|
case 125: ptr = blur_func29; break;
|
|
case 126: ptr = blur_func30; break;
|
|
case 127: ptr = blur_func31; break;
|
|
default: ptr=NULL;
|
|
}
|
|
(*ptr)(&sum, &array[x2][y2]);
|
|
}
|
|
array2[x][y] = sum;
|
|
}
|
|
}
|
|
#endif //ifdef BLUR_FNPTR_SWITCH
|
|
|
|
#ifdef BLUR_FNPTR_TABLE
|
|
void blur_func00 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func01 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func02 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func03 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func04 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func05 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func06 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func07 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func08 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func09 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func10 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func11 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func12 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func13 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func14 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func15 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func16 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func17 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func18 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func19 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func20 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func21 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func22 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func23 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func24 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func25 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func26 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func27 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func28 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func29 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func30 (int *sum, int *data) { *sum += *data; }
|
|
void blur_func31 (int *sum, int *data) { *sum += *data; }
|
|
|
|
typedef void (*blurfn)(int *sum, int *data);
|
|
blurfn fntable[MAX] = {
|
|
blur_func00,
|
|
blur_func01,
|
|
blur_func02,
|
|
blur_func03,
|
|
blur_func04,
|
|
blur_func05,
|
|
blur_func06,
|
|
blur_func07,
|
|
blur_func08,
|
|
blur_func09,
|
|
blur_func10,
|
|
blur_func11,
|
|
blur_func12,
|
|
blur_func13,
|
|
blur_func14,
|
|
blur_func15,
|
|
blur_func16,
|
|
blur_func17,
|
|
blur_func18,
|
|
blur_func19,
|
|
blur_func20,
|
|
blur_func21,
|
|
blur_func22,
|
|
blur_func23,
|
|
blur_func24,
|
|
blur_func25,
|
|
blur_func26,
|
|
blur_func27,
|
|
blur_func28,
|
|
blur_func29,
|
|
blur_func30,
|
|
blur_func31,
|
|
blur_func00,
|
|
blur_func01,
|
|
blur_func02,
|
|
blur_func03,
|
|
blur_func04,
|
|
blur_func05,
|
|
blur_func06,
|
|
blur_func07,
|
|
blur_func08,
|
|
blur_func09,
|
|
blur_func10,
|
|
blur_func11,
|
|
blur_func12,
|
|
blur_func13,
|
|
blur_func14,
|
|
blur_func15,
|
|
blur_func16,
|
|
blur_func17,
|
|
blur_func18,
|
|
blur_func19,
|
|
blur_func20,
|
|
blur_func21,
|
|
blur_func22,
|
|
blur_func23,
|
|
blur_func24,
|
|
blur_func25,
|
|
blur_func26,
|
|
blur_func27,
|
|
blur_func28,
|
|
blur_func29,
|
|
blur_func30,
|
|
blur_func31,
|
|
blur_func00,
|
|
blur_func01,
|
|
blur_func02,
|
|
blur_func03,
|
|
blur_func04,
|
|
blur_func05,
|
|
blur_func06,
|
|
blur_func07,
|
|
blur_func08,
|
|
blur_func09,
|
|
blur_func10,
|
|
blur_func11,
|
|
blur_func12,
|
|
blur_func13,
|
|
blur_func14,
|
|
blur_func15,
|
|
blur_func16,
|
|
blur_func17,
|
|
blur_func18,
|
|
blur_func19,
|
|
blur_func20,
|
|
blur_func21,
|
|
blur_func22,
|
|
blur_func23,
|
|
blur_func24,
|
|
blur_func25,
|
|
blur_func26,
|
|
blur_func27,
|
|
blur_func28,
|
|
blur_func29,
|
|
blur_func30,
|
|
blur_func31,
|
|
blur_func00,
|
|
blur_func01,
|
|
blur_func02,
|
|
blur_func03,
|
|
blur_func04,
|
|
blur_func05,
|
|
blur_func06,
|
|
blur_func07,
|
|
blur_func08,
|
|
blur_func09,
|
|
blur_func10,
|
|
blur_func11,
|
|
blur_func12,
|
|
blur_func13,
|
|
blur_func14,
|
|
blur_func15,
|
|
blur_func16,
|
|
blur_func17,
|
|
blur_func18,
|
|
blur_func19,
|
|
blur_func20,
|
|
blur_func21,
|
|
blur_func22,
|
|
blur_func23,
|
|
blur_func24,
|
|
blur_func25,
|
|
blur_func26,
|
|
blur_func27,
|
|
blur_func28,
|
|
blur_func29,
|
|
blur_func30,
|
|
blur_func31,
|
|
};
|
|
|
|
void blur_fnptr_table()
|
|
{
|
|
int sum;
|
|
int x,y,x2,y2;
|
|
void (*ptr)(int *sum, int *data) = NULL;
|
|
for (x=1; x<MAX-1; x++)
|
|
for (y=1; y<MAX-1; y++)
|
|
{
|
|
sum = 0;
|
|
for (x2=x-BLUR_WINDOW_HALF; x2<=x+BLUR_WINDOW_HALF; x2++)
|
|
for (y2=y-BLUR_WINDOW_HALF; y2<=y+BLUR_WINDOW_HALF; y2++) {
|
|
ptr = fntable[y2];
|
|
(*ptr)(&sum, &array[x2][y2]);
|
|
}
|
|
array2[x][y] = sum;
|
|
}
|
|
}
|
|
#endif // ifdef BLUR_FNPTR_TABLE
|
|
|
|
void fill_array()
|
|
{
|
|
int x,y;
|
|
for (x=0; x<MAX; x++)
|
|
for (y=0; y<MAX; y++)
|
|
array[x][y] = (x*17+y*31)%29;
|
|
}
|
|
|
|
void dump_array (FILE *fp, int ptr[MAX][MAX])
|
|
{
|
|
int x,y;
|
|
for (x=0; x<MAX; x++) {
|
|
for (y=0; y<MAX; y++) {
|
|
fprintf (fp, "%3d ", ptr[x][y]);
|
|
}
|
|
fprintf (fp, "\n");
|
|
}
|
|
}
|
|
|
|
struct timeval start, stop;
|
|
#define start_timer() gettimeofday (&start, NULL);
|
|
#define stop_timer() gettimeofday (&stop, NULL);
|
|
|
|
void report_time (FILE *fp, int iters)
|
|
{
|
|
int usec_duration =
|
|
(stop.tv_sec*1000000 + stop.tv_usec)
|
|
- (start.tv_sec*1000000 + start.tv_usec);
|
|
double sec = (double)usec_duration / 1.0e3;
|
|
double sec_per_iter = sec / (double)iters;
|
|
fprintf (fp, "Total time elapsed = %f msec\n", sec);
|
|
fprintf (fp, "Iterations = %d\n", iters);
|
|
fprintf (fp, "Time per iteration = %f msec\n", sec_per_iter);
|
|
}
|
|
|
|
int main (int argc, char *argv[])
|
|
{
|
|
int i;
|
|
int times = 0;
|
|
FILE *out;
|
|
if (argc>1) {
|
|
assert (sscanf (argv[1], "%d", ×) == 1);
|
|
} else {
|
|
times = DEFAULT_TIMES;
|
|
}
|
|
fill_array();
|
|
//dump_array (stderr);
|
|
start_timer();
|
|
for (i=0; i<times; i++) {
|
|
#if defined BLUR_UNROLL_INNER
|
|
blur_unroll_inner ();
|
|
#elif defined BLUR_USE_FUNCTION_CALL
|
|
blur_funcall ();
|
|
#elif defined BLUR_USE_SWITCH
|
|
blur_switch();
|
|
/// With BLUR_USE_SWITCH, don't expect the result to match. I had to make
|
|
//all the cases in the switch do something slightly different, so that they
|
|
//weren't all merged together by a compiler optimization.\n");
|
|
#elif defined BLUR_USE_SWITCH_CALL
|
|
blur_switch_call();
|
|
#elif defined BLUR_FNPTR_SWITCH
|
|
bluf_fnptr_switch();
|
|
#elif defined BLUR_FNPTR_TABLE
|
|
blur_fnptr_table();
|
|
#else
|
|
blur_simple();
|
|
#endif
|
|
}
|
|
stop_timer();
|
|
report_time (stdout, times);
|
|
//fprintf (stderr, "-----------------------------------\n");
|
|
out = fopen ("blur.out", "w");
|
|
assert (out != NULL);
|
|
dump_array (out, array2);
|
|
fclose (out);
|
|
return 0;
|
|
}
|