23e1766ad07490b3226e937aec0008fd6c15b8fb
[ffmpeg.git] / tests / dnn / dnn-layer-mathunary-test.c
1 /*
2  * Copyright (c) 2020
3  *
4  * This file is part of FFmpeg.
5  *
6  * FFmpeg is free software; you can redistribute it and/or
7  * modify it under the terms of the GNU Lesser General Public
8  * License as published by the Free Software Foundation; either
9  * version 2.1 of the License, or (at your option) any later version.
10  *
11  * FFmpeg is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
14  * Lesser General Public License for more details.
15  *
16  * You should have received a copy of the GNU Lesser General Public
17  * License along with FFmpeg; if not, write to the Free Software
18  * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
19  */
20
21 #include <stdio.h>
22 #include <string.h>
23 #include <math.h>
24 #include "libavfilter/dnn/dnn_backend_native_layer_mathunary.h"
25 #include "libavutil/avassert.h"
26
27 #define EPS 0.00001
28
29 static float get_expected(float f, DNNMathUnaryOperation op)
30 {
31     switch (op)
32     {
33     case DMUO_ABS:
34         return (f >= 0) ? f : -f;
35     case DMUO_SIN:
36         return sin(f);
37     case DMUO_COS:
38         return cos(f);
39     default:
40         av_assert0(!"not supported yet");
41         return 0.f;
42     }
43 }
44
45 static int test(DNNMathUnaryOperation op)
46 {
47     DnnLayerMathUnaryParams params;
48     DnnOperand operands[2];
49     int32_t input_indexes[1];
50     float input[1*1*2*3] = {
51         -3, 2.5, 2, -2.1, 7.8, 100};
52     float *output;
53
54     params.un_op = op;
55
56     operands[0].data = input;
57     operands[0].dims[0] = 1;
58     operands[0].dims[1] = 1;
59     operands[0].dims[2] = 2;
60     operands[0].dims[3] = 3;
61     operands[1].data = NULL;
62
63     input_indexes[0] = 0;
64     dnn_execute_layer_math_unary(operands, input_indexes, 1, &params);
65
66     output = operands[1].data;
67     for (int i = 0; i < sizeof(input) / sizeof(float); ++i) {
68         float expected_output = get_expected(input[i], op);
69         if(fabs(output[i] - expected_output) > EPS) {
70             printf("at index %d, output: %f, expected_output: %f\n", i, output[i], expected_output);
71             av_freep(&output);
72             return 1;
73         }
74     }
75
76     av_freep(&output);
77     return 0;
78 }
79
80 int main(int agrc, char **argv)
81 {
82     if (test(DMUO_ABS))
83         return 1;
84     if (test(DMUO_SIN))
85         return 1;
86     if (test(DMUO_COS))
87         return 1;
88     return 0;
89 }