I have a sample of 100 particles with their 3D positions and velocities. I want to decompose the velocity vector field to its curl free and divergence free components.Is it a feasible to do so for the sampled data using Python, Matlab or Mathematica?? 
How to compute the Helmholtz decomposition of 2D and 3D vector fields?
matlab;mathematica;vector