parallelization not yet working, use set type of particle to the corresponding thread number, paraview helps for debuggin