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