This procedure computes the Helmholz decomposition of a 3D vector field into its divergence and curl free components.
The input vector is projected into fourier space, and decomposed into modes parallel and perpendicular to k. These are then projected back into the original space, and are equal to the curl- and divergence-free vectors, respectively.
- x in required
The x component of the vector. A 3D array
- y in required
The y component of the vector. A 3D array
- z in required
The z component of the vector. A 3D array
|Modifcation date:||Fri Jul 30 23:11:09 2010|