oxiphysics_geometry/csg/
sdfbox_traits.rs1use super::functions::ImplicitSurface;
12use super::functions::{length, normalize};
13use super::types::SdfBox;
14
15impl ImplicitSurface for SdfBox {
16 fn sdf(&self, p: [f64; 3]) -> f64 {
17 let d = [
18 (p[0] - self.center[0]).abs() - self.half_extents[0],
19 (p[1] - self.center[1]).abs() - self.half_extents[1],
20 (p[2] - self.center[2]).abs() - self.half_extents[2],
21 ];
22 let exterior = length([d[0].max(0.0), d[1].max(0.0), d[2].max(0.0)]);
23 let interior = d[0].max(d[1]).max(d[2]).min(0.0);
24 exterior + interior
25 }
26 fn gradient(&self, p: [f64; 3]) -> [f64; 3] {
27 const EPS: f64 = 1e-5;
28 let dx = self.sdf([p[0] + EPS, p[1], p[2]]) - self.sdf([p[0] - EPS, p[1], p[2]]);
29 let dy = self.sdf([p[0], p[1] + EPS, p[2]]) - self.sdf([p[0], p[1] - EPS, p[2]]);
30 let dz = self.sdf([p[0], p[1], p[2] + EPS]) - self.sdf([p[0], p[1], p[2] - EPS]);
31 normalize([dx, dy, dz])
32 }
33}