1*437bfbebSnyanmisaka /*
2*437bfbebSnyanmisaka * Copyright 2020 Rockchip Electronics Co. LTD
3*437bfbebSnyanmisaka *
4*437bfbebSnyanmisaka * Licensed under the Apache License, Version 2.0 (the "License");
5*437bfbebSnyanmisaka * you may not use this file except in compliance with the License.
6*437bfbebSnyanmisaka * You may obtain a copy of the License at
7*437bfbebSnyanmisaka *
8*437bfbebSnyanmisaka * http://www.apache.org/licenses/LICENSE-2.0
9*437bfbebSnyanmisaka *
10*437bfbebSnyanmisaka * Unless required by applicable law or agreed to in writing, software
11*437bfbebSnyanmisaka * distributed under the License is distributed on an "AS IS" BASIS,
12*437bfbebSnyanmisaka * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13*437bfbebSnyanmisaka * See the License for the specific language governing permissions and
14*437bfbebSnyanmisaka * limitations under the License.
15*437bfbebSnyanmisaka */
16*437bfbebSnyanmisaka
17*437bfbebSnyanmisaka #define MODULE_TAG "iep2"
18*437bfbebSnyanmisaka
19*437bfbebSnyanmisaka #include <stdio.h>
20*437bfbebSnyanmisaka #include <string.h>
21*437bfbebSnyanmisaka #include <stdlib.h>
22*437bfbebSnyanmisaka
23*437bfbebSnyanmisaka #include "mpp_common.h"
24*437bfbebSnyanmisaka
25*437bfbebSnyanmisaka #include "iep2_api.h"
26*437bfbebSnyanmisaka #include "iep2_gmv.h"
27*437bfbebSnyanmisaka
iep2_sort(uint32_t bin[],int map[],int size)28*437bfbebSnyanmisaka static void iep2_sort(uint32_t bin[], int map[], int size)
29*437bfbebSnyanmisaka {
30*437bfbebSnyanmisaka int i, m, n;
31*437bfbebSnyanmisaka uint32_t *dat = malloc(size * sizeof(uint32_t));
32*437bfbebSnyanmisaka
33*437bfbebSnyanmisaka for (i = 0; i < size; ++i) {
34*437bfbebSnyanmisaka map[i] = i;
35*437bfbebSnyanmisaka dat[i] = bin[i];
36*437bfbebSnyanmisaka }
37*437bfbebSnyanmisaka
38*437bfbebSnyanmisaka for (m = 0; m < size; ++m) {
39*437bfbebSnyanmisaka int max = m;
40*437bfbebSnyanmisaka uint32_t temp;
41*437bfbebSnyanmisaka int p;
42*437bfbebSnyanmisaka
43*437bfbebSnyanmisaka for (n = m + 1; n < size; ++n) if (dat[n] > dat[max]) max = n;
44*437bfbebSnyanmisaka // Put found minimum element in its place
45*437bfbebSnyanmisaka temp = dat[m];
46*437bfbebSnyanmisaka p = map[m];
47*437bfbebSnyanmisaka
48*437bfbebSnyanmisaka map[m] = map[max];
49*437bfbebSnyanmisaka map[max] = p;
50*437bfbebSnyanmisaka dat[m] = dat[max];
51*437bfbebSnyanmisaka dat[max] = temp;
52*437bfbebSnyanmisaka }
53*437bfbebSnyanmisaka
54*437bfbebSnyanmisaka free(dat);
55*437bfbebSnyanmisaka }
56*437bfbebSnyanmisaka
iep2_is_subt_mv(int mv,struct mv_list * mv_ls)57*437bfbebSnyanmisaka static int iep2_is_subt_mv(int mv, struct mv_list *mv_ls)
58*437bfbebSnyanmisaka {
59*437bfbebSnyanmisaka int i;
60*437bfbebSnyanmisaka
61*437bfbebSnyanmisaka for (i = 0; i < mv_ls->idx; ++i) {
62*437bfbebSnyanmisaka if (RKABS(mv_ls->mv[i] - (mv * 4)) < 3)
63*437bfbebSnyanmisaka return 1;
64*437bfbebSnyanmisaka }
65*437bfbebSnyanmisaka
66*437bfbebSnyanmisaka return 0;
67*437bfbebSnyanmisaka }
68*437bfbebSnyanmisaka
iep2_update_gmv(struct iep2_api_ctx * ctx,struct mv_list * mv_ls)69*437bfbebSnyanmisaka void iep2_update_gmv(struct iep2_api_ctx *ctx, struct mv_list *mv_ls)
70*437bfbebSnyanmisaka {
71*437bfbebSnyanmisaka int rows = ctx->params.tile_rows;
72*437bfbebSnyanmisaka int cols = ctx->params.tile_cols;
73*437bfbebSnyanmisaka uint32_t *bin = ctx->output.mv_hist;
74*437bfbebSnyanmisaka int lbin = MPP_ARRAY_ELEMS(ctx->output.mv_hist);
75*437bfbebSnyanmisaka int i;
76*437bfbebSnyanmisaka
77*437bfbebSnyanmisaka int map[MPP_ARRAY_ELEMS(ctx->output.mv_hist)];
78*437bfbebSnyanmisaka
79*437bfbebSnyanmisaka uint32_t r = 6;
80*437bfbebSnyanmisaka
81*437bfbebSnyanmisaka // print mvc histogram of current motion estimation.
82*437bfbebSnyanmisaka for (i = 0; i < lbin; ++i) {
83*437bfbebSnyanmisaka if (bin[i] == 0)
84*437bfbebSnyanmisaka continue;
85*437bfbebSnyanmisaka iep_dbg_trace("mv(%d) %d\n", i - MVL, bin[i]);
86*437bfbebSnyanmisaka }
87*437bfbebSnyanmisaka
88*437bfbebSnyanmisaka bin[MVL] = 0; // disable 0 mv
89*437bfbebSnyanmisaka
90*437bfbebSnyanmisaka // update motion vector candidates
91*437bfbebSnyanmisaka iep2_sort(bin, map, lbin);
92*437bfbebSnyanmisaka
93*437bfbebSnyanmisaka iep_dbg_trace("sort map:\n");
94*437bfbebSnyanmisaka if (0) {
95*437bfbebSnyanmisaka for (i = 0; i < lbin; ++i) {
96*437bfbebSnyanmisaka fprintf(stderr, "%d ", map[i]);
97*437bfbebSnyanmisaka }
98*437bfbebSnyanmisaka fprintf(stderr, "\n");
99*437bfbebSnyanmisaka }
100*437bfbebSnyanmisaka
101*437bfbebSnyanmisaka memset(ctx->params.mv_tru_list, 0, sizeof(ctx->params.mv_tru_list));
102*437bfbebSnyanmisaka memset(ctx->params.mv_tru_vld, 0, sizeof(ctx->params.mv_tru_vld));
103*437bfbebSnyanmisaka
104*437bfbebSnyanmisaka // Get top 8 candidates of current motion estimation.
105*437bfbebSnyanmisaka for (i = 0; i < 8; ++i) {
106*437bfbebSnyanmisaka int8_t x = map[i] - MVL;
107*437bfbebSnyanmisaka
108*437bfbebSnyanmisaka if (bin[map[i]] > r * ((rows * cols) >> 7) ||
109*437bfbebSnyanmisaka iep2_is_subt_mv(x, mv_ls)) {
110*437bfbebSnyanmisaka
111*437bfbebSnyanmisaka // 1 bit at low endian for mv valid check
112*437bfbebSnyanmisaka ctx->params.mv_tru_list[i] = x;
113*437bfbebSnyanmisaka ctx->params.mv_tru_vld[i] = 1;
114*437bfbebSnyanmisaka } else {
115*437bfbebSnyanmisaka if (i == 0) {
116*437bfbebSnyanmisaka ctx->params.mv_tru_list[0] = 0;
117*437bfbebSnyanmisaka ctx->params.mv_tru_vld[0] = 1;
118*437bfbebSnyanmisaka }
119*437bfbebSnyanmisaka break;
120*437bfbebSnyanmisaka }
121*437bfbebSnyanmisaka }
122*437bfbebSnyanmisaka
123*437bfbebSnyanmisaka for (i = 0; i < 8; ++i)
124*437bfbebSnyanmisaka iep_dbg_trace("new mv candidates list[%d] (%d,%d) %d\n",
125*437bfbebSnyanmisaka i, ctx->params.mv_tru_list[i], 0, ctx->params.mv_tru_vld[i]);
126*437bfbebSnyanmisaka }
127*437bfbebSnyanmisaka
128