fast-ssim2 0.8.2

Fast SSIMULACRA2 image quality metric with SIMD acceleration
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128

//! Benchmark comparing blur implementations and full SSIMULACRA2 pipeline
//!
//! Run with:
//!   cargo run --release --example benchmark_simd

#![allow(deprecated)]

use std::time::Instant;

use fast_ssim2::{Blur, SimdImpl, Ssimulacra2Config, compute_frame_ssimulacra2_with_config};
use yuvxyb::{ColorPrimaries, Rgb, TransferCharacteristic};

fn create_test_image(width: usize, height: usize, seed: u64) -> Rgb {
    let mut state = seed;
    let data: Vec<[f32; 3]> = (0..width * height)
        .map(|_| {
            state = state.wrapping_mul(6364136223846793005).wrapping_add(1);
            let r = ((state >> 33) & 0xFF) as f32 / 255.0;
            state = state.wrapping_mul(6364136223846793005).wrapping_add(1);
            let g = ((state >> 33) & 0xFF) as f32 / 255.0;
            state = state.wrapping_mul(6364136223846793005).wrapping_add(1);
            let b = ((state >> 33) & 0xFF) as f32 / 255.0;
            [r, g, b]
        })
        .collect();

    let nz_width = std::num::NonZeroUsize::new(width).unwrap();
    let nz_height = std::num::NonZeroUsize::new(height).unwrap();
    Rgb::new(
        data,
        nz_width,
        nz_height,
        TransferCharacteristic::SRGB,
        ColorPrimaries::BT709,
    )
    .unwrap()
}

fn benchmark_blur(width: usize, height: usize, impl_type: SimdImpl, iterations: usize) -> f64 {
    // Create test plane
    let plane: Vec<f32> = (0..width * height)
        .map(|i| i as f32 / (width * height) as f32)
        .collect();

    let mut blur = Blur::with_simd_impl(width, height, impl_type);
    let img = [plane.clone(), plane.clone(), plane.clone()];

    // Warmup
    for _ in 0..5 {
        let _ = blur.blur(&img);
    }

    // Timed runs
    let start = Instant::now();
    for _ in 0..iterations {
        let _ = blur.blur(&img);
    }
    let elapsed = start.elapsed();

    elapsed.as_secs_f64() * 1000.0 / iterations as f64
}

fn benchmark_full_ssimulacra2(
    width: usize,
    height: usize,
    config: Ssimulacra2Config,
    iterations: usize,
) -> f64 {
    let source = create_test_image(width, height, 12345);
    let distorted = create_test_image(width, height, 67890);

    // Warmup
    for _ in 0..3 {
        let _ = compute_frame_ssimulacra2_with_config(source.clone(), distorted.clone(), config);
    }

    // Timed runs
    let start = Instant::now();
    for _ in 0..iterations {
        let _ = compute_frame_ssimulacra2_with_config(source.clone(), distorted.clone(), config);
    }
    let elapsed = start.elapsed();

    elapsed.as_secs_f64() * 1000.0 / iterations as f64
}

fn main() {
    println!("SSIMULACRA2 Implementation Benchmark");
    println!("=====================================\n");

    // Test different image sizes
    let sizes = [
        (512, 512, "512x512", 100),
        (1024, 1024, "1024x1024", 50),
        (1920, 1080, "1920x1080 (FHD)", 30),
        (2560, 1440, "2560x1440 (QHD)", 20),
        (3840, 2160, "3840x2160 (4K)", 10),
        (4096, 2160, "4096x2160 (DCI 4K)", 10),
    ];

    // Blur-only benchmarks
    println!("Blur-only benchmark (3 planes):");
    println!("{:20} {:>12} {:>12}", "Size", "Scalar", "SIMD");
    println!("{:-<48}", "");

    for (width, height, name, iters) in sizes.iter() {
        let scalar_ms = benchmark_blur(*width, *height, SimdImpl::Scalar, *iters);
        let simd_ms = benchmark_blur(*width, *height, SimdImpl::Simd, *iters);

        println!("{:20} {:>12.3} {:>12.3}", name, scalar_ms, simd_ms);
    }

    println!("\nFull SSIMULACRA2 benchmark:");
    println!("{:20} {:>12} {:>12}", "Size", "Scalar", "SIMD");
    println!("{:-<48}", "");

    for (width, height, name, iters) in sizes.iter() {
        let iters = iters / 2;

        let scalar_ms =
            benchmark_full_ssimulacra2(*width, *height, Ssimulacra2Config::scalar(), iters);
        let simd_ms = benchmark_full_ssimulacra2(*width, *height, Ssimulacra2Config::simd(), iters);

        println!("{:20} {:>12.3} {:>12.3}", name, scalar_ms, simd_ms);
    }

    println!("\nDone.");
}