Nonsensical magnitudes on number densities

[sigvaldm]

I made a slice in ParaView of number densities from the punc++/test simulation and it ranged from 0 to 44. Certainly this is not physical, since the background number density is 1e10. I believe this must be the density of simulation particles then. It would be better to output the physical number density, which is obtained by multiplying by the statistical weight (i.e. ratio of physical to simulation particles). I am not entirely sure how to do this in a good way. We do not store the statistical weight, and it might be different for different species. For sure, it's just one number, so it could be saved in Species during its construction.

[diakodarian]

Yes, you're right. What is computed by density_cg1() and density_dg0() functions is the density of simulation particles. I have now modified these functions so that they compute the physical number densities instead.

What I have done is to add a new member toSpecies class named weight, which stores the statistical weight for each species. The species vector is passed as argument to density_cg1() and density_dg0(), and the statistical weight for each species is taken into account when the number densities are computed.

[sigvaldm]

Having a member named weight makes sense, but what I do not understand, is how you would apply the right weight for each particle? Do you compare the charge and mass of each particle with the species? I would assume that is terribly slow.