fast image operations
| -rw-r--r-- | src/drawing/tri.rs | 23 | ||||
| -rw-r--r-- | src/math.rs | 2 |
2 files changed, 10 insertions, 15 deletions
diff --git a/src/drawing/tri.rs b/src/drawing/tri.rs index ca44509..1333abe 100644 --- a/src/drawing/tri.rs +++ b/src/drawing/tri.rs @@ -1,9 +1,7 @@ //! trongle drawing - use crate::math::madd; -use std::cmp::{max, min}; - use crate::Image; +use std::cmp::{max, min}; impl<T: AsMut<[u8]> + AsRef<[u8]>, const CHANNELS: usize> Image<T, CHANNELS> { /// Draw a (filled) triangle @@ -17,31 +15,26 @@ impl<T: AsMut<[u8]> + AsRef<[u8]>, const CHANNELS: usize> Image<T, CHANNELS> { /// (1.0, 8.0), // point c /// [255] // white /// ); - /// # assert_eq!(a.buffer(), b"\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\xff\xff\x00\x00\x00\x00\x00\x00\x00\x00\xff\xff\xff\x00\x00\x00\x00\x00\x00\x00\xff\xff\xff\xff\x00\x00\x00\x00\x00\xff\xff\xff\xff\xff\xff\x00\x00\x00\x00\xff\xff\xff\xff\xff\xff\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00"); + /// # assert_eq!(a.buffer(), b"\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\xff\x00\x00\x00\x00\x00\x00\x00\x00\x00\xff\xff\x00\x00\x00\x00\x00\x00\x00\x00\xff\xff\xff\x00\x00\x00\x00\x00\x00\xff\xff\xff\xff\xff\x00\x00\x00\x00\x00\xff\xff\xff\xff\xff\xff\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00\x00"); /// ``` pub fn tri( &mut self, - (x1, y1): (f32, f32), (x2, y2): (f32, f32), + (x1, y1): (f32, f32), (x3, y3): (f32, f32), c: [u8; CHANNELS], ) { - // TODO optimize let ymin = max(y1.min(y2).min(y3) as u32, 0); let ymax = min(y1.max(y2).max(y3) as u32, self.height()); let xmin = max(x1.min(x2).min(x3) as u32, 0); let xmax = min(x1.max(x2).max(x3) as u32, self.width()); for y in ymin..ymax { for x in xmin..xmax { - let s = madd(x1 - x3, y as f32 - y3, -(y1 - y2) * (x as f32 - x3)); - let t = madd(x2 - x1, y as f32 - y1, -(y2 - y1) * (x as f32 - x1)); - - if (s < 0.0) != (t < 0.0) && s != 0.0 && t != 0.0 { - continue; - } - - let d = madd(x3 - x2, y as f32 - y2, -(y3 - y2) * (x as f32 - x2)); - if d == 0.0 || (d < 0.0) == (s + t <= 0.0) { + // algorithm from https://web.archive.org/web/20050408192410/http://sw-shader.sourceforge.net/rasterizer.html, but im too dumb to implement the faster ones + if madd(x1 - x2, y as f32 - y1, -(y1 - y2) * (x as f32 - x1)) > 0. + && madd(x2 - x3, y as f32 - y2, -(y2 - y3) * (x as f32 - x2)) > 0. + && madd(x3 - x1, y as f32 - y3, -(y3 - y1) * (x as f32 - x3)) > 0. + { // SAFETY: x, y are bounded unsafe { self.set_pixel(x, y, c) }; } diff --git a/src/math.rs b/src/math.rs index 212798c..6ff5572 100644 --- a/src/math.rs +++ b/src/math.rs @@ -1,4 +1,6 @@ +//! utility math /// Calculates `a * b + c`, with hardware support if possible. +#[allow(clippy::suboptimal_flops)] pub fn madd(a: f32, b: f32, c: f32) -> f32 { if cfg!(target_feature = "fma") { a.mul_add(b, c) |