slosh3d 0.1.0

Cross-platform GPU 3D Material Point Method implementation.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65

use crate::solver::{GpuRigidParticles, GpuSampleIds};
use nexus::dynamics::GpuBodySet;
use nexus::math::{GpuSim, Point};
use slang_hal::backend::Backend;
use slang_hal::function::GpuFunction;
use slang_hal::{Shader, ShaderArgs};
use stensor::tensor::GpuTensor;

#[derive(Shader)]
#[shader(module = "slosh::solver::rigid_particle_update")]
pub struct WgRigidParticleUpdate<B: Backend> {
    pub transform_sample_points: GpuFunction<B>,
    pub transform_shape_points: GpuFunction<B>,
}

#[derive(ShaderArgs)]
struct RigidParticleUpdateArgs<'a, B: Backend> {
    vertex_collider_ids: Option<&'a GpuTensor<u32, B>>,
    rigid_particle_indices: Option<&'a GpuTensor<GpuSampleIds, B>>,
    poses: &'a GpuTensor<GpuSim, B>,
    local_pts: &'a GpuTensor<Point<f32>, B>,
    world_pts: &'a GpuTensor<Point<f32>, B>,
}

impl<B: Backend> WgRigidParticleUpdate<B> {
    pub fn launch(
        &self,
        backend: &B,
        pass: &mut B::Pass,
        bodies: &GpuBodySet<B>,
        rigid_particles: &GpuRigidParticles<B>,
    ) -> Result<(), B::Error> {
        if rigid_particles.sample_points.is_empty() {
            return Ok(());
        }

        let args = RigidParticleUpdateArgs {
            vertex_collider_ids: None,
            rigid_particle_indices: Some(&rigid_particles.sample_ids),
            poses: bodies.poses(),
            local_pts: &rigid_particles.local_sample_points,
            world_pts: &rigid_particles.sample_points,
        };
        self.transform_sample_points.launch(
            backend,
            pass,
            &args,
            [rigid_particles.local_sample_points.len() as u32, 1, 1],
        )?;

        let args = RigidParticleUpdateArgs {
            vertex_collider_ids: Some(bodies.shapes_vertex_collider_id()),
            rigid_particle_indices: None,
            poses: bodies.poses(),
            local_pts: bodies.shapes_local_vertex_buffers(),
            world_pts: bodies.shapes_vertex_buffers(),
        };
        self.transform_shape_points.launch(
            backend,
            pass,
            &args,
            [bodies.shapes_vertex_buffers().len() as u32, 1, 1],
        )
    }
}