splashsurf_lib 0.14.1

Library for surface reconstruction of SPH particle data
Documentation 
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
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126

use nalgebra::Vector3;
use splashsurf_lib::marching_cubes::check_mesh_consistency;
use splashsurf_lib::{
    Aabb3d, GridDecompositionParameters, Parameters, Real, SpatialDecomposition,
    reconstruct_surface,
};

enum Strategy {
    Global,
    SubdomainGrid,
}

fn params_with_aabb<R: Real>(
    particle_radius: R,
    compact_support_radius: R,
    cube_size: R,
    iso_surface_threshold: R,
    domain_aabb: Option<Aabb3d<R>>,
    strategy: Strategy,
) -> Parameters<R> {
    let compact_support_radius = particle_radius * compact_support_radius;
    let cube_size = particle_radius * cube_size;

    let mut parameters = Parameters {
        particle_radius,
        rest_density: R::from_float(1000.0),
        compact_support_radius,
        cube_size,
        iso_surface_threshold,
        particle_aabb: domain_aabb,
        enable_multi_threading: false,
        enable_simd: false,
        spatial_decomposition: SpatialDecomposition::None,
        global_neighborhood_list: false,
    };

    match strategy {
        Strategy::Global => {}
        Strategy::SubdomainGrid => {
            parameters.spatial_decomposition =
                SpatialDecomposition::UniformGrid(GridDecompositionParameters {
                    subdomain_num_cubes_per_dim: 64,
                    auto_disable: false,
                })
        }
    }

    parameters
}

fn params<R: Real>(
    particle_radius: R,
    compact_support_radius: R,
    cube_size: R,
    iso_surface_threshold: R,
    strategy: Strategy,
) -> Parameters<R> {
    params_with_aabb(
        particle_radius,
        compact_support_radius,
        cube_size,
        iso_surface_threshold,
        None,
        strategy,
    )
}

/// This tests ensures that a surface is reconstructed properly if a single edge is above the threshold
/// on one side but outside the compact support on the other side.
#[test]
fn test_edge_above_threshold_to_outside_of_compact_support_global() {
    let particle_positions: Vec<Vector3<f32>> = vec![Vector3::new(0.01, 0.0, 0.0)];
    let parameters = params(1.0, 1.0, 1.0, 0.1, Strategy::Global);
    let reconstruction =
        reconstruct_surface::<i64, _>(particle_positions.as_slice(), &parameters).unwrap();

    println!(
        "Reconstructed mesh from {} particles has {} triangles.",
        particle_positions.len(),
        reconstruction.mesh.triangles.len()
    );

    assert_eq!(reconstruction.mesh.vertices.len(), 6, "Number of vertices");
    assert_eq!(
        reconstruction.mesh.triangles.len(),
        8,
        "Number of triangles"
    );

    // Ensure that the mesh does not have a boundary
    if let Err(e) =
        check_mesh_consistency(&reconstruction.grid, &reconstruction.mesh, true, true, true)
    {
        eprintln!("{}", e);
        panic!("Mesh contains topological/manifold errors");
    }
}

#[test]
fn test_edge_above_threshold_to_outside_of_compact_support_subdomains() {
    let particle_positions: Vec<Vector3<f32>> = vec![Vector3::new(0.01, 0.0, 0.0)];
    let parameters = params(1.0, 1.0, 1.0, 0.1, Strategy::SubdomainGrid);
    let reconstruction =
        reconstruct_surface::<i64, _>(particle_positions.as_slice(), &parameters).unwrap();

    println!(
        "Reconstructed mesh from {} particles has {} triangles.",
        particle_positions.len(),
        reconstruction.mesh.triangles.len()
    );

    assert_eq!(reconstruction.mesh.vertices.len(), 6, "Number of vertices");
    assert_eq!(
        reconstruction.mesh.triangles.len(),
        8,
        "Number of triangles"
    );

    // Ensure that the mesh does not have a boundary
    if let Err(e) =
        check_mesh_consistency(&reconstruction.grid, &reconstruction.mesh, true, true, true)
    {
        eprintln!("{}", e);
        panic!("Mesh contains topological/manifold errors");
    }
}