git.sesse.net Git - movit/blob - test_util.cpp

   1 #include <assert.h>
   2 #include <math.h>
   3 #include <stdio.h>
   4 #include <algorithm>
   5 #include <epoxy/gl.h>
   6 #include <gtest/gtest.h>
   7 #include <gtest/gtest-message.h>
   8
   9 #include "flat_input.h"
  10 #include "init.h"
  11 #include "resource_pool.h"
  12 #include "test_util.h"
  13 #include "util.h"
  14
  15 using namespace std;
  16
  17 namespace movit {
  18
  19 class Input;
  20
  21 namespace {
  22
  23 // Not thread-safe, but this isn't a big problem for testing.
  24 ResourcePool *get_static_pool()
  25 {
  26         static ResourcePool *resource_pool = NULL;
  27         if (!resource_pool) {
  28                 resource_pool = new ResourcePool();
  29         }
  30         return resource_pool;
  31 }
  32
  33 // Flip upside-down to compensate for different origin.
  34 template<class T>
  35 void vertical_flip(T *data, unsigned width, unsigned height)
  36 {
  37         for (unsigned y = 0; y < height / 2; ++y) {
  38                 unsigned flip_y = height - y - 1;
  39                 for (unsigned x = 0; x < width; ++x) {
  40                         swap(data[y * width + x], data[flip_y * width + x]);
  41                 }
  42         }
  43 }
  44
  45 }  // namespace
  46
  47 EffectChainTester::EffectChainTester(const float *data, unsigned width, unsigned height,
  48                                      MovitPixelFormat pixel_format, Colorspace color_space, GammaCurve gamma_curve,
  49                                      GLenum framebuffer_format)
  50         : chain(width, height, get_static_pool()), width(width), height(height), finalized(false)
  51 {
  52         CHECK(init_movit(".", MOVIT_DEBUG_OFF));
  53
  54         if (data != NULL) {
  55                 add_input(data, pixel_format, color_space, gamma_curve);
  56         }
  57
  58         glGenTextures(1, &texnum);
  59         check_error();
  60         glBindTexture(GL_TEXTURE_2D, texnum);
  61         check_error();
  62         glTexImage2D(GL_TEXTURE_2D, 0, framebuffer_format, width, height, 0, GL_RGBA, GL_UNSIGNED_BYTE, NULL);
  63         check_error();
  64
  65         glGenFramebuffers(1, &fbo);
  66         check_error();
  67         glBindFramebuffer(GL_FRAMEBUFFER, fbo);
  68         check_error();
  69         glFramebufferTexture2D(
  70                 GL_FRAMEBUFFER,
  71                 GL_COLOR_ATTACHMENT0,
  72                 GL_TEXTURE_2D,
  73                 texnum,
  74                 0);
  75         check_error();
  76         glBindFramebuffer(GL_FRAMEBUFFER, 0);
  77         check_error();
  78 }
  79
  80 EffectChainTester::~EffectChainTester()
  81 {
  82         glDeleteFramebuffers(1, &fbo);
  83         check_error();
  84         glDeleteTextures(1, &texnum);
  85         check_error();
  86 }
  87
  88 Input *EffectChainTester::add_input(const float *data, MovitPixelFormat pixel_format, Colorspace color_space, GammaCurve gamma_curve, int input_width, int input_height)
  89 {
  90         ImageFormat format;
  91         format.color_space = color_space;
  92         format.gamma_curve = gamma_curve;
  93
  94         if (input_width == -1) {
  95                 input_width = width;
  96         }
  97         if (input_height == -1) {
  98                 input_height = height;
  99         }
 100
 101         FlatInput *input = new FlatInput(format, pixel_format, GL_FLOAT, input_width, input_height);
 102         input->set_pixel_data(data);
 103         chain.add_input(input);
 104         return input;
 105 }
 106
 107 Input *EffectChainTester::add_input(const unsigned char *data, MovitPixelFormat pixel_format, Colorspace color_space, GammaCurve gamma_curve, int input_width, int input_height)
 108 {
 109         ImageFormat format;
 110         format.color_space = color_space;
 111         format.gamma_curve = gamma_curve;
 112
 113         if (input_width == -1) {
 114                 input_width = width;
 115         }
 116         if (input_height == -1) {
 117                 input_height = height;
 118         }
 119
 120         FlatInput *input = new FlatInput(format, pixel_format, GL_UNSIGNED_BYTE, input_width, input_height);
 121         input->set_pixel_data(data);
 122         chain.add_input(input);
 123         return input;
 124 }
 125
 126 void EffectChainTester::run(float *out_data, GLenum format, Colorspace color_space, GammaCurve gamma_curve, OutputAlphaFormat alpha_format)
 127 {
 128         if (!finalized) {
 129                 finalize_chain(color_space, gamma_curve, alpha_format);
 130         }
 131
 132         chain.render_to_fbo(fbo, width, height);
 133
 134         glBindFramebuffer(GL_FRAMEBUFFER, fbo);
 135         glReadPixels(0, 0, width, height, format, GL_FLOAT, out_data);
 136
 137         if (format == GL_RGBA) {
 138                 width *= 4;
 139         }
 140
 141         vertical_flip(out_data, width, height);
 142 }
 143
 144 void EffectChainTester::run(unsigned char *out_data, GLenum format, Colorspace color_space, GammaCurve gamma_curve, OutputAlphaFormat alpha_format)
 145 {
 146         if (!finalized) {
 147                 finalize_chain(color_space, gamma_curve, alpha_format);
 148         }
 149
 150         chain.render_to_fbo(fbo, width, height);
 151
 152         glBindFramebuffer(GL_FRAMEBUFFER, fbo);
 153         glReadPixels(0, 0, width, height, format, GL_UNSIGNED_BYTE, out_data);
 154
 155         if (format == GL_RGBA) {
 156                 width *= 4;
 157         }
 158
 159         vertical_flip(out_data, width, height);
 160 }
 161
 162 void EffectChainTester::finalize_chain(Colorspace color_space, GammaCurve gamma_curve, OutputAlphaFormat alpha_format)
 163 {
 164         assert(!finalized);
 165         ImageFormat image_format;
 166         image_format.color_space = color_space;
 167         image_format.gamma_curve = gamma_curve;
 168         chain.add_output(image_format, alpha_format);
 169         chain.finalize();
 170         finalized = true;
 171 }
 172
 173 void expect_equal(const float *ref, const float *result, unsigned width, unsigned height, float largest_difference_limit, float rms_limit)
 174 {
 175         float largest_difference = -1.0f;
 176         float squared_difference = 0.0f;
 177         int largest_diff_x = -1, largest_diff_y = -1;
 178
 179         for (unsigned y = 0; y < height; ++y) {
 180                 for (unsigned x = 0; x < width; ++x) {
 181                         float diff = ref[y * width + x] - result[y * width + x];
 182                         if (fabs(diff) > largest_difference) {
 183                                 largest_difference = fabs(diff);
 184                                 largest_diff_x = x;
 185                                 largest_diff_y = y;
 186                         }
 187                         squared_difference += diff * diff;
 188                 }
 189         }
 190
 191         EXPECT_LT(largest_difference, largest_difference_limit)
 192                 << "Largest difference is in x=" << largest_diff_x << ", y=" << largest_diff_y << ":\n"
 193                 << "Reference: " << ref[largest_diff_y * width + largest_diff_x] << "\n"
 194                 << "Result:    " << result[largest_diff_y * width + largest_diff_x];
 195
 196         float rms = sqrt(squared_difference) / (width * height);
 197         EXPECT_LT(rms, rms_limit);
 198
 199         if (largest_difference >= largest_difference_limit || rms >= rms_limit) {
 200                 fprintf(stderr, "Dumping matrices for easier debugging, since at least one test failed.\n");
 201
 202                 fprintf(stderr, "Reference:\n");
 203                 for (unsigned y = 0; y < height; ++y) {
 204                         for (unsigned x = 0; x < width; ++x) {
 205                                 fprintf(stderr, "%7.4f ", ref[y * width + x]);
 206                         }
 207                         fprintf(stderr, "\n");
 208                 }
 209
 210                 fprintf(stderr, "\nResult:\n");
 211                 for (unsigned y = 0; y < height; ++y) {
 212                         for (unsigned x = 0; x < width; ++x) {
 213                                 fprintf(stderr, "%7.4f ", result[y * width + x]);
 214                         }
 215                         fprintf(stderr, "\n");
 216                 }
 217         }
 218 }
 219
 220 void expect_equal(const unsigned char *ref, const unsigned char *result, unsigned width, unsigned height, unsigned largest_difference_limit, float rms_limit)
 221 {
 222         assert(width > 0);
 223         assert(height > 0);
 224
 225         float *ref_float = new float[width * height];
 226         float *result_float = new float[width * height];
 227
 228         for (unsigned y = 0; y < height; ++y) {
 229                 for (unsigned x = 0; x < width; ++x) {
 230                         ref_float[y * width + x] = ref[y * width + x];
 231                         result_float[y * width + x] = result[y * width + x];
 232                 }
 233         }
 234
 235         expect_equal(ref_float, result_float, width, height, largest_difference_limit, rms_limit);
 236
 237         delete[] ref_float;
 238         delete[] result_float;
 239 }
 240
 241 void test_accuracy(const float *expected, const float *result, unsigned num_values, double absolute_error_limit, double relative_error_limit, double local_relative_error_limit, double rms_limit)
 242 {
 243         double squared_difference = 0.0;
 244         for (unsigned i = 0; i < num_values; ++i) {
 245                 double absolute_error = fabs(expected[i] - result[i]);
 246                 squared_difference += absolute_error * absolute_error;
 247                 EXPECT_LT(absolute_error, absolute_error_limit);
 248
 249                 if (expected[i] > 0.0) {
 250                         double relative_error = fabs(absolute_error / expected[i]);
 251
 252                         EXPECT_LT(relative_error, relative_error_limit);
 253                 }
 254                 if (i < num_values - 1) {
 255                         double delta = expected[i + 1] - expected[i];
 256                         double local_relative_error = fabs(absolute_error / delta);
 257                         EXPECT_LT(local_relative_error, local_relative_error_limit);
 258                 }
 259         }
 260         double rms = sqrt(squared_difference) / num_values;
 261         EXPECT_LT(rms, rms_limit);
 262 }
 263
 264 }  // namespace movit