can we just take this one: > for (Particles::iterator it = ps.begin(); it != ps.end();it++){ > core::XYZ d(*it); > d.set_coordinates(t.transform(d.get_coordinates()); > } and call it transform, something like: void transform(XYZs, Transformation3D)