oxiphysics_geometry/csg/
csgsmoothdifference_traits.rs1use super::functions::ImplicitSurface;
12#[allow(unused_imports)]
13use super::functions::*;
14use super::functions::{clamp, mix, normalize};
15use super::types::CsgSmoothDifference;
16
17impl ImplicitSurface for CsgSmoothDifference {
18 fn sdf(&self, p: [f64; 3]) -> f64 {
19 let da = self.a.sdf(p);
20 let db = -self.b.sdf(p);
21 let h = clamp(0.5 - 0.5 * (db - da) / self.k, 0.0, 1.0);
22 mix(db, da, h) + self.k * h * (1.0 - h)
23 }
24 fn gradient(&self, p: [f64; 3]) -> [f64; 3] {
25 const EPS: f64 = 1e-5;
26 let dx = self.sdf([p[0] + EPS, p[1], p[2]]) - self.sdf([p[0] - EPS, p[1], p[2]]);
27 let dy = self.sdf([p[0], p[1] + EPS, p[2]]) - self.sdf([p[0], p[1] - EPS, p[2]]);
28 let dz = self.sdf([p[0], p[1], p[2] + EPS]) - self.sdf([p[0], p[1], p[2] - EPS]);
29 normalize([dx, dy, dz])
30 }
31}