Thanks for using Compiler Explorer
Sponsors
Jakt
C++
Ada
Analysis
Android Java
Android Kotlin
Assembly
C
C3
Carbon
C++ (Circle)
CIRCT
Clean
CMake
CMakeScript
COBOL
C++ for OpenCL
MLIR
Cppx
Cppx-Blue
Cppx-Gold
Cpp2-cppfront
Crystal
C#
CUDA C++
D
Dart
Elixir
Erlang
Fortran
F#
Go
Haskell
HLSL
Hook
Hylo
ispc
Java
Julia
Kotlin
LLVM IR
LLVM MIR
Modula-2
Nim
Objective-C
Objective-C++
OCaml
OpenCL C
Pascal
Pony
Python
Racket
Ruby
Rust
Snowball
Scala
Solidity
Spice
Swift
LLVM TableGen
Toit
TypeScript Native
V
Vala
Visual Basic
Zig
Javascript
GIMPLE
rust source #1
Output
Compile to binary object
Link to binary
Execute the code
Intel asm syntax
Demangle identifiers
Verbose demangling
Filters
Unused labels
Library functions
Directives
Comments
Horizontal whitespace
Debug intrinsics
Compiler
BPF gcc (trunk)
mrustc (master)
rustc 1.0.0
rustc 1.1.0
rustc 1.10.0
rustc 1.11.0
rustc 1.12.0
rustc 1.13.0
rustc 1.14.0
rustc 1.15.1
rustc 1.16.0
rustc 1.17.0
rustc 1.18.0
rustc 1.19.0
rustc 1.2.0
rustc 1.20.0
rustc 1.21.0
rustc 1.22.0
rustc 1.23.0
rustc 1.24.0
rustc 1.25.0
rustc 1.26.0
rustc 1.27.0
rustc 1.27.1
rustc 1.28.0
rustc 1.29.0
rustc 1.3.0
rustc 1.30.0
rustc 1.31.0
rustc 1.32.0
rustc 1.33.0
rustc 1.34.0
rustc 1.35.0
rustc 1.36.0
rustc 1.37.0
rustc 1.38.0
rustc 1.39.0
rustc 1.4.0
rustc 1.40.0
rustc 1.41.0
rustc 1.42.0
rustc 1.43.0
rustc 1.44.0
rustc 1.45.0
rustc 1.45.2
rustc 1.46.0
rustc 1.47.0
rustc 1.48.0
rustc 1.49.0
rustc 1.5.0
rustc 1.50.0
rustc 1.51.0
rustc 1.52.0
rustc 1.53.0
rustc 1.54.0
rustc 1.55.0
rustc 1.56.0
rustc 1.57.0
rustc 1.58.0
rustc 1.59.0
rustc 1.6.0
rustc 1.60.0
rustc 1.61.0
rustc 1.62.0
rustc 1.63.0
rustc 1.64.0
rustc 1.65.0
rustc 1.66.0
rustc 1.67.0
rustc 1.68.0
rustc 1.69.0
rustc 1.7.0
rustc 1.70.0
rustc 1.71.0
rustc 1.72.0
rustc 1.73.0
rustc 1.74.0
rustc 1.75.0
rustc 1.76.0
rustc 1.77.0
rustc 1.8.0
rustc 1.9.0
rustc beta
rustc nightly
rustc-cg-gcc (master)
x86-64 GCCRS (GCC master)
x86-64 GCCRS (GCCRS master)
Options
Source code
use std::ops::{Add, Index, IndexMut, Sub}; pub struct SimpleMat<T> { rows: usize, cols: usize, data: Vec<T>, } impl<T: Copy> SimpleMat<T> { pub fn new(rows: usize, cols: usize, default_value: T) -> SimpleMat<T> { SimpleMat { cols, rows, data: vec![default_value; cols * rows] } } unsafe fn get_unchecked(&self, i: usize, j: usize) -> &T { self.data.get_unchecked(self.calc_index_unchecked(i, j)) } unsafe fn get_unchecked_mut(&mut self, i: usize, j: usize) -> &mut T { let idx = self.calc_index_unchecked(i, j); self.data.get_unchecked_mut(idx) } } impl<T> SimpleMat<T> { fn calc_index(&self, i: usize, j: usize) -> usize { assert!(i < self.rows); assert!(j < self.cols); self.calc_index_unchecked(i, j) } fn calc_index_unchecked(&self, i: usize, j: usize) -> usize { i * self.cols + j } } impl<T> Index<(usize, usize)> for SimpleMat<T> { type Output = T; fn index(&self, (i, j): (usize, usize)) -> &Self::Output { &self.data[self.calc_index(i, j)] } } impl<T> IndexMut<(usize, usize)> for SimpleMat<T> { fn index_mut(&mut self, (i, j): (usize, usize)) -> &mut Self::Output { let idx = self.calc_index(i, j); &mut self.data[idx] } } // fn local_sum_for_integral<T: Copy + Add + Sub>(im_sum: &SimpleMat<T>, x1: usize, x2: usize, y1: usize, y2: usize) -> T { // (im_sum[(y2, x2)] + im_sum[(y1, x1)]) - (im_sum[(y1, x2)] + im_sum[(y2, x1)]) // } unsafe fn local_sum_for_integral_unchecked<T: Copy + Add<T, Output=T> + Sub<T, Output=T>>(im_sum: &SimpleMat<T>, x1: usize, x2: usize, y1: usize, y2: usize) -> T { (*im_sum.get_unchecked(y2, x2) + *im_sum.get_unchecked(y1, x1)) - (*im_sum.get_unchecked(y1, x2) + *im_sum.get_unchecked(y2, x1)) } pub fn calc_mean_and_std_given_sum_and_sum_sq(im_sum: &SimpleMat<i32>, im_sum_sq: &SimpleMat<f64>, sqrt_table: &[u8; 65536], win_x: usize, win_y: usize) -> (SimpleMat<u8>, SimpleMat<u8>, u8) { let mut im_mean: SimpleMat<u8> = SimpleMat::new(im_sum.rows, im_sum.cols, 0); let mut im_std: SimpleMat<u8> = SimpleMat::new(im_sum.rows, im_sum.cols, 0); let std_max: u8 = 0; assert!((win_x % 2 == 0) && (win_y % 2 == 0)); let (win_x_half, win_y_half) = (win_x / 2, win_y / 2); let win_area = win_x * win_y; for y in win_y_half..im_sum.rows - win_y_half { for x in win_x_half..im_sum.cols - win_x_half { unsafe { let (x1, x2) = (x - win_x_half, x + win_x_half); let (y1, y2) = (y - win_y_half, y + win_y_half); let local_sum = local_sum_for_integral_unchecked(im_sum, x1, x2, y1, y2); let local_sum_sq = local_sum_for_integral_unchecked(im_sum_sq, x1, x2, y1, y2); let local_mean = local_sum / win_area as i32; let local_variance = ((local_sum_sq / win_area as f64) as i32 - local_mean * local_mean) as usize; // use a table, instead of directly calculating the sqrt let local_std = sqrt_table[local_variance]; // im_mean[(y, x)] = local_mean as u8; // im_std[(y, x)] = local_std; *im_mean.get_unchecked_mut(y, x) = local_mean as u8; *im_std.get_unchecked_mut(y, x) = local_std; } } } (im_mean, im_std, std_max) }
Become a Patron
Sponsor on GitHub
Donate via PayPal
Source on GitHub
Mailing list
Installed libraries
Wiki
Report an issue
How it works
Contact the author
CE on Mastodon
About the author
Statistics
Changelog
Version tree