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