From 38d2ea643bb0308c020690b5ef6df26ad41d8dc0 Mon Sep 17 00:00:00 2001 From: Pedro Gonnet <pedro.gonnet@durham.ac.uk> Date: Sat, 13 Apr 2013 16:44:04 +0000 Subject: [PATCH] found sneaky little bug... Former-commit-id: 1bfe0180f5cd559a823a106d1701fe9d52c03803 --- src/runner_iact.h | 9 +-------- 1 file changed, 1 insertion(+), 8 deletions(-) diff --git a/src/runner_iact.h b/src/runner_iact.h index 551e0bd9c6..f3e0af15d6 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 - } /** -- GitLab