michael@0: /* michael@0: * Copyright 2011 The LibYuv Project Authors. All rights reserved. michael@0: * michael@0: * Use of this source code is governed by a BSD-style license michael@0: * that can be found in the LICENSE file in the root of the source michael@0: * tree. An additional intellectual property rights grant can be found michael@0: * in the file PATENTS. All contributing project authors may michael@0: * be found in the AUTHORS file in the root of the source tree. michael@0: */ michael@0: michael@0: #include michael@0: #include michael@0: michael@0: #include "libyuv/cpu_id.h" michael@0: #include "libyuv/scale.h" michael@0: #include "../unit_test/unit_test.h" michael@0: michael@0: namespace libyuv { michael@0: michael@0: // Test scaling with C vs Opt and return maximum pixel difference. 0 = exact. michael@0: static int TestFilter(int src_width, int src_height, michael@0: int dst_width, int dst_height, michael@0: FilterMode f, int benchmark_iterations) { michael@0: int i, j; michael@0: const int b = 128; michael@0: int src_width_uv = (Abs(src_width) + 1) >> 1; michael@0: int src_height_uv = (Abs(src_height) + 1) >> 1; michael@0: michael@0: int src_y_plane_size = (Abs(src_width) + b * 2) * (Abs(src_height) + b * 2); michael@0: int src_uv_plane_size = (src_width_uv + b * 2) * (src_height_uv + b * 2); michael@0: michael@0: int src_stride_y = b * 2 + Abs(src_width); michael@0: int src_stride_uv = b * 2 + src_width_uv; michael@0: michael@0: align_buffer_page_end(src_y, src_y_plane_size) michael@0: align_buffer_page_end(src_u, src_uv_plane_size) michael@0: align_buffer_page_end(src_v, src_uv_plane_size) michael@0: srandom(time(NULL)); michael@0: MemRandomize(src_y, src_y_plane_size); michael@0: MemRandomize(src_u, src_uv_plane_size); michael@0: MemRandomize(src_v, src_uv_plane_size); michael@0: michael@0: int dst_width_uv = (dst_width + 1) >> 1; michael@0: int dst_height_uv = (dst_height + 1) >> 1; michael@0: michael@0: int dst_y_plane_size = (dst_width + b * 2) * (dst_height + b * 2); michael@0: int dst_uv_plane_size = (dst_width_uv + b * 2) * (dst_height_uv + b * 2); michael@0: michael@0: int dst_stride_y = b * 2 + dst_width; michael@0: int dst_stride_uv = b * 2 + dst_width_uv; michael@0: michael@0: align_buffer_page_end(dst_y_c, dst_y_plane_size) michael@0: align_buffer_page_end(dst_u_c, dst_uv_plane_size) michael@0: align_buffer_page_end(dst_v_c, dst_uv_plane_size) michael@0: align_buffer_page_end(dst_y_opt, dst_y_plane_size) michael@0: align_buffer_page_end(dst_u_opt, dst_uv_plane_size) michael@0: align_buffer_page_end(dst_v_opt, dst_uv_plane_size) michael@0: michael@0: michael@0: MaskCpuFlags(0); // Disable all CPU optimization. michael@0: double c_time = get_time(); michael@0: I420Scale(src_y + (src_stride_y * b) + b, src_stride_y, michael@0: src_u + (src_stride_uv * b) + b, src_stride_uv, michael@0: src_v + (src_stride_uv * b) + b, src_stride_uv, michael@0: src_width, src_height, michael@0: dst_y_c + (dst_stride_y * b) + b, dst_stride_y, michael@0: dst_u_c + (dst_stride_uv * b) + b, dst_stride_uv, michael@0: dst_v_c + (dst_stride_uv * b) + b, dst_stride_uv, michael@0: dst_width, dst_height, f); michael@0: c_time = (get_time() - c_time); michael@0: michael@0: MaskCpuFlags(-1); // Enable all CPU optimization. michael@0: double opt_time = get_time(); michael@0: for (i = 0; i < benchmark_iterations; ++i) { michael@0: I420Scale(src_y + (src_stride_y * b) + b, src_stride_y, michael@0: src_u + (src_stride_uv * b) + b, src_stride_uv, michael@0: src_v + (src_stride_uv * b) + b, src_stride_uv, michael@0: src_width, src_height, michael@0: dst_y_opt + (dst_stride_y * b) + b, dst_stride_y, michael@0: dst_u_opt + (dst_stride_uv * b) + b, dst_stride_uv, michael@0: dst_v_opt + (dst_stride_uv * b) + b, dst_stride_uv, michael@0: dst_width, dst_height, f); michael@0: } michael@0: opt_time = (get_time() - opt_time) / benchmark_iterations; michael@0: // Report performance of C vs OPT michael@0: printf("filter %d - %8d us C - %8d us OPT\n", michael@0: f, michael@0: static_cast(c_time * 1e6), michael@0: static_cast(opt_time * 1e6)); michael@0: michael@0: // C version may be a little off from the optimized. Order of michael@0: // operations may introduce rounding somewhere. So do a difference michael@0: // of the buffers and look to see that the max difference isn't michael@0: // over 2. michael@0: int max_diff = 0; michael@0: for (i = b; i < (dst_height + b); ++i) { michael@0: for (j = b; j < (dst_width + b); ++j) { michael@0: int abs_diff = Abs(dst_y_c[(i * dst_stride_y) + j] - michael@0: dst_y_opt[(i * dst_stride_y) + j]); michael@0: if (abs_diff > max_diff) { michael@0: max_diff = abs_diff; michael@0: } michael@0: } michael@0: } michael@0: michael@0: for (i = b; i < (dst_height_uv + b); ++i) { michael@0: for (j = b; j < (dst_width_uv + b); ++j) { michael@0: int abs_diff = Abs(dst_u_c[(i * dst_stride_uv) + j] - michael@0: dst_u_opt[(i * dst_stride_uv) + j]); michael@0: if (abs_diff > max_diff) { michael@0: max_diff = abs_diff; michael@0: } michael@0: abs_diff = Abs(dst_v_c[(i * dst_stride_uv) + j] - michael@0: dst_v_opt[(i * dst_stride_uv) + j]); michael@0: if (abs_diff > max_diff) { michael@0: max_diff = abs_diff; michael@0: } michael@0: } michael@0: } michael@0: michael@0: free_aligned_buffer_page_end(dst_y_c) michael@0: free_aligned_buffer_page_end(dst_u_c) michael@0: free_aligned_buffer_page_end(dst_v_c) michael@0: free_aligned_buffer_page_end(dst_y_opt) michael@0: free_aligned_buffer_page_end(dst_u_opt) michael@0: free_aligned_buffer_page_end(dst_v_opt) michael@0: michael@0: free_aligned_buffer_page_end(src_y) michael@0: free_aligned_buffer_page_end(src_u) michael@0: free_aligned_buffer_page_end(src_v) michael@0: michael@0: return max_diff; michael@0: } michael@0: michael@0: #define TEST_FACTOR1(name, filter, hfactor, vfactor, max_diff) \ michael@0: TEST_F(libyuvTest, ScaleDownBy##name##_##filter) { \ michael@0: int diff = TestFilter(benchmark_width_, benchmark_height_, \ michael@0: Abs(benchmark_width_) * hfactor, \ michael@0: Abs(benchmark_height_) * vfactor, \ michael@0: kFilter##filter, benchmark_iterations_); \ michael@0: EXPECT_LE(diff, max_diff); \ michael@0: } michael@0: michael@0: // Test a scale factor with all 4 filters. Expect unfiltered to be exact, but michael@0: // filtering is different fixed point implementations for SSSE3, Neon and C. michael@0: #define TEST_FACTOR(name, hfactor, vfactor) \ michael@0: TEST_FACTOR1(name, None, hfactor, vfactor, 0) \ michael@0: TEST_FACTOR1(name, Linear, hfactor, vfactor, 3) \ michael@0: TEST_FACTOR1(name, Bilinear, hfactor, vfactor, 3) \ michael@0: TEST_FACTOR1(name, Box, hfactor, vfactor, 3) \ michael@0: michael@0: TEST_FACTOR(2, 1 / 2, 1 / 2) michael@0: TEST_FACTOR(4, 1 / 4, 1 / 4) michael@0: TEST_FACTOR(8, 1 / 8, 1 / 8) michael@0: TEST_FACTOR(3by4, 3 / 4, 3 / 4) michael@0: #undef TEST_FACTOR1 michael@0: #undef TEST_FACTOR michael@0: michael@0: #define TEST_SCALETO1(name, width, height, filter, max_diff) \ michael@0: TEST_F(libyuvTest, name##To##width##x##height##_##filter) { \ michael@0: int diff = TestFilter(benchmark_width_, benchmark_height_, \ michael@0: width, height, \ michael@0: kFilter##filter, benchmark_iterations_); \ michael@0: EXPECT_LE(diff, max_diff); \ michael@0: } \ michael@0: TEST_F(libyuvTest, name##From##width##x##height##_##filter) { \ michael@0: int diff = TestFilter(width, height, \ michael@0: Abs(benchmark_width_), Abs(benchmark_height_), \ michael@0: kFilter##filter, benchmark_iterations_); \ michael@0: EXPECT_LE(diff, max_diff); \ michael@0: } michael@0: michael@0: // Test scale to a specified size with all 4 filters. michael@0: #define TEST_SCALETO(name, width, height) \ michael@0: TEST_SCALETO1(name, width, height, None, 0) \ michael@0: TEST_SCALETO1(name, width, height, Linear, 3) \ michael@0: TEST_SCALETO1(name, width, height, Bilinear, 3) \ michael@0: TEST_SCALETO1(name, width, height, Box, 3) michael@0: michael@0: TEST_SCALETO(Scale, 1, 1) michael@0: TEST_SCALETO(Scale, 320, 240) michael@0: TEST_SCALETO(Scale, 352, 288) michael@0: TEST_SCALETO(Scale, 640, 360) michael@0: TEST_SCALETO(Scale, 1280, 720) michael@0: #undef TEST_SCALETO1 michael@0: #undef TEST_SCALETO michael@0: michael@0: } // namespace libyuv