oxiphysics_gpu/kernels/rigid/
integratepositionkernel_traits.rs1use crate::compute::ComputeKernel;
11
12use super::types::IntegratePositionKernel;
13
14impl ComputeKernel for IntegratePositionKernel {
15 fn name(&self) -> &str {
16 "IntegratePositionKernel"
17 }
18 fn execute(&self, inputs: &[&[f64]], outputs: &mut [Vec<f64>], work_size: usize) {
19 if inputs.len() < 3 || outputs.is_empty() {
20 return;
21 }
22 let pos = inputs[0];
23 let vel = inputs[1];
24 let dt = inputs[2][0];
25 let n = work_size;
26 let mut new_pos = vec![0.0; n * 3];
27 for i in 0..n {
28 new_pos[i * 3] = pos[i * 3] + vel[i * 3] * dt;
29 new_pos[i * 3 + 1] = pos[i * 3 + 1] + vel[i * 3 + 1] * dt;
30 new_pos[i * 3 + 2] = pos[i * 3 + 2] + vel[i * 3 + 2] * dt;
31 }
32 outputs[0] = new_pos;
33 }
34}