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