Aggregation solver volume conservation

[TiPlath]

Dear Dyssol developers,

I recently tried to write a script to compare the cell average technique with my quadrature method of moments and I tripped over something which I wanted to share with you. It seems like the cell average technique in the way it is implemented in Dyssol does not conserve the volume of the density function, see my attached MATLAB scripts. (I have this somewhere also in C++, but cannot find it now, sorry. Will update this as soon as I find it).
For growth the cell average works well, but for aggregation the first moment (M1) deviates from an analytical solution for constant aggregation. I tried to fix this, but had zero success. So I decided to just flag this issue here to let you know in case you did not. I would be highly interested in the solution of this if you manage to find it. I do not exclude the possibility that there is an error on my part and I am just too stupid, feel free to hint me to this as well ;)

CellAverageTechnique.zip

kind regards
Timo