diff --git a/src/runner_iact.h b/src/runner_iact.h index 551e0bd9c61126b562de4c35a675441c2221df21..f3e0af15d61cfdc70e09e1f6fac11ecda5af6da6 100644 --- a/src/runner_iact.h +++ b/src/runner_iact.h @@ -91,19 +91,12 @@ __attribute__ ((always_inline)) INLINE static void runner_iact_density ( float r pj->density.wcount_dh -= xj * h_inv * wj_dx * ( 4.0f * M_PI / 3.0f * kernel_igamma3 ); // pj->icount += 1; - pj->density.div_v = pi->mass * dvdr * wj_dx; + pj->density.div_v += pi->mass * dvdr * wj_dx; for ( k = 0 ; k < 3 ; k++ ) pj->density.curl_v[k] += pi->mass * curlvr[k] * wj_dx; } - #ifdef HIST - if ( hi > hj ) - runner_hist_hit( hi / hj ); - else - runner_hist_hit( hj / hi ); - #endif - } /**