Bug in compute nparticles/tracer/region

Submitted by Christoph on Fri, 01/29/2016 - 12:47

Hi,

It seems to me that the mass calculation of marked particles with the compute nparticles/Tracer/Region is not working correct.
The reason therefore may be that the marker of the Particle will be reset by the first call of compute_vector_eval.
As the Method compute_vector calls the compute_vector_eval twice, it have to be ensured that the marker will not be reset by the first call.

The following modification worked for me:
compute_nparticles_tracer_region.cpp 253
replace if(reset_marker_)
by if(countMass && reset_marker_)

BR,
Christoph C.

ckloss's picture

ckloss | Mon, 02/08/2016 - 12:54

Hi Christoph,

thanks for the note! We'll look into it before the next release!

Best wishes
Christoph

ckloss's picture

ckloss | Fri, 05/13/2016 - 15:08

Hi Christoph,

thanks for the note! Will be fixed in the next release!

best wishes
Christoph