hevc: Add names for reserved NAL unit types
[libav.git] / tests / rotozoom.c
CommitLineData
66961636 1/*
58c42af7 2 * Generate a synthetic YUV video sequence suitable for codec testing.
01d193aa
DB
3 *
4 * copyright (c) Sebastien Bechet <s.bechet@av7.net>
5 *
2912e87a 6 * This file is part of Libav.
01d193aa 7 *
2912e87a 8 * Libav is free software; you can redistribute it and/or
01d193aa
DB
9 * modify it under the terms of the GNU Lesser General Public
10 * License as published by the Free Software Foundation; either
11 * version 2.1 of the License, or (at your option) any later version.
12 *
2912e87a 13 * Libav is distributed in the hope that it will be useful,
01d193aa
DB
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 * Lesser General Public License for more details.
17 *
18 * You should have received a copy of the GNU Lesser General Public
2912e87a 19 * License along with Libav; if not, write to the Free Software
01d193aa 20 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
66961636 21 */
01d193aa 22
66961636
MN
23#include <stdlib.h>
24#include <stdio.h>
6afd3b92
MN
25#include <inttypes.h>
26
f0ccd53a
DB
27#include "utils.c"
28
f628559d 29#define FIXP (1 << 16)
b481bbc3 30#define MY_PI 205887 // (M_PI * FIX)
6afd3b92 31
f628559d
DB
32static int64_t int_pow(int64_t a, int p)
33{
34 int64_t v = FIXP;
6afd3b92 35
f628559d
DB
36 for (; p; p--) {
37 v *= a;
38 v /= FIXP;
6afd3b92
MN
39 }
40
41 return v;
42}
43
f628559d
DB
44static int64_t int_sin(int64_t a)
45{
46 if (a < 0)
47 a = MY_PI - a; // 0..inf
48 a %= 2 * MY_PI; // 0..2PI
6afd3b92 49
f628559d
DB
50 if (a >= MY_PI * 3 / 2)
51 a -= 2 * MY_PI; // -PI / 2 .. 3PI / 2
b481bbc3 52 if (a >= MY_PI / 2)
f628559d 53 a = MY_PI - a; // -PI / 2 .. PI / 2
115329f1 54
f628559d 55 return a - int_pow(a, 3) / 6 + int_pow(a, 5) / 120 - int_pow(a, 7) / 5040;
6afd3b92 56}
66961636 57
814208a7
DB
58static unsigned char tab_r[256 * 256];
59static unsigned char tab_g[256 * 256];
60static unsigned char tab_b[256 * 256];
66961636 61
814208a7
DB
62static int h_cos[360];
63static int h_sin[360];
66961636 64
f628559d
DB
65static int ipol(uint8_t *src, int x, int y)
66{
67 int int_x = x >> 16;
68 int int_y = y >> 16;
69 int frac_x = x & 0xFFFF;
70 int frac_y = y & 0xFFFF;
71 int s00 = src[( int_x & 255) + 256 * ( int_y & 255)];
72 int s01 = src[((int_x + 1) & 255) + 256 * ( int_y & 255)];
73 int s10 = src[( int_x & 255) + 256 * ((int_y + 1) & 255)];
74 int s11 = src[((int_x + 1) & 255) + 256 * ((int_y + 1) & 255)];
75 int s0 = (((1 << 16) - frac_x) * s00 + frac_x * s01) >> 8;
76 int s1 = (((1 << 16) - frac_x) * s10 + frac_x * s11) >> 8;
77
78 return (((1 << 16) - frac_y) * s0 + frac_y * s1) >> 24;
2add6b5d
MN
79}
80
6f2162b4 81static void gen_image(int num, int w, int h)
66961636 82{
b481bbc3
GN
83 const int c = h_cos[num % 360];
84 const int s = h_sin[num % 360];
115329f1 85
f628559d
DB
86 const int xi = -(w / 2) * c;
87 const int yi = (w / 2) * s;
115329f1 88
f628559d
DB
89 const int xj = -(h / 2) * s;
90 const int yj = -(h / 2) * c;
91 int i, j;
66961636 92
f628559d
DB
93 int x, y;
94 int xprime = xj;
95 int yprime = yj;
66961636 96
f628559d 97 for (j = 0; j < h; j++) {
b481bbc3 98 x = xprime + xi + FIXP * w / 2;
f628559d 99 xprime += s;
66961636 100
b481bbc3 101 y = yprime + yi + FIXP * h / 2;
f628559d 102 yprime += c;
66961636 103
b481bbc3 104 for (i = 0; i < w; i++) {
f628559d
DB
105 x += c;
106 y -= s;
b481bbc3
GN
107 put_pixel(i, j,
108 ipol(tab_r, x, y),
109 ipol(tab_g, x, y),
110 ipol(tab_b, x, y));
f628559d 111 }
66961636 112 }
66961636
MN
113}
114
6afd3b92
MN
115#define W 256
116#define H 256
117
2131e859 118static int init_demo(const char *filename)
f628559d
DB
119{
120 int i, j;
121 int h;
122 int radian;
123 char line[3 * W];
124
e1e0ca70 125 FILE *input_file;
f628559d 126
e1e0ca70
DB
127 input_file = fopen(filename, "rb");
128 if (!input_file) {
f628559d 129 perror(filename);
2131e859 130 return 1;
f628559d
DB
131 }
132
e1e0ca70 133 if (fread(line, 1, 15, input_file) != 15)
cbb0930f 134 return 1;
f628559d 135 for (i = 0; i < H; i++) {
e1e0ca70 136 if (fread(line, 1, 3 * W, input_file) != 3 * W)
cbb0930f 137 return 1;
f628559d
DB
138 for (j = 0; j < W; j++) {
139 tab_r[W * i + j] = line[3 * j ];
140 tab_g[W * i + j] = line[3 * j + 1];
141 tab_b[W * i + j] = line[3 * j + 2];
142 }
143 }
e1e0ca70 144 fclose(input_file);
f628559d
DB
145
146 /* tables sin/cos */
147 for (i = 0; i < 360; i++) {
b481bbc3
GN
148 radian = 2 * i * MY_PI / 360;
149 h = 2 * FIXP + int_sin(radian);
e9c10459
DB
150 h_cos[i] = h * int_sin(radian + MY_PI / 2) / 2 / FIXP;
151 h_sin[i] = h * int_sin(radian) / 2 / FIXP;
66961636 152 }
2131e859 153
b481bbc3 154 return 0;
66961636
MN
155}
156
157int main(int argc, char **argv)
158{
159 int w, h, i;
160 char buf[1024];
d0e9415d 161 int isdir = 0;
66961636 162
d0e9415d
MR
163 if (argc != 3) {
164 printf("usage: %s image.pnm file|dir\n"
66961636 165 "generate a test video stream\n", argv[0]);
771339ca 166 return 1;
66961636
MN
167 }
168
d0e9415d
MR
169 if (!freopen(argv[2], "wb", stdout))
170 isdir = 1;
7e5880e0 171
66961636
MN
172 w = DEFAULT_WIDTH;
173 h = DEFAULT_HEIGHT;
174
175 rgb_tab = malloc(w * h * 3);
f628559d
DB
176 wrap = w * 3;
177 width = w;
178 height = h;
66961636 179
7e5880e0 180 if (init_demo(argv[1]))
2131e859 181 return 1;
66961636 182
f628559d 183 for (i = 0; i < DEFAULT_NB_PICT; i++) {
66961636 184 gen_image(i, w, h);
d0e9415d 185 if (isdir) {
7e5880e0
MR
186 snprintf(buf, sizeof(buf), "%s%02d.pgm", argv[2], i);
187 pgmyuv_save(buf, w, h, rgb_tab);
188 } else {
189 pgmyuv_save(NULL, w, h, rgb_tab);
190 }
66961636 191 }
115329f1 192
66961636
MN
193 free(rgb_tab);
194 return 0;
195}