Abstraction of sphinxsys range for parallel computing

[Xiangyu-Hu]

I am thinking that we can arrange all the loop ranges in sphinxsys together. Depending on execution policy, they will give the right way of looping. One important part of this to define right upper bound of the range. The other is the mapping fron the compute item to particle index. For device computing, this information will be on device too.