|
AQUAgpusph 5.0.4
|
Euler-XYZ based transformation script. More...
#include "resources/Scripts/types/types.h"
Functions | |
| __kernel void | entry (const __global uint *iset, const __global int *imove, __global vec *r, __global vec *normal, __global vec *tangent, usize N, unsigned int motion_iset, vec motion_r, vec4 motion_a) |
| Transform the boundaries applying Euler-XYZ motion. | |
Euler-XYZ based transformation script.
| __kernel void entry | ( | const __global uint * | iset, |
| const __global int * | imove, | ||
| __global vec * | r, | ||
| __global vec * | normal, | ||
| __global vec * | tangent, | ||
| usize | N, | ||
| unsigned int | motion_iset, | ||
| vec | motion_r, | ||
| vec4 | motion_a ) |
Transform the boundaries applying Euler-XYZ motion.
Just the position and the normal of the particle are modified, but not the velocity which is changed by MotionVelocity.cl.
In Euler-XYZ the following transformation is applied to a particle \( a \):
\[ R_z \cdot R_y \cdot R_x \cdot \mathbf{x_a} + \mathbf{cor}, \]
where \( \mathbf{cor} \) is the position of the center of rotation (global translations), and \( R_x, R_y, R_z \) are the rotation matrices:
\[ R_x = \left[ \begin{matrix} 1 & 0 & 0 \\ 0 & \mathrm{cos}(\phi) & -\mathrm{sin}(\phi) \\ 0 & \mathrm{sin}(\phi) & \mathrm{cos}(\phi) \\ \end{matrix} \right], \]
\[ R_y = \left[ \begin{matrix} \mathrm{cos}(\theta) & 0 & \mathrm{sin}(\theta) \\ 0 & 1 & 0 \\ -\mathrm{sin}(\theta) & 0 & \mathrm{cos}(\theta) \\ \end{matrix} \right], \]
\[ R_z = \left[ \begin{matrix} \mathrm{cos}(\psi) & -\mathrm{sin}(\psi) & 0 \\ \mathrm{sin}(\psi) & \mathrm{cos}(\psi) & 0 \\ 0 & 0 & 1 \\ \end{matrix} \right]. \]
| iset | Set of particles index. |
| imove | Moving flags.
|
| r | Position \( \mathbf{r} \). |
| normal | Normal \( \mathbf{n} \). |
| tangent | Tangent \( \mathbf{t} \). |
| N | Number of particles. |
| motion_iset | Set of particles affected. |
| motion_r | Center of rotation. |
| motion_a | Rotation angles \( \phi, \theta, \psi \). |