goad 1.1.9

Compute the single scattering properties of particles much larger than the wavelength of light with geometric optics and aperture diffraction theory.
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
127

use goad::clip::Clipping;
use goad::geom::{Geom, Face};
use nalgebra::{Point3, Vector3};

fn main() {
    println!("=== projection-debug (hex.obj, face 4, projection (0,0,-1)) ===");
    test_projection_debug();
    
    println!("\n=== projection1 (concave1.obj, face 4, projection (-0.3,0,-1)) ===");
    test_projection1();
    
    println!("\n=== projection2 (cube_inside_ico.obj) ===");
    test_projection2();
    
    println!("\n=== projection_multi (multiple.obj, shape[0].face[5], projection (-1,0,0)) ===");
    test_projection_multi();
    
    println!("\n=== clip_test (clip_test.obj) ===");
    test_clip_test();
    
    println!("\n=== remainder (multiple.obj, custom rect, projection (0,0,-1)) ===");
    test_remainder();
}

fn test_projection_debug() {
    let geoms = Geom::load("./examples/data/hex.obj").unwrap();
    let mut geom = geoms[0].clone();
    let mut clip = geom.shapes[0].faces.remove(4);
    let projection = Vector3::new(0.0, 0.0, -1.0);
    
    let mut clipping = Clipping::new(&mut geom, &mut clip, &projection);
    let _ = clipping.clip(0.01);
    
    print_vertices(&clipping);
}

fn test_projection1() {
    let geoms = Geom::load("./examples/data/concave1.obj").unwrap();
    let mut geom = geoms[0].clone();
    let mut clip = geom.shapes[0].faces.remove(4);
    let projection = Vector3::new(-0.3, 0.0, -1.0);
    
    let mut clipping = Clipping::new(&mut geom, &mut clip, &projection);
    let _ = clipping.clip(0.01);
    
    print_vertices(&clipping);
}

fn test_projection2() {
    let geoms = Geom::load("./examples/data/cube_inside_ico.obj").unwrap();
    let mut geom = geoms[0].clone();
    let mut clip = geom.shapes[0].faces.remove(5);
    let projection = Vector3::new(-0.2, 0.0, -1.0);
    
    let mut clipping = Clipping::new(&mut geom, &mut clip, &projection);
    let _ = clipping.clip(0.01);
    
    print_vertices(&clipping);
}

fn test_projection_multi() {
    let geoms = Geom::load("./examples/data/multiple.obj").unwrap();
    let mut geom = geoms[0].clone();
    let mut clip = geom.shapes[0].faces.remove(5);
    let projection = Vector3::new(-1.0, 0.0, 0.0);
    
    let mut clipping = Clipping::new(&mut geom, &mut clip, &projection);
    let _ = clipping.clip(0.01);
    
    print_vertices(&clipping);
}

fn test_clip_test() {
    let geoms = Geom::load("./examples/data/clip_test.obj").unwrap();
    let mut geom = geoms[0].clone();
    let mut clip = geom.shapes[1].faces.remove(1);
    let projection = Vector3::new(1.0, 1.0, 0.0);
    
    let mut clipping = Clipping::new(&mut geom, &mut clip, &projection);
    let _ = clipping.clip(0.01);
    
    print_vertices(&clipping);
}

fn test_remainder() {
    let geoms = Geom::load("./examples/data/multiple.obj").unwrap();
    let mut geom = geoms[0].clone();
    let projection = Vector3::new(0.0, 0.0, -1.0);
    
    let mut clip_vertices = vec![
        Point3::new(-19.0, 3.0, 10.0),
        Point3::new(-19.0, -3.0, 10.0),
        Point3::new(10.0, -3.0, 10.0),
        Point3::new(10.0, 3.0, 10.0),
    ];
    clip_vertices.reverse();
    let mut clip = Face::new_simple(clip_vertices, None, None).unwrap();
    
    let mut clipping = Clipping::new(&mut geom, &mut clip, &projection);
    let _ = clipping.clip(0.01);
    
    print_vertices(&clipping);
}

fn print_vertices(clipping: &Clipping) {
    for (i, intsn) in clipping.intersections.iter().enumerate() {
        let data = intsn.data();
        println!("intersection[{}]:", i);
        println!("  midpoint: ({:.6}, {:.6}, {:.6})", data.midpoint.x, data.midpoint.y, data.midpoint.z);
        println!("  normal: ({:.6}, {:.6}, {:.6})", data.normal.x, data.normal.y, data.normal.z);
        println!("  num_vertices: {}", data.num_vertices);
        println!("  vertices:");
        for (j, v) in data.exterior.iter().enumerate() {
            println!("    [{}]: ({:.6}, {:.6}, {:.6})", j, v.x, v.y, v.z);
        }
    }
    for (i, rem) in clipping.remaining.iter().enumerate() {
        let data = rem.data();
        println!("remaining[{}]:", i);
        println!("  midpoint: ({:.6}, {:.6}, {:.6})", data.midpoint.x, data.midpoint.y, data.midpoint.z);
        println!("  num_vertices: {}", data.num_vertices);
        println!("  vertices:");
        for (j, v) in data.exterior.iter().enumerate() {
            println!("    [{}]: ({:.6}, {:.6}, {:.6})", j, v.x, v.y, v.z);
        }
    }
}