Store velocities in set_rsd

[florpi]

It'd be great if the 3-D velocity field could be stored together with the real and redshift space velocity mesh. I'm afraid I don't really understand the code well enough as to extract it (and I guess only the line of sight is being computed here), but that functionality would be useful to generate mocks with velocities

[adematti]

Hi, sorry for the delay. You mean, in LagrangianLinearMock.set_rsd, here https://github.com/adematti/mockfactory/blob/116610f940f9f553a57daf53c31db8624180c85e/mockfactory/lagrangian_mock.py#L133?
Zeldovich displacements are stored as the disps attribute; I think the conversion to (physical) velocities is just 100 * f * a * E(z) * disps, with E(z) = H(z)/H(0).

[florpi]

Thanks for your answer! I was looking at the Eulerian mock, and it doesn't seem to have a disps attribute?