dnn-layer-mathbinary-test.c 5.24 KB
Newer Older
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
/*
 * Copyright (c) 2020
 *
 * This file is part of FFmpeg.
 *
 * FFmpeg is free software; you can redistribute it and/or
 * modify it under the terms of the GNU Lesser General Public
 * License as published by the Free Software Foundation; either
 * version 2.1 of the License, or (at your option) any later version.
 *
 * FFmpeg is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * Lesser General Public License for more details.
 *
 * You should have received a copy of the GNU Lesser General Public
 * License along with FFmpeg; if not, write to the Free Software
 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
 */

#include <stdio.h>
#include <string.h>
#include <math.h>
#include "libavfilter/dnn/dnn_backend_native_layer_mathbinary.h"
25
#include "libavutil/avassert.h"
26

27
#define EPSON 0.00005
28

29 30 31 32 33 34 35 36
static float get_expected(float f1, float f2, DNNMathBinaryOperation op)
{
    switch (op)
    {
    case DMBO_SUB:
        return f1 - f2;
    case DMBO_ADD:
        return f1 + f2;
37 38
    case DMBO_MUL:
        return f1 * f2;
39 40
    case DMBO_REALDIV:
        return f1 / f2;
41 42 43 44 45 46 47
    default:
        av_assert0(!"not supported yet");
        return 0.f;
    }
}

static int test_broadcast_input0(DNNMathBinaryOperation op)
48 49 50 51 52 53 54 55 56
{
    DnnLayerMathBinaryParams params;
    DnnOperand operands[2];
    int32_t input_indexes[1];
    float input[1*1*2*3] = {
        -3, 2.5, 2, -2.1, 7.8, 100
    };
    float *output;

57
    params.bin_op = op;
58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73
    params.input0_broadcast = 1;
    params.input1_broadcast = 0;
    params.v = 7.28;

    operands[0].data = input;
    operands[0].dims[0] = 1;
    operands[0].dims[1] = 1;
    operands[0].dims[2] = 2;
    operands[0].dims[3] = 3;
    operands[1].data = NULL;

    input_indexes[0] = 0;
    dnn_execute_layer_math_binary(operands, input_indexes, 1, &params);

    output = operands[1].data;
    for (int i = 0; i < sizeof(input) / sizeof(float); i++) {
74
        float expected_output = get_expected(params.v, input[i], op);
75
        if (fabs(output[i] - expected_output) > EPSON) {
76 77
            printf("op %d, at index %d, output: %f, expected_output: %f (%s:%d)\n",
                    op, i, output[i], expected_output, __FILE__, __LINE__);
78 79 80 81 82 83 84 85 86
            av_freep(&output);
            return 1;
        }
    }

    av_freep(&output);
    return 0;
}

87
static int test_broadcast_input1(DNNMathBinaryOperation op)
88 89 90 91 92 93 94 95 96
{
    DnnLayerMathBinaryParams params;
    DnnOperand operands[2];
    int32_t input_indexes[1];
    float input[1*1*2*3] = {
        -3, 2.5, 2, -2.1, 7.8, 100
    };
    float *output;

97
    params.bin_op = op;
98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113
    params.input0_broadcast = 0;
    params.input1_broadcast = 1;
    params.v = 7.28;

    operands[0].data = input;
    operands[0].dims[0] = 1;
    operands[0].dims[1] = 1;
    operands[0].dims[2] = 2;
    operands[0].dims[3] = 3;
    operands[1].data = NULL;

    input_indexes[0] = 0;
    dnn_execute_layer_math_binary(operands, input_indexes, 1, &params);

    output = operands[1].data;
    for (int i = 0; i < sizeof(input) / sizeof(float); i++) {
114
        float expected_output = get_expected(input[i], params.v, op);
115
        if (fabs(output[i] - expected_output) > EPSON) {
116 117
            printf("op %d, at index %d, output: %f, expected_output: %f (%s:%d)\n",
                    op, i, output[i], expected_output, __FILE__, __LINE__);
118 119 120 121 122 123 124 125 126
            av_freep(&output);
            return 1;
        }
    }

    av_freep(&output);
    return 0;
}

127
static int test_no_broadcast(DNNMathBinaryOperation op)
128 129 130 131 132 133 134 135 136 137 138 139
{
    DnnLayerMathBinaryParams params;
    DnnOperand operands[3];
    int32_t input_indexes[2];
    float input0[1*1*2*3] = {
        -3, 2.5, 2, -2.1, 7.8, 100
    };
    float input1[1*1*2*3] = {
        -1, 2, 3, -21, 8, 10.0
    };
    float *output;

140
    params.bin_op = op;
141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161
    params.input0_broadcast = 0;
    params.input1_broadcast = 0;

    operands[0].data = input0;
    operands[0].dims[0] = 1;
    operands[0].dims[1] = 1;
    operands[0].dims[2] = 2;
    operands[0].dims[3] = 3;
    operands[1].data = input1;
    operands[1].dims[0] = 1;
    operands[1].dims[1] = 1;
    operands[1].dims[2] = 2;
    operands[1].dims[3] = 3;
    operands[2].data = NULL;

    input_indexes[0] = 0;
    input_indexes[1] = 1;
    dnn_execute_layer_math_binary(operands, input_indexes, 2, &params);

    output = operands[2].data;
    for (int i = 0; i < sizeof(input0) / sizeof(float); i++) {
162
        float expected_output = get_expected(input0[i], input1[i], op);
163
        if (fabs(output[i] - expected_output) > EPSON) {
164 165
            printf("op %d, at index %d, output: %f, expected_output: %f (%s:%d)\n",
                    op, i, output[i], expected_output, __FILE__, __LINE__);
166 167 168 169 170 171 172 173 174
            av_freep(&output);
            return 1;
        }
    }

    av_freep(&output);
    return 0;
}

175
static int test(DNNMathBinaryOperation op)
176
{
177
    if (test_broadcast_input0(op))
178 179
        return 1;

180
    if (test_broadcast_input1(op))
181 182
        return 1;

183
    if (test_no_broadcast(op))
184 185 186 187 188 189 190
        return 1;

    return 0;
}

int main(int argc, char **argv)
{
191 192 193 194
    if (test(DMBO_SUB))
        return 1;

    if (test(DMBO_ADD))
195 196
        return 1;

197 198 199
    if (test(DMBO_MUL))
        return 1;

200 201 202
    if (test(DMBO_REALDIV))
        return 1;

203 204
    return 0;
}