imlet 0.2.0

A lightweight engine for implicit modeling.
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

use imlet::types::{
    computation::{
        data::sampler::{DenseSampler, Sampler},
        functions::XYZValue,
        model::ImplicitModel,
        operations::math::{Divide, Multiply, Subtract},
        ModelError,
    },
    geometry::{BoundingBox, Vec3},
};
use imlet::utils;

/// Create a new implicit model representing a hyperbolic paraboloid surface.
fn create_hyperbolic_paraboliod(a: f32, b: f32) -> Result<ImplicitModel<f32>, ModelError> {
    let mut model = ImplicitModel::new();

    // Coordinate values
    let x_coord = model.add_function("x_coord", XYZValue::x()).unwrap();
    let y_coord = model.add_function("y_coord", XYZValue::y()).unwrap();
    let z_coord = model.add_function("z_coord", XYZValue::z()).unwrap();

    // Constants
    let a_sqrt = model.add_constant("a", a.powi(2)).unwrap();
    let b_sqrt = model.add_constant("b", b.powi(2)).unwrap();

    // X and Y terms
    let x_sqrt = model.add_operation("x_sqrt", Multiply::new(), Some(&[&x_coord, &x_coord]))?;
    let y_sqrt = model.add_operation("y_sqrt", Multiply::new(), Some(&[&y_coord, &y_coord]))?;

    let x_term = model.add_operation("x_term", Divide::new(), Some(&[&x_sqrt, &a_sqrt]))?;
    let y_term = model.add_operation("y_term", Divide::new(), Some(&[&y_sqrt, &b_sqrt]))?;

    // Final iso value
    let sub_1 = model.add_operation("sub_1", Subtract::new(), Some(&[&y_term, &x_term]))?;

    let _ = model.add_operation("sub_2", Subtract::new(), Some(&[&z_coord, &sub_1]))?;

    Ok(model)
}

pub fn main() {
    utils::logging::init_info();

    let size: f32 = 5.0;
    let cell_size = 0.1;
    let bounds = BoundingBox::new(Vec3::new(-size, -size, -size), Vec3::new(size, size, size));

    let a: f32 = 2.0;
    let b: f32 = 3.0;
    let model = create_hyperbolic_paraboliod(a, b).unwrap();

    let mut sampler = DenseSampler::builder()
        .with_bounds(bounds)
        .with_cell_size(cell_size)
        .build()
        .unwrap();

    sampler.sample_field(&model).expect("Sampling should work.");

    let mesh = sampler
        .iso_surface(0.0)
        .expect("Extracting iso-surface should work.");

    utils::io::write_obj_file(&mesh, "gyroid_example").unwrap();

    #[cfg(feature = "viewer")]
    {
        imlet::viewer::show_mesh(&mesh, Some(bounds));
    }
}