/beegfs/hackathon/users/eoseret/qaas_runs_test/178-654-5529/intel/HACCmk/build/HACCmk/src/Step10_orig.c: 10 - 41
--------------------------------------------------------------------------------

10: {
[...]
19:     for ( j = 0; j < count1; j++ ) 
20:     {
21:         dxc = xx1[j] - xxi;
22:         dyc = yy1[j] - yyi;
23:         dzc = zz1[j] - zzi;
24:   
25:         r2 = dxc * dxc + dyc * dyc + dzc * dzc;
26:        
27:         m = ( r2 < fsrrmax2 ) ? mass1[j] : 0.0f;
28: 
29:         f =  pow( r2 + mp_rsm2, -1.5 ) - ( ma0 + r2*(ma1 + r2*(ma2 + r2*(ma3 + r2*(ma4 + r2*ma5)))));
30:         
31:         f = ( r2 > 0.0f ) ? m * f : 0.0f;
32: 
33:         xi = xi + f * dxc;
34:         yi = yi + f * dyc;
35:         zi = zi + f * dzc;
36:     }
37: 
38:     *dxi = xi;
39:     *dyi = yi;
40:     *dzi = zi;
41: }
