50 #define N      4000000  macro
64 double a[N + OFFSET], b[N + OFFSET], c[N + OFFSET];
73 …uble bytes[4]   = {2 * sizeof(double) * N, 2 * sizeof(double) * N, 3 * sizeof(double) * N, 3 * siz…  in main()
93   PetscPrintf(PETSC_COMM_WORLD, "Array size = %d, Offset = %d\n", N, OFFSET);  in main()
94 …PetscPrintf(PETSC_COMM_WORLD, "Total memory required = %.1f MB per process.\n", (3 * N * BytesPerW…  in main()
102     posix_memalign((void **)&a, 64, N * sizeof(double));  in main()
103     posix_memalign((void **)&b, 64, N * sizeof(double));  in main()
104     posix_memalign((void **)&c, 64, N * sizeof(double));  in main()
106     a = malloc(N * sizeof(double));  in main()
107     b = malloc(N * sizeof(double));  in main()
108     c = malloc(N * sizeof(double));  in main()
111     a = numa_alloc_onnode(N * sizeof(double), node);  in main()
112     b = numa_alloc_onnode(N * sizeof(double), node);  in main()
113     c = numa_alloc_onnode(N * sizeof(double), node);  in main()
118   for (j = 0; j < N; j++) {  in main()
124   for (j = 0; j < N; j++) a[j] = 1.0;  in main()
125   for (j = 0; j < N; j++) b[j] = 2.0;  in main()
126   for (j = 0; j < N; j++) c[j] = 0.0;  in main()
135   for (j = 0; j < N; j++) a[j] = 2.0E0 * a[j];  in main()
159     memcpy(c, a, N * sizeof(double));  in main()
161     for (j = 0; j < N; j += 8) {  in main()
171     for (j = 0; j < N; j++) c[j] = a[j];  in main()
182       for (j = 0; j < N; j += 8) {  in main()
193     for (j = 0; j < N; j++) b[j] = scalar * c[j];  in main()
203       for (j = 0; j < N; j += 8) {  in main()
215     for (j = 0; j < N; j++) c[j] = a[j] + b[j];  in main()
226       for (j = 0; j < N; j += 8) {  in main()
238     for (j = 0; j < N; j++) a[j] = b[j] + scalar * c[j];  in main()