oxiphysics_geometry/implicit_geometry/
sdftriangularprism_traits.rs1use super::functions::*;
12use super::types::SdfTriangularPrism;
13
14impl Sdf for SdfTriangularPrism {
15 fn dist(&self, p: [f64; 3]) -> f64 {
16 let h = [self.side, self.half_height];
17 let q = [p[0].abs(), p[1], p[2].abs()];
18 let k = 3.0_f64.sqrt();
19 let qx = q[0] - clamp_f(q[0] - k * q[2], 0.0, 1.0) * h[0] / 2.0;
20 let qz = q[2] - clamp_f(q[0] - k * q[2], 0.0, 1.0) * h[0] / (2.0 * k);
21 let d1 = len2([qx, qz])
22 * if q[0] - h[0] / 2.0 < 0.0 && k * q[2] - q[0] < 0.0 {
23 -1.0
24 } else {
25 1.0
26 };
27 let d2 = q[1] - h[1];
28 -((-d1).max(-d2)).min(0.0) + len2([d1.max(0.0), d2.max(0.0)])
29 }
30}