fast image operations
Diffstat (limited to 'src/overlay.rs')
-rw-r--r--src/overlay.rs82
1 files changed, 21 insertions, 61 deletions
diff --git a/src/overlay.rs b/src/overlay.rs
index 3dbc45e..d1211a9 100644
--- a/src/overlay.rs
+++ b/src/overlay.rs
@@ -1,10 +1,8 @@
-//! Handles image overlay
use super::{assert_unchecked, really_unsafe_index, Image};
use std::simd::SimdInt;
use std::simd::SimdPartialOrd;
use std::simd::{simd_swizzle, Simd};
-/// Trait for layering a image ontop of another, with a offset to the second image.
pub trait OverlayAt<W> {
/// Overlay with => self at coordinates x, y, without blending
/// # Safety
@@ -12,8 +10,7 @@ pub trait OverlayAt<W> {
/// UB if x, y is out of bounds
unsafe fn overlay_at(&mut self, with: &W, x: u32, y: u32) -> &mut Self;
}
-/// Trait for layering images ontop of each other.
-/// Think `magick a b -layers flatten a`
+
pub trait Overlay<W> {
/// Overlay with => self (does not blend)
/// # Safety
@@ -23,20 +20,15 @@ pub trait Overlay<W> {
}
#[inline]
-/// SIMD accelerated rgba => rgb overlay.
-///
-/// See [blit](https://en.wikipedia.org/wiki/Bit_blit)
-///
-/// # Safety
-/// - UB if rgb.len() % 3 != 0
-/// - UB if rgba.len() % 4 != 0
unsafe fn blit(rgb: &mut [u8], rgba: &[u8]) {
+ const LAST4: Simd<u8, 16> = Simd::from_array([
+ 255, 255, 255, 255, 255, 255, 255, 255, 255, 255, 255, 255, 0, 0, 0, 0,
+ ]);
+
let mut srci = 0;
let mut dsti = 0;
while dsti + 16 <= rgb.len() {
- // SAFETY: i think it ok
let old: Simd<u8, 16> = Simd::from_slice(unsafe { rgb.get_unchecked(dsti..dsti + 16) });
- // SAFETY: definetly ok
let new: Simd<u8, 16> = Simd::from_slice(unsafe { rgba.get_unchecked(srci..srci + 16) });
let threshold = new.simd_ge(Simd::splat(128)).to_int().cast::<u8>();
@@ -44,13 +36,10 @@ unsafe fn blit(rgb: &mut [u8], rgba: &[u8]) {
threshold,
[3, 3, 3, 7, 7, 7, 11, 11, 11, 15, 15, 15, 0, 0, 0, 0]
);
- mask &= Simd::from_array([
- 255, 255, 255, 255, 255, 255, 255, 255, 255, 255, 255, 255, 0, 0, 0, 0,
- ]);
+ mask &= LAST4;
let new_rgb = simd_swizzle!(new, [0, 1, 2, 4, 5, 6, 8, 9, 10, 12, 13, 14, 0, 0, 0, 0]);
let blended = (new_rgb & mask) | (old & !mask);
- // SAFETY: 4 * 4 == 16, so in bounds
blended.copy_to_slice(unsafe { rgb.get_unchecked_mut(dsti..dsti + 16) });
srci += 16;
@@ -58,13 +47,10 @@ unsafe fn blit(rgb: &mut [u8], rgba: &[u8]) {
}
while dsti + 3 <= rgb.len() {
- // SAFETY: caller gurantees slice is big enough
if unsafe { *rgba.get_unchecked(srci + 3) } >= 128 {
- // SAFETY: slice is big enough!
- let src = unsafe { rgba.get_unchecked(srci..=srci + 2) };
- // SAFETY: i hear it bound
- let end = unsafe { rgb.get_unchecked_mut(dsti..=dsti + 2) };
- end.copy_from_slice(src);
+ let src = unsafe { rgba.get_unchecked(srci..srci + 3) };
+ let end = unsafe { rgb.get_unchecked_mut(dsti..dsti + 3) };
+ unsafe { std::ptr::copy_nonoverlapping(src.as_ptr(), end.as_mut_ptr(), 3) };
}
srci += 4;
@@ -79,9 +65,16 @@ impl Overlay<Image<&[u8], 4>> for Image<&mut [u8], 4> {
debug_assert!(self.height() == with.height());
for (i, other_pixels) in with.chunked().enumerate() {
if other_pixels[3] >= 128 {
- // SAFETY: outside are bounds of index from slice
- let own_pixels = unsafe { self.buffer.get_unchecked_mut(i * 4..i * 4 + 4) };
- own_pixels.copy_from_slice(other_pixels);
+ let idx_begin = unsafe { i.unchecked_mul(4) };
+ let idx_end = unsafe { idx_begin.unchecked_add(4) };
+ let own_pixels = unsafe { self.buffer.get_unchecked_mut(idx_begin..idx_end) };
+ unsafe {
+ std::ptr::copy_nonoverlapping(
+ other_pixels.as_ptr(),
+ own_pixels.as_mut_ptr(),
+ 4,
+ );
+ };
}
}
self
@@ -91,9 +84,9 @@ impl Overlay<Image<&[u8], 4>> for Image<&mut [u8], 4> {
impl OverlayAt<Image<&[u8], 4>> for Image<&mut [u8], 3> {
#[inline]
unsafe fn overlay_at(&mut self, with: &Image<&[u8], 4>, x: u32, y: u32) -> &mut Self {
- // SAFETY: caller upholds this
+ // SAFETY: caller upholds these
unsafe { assert_unchecked!(x + with.width() <= self.width()) };
- debug_assert!(y + with.height() <= self.height());
+ unsafe { assert_unchecked!(y + with.height() <= self.height()) };
for j in 0..with.height() {
let i_x = j as usize * with.width() as usize * 4
..(j as usize + 1) * with.width() as usize * 4;
@@ -102,11 +95,8 @@ impl OverlayAt<Image<&[u8], 4>> for Image<&mut [u8], 3> {
+ x as usize
+ with.width() as usize)
* 3;
- // SAFETY: index is in bounds
let rgb = unsafe { self.buffer.get_unchecked_mut(o_x) };
- // SAFETY: bounds are outside index
let rgba = unsafe { with.buffer.get_unchecked(i_x) };
- // SAFETY: arguments are 🟢
unsafe { blit(rgb, rgba) }
}
self
@@ -114,16 +104,8 @@ impl OverlayAt<Image<&[u8], 4>> for Image<&mut [u8], 3> {
}
impl OverlayAt<Image<&[u8], 3>> for Image<&mut [u8], 3> {
- /// Overlay a RGB image(with) => self at coordinates x, y.
- /// As this is a `RGBxRGB` operation, blending is unnecessary,
- /// and this is simply a copy.
- ///
- /// # Safety
- ///
- /// UB if x, y is out of bounds
#[inline]
unsafe fn overlay_at(&mut self, with: &Image<&[u8], 3>, x: u32, y: u32) -> &mut Self {
- /// helper macro for defining rgb=>rgb overlays. allows unrolling
macro_rules! o3x3 {
($n:expr) => {{
for j in 0..($n as usize) {
@@ -131,12 +113,7 @@ impl OverlayAt<Image<&[u8], 3>> for Image<&mut [u8], 3> {
let o_x = ((j + y as usize) * self.width() as usize + x as usize) * 3
..((j + y as usize) * self.width() as usize + x as usize + ($n as usize))
* 3;
- // <= because ".." range
- debug_assert!(o_x.end <= self.buffer().len());
- debug_assert!(i_x.end <= with.buffer().len());
- // SAFETY: bounds are ✅
let a = unsafe { self.buffer.get_unchecked_mut(o_x) };
- // SAFETY: we are in ⬜!
let b = unsafe { with.buffer.get_unchecked(i_x) };
a.copy_from_slice(b);
}
@@ -162,13 +139,11 @@ impl Overlay<Image<&[u8], 4>> for Image<&mut [u8], 3> {
.chunks_exact(with.width() as usize * 4)
.enumerate()
{
- // SAFETY: all the bounds are good
let rgb = unsafe {
self.buffer.get_unchecked_mut(
i * with.width() as usize * 3..(i + 1) * with.width() as usize * 3,
)
};
- // SAFETY: we have the rgb and rgba arguments right
unsafe { blit(rgb, chunk) };
}
self
@@ -177,30 +152,15 @@ impl Overlay<Image<&[u8], 4>> for Image<&mut [u8], 3> {
impl OverlayAt<Image<&[u8], 4>> for Image<&mut [u8], 4> {
#[inline]
- /// Overlay with => self at coordinates x, y, without blending
- ///
- /// # Safety
- /// - UB if x, y is out of bounds
- /// - UB if x + with.width() > [`u32::MAX`]
- /// - UB if y + with.height() > [`u32::MAX`]
unsafe fn overlay_at(&mut self, with: &Image<&[u8], 4>, x: u32, y: u32) -> &mut Self {
for j in 0..with.height() {
for i in 0..with.width() {
- // SAFETY: i, j is in bounds.
let index = unsafe { really_unsafe_index(i, j, with.width()) };
- // SAFETY: using .pixel() results in horrible asm (+5k ns/iter)
let their_px = unsafe { with.buffer.get_unchecked(index * 4..index * 4 + 4) };
- // SAFETY: must be sized right
if unsafe { *their_px.get_unchecked(3) } >= 128 {
- // SAFETY:
- // they said it cant go over.
- // i dont know why, but this has performance importanceâ„¢
let x = unsafe { i.unchecked_add(x) };
- // SAFETY: caller gurantees this cannot overflow.
let y = unsafe { j.unchecked_add(y) };
- // SAFETY: compute the offset index.
let index = unsafe { really_unsafe_index(x, y, self.width()) };
- // SAFETY: if everything else goes well, this is fine
let our_px = unsafe { self.buffer.get_unchecked_mut(index * 4..index * 4 + 4) };
our_px.copy_from_slice(their_px);
}