/beegfs/hackathon/users/eoseret/qaas_runs/171-094-7986/intel/HACCmk/build/HACCmk/src/main.c: 50 - 191
--------------------------------------------------------------------------------

50: {
[...]
65:   rank = 0;
66:   nprocs = 1;
67: 
68:   //count = 327;
69:   count = 3000;
70: 
71:   if ( rank == 0 ) 
72:   {  
73:       printf( "count is set %d\n", count );
74:       printf( "Total MPI ranks %d\n", nprocs );
75:   } 
76: 
77:   if (argc == 2 && strncmp(argv[1], "-s", 2) == 0)
78:       NN = 15000;
79: 
80:   printf( "N is set %ld\n", NN );
81: 
82: #pragma omp parallel
[...]
97:   for ( n = 400; n < NN; n = n + 20 ) 
98:   {
99:       /* Initial data preparation */
100:       fcoeff = 0.23f;  
101:       fsrrmax2 = 0.5f; 
102:       mp_rsm2 = 0.03f;
103:       dx1 = 1.0f/(float)n;
104:       dy1 = 2.0f/(float)n;
105:       dz1 = 3.0f/(float)n;
106:       xx[0] = 0.f;
107:       yy[0] = 0.f;
108:       zz[0] = 0.f;
109:       mass[0] = 2.f;
110:       
111:       for ( i = 1; i < n; i++ )
112:       {
113:           xx[i] = xx[i-1] + dx1;
114:           yy[i] = yy[i-1] + dy1;
115:           zz[i] = zz[i-1] + dz1;
116:           mass[i] = (float)i * 0.01f + xx[i];
117:       }
118:     
119:       for ( i = 0; i < n; i++ )
120:       {
121:           vx1[i] = 0.f;
122:           vy1[i] = 0.f;
123:           vz1[i] = 0.f;
[...]
136:       t1 = mysecond();
137: #endif
138: 
139:     #pragma omp parallel for private( dx1, dy1, dz1 )
[...]
152:       t2 = mysecond();
[...]
166:       t3 = (t2 - t1) * 1e6;
167: #endif 
168: 
169:       elapsed = elapsed + t3;
[...]
183:   if ( rank == 0 )
184:   {
185:       printf( "\nKernel elapsed time, s: %18.8lf\n", elapsed*1e-6 );
[...]
191:   return 0;
