Hi @achodankar, I could output the density and the lattice point’s position using the density functor. But, I am having a problem with the input and output array. Things go fine if I don’t use parallelism but as I increase the ranks I do not get correct values for the output. Can you help me out? It seems that I am not providing correct input, as with more ranks the cuboid number becomes a member.