/src/aom/av1/encoder/model_rd.h

Source (jump to first uncovered line)
/*
 * Copyright (c) 2020, Alliance for Open Media. All rights reserved.
 *
 * This source code is subject to the terms of the BSD 2 Clause License and
 * the Alliance for Open Media Patent License 1.0. If the BSD 2 Clause License
 * was not distributed with this source code in the LICENSE file, you can
 * obtain it at www.aomedia.org/license/software. If the Alliance for Open
 * Media Patent License 1.0 was not distributed with this source code in the
 * PATENTS file, you can obtain it at www.aomedia.org/license/patent.
 */

#ifndef AOM_AV1_ENCODER_MODEL_RD_H_
#define AOM_AV1_ENCODER_MODEL_RD_H_

#include "aom/aom_integer.h"
#include "av1/encoder/block.h"
#include "av1/encoder/encoder.h"
#include "av1/encoder/pustats.h"
#include "av1/encoder/rdopt_utils.h"
#include "config/aom_dsp_rtcd.h"

#ifdef __cplusplus
extern "C" {
#endif

// 0: Legacy model
// 1: Curve fit model
// 2: Surface fit model
// 3: DNN regression model
// 4: Full rd model
#define MODELRD_TYPE_INTERP_FILTER 1
#define MODELRD_TYPE_TX_SEARCH_PRUNE 1
#define MODELRD_TYPE_MASKED_COMPOUND 1
#define MODELRD_TYPE_INTERINTRA 1
#define MODELRD_TYPE_INTRA 1
#define MODELRD_TYPE_MOTION_MODE_RD 1

typedef void (*model_rd_for_sb_type)(
    const AV1_COMP *const cpi, BLOCK_SIZE bsize, MACROBLOCK *x, MACROBLOCKD *xd,
    int plane_from, int plane_to, int *out_rate_sum, int64_t *out_dist_sum,
    uint8_t *skip_txfm_sb, int64_t *skip_sse_sb, int *plane_rate,
    int64_t *plane_sse, int64_t *plane_dist);
typedef void (*model_rd_from_sse_type)(const AV1_COMP *const cpi,
                                       const MACROBLOCK *const x,
                                       BLOCK_SIZE plane_bsize, int plane,
                                       int64_t sse, int num_samples, int *rate,
                                       int64_t *dist);

static int64_t calculate_sse(MACROBLOCKD *const xd,
                             const struct macroblock_plane *p,
                             struct macroblockd_plane *pd, const int bw,
                             const int bh) {
  int64_t sse = 0;
  const int shift = xd->bd - 8;
#if CONFIG_AV1_HIGHBITDEPTH
  if (is_cur_buf_hbd(xd)) {
    sse = aom_highbd_sse(p->src.buf, p->src.stride, pd->dst.buf, pd->dst.stride,
                         bw, bh);
  } else {
    sse =
        aom_sse(p->src.buf, p->src.stride, pd->dst.buf, pd->dst.stride, bw, bh);
  }
#else
  sse = aom_sse(p->src.buf, p->src.stride, pd->dst.buf, pd->dst.stride, bw, bh);
#endif
  sse = ROUND_POWER_OF_TWO(sse, shift * 2);
  return sse;
}

static inline int64_t compute_sse_plane(MACROBLOCK *x, MACROBLOCKD *xd,
                                        int plane, const BLOCK_SIZE bsize) {
  struct macroblockd_plane *const pd = &xd->plane[plane];
  const BLOCK_SIZE plane_bsize =
      get_plane_block_size(bsize, pd->subsampling_x, pd->subsampling_y);
  int bw, bh;
  const struct macroblock_plane *const p = &x->plane[plane];
  get_txb_dimensions(xd, plane, plane_bsize, 0, 0, plane_bsize, NULL, NULL, &bw,
                     &bh);

  int64_t sse = calculate_sse(xd, p, pd, bw, bh);

  return sse;
}

static inline void model_rd_from_sse(const AV1_COMP *const cpi,
                                     const MACROBLOCK *const x,
                                     BLOCK_SIZE plane_bsize, int plane,
                                     int64_t sse, int num_samples, int *rate,
                                     int64_t *dist) {
  (void)num_samples;
  const MACROBLOCKD *const xd = &x->e_mbd;
  const struct macroblock_plane *const p = &x->plane[plane];
  const int dequant_shift = (is_cur_buf_hbd(xd)) ? xd->bd - 5 : 3;

  // Fast approximate the modelling function.
  if (cpi->sf.rd_sf.simple_model_rd_from_var) {
    const int64_t square_error = sse;
    int quantizer = p->dequant_QTX[1] >> dequant_shift;
    if (quantizer < 120)
      *rate = (int)AOMMIN(
          (square_error * (280 - quantizer)) >> (16 - AV1_PROB_COST_SHIFT),
          INT_MAX);
    else
      *rate = 0;
    assert(*rate >= 0);
    *dist = (square_error * quantizer) >> 8;
  } else {
    av1_model_rd_from_var_lapndz(sse, num_pels_log2_lookup[plane_bsize],
                                 p->dequant_QTX[1] >> dequant_shift, rate,
                                 dist);
  }
  *dist <<= 4;
}

// Fits a curve for rate and distortion using as feature:
// log2(sse_norm/qstep^2)
static inline void model_rd_with_curvfit(const AV1_COMP *const cpi,
                                         const MACROBLOCK *const x,
                                         BLOCK_SIZE plane_bsize, int plane,
                                         int64_t sse, int num_samples,
                                         int *rate, int64_t *dist) {
  (void)cpi;
  (void)plane_bsize;
  const MACROBLOCKD *const xd = &x->e_mbd;
  const struct macroblock_plane *const p = &x->plane[plane];
  const int dequant_shift = (is_cur_buf_hbd(xd)) ? xd->bd - 5 : 3;
  const int qstep = AOMMAX(p->dequant_QTX[1] >> dequant_shift, 1);

  if (sse == 0) {
    if (rate) *rate = 0;
    if (dist) *dist = 0;
    return;
  }
  const double sse_norm = (double)sse / num_samples;
  const double qstepsqr = (double)qstep * qstep;
  const double xqr = log2(sse_norm / qstepsqr);
  double rate_f, dist_by_sse_norm_f;
  av1_model_rd_curvfit(plane_bsize, sse_norm, xqr, &rate_f,
                       &dist_by_sse_norm_f);

  const double dist_f = dist_by_sse_norm_f * sse_norm;
  int rate_i = (int)(AOMMAX(0.0, rate_f * num_samples) + 0.5);
  int64_t dist_i = (int64_t)(AOMMAX(0.0, dist_f * num_samples) + 0.5);

  // Check if skip is better
  if (rate_i == 0) {
    dist_i = sse << 4;
  } else if (RDCOST(x->rdmult, rate_i, dist_i) >=
             RDCOST(x->rdmult, 0, sse << 4)) {
    rate_i = 0;
    dist_i = sse << 4;
  }

  if (rate) *rate = rate_i;
  if (dist) *dist = dist_i;
}

static inline void model_rd_for_sb(const AV1_COMP *const cpi, BLOCK_SIZE bsize,
                                   MACROBLOCK *x, MACROBLOCKD *xd,
                                   int plane_from, int plane_to,
                                   int *out_rate_sum, int64_t *out_dist_sum,
                                   uint8_t *skip_txfm_sb, int64_t *skip_sse_sb,
                                   int *plane_rate, int64_t *plane_sse,
                                   int64_t *plane_dist) {
  // Note our transform coeffs are 8 times an orthogonal transform.
  // Hence quantizer step is also 8 times. To get effective quantizer
  // we need to divide by 8 before sending to modeling function.
  int plane;
  const int ref = xd->mi[0]->ref_frame[0];

  int64_t rate_sum = 0;
  int64_t dist_sum = 0;
  int64_t total_sse = 0;

  assert(bsize < BLOCK_SIZES_ALL);

  for (plane = plane_from; plane <= plane_to; ++plane) {
    if (plane && !xd->is_chroma_ref) break;
    struct macroblock_plane *const p = &x->plane[plane];
    struct macroblockd_plane *const pd = &xd->plane[plane];
    const BLOCK_SIZE plane_bsize =
        get_plane_block_size(bsize, pd->subsampling_x, pd->subsampling_y);
    assert(plane_bsize < BLOCK_SIZES_ALL);
    const int bw = block_size_wide[plane_bsize];
    const int bh = block_size_high[plane_bsize];
    int64_t sse;
    int rate;
    int64_t dist;

    sse = calculate_sse(xd, p, pd, bw, bh);

    model_rd_from_sse(cpi, x, plane_bsize, plane, sse, bw * bh, &rate, &dist);

    if (plane == 0) x->pred_sse[ref] = (unsigned int)AOMMIN(sse, UINT_MAX);

    total_sse += sse;
    rate_sum += rate;
    dist_sum += dist;
    if (plane_rate) plane_rate[plane] = rate;
    if (plane_sse) plane_sse[plane] = sse;
    if (plane_dist) plane_dist[plane] = dist;
    assert(rate_sum >= 0);
  }

  if (skip_txfm_sb) *skip_txfm_sb = total_sse == 0;
  if (skip_sse_sb) *skip_sse_sb = total_sse << 4;
  rate_sum = AOMMIN(rate_sum, INT_MAX);
  *out_rate_sum = (int)rate_sum;
  *out_dist_sum = dist_sum;
}

static inline void model_rd_for_sb_with_curvfit(
    const AV1_COMP *const cpi, BLOCK_SIZE bsize, MACROBLOCK *x, MACROBLOCKD *xd,
    int plane_from, int plane_to, int *out_rate_sum, int64_t *out_dist_sum,
    uint8_t *skip_txfm_sb, int64_t *skip_sse_sb, int *plane_rate,
    int64_t *plane_sse, int64_t *plane_dist) {
  // Note our transform coeffs are 8 times an orthogonal transform.
  // Hence quantizer step is also 8 times. To get effective quantizer
  // we need to divide by 8 before sending to modeling function.
  const int ref = xd->mi[0]->ref_frame[0];

  int64_t rate_sum = 0;
  int64_t dist_sum = 0;
  int64_t total_sse = 0;

  for (int plane = plane_from; plane <= plane_to; ++plane) {
    if (plane && !xd->is_chroma_ref) break;
    struct macroblockd_plane *const pd = &xd->plane[plane];
    const BLOCK_SIZE plane_bsize =
        get_plane_block_size(bsize, pd->subsampling_x, pd->subsampling_y);
    int64_t dist, sse;
    int rate;
    int bw, bh;
    const struct macroblock_plane *const p = &x->plane[plane];
    get_txb_dimensions(xd, plane, plane_bsize, 0, 0, plane_bsize, NULL, NULL,
                       &bw, &bh);

    sse = calculate_sse(xd, p, pd, bw, bh);
    model_rd_with_curvfit(cpi, x, plane_bsize, plane, sse, bw * bh, &rate,
                          &dist);

    if (plane == 0) x->pred_sse[ref] = (unsigned int)AOMMIN(sse, UINT_MAX);

    total_sse += sse;
    rate_sum += rate;
    dist_sum += dist;

    if (plane_rate) plane_rate[plane] = rate;
    if (plane_sse) plane_sse[plane] = sse;
    if (plane_dist) plane_dist[plane] = dist;
  }

  if (skip_txfm_sb) *skip_txfm_sb = rate_sum == 0;
  if (skip_sse_sb) *skip_sse_sb = total_sse << 4;
  *out_rate_sum = (int)rate_sum;
  *out_dist_sum = dist_sum;
}

enum { MODELRD_LEGACY, MODELRD_CURVFIT, MODELRD_TYPES } UENUM1BYTE(ModelRdType);

static const model_rd_for_sb_type model_rd_sb_fn[MODELRD_TYPES] = {
  model_rd_for_sb, model_rd_for_sb_with_curvfit
};

static const model_rd_from_sse_type model_rd_sse_fn[MODELRD_TYPES] = {
  model_rd_from_sse, model_rd_with_curvfit
};

#ifdef __cplusplus
}  // extern "C"
#endif
#endif  // AOM_AV1_ENCODER_MODEL_RD_H_

Coverage Report

Created: 2025-06-22 08:04

Line	Count	Source (jump to first uncovered line)
1		/*
2		* Copyright (c) 2020, Alliance for Open Media. All rights reserved.
3		*
4		* This source code is subject to the terms of the BSD 2 Clause License and
5		* the Alliance for Open Media Patent License 1.0. If the BSD 2 Clause License
6		* was not distributed with this source code in the LICENSE file, you can
7		* obtain it at www.aomedia.org/license/software. If the Alliance for Open
8		* Media Patent License 1.0 was not distributed with this source code in the
9		* PATENTS file, you can obtain it at www.aomedia.org/license/patent.
10		*/
11
12		#ifndef AOM_AV1_ENCODER_MODEL_RD_H_
13		#define AOM_AV1_ENCODER_MODEL_RD_H_
14
15		#include "aom/aom_integer.h"
16		#include "av1/encoder/block.h"
17		#include "av1/encoder/encoder.h"
18		#include "av1/encoder/pustats.h"
19		#include "av1/encoder/rdopt_utils.h"
20		#include "config/aom_dsp_rtcd.h"
21
22		#ifdef __cplusplus
23		extern "C" {
24		#endif
25
26		// 0: Legacy model
27		// 1: Curve fit model
28		// 2: Surface fit model
29		// 3: DNN regression model
30		// 4: Full rd model
31	0	#define MODELRD_TYPE_INTERP_FILTER 1
32	0	#define MODELRD_TYPE_TX_SEARCH_PRUNE 1
33	0	#define MODELRD_TYPE_MASKED_COMPOUND 1
34	0	#define MODELRD_TYPE_INTERINTRA 1
35		#define MODELRD_TYPE_INTRA 1
36	0	#define MODELRD_TYPE_MOTION_MODE_RD 1
37
38		typedef void (*model_rd_for_sb_type)(
39		const AV1_COMP const cpi, BLOCK_SIZE bsize, MACROBLOCK x, MACROBLOCKD *xd,
40		int plane_from, int plane_to, int out_rate_sum, int64_t out_dist_sum,
41		uint8_t skip_txfm_sb, int64_t skip_sse_sb, int *plane_rate,
42		int64_t plane_sse, int64_t plane_dist);
43		typedef void (model_rd_from_sse_type)(const AV1_COMP const cpi,
44		const MACROBLOCK *const x,
45		BLOCK_SIZE plane_bsize, int plane,
46		int64_t sse, int num_samples, int *rate,
47		int64_t *dist);
48
49		static int64_t calculate_sse(MACROBLOCKD *const xd,
50		const struct macroblock_plane *p,
51		struct macroblockd_plane *pd, const int bw,
52	0	const int bh) {
53	0	int64_t sse = 0;
54	0	const int shift = xd->bd - 8;
55	0	#if CONFIG_AV1_HIGHBITDEPTH
56	0	if (is_cur_buf_hbd(xd)) {
57	0	sse = aom_highbd_sse(p->src.buf, p->src.stride, pd->dst.buf, pd->dst.stride,
58	0	bw, bh);
59	0	} else {
60	0	sse =
61	0	aom_sse(p->src.buf, p->src.stride, pd->dst.buf, pd->dst.stride, bw, bh);
62	0	}
63		#else
64		sse = aom_sse(p->src.buf, p->src.stride, pd->dst.buf, pd->dst.stride, bw, bh);
65		#endif
66	0	sse = ROUND_POWER_OF_TWO(sse, shift * 2);
67	0	return sse;
68	0	} Unexecuted instantiation: av1_cx_iface.c:calculate_sse Unexecuted instantiation: allintra_vis.c:calculate_sse Unexecuted instantiation: encodeframe.c:calculate_sse Unexecuted instantiation: encoder.c:calculate_sse Unexecuted instantiation: encoder_utils.c:calculate_sse Unexecuted instantiation: ethread.c:calculate_sse Unexecuted instantiation: palette.c:calculate_sse Unexecuted instantiation: partition_search.c:calculate_sse Unexecuted instantiation: rdopt.c:calculate_sse Unexecuted instantiation: nonrd_pickmode.c:calculate_sse Unexecuted instantiation: superres_scale.c:calculate_sse Unexecuted instantiation: svc_layercontext.c:calculate_sse Unexecuted instantiation: temporal_filter.c:calculate_sse Unexecuted instantiation: tx_search.c:calculate_sse Unexecuted instantiation: intra_mode_search.c:calculate_sse Unexecuted instantiation: compound_type.c:calculate_sse Unexecuted instantiation: encode_strategy.c:calculate_sse Unexecuted instantiation: interp_search.c:calculate_sse
69
70		static inline int64_t compute_sse_plane(MACROBLOCK x, MACROBLOCKD xd,
71	0	int plane, const BLOCK_SIZE bsize) {
72	0	struct macroblockd_plane *const pd = &xd->plane[plane];
73	0	const BLOCK_SIZE plane_bsize =
74	0	get_plane_block_size(bsize, pd->subsampling_x, pd->subsampling_y);
75	0	int bw, bh;
76	0	const struct macroblock_plane *const p = &x->plane[plane];
77	0	get_txb_dimensions(xd, plane, plane_bsize, 0, 0, plane_bsize, NULL, NULL, &bw,
78	0	&bh);
79
80	0	int64_t sse = calculate_sse(xd, p, pd, bw, bh);
81
82	0	return sse;
83	0	} Unexecuted instantiation: av1_cx_iface.c:compute_sse_plane Unexecuted instantiation: allintra_vis.c:compute_sse_plane Unexecuted instantiation: encodeframe.c:compute_sse_plane Unexecuted instantiation: encoder.c:compute_sse_plane Unexecuted instantiation: encoder_utils.c:compute_sse_plane Unexecuted instantiation: ethread.c:compute_sse_plane Unexecuted instantiation: palette.c:compute_sse_plane Unexecuted instantiation: partition_search.c:compute_sse_plane Unexecuted instantiation: rdopt.c:compute_sse_plane Unexecuted instantiation: nonrd_pickmode.c:compute_sse_plane Unexecuted instantiation: superres_scale.c:compute_sse_plane Unexecuted instantiation: svc_layercontext.c:compute_sse_plane Unexecuted instantiation: temporal_filter.c:compute_sse_plane Unexecuted instantiation: tx_search.c:compute_sse_plane Unexecuted instantiation: intra_mode_search.c:compute_sse_plane Unexecuted instantiation: compound_type.c:compute_sse_plane Unexecuted instantiation: encode_strategy.c:compute_sse_plane Unexecuted instantiation: interp_search.c:compute_sse_plane
84
85		static inline void model_rd_from_sse(const AV1_COMP *const cpi,
86		const MACROBLOCK *const x,
87		BLOCK_SIZE plane_bsize, int plane,
88		int64_t sse, int num_samples, int *rate,
89	0	int64_t *dist) {
90	0	(void)num_samples;
91	0	const MACROBLOCKD *const xd = &x->e_mbd;
92	0	const struct macroblock_plane *const p = &x->plane[plane];
93	0	const int dequant_shift = (is_cur_buf_hbd(xd)) ? xd->bd - 5 : 3;
94
95		// Fast approximate the modelling function.
96	0	if (cpi->sf.rd_sf.simple_model_rd_from_var) {
97	0	const int64_t square_error = sse;
98	0	int quantizer = p->dequant_QTX[1] >> dequant_shift;
99	0	if (quantizer < 120)
100	0	*rate = (int)AOMMIN(
101	0	(square_error * (280 - quantizer)) >> (16 - AV1_PROB_COST_SHIFT),
102	0	INT_MAX);
103	0	else
104	0	*rate = 0;
105	0	assert(*rate >= 0);
106	0	dist = (square_error quantizer) >> 8;
107	0	} else {
108	0	av1_model_rd_from_var_lapndz(sse, num_pels_log2_lookup[plane_bsize],
109	0	p->dequant_QTX[1] >> dequant_shift, rate,
110	0	dist);
111	0	}
112	0	*dist <<= 4;
113	0	} Unexecuted instantiation: av1_cx_iface.c:model_rd_from_sse Unexecuted instantiation: allintra_vis.c:model_rd_from_sse Unexecuted instantiation: encodeframe.c:model_rd_from_sse Unexecuted instantiation: encoder.c:model_rd_from_sse Unexecuted instantiation: encoder_utils.c:model_rd_from_sse Unexecuted instantiation: ethread.c:model_rd_from_sse Unexecuted instantiation: palette.c:model_rd_from_sse Unexecuted instantiation: partition_search.c:model_rd_from_sse Unexecuted instantiation: rdopt.c:model_rd_from_sse Unexecuted instantiation: nonrd_pickmode.c:model_rd_from_sse Unexecuted instantiation: superres_scale.c:model_rd_from_sse Unexecuted instantiation: svc_layercontext.c:model_rd_from_sse Unexecuted instantiation: temporal_filter.c:model_rd_from_sse Unexecuted instantiation: tx_search.c:model_rd_from_sse Unexecuted instantiation: intra_mode_search.c:model_rd_from_sse Unexecuted instantiation: compound_type.c:model_rd_from_sse Unexecuted instantiation: encode_strategy.c:model_rd_from_sse Unexecuted instantiation: interp_search.c:model_rd_from_sse
114
115		// Fits a curve for rate and distortion using as feature:
116		// log2(sse_norm/qstep^2)
117		static inline void model_rd_with_curvfit(const AV1_COMP *const cpi,
118		const MACROBLOCK *const x,
119		BLOCK_SIZE plane_bsize, int plane,
120		int64_t sse, int num_samples,
121	0	int rate, int64_t dist) {
122	0	(void)cpi;
123	0	(void)plane_bsize;
124	0	const MACROBLOCKD *const xd = &x->e_mbd;
125	0	const struct macroblock_plane *const p = &x->plane[plane];
126	0	const int dequant_shift = (is_cur_buf_hbd(xd)) ? xd->bd - 5 : 3;
127	0	const int qstep = AOMMAX(p->dequant_QTX[1] >> dequant_shift, 1);
128
129	0	if (sse == 0) {
130	0	if (rate) *rate = 0;
131	0	if (dist) *dist = 0;
132	0	return;
133	0	}
134	0	const double sse_norm = (double)sse / num_samples;
135	0	const double qstepsqr = (double)qstep * qstep;
136	0	const double xqr = log2(sse_norm / qstepsqr);
137	0	double rate_f, dist_by_sse_norm_f;
138	0	av1_model_rd_curvfit(plane_bsize, sse_norm, xqr, &rate_f,
139	0	&dist_by_sse_norm_f);
140
141	0	const double dist_f = dist_by_sse_norm_f * sse_norm;
142	0	int rate_i = (int)(AOMMAX(0.0, rate_f * num_samples) + 0.5);
143	0	int64_t dist_i = (int64_t)(AOMMAX(0.0, dist_f * num_samples) + 0.5);
144
145		// Check if skip is better
146	0	if (rate_i == 0) {
147	0	dist_i = sse << 4;
148	0	} else if (RDCOST(x->rdmult, rate_i, dist_i) >=
149	0	RDCOST(x->rdmult, 0, sse << 4)) {
150	0	rate_i = 0;
151	0	dist_i = sse << 4;
152	0	}
153
154	0	if (rate) *rate = rate_i;
155	0	if (dist) *dist = dist_i;
156	0	} Unexecuted instantiation: av1_cx_iface.c:model_rd_with_curvfit Unexecuted instantiation: allintra_vis.c:model_rd_with_curvfit Unexecuted instantiation: encodeframe.c:model_rd_with_curvfit Unexecuted instantiation: encoder.c:model_rd_with_curvfit Unexecuted instantiation: encoder_utils.c:model_rd_with_curvfit Unexecuted instantiation: ethread.c:model_rd_with_curvfit Unexecuted instantiation: palette.c:model_rd_with_curvfit Unexecuted instantiation: partition_search.c:model_rd_with_curvfit Unexecuted instantiation: rdopt.c:model_rd_with_curvfit Unexecuted instantiation: nonrd_pickmode.c:model_rd_with_curvfit Unexecuted instantiation: superres_scale.c:model_rd_with_curvfit Unexecuted instantiation: svc_layercontext.c:model_rd_with_curvfit Unexecuted instantiation: temporal_filter.c:model_rd_with_curvfit Unexecuted instantiation: tx_search.c:model_rd_with_curvfit Unexecuted instantiation: intra_mode_search.c:model_rd_with_curvfit Unexecuted instantiation: compound_type.c:model_rd_with_curvfit Unexecuted instantiation: encode_strategy.c:model_rd_with_curvfit Unexecuted instantiation: interp_search.c:model_rd_with_curvfit
157
158		static inline void model_rd_for_sb(const AV1_COMP *const cpi, BLOCK_SIZE bsize,
159		MACROBLOCK x, MACROBLOCKD xd,
160		int plane_from, int plane_to,
161		int out_rate_sum, int64_t out_dist_sum,
162		uint8_t skip_txfm_sb, int64_t skip_sse_sb,
163		int plane_rate, int64_t plane_sse,
164	0	int64_t *plane_dist) {
165		// Note our transform coeffs are 8 times an orthogonal transform.
166		// Hence quantizer step is also 8 times. To get effective quantizer
167		// we need to divide by 8 before sending to modeling function.
168	0	int plane;
169	0	const int ref = xd->mi[0]->ref_frame[0];
170
171	0	int64_t rate_sum = 0;
172	0	int64_t dist_sum = 0;
173	0	int64_t total_sse = 0;
174
175	0	assert(bsize < BLOCK_SIZES_ALL);
176
177	0	for (plane = plane_from; plane <= plane_to; ++plane) {
178	0	if (plane && !xd->is_chroma_ref) break;
179	0	struct macroblock_plane *const p = &x->plane[plane];
180	0	struct macroblockd_plane *const pd = &xd->plane[plane];
181	0	const BLOCK_SIZE plane_bsize =
182	0	get_plane_block_size(bsize, pd->subsampling_x, pd->subsampling_y);
183	0	assert(plane_bsize < BLOCK_SIZES_ALL);
184	0	const int bw = block_size_wide[plane_bsize];
185	0	const int bh = block_size_high[plane_bsize];
186	0	int64_t sse;
187	0	int rate;
188	0	int64_t dist;
189
190	0	sse = calculate_sse(xd, p, pd, bw, bh);
191
192	0	model_rd_from_sse(cpi, x, plane_bsize, plane, sse, bw * bh, &rate, &dist);
193
194	0	if (plane == 0) x->pred_sse[ref] = (unsigned int)AOMMIN(sse, UINT_MAX);
195
196	0	total_sse += sse;
197	0	rate_sum += rate;
198	0	dist_sum += dist;
199	0	if (plane_rate) plane_rate[plane] = rate;
200	0	if (plane_sse) plane_sse[plane] = sse;
201	0	if (plane_dist) plane_dist[plane] = dist;
202	0	assert(rate_sum >= 0);
203	0	}
204
205	0	if (skip_txfm_sb) *skip_txfm_sb = total_sse == 0;
206	0	if (skip_sse_sb) *skip_sse_sb = total_sse << 4;
207	0	rate_sum = AOMMIN(rate_sum, INT_MAX);
208	0	*out_rate_sum = (int)rate_sum;
209	0	*out_dist_sum = dist_sum;
210	0	} Unexecuted instantiation: av1_cx_iface.c:model_rd_for_sb Unexecuted instantiation: allintra_vis.c:model_rd_for_sb Unexecuted instantiation: encodeframe.c:model_rd_for_sb Unexecuted instantiation: encoder.c:model_rd_for_sb Unexecuted instantiation: encoder_utils.c:model_rd_for_sb Unexecuted instantiation: ethread.c:model_rd_for_sb Unexecuted instantiation: palette.c:model_rd_for_sb Unexecuted instantiation: partition_search.c:model_rd_for_sb Unexecuted instantiation: rdopt.c:model_rd_for_sb Unexecuted instantiation: nonrd_pickmode.c:model_rd_for_sb Unexecuted instantiation: superres_scale.c:model_rd_for_sb Unexecuted instantiation: svc_layercontext.c:model_rd_for_sb Unexecuted instantiation: temporal_filter.c:model_rd_for_sb Unexecuted instantiation: tx_search.c:model_rd_for_sb Unexecuted instantiation: intra_mode_search.c:model_rd_for_sb Unexecuted instantiation: compound_type.c:model_rd_for_sb Unexecuted instantiation: encode_strategy.c:model_rd_for_sb Unexecuted instantiation: interp_search.c:model_rd_for_sb
211
212		static inline void model_rd_for_sb_with_curvfit(
213		const AV1_COMP const cpi, BLOCK_SIZE bsize, MACROBLOCK x, MACROBLOCKD *xd,
214		int plane_from, int plane_to, int out_rate_sum, int64_t out_dist_sum,
215		uint8_t skip_txfm_sb, int64_t skip_sse_sb, int *plane_rate,
216	0	int64_t plane_sse, int64_t plane_dist) {
217		// Note our transform coeffs are 8 times an orthogonal transform.
218		// Hence quantizer step is also 8 times. To get effective quantizer
219		// we need to divide by 8 before sending to modeling function.
220	0	const int ref = xd->mi[0]->ref_frame[0];
221
222	0	int64_t rate_sum = 0;
223	0	int64_t dist_sum = 0;
224	0	int64_t total_sse = 0;
225
226	0	for (int plane = plane_from; plane <= plane_to; ++plane) {
227	0	if (plane && !xd->is_chroma_ref) break;
228	0	struct macroblockd_plane *const pd = &xd->plane[plane];
229	0	const BLOCK_SIZE plane_bsize =
230	0	get_plane_block_size(bsize, pd->subsampling_x, pd->subsampling_y);
231	0	int64_t dist, sse;
232	0	int rate;
233	0	int bw, bh;
234	0	const struct macroblock_plane *const p = &x->plane[plane];
235	0	get_txb_dimensions(xd, plane, plane_bsize, 0, 0, plane_bsize, NULL, NULL,
236	0	&bw, &bh);
237
238	0	sse = calculate_sse(xd, p, pd, bw, bh);
239	0	model_rd_with_curvfit(cpi, x, plane_bsize, plane, sse, bw * bh, &rate,
240	0	&dist);
241
242	0	if (plane == 0) x->pred_sse[ref] = (unsigned int)AOMMIN(sse, UINT_MAX);
243
244	0	total_sse += sse;
245	0	rate_sum += rate;
246	0	dist_sum += dist;
247
248	0	if (plane_rate) plane_rate[plane] = rate;
249	0	if (plane_sse) plane_sse[plane] = sse;
250	0	if (plane_dist) plane_dist[plane] = dist;
251	0	}
252
253	0	if (skip_txfm_sb) *skip_txfm_sb = rate_sum == 0;
254	0	if (skip_sse_sb) *skip_sse_sb = total_sse << 4;
255	0	*out_rate_sum = (int)rate_sum;
256	0	*out_dist_sum = dist_sum;
257	0	} Unexecuted instantiation: av1_cx_iface.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: allintra_vis.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: encodeframe.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: encoder.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: encoder_utils.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: ethread.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: palette.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: partition_search.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: rdopt.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: nonrd_pickmode.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: superres_scale.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: svc_layercontext.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: temporal_filter.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: tx_search.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: intra_mode_search.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: compound_type.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: encode_strategy.c:model_rd_for_sb_with_curvfit Unexecuted instantiation: interp_search.c:model_rd_for_sb_with_curvfit
258
259		enum { MODELRD_LEGACY, MODELRD_CURVFIT, MODELRD_TYPES } UENUM1BYTE(ModelRdType);
260
261		static const model_rd_for_sb_type model_rd_sb_fn[MODELRD_TYPES] = {
262		model_rd_for_sb, model_rd_for_sb_with_curvfit
263		};
264
265		static const model_rd_from_sse_type model_rd_sse_fn[MODELRD_TYPES] = {
266		model_rd_from_sse, model_rd_with_curvfit
267		};
268
269		#ifdef __cplusplus
270		} // extern "C"
271		#endif
272		#endif // AOM_AV1_ENCODER_MODEL_RD_H_