| /** |
| * dynprog.c: This file is part of the PolyBench/C 3.2 test suite. |
| * |
| * |
| * Contact: Louis-Noel Pouchet <pouchet@cse.ohio-state.edu> |
| * Web address: http://polybench.sourceforge.net |
| */ |
| #include <stdio.h> |
| #include <unistd.h> |
| #include <string.h> |
| #include <math.h> |
| |
| /* Include polybench common header. */ |
| #include <polybench.h> |
| |
| /* Include benchmark-specific header. */ |
| /* Default data type is int, default size is 50. */ |
| #include "dynprog.h" |
| |
| |
| /* Array initialization. */ |
| static |
| void init_array(int length, |
| DATA_TYPE POLYBENCH_2D(c,LENGTH,LENGTH,length,length), |
| DATA_TYPE POLYBENCH_2D(W,LENGTH,LENGTH,length,length)) |
| { |
| #pragma STDC FP_CONTRACT OFF |
| int i, j; |
| for (i = 0; i < length; i++) |
| for (j = 0; j < length; j++) { |
| c[i][j] = i*j % 2; |
| W[i][j] = ((DATA_TYPE) i-j) / length; |
| } |
| } |
| |
| |
| /* DCE code. Must scan the entire live-out data. |
| Can be used also to check the correctness of the output. */ |
| static |
| void print_array(DATA_TYPE out) |
| { |
| fprintf (stderr, DATA_PRINTF_MODIFIER, out); |
| fprintf (stderr, "\n"); |
| } |
| |
| |
| /* Main computational kernel. The whole function will be timed, |
| including the call and return. */ |
| static |
| void kernel_dynprog(int tsteps, int length, |
| DATA_TYPE POLYBENCH_2D(c,LENGTH,LENGTH,length,length), |
| DATA_TYPE POLYBENCH_2D(W,LENGTH,LENGTH,length,length), |
| DATA_TYPE POLYBENCH_3D(sum_c,LENGTH,LENGTH,LENGTH,length,length,length), |
| DATA_TYPE *out) |
| { |
| int iter, i, j, k; |
| |
| DATA_TYPE out_l = 0; |
| |
| #pragma scop |
| for (iter = 0; iter < _PB_TSTEPS; iter++) |
| { |
| for (i = 0; i <= _PB_LENGTH - 1; i++) |
| for (j = 0; j <= _PB_LENGTH - 1; j++) |
| c[i][j] = 0; |
| |
| for (i = 0; i <= _PB_LENGTH - 2; i++) |
| { |
| for (j = i + 1; j <= _PB_LENGTH - 1; j++) |
| { |
| sum_c[i][j][i] = 0; |
| for (k = i + 1; k <= j-1; k++) |
| sum_c[i][j][k] = sum_c[i][j][k - 1] + c[i][k] + c[k][j]; |
| c[i][j] = sum_c[i][j][j-1] + W[i][j]; |
| } |
| } |
| out_l += c[0][_PB_LENGTH - 1]; |
| } |
| #pragma endscop |
| |
| *out = out_l; |
| } |
| |
| #if !FMA_DISABLED |
| // NOTE: FMA_DISABLED is true for targets where FMA contraction causes |
| // discrepancies which cause the accuracy checks to fail. |
| // In this case, the test runs with the option -ffp-contract=off |
| static void |
| kernel_dynprog_StrictFP(int tsteps, int length, |
| DATA_TYPE POLYBENCH_2D(c,LENGTH,LENGTH,length,length), |
| DATA_TYPE POLYBENCH_2D(W,LENGTH,LENGTH,length,length), |
| DATA_TYPE POLYBENCH_3D(sum_c,LENGTH,LENGTH,LENGTH,length,length,length), |
| DATA_TYPE *out) |
| { |
| #pragma STDC FP_CONTRACT OFF |
| int iter, i, j, k; |
| |
| DATA_TYPE out_l = 0; |
| |
| for (iter = 0; iter < _PB_TSTEPS; iter++) |
| { |
| for (i = 0; i <= _PB_LENGTH - 1; i++) |
| for (j = 0; j <= _PB_LENGTH - 1; j++) |
| c[i][j] = 0; |
| |
| for (i = 0; i <= _PB_LENGTH - 2; i++) |
| { |
| for (j = i + 1; j <= _PB_LENGTH - 1; j++) |
| { |
| sum_c[i][j][i] = 0; |
| for (k = i + 1; k <= j-1; k++) |
| sum_c[i][j][k] = sum_c[i][j][k - 1] + c[i][k] + c[k][j]; |
| c[i][j] = sum_c[i][j][j-1] + W[i][j]; |
| } |
| } |
| out_l += c[0][_PB_LENGTH - 1]; |
| } |
| |
| *out = out_l; |
| } |
| |
| /* Return 0 when one of the elements of arrays A and B do not match within the |
| allowed FP_ABSTOLERANCE. Return 1 when all elements match. */ |
| static int |
| check_FP(DATA_TYPE A, |
| DATA_TYPE B) { |
| double AbsTolerance = FP_ABSTOLERANCE; |
| double V1 = A; |
| double V2 = B; |
| double Diff = fabs(V1 - V2); |
| if (Diff > AbsTolerance) { |
| fprintf(stderr, "A = %lf and B = %lf differ more than" |
| " FP_ABSTOLERANCE = %lf\n", V1, V2, AbsTolerance); |
| return 0; |
| } |
| return 1; |
| } |
| #endif |
| |
| int main(int argc, char** argv) |
| { |
| /* Retrieve problem size. */ |
| int length = LENGTH; |
| int tsteps = TSTEPS; |
| |
| /* Variable declaration/allocation. */ |
| DATA_TYPE out; |
| #if !FMA_DISABLED |
| DATA_TYPE out_StrictFP; |
| #endif |
| POLYBENCH_3D_ARRAY_DECL(sum_c,DATA_TYPE,LENGTH,LENGTH,LENGTH,length,length,length); |
| POLYBENCH_2D_ARRAY_DECL(c,DATA_TYPE,LENGTH,LENGTH,length,length); |
| POLYBENCH_2D_ARRAY_DECL(W,DATA_TYPE,LENGTH,LENGTH,length,length); |
| |
| /* Initialize array(s). */ |
| init_array (length, POLYBENCH_ARRAY(c), POLYBENCH_ARRAY(W)); |
| |
| /* Start timer. */ |
| polybench_start_instruments; |
| |
| /* Run kernel. */ |
| kernel_dynprog (tsteps, length, |
| POLYBENCH_ARRAY(c), |
| POLYBENCH_ARRAY(W), |
| POLYBENCH_ARRAY(sum_c), |
| &out); |
| |
| /* Stop and print timer. */ |
| polybench_stop_instruments; |
| polybench_print_instruments; |
| |
| #if FMA_DISABLED |
| /* Prevent dead-code elimination. All live-out data must be printed |
| by the function call in argument. */ |
| polybench_prevent_dce(print_array(out)); |
| #else |
| init_array (length, POLYBENCH_ARRAY(c), POLYBENCH_ARRAY(W)); |
| kernel_dynprog (tsteps, length, |
| POLYBENCH_ARRAY(c), |
| POLYBENCH_ARRAY(W), |
| POLYBENCH_ARRAY(sum_c), |
| &out_StrictFP); |
| |
| if (!check_FP(out, out_StrictFP)) |
| return 1; |
| |
| /* Prevent dead-code elimination. All live-out data must be printed |
| by the function call in argument. */ |
| polybench_prevent_dce(print_array(out_StrictFP)); |
| #endif |
| |
| /* Be clean. */ |
| POLYBENCH_FREE_ARRAY(sum_c); |
| POLYBENCH_FREE_ARRAY(c); |
| POLYBENCH_FREE_ARRAY(W); |
| |
| return 0; |
| } |