func DomainDecompose(d particles.Domainer, s *structvec.StructVec) {
pp := GetArray(s)
localndx, mpirank := d.Domain(pp)
s.RestoreArray()
s.Scatter(localndx, mpirank)
}
func GetArray(s *structvec.StructVec) Arr {
return Arr(s.GetArray().([]One))
}