Skip to content
Snippets Groups Projects
Commit 38d2ea64 authored by Pedro Gonnet's avatar Pedro Gonnet
Browse files

found sneaky little bug... 

Former-commit-id: 1bfe0180f5cd559a823a106d1701fe9d52c03803
parent 44490d9b
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 1 addition and 8 deletions
src/runner_iact.h
+ 1
− 8
View file @ 38d2ea64
...@@ -91,19 +91,12 @@ __attribute__ ((always_inline)) INLINE static void runner_iact_density ( float r ...@@ -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->density.wcount_dh -= xj * h_inv * wj_dx * ( 4.0f * M_PI / 3.0f * kernel_igamma3 );
// pj->icount += 1; // 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++ ) for ( k = 0 ; k < 3 ; k++ )
pj->density.curl_v[k] += pi->mass * curlvr[k] * wj_dx; 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
} }
/** /**
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment