31 for (i = 0; i < n; i++)
32 for (j = 0; j < n; j++)
50 for (i = 0; i < n; i++)
51 for (j = 0; j < n; j++) {
67 void kernel_adi(
int tsteps,
int n,
86 mul1 = B1 * DT / (DX * DX);
87 mul2 = B2 * DT / (DY * DY);
98 for (i=1; i<
_PB_N-1; i++) {
102 for (j=1; j<_PB_N-1; j++) {
103 p[i][j] = -c / (a*p[i][j-1]+b);
104 q[i][j] = (-d*u[j][i-1]+(
SCALAR_VAL(1.0)+
SCALAR_VAL(2.0)*d)*u[j][i] - f*u[j][i+1]-a*q[i][j-1])/(a*p[i][j-1]+b);
108 for (j=_PB_N-2; j>=1; j--) {
109 v[j][i] = p[i][j] * v[j+1][i] + q[i][j];
113 for (i=1; i<_PB_N-1; i++) {
117 for (j=1; j<_PB_N-1; j++) {
118 p[i][j] = -f / (d*p[i][j-1]+e);
119 q[i][j] = (-a*v[i-1][j]+(
SCALAR_VAL(1.0)+
SCALAR_VAL(2.0)*a)*v[i][j] - c*v[i+1][j]-d*q[i][j-1])/(d*p[i][j-1]+e);
122 for (j=_PB_N-2; j>=1; j--) {
123 u[i][j] = p[i][j] * u[i][j+1] + q[i][j];
131 int main(
int argc,
char** argv)
#define POLYBENCH_ARRAY(x)
#define POLYBENCH_DUMP_BEGIN(s)
#define POLYBENCH_FREE_ARRAY(x)
#define POLYBENCH_2D(var, dim1, dim2, ddim1, ddim2)
int main(int argc, char **argv)
#define POLYBENCH_DUMP_START
#define POLYBENCH_2D_ARRAY_DECL(var, type, dim1, dim2, ddim1, ddim2)
__attribute__((noinline))
Convert the given fixedpt number to a decimal string.
#define DATA_PRINTF_MODIFIER
#define polybench_prevent_dce(func)
#define POLYBENCH_DUMP_TARGET
#define POLYBENCH_DUMP_END(s)
static void print_array(int n, DATA_TYPE POLYBENCH_2D(u, N, N, n, n))
#define POLYBENCH_DUMP_FINISH
static void init_array(int n, DATA_TYPE POLYBENCH_2D(u, N, N, n, n))
This version is stamped on May 10, 2016.
#define polybench_stop_instruments
#define polybench_print_instruments
#define polybench_start_instruments