michael@0: /* michael@0: * Copyright (c) 2010 The WebM project authors. All Rights Reserved. michael@0: * michael@0: * Use of this source code is governed by a BSD-style license michael@0: * that can be found in the LICENSE file in the root of the source michael@0: * tree. An additional intellectual property rights grant can be found michael@0: * in the file PATENTS. All contributing project authors may michael@0: * be found in the AUTHORS file in the root of the source tree. michael@0: */ michael@0: michael@0: michael@0: #include michael@0: #include michael@0: #include michael@0: #include michael@0: #include michael@0: michael@0: #include "math.h" michael@0: #include "vp8/common/common.h" michael@0: #include "ratectrl.h" michael@0: #include "vp8/common/entropymode.h" michael@0: #include "vpx_mem/vpx_mem.h" michael@0: #include "vp8/common/systemdependent.h" michael@0: #include "encodemv.h" michael@0: michael@0: michael@0: #define MIN_BPB_FACTOR 0.01 michael@0: #define MAX_BPB_FACTOR 50 michael@0: michael@0: extern const MB_PREDICTION_MODE vp8_mode_order[MAX_MODES]; michael@0: michael@0: michael@0: michael@0: #ifdef MODE_STATS michael@0: extern int y_modes[5]; michael@0: extern int uv_modes[4]; michael@0: extern int b_modes[10]; michael@0: michael@0: extern int inter_y_modes[10]; michael@0: extern int inter_uv_modes[4]; michael@0: extern int inter_b_modes[10]; michael@0: #endif michael@0: michael@0: /* Bits Per MB at different Q (Multiplied by 512) */ michael@0: #define BPER_MB_NORMBITS 9 michael@0: michael@0: /* Work in progress recalibration of baseline rate tables based on michael@0: * the assumption that bits per mb is inversely proportional to the michael@0: * quantizer value. michael@0: */ michael@0: const int vp8_bits_per_mb[2][QINDEX_RANGE] = michael@0: { michael@0: /* Intra case 450000/Qintra */ michael@0: { michael@0: 1125000,900000, 750000, 642857, 562500, 500000, 450000, 450000, michael@0: 409090, 375000, 346153, 321428, 300000, 281250, 264705, 264705, michael@0: 250000, 236842, 225000, 225000, 214285, 214285, 204545, 204545, michael@0: 195652, 195652, 187500, 180000, 180000, 173076, 166666, 160714, michael@0: 155172, 150000, 145161, 140625, 136363, 132352, 128571, 125000, michael@0: 121621, 121621, 118421, 115384, 112500, 109756, 107142, 104651, michael@0: 102272, 100000, 97826, 97826, 95744, 93750, 91836, 90000, michael@0: 88235, 86538, 84905, 83333, 81818, 80357, 78947, 77586, michael@0: 76271, 75000, 73770, 72580, 71428, 70312, 69230, 68181, michael@0: 67164, 66176, 65217, 64285, 63380, 62500, 61643, 60810, michael@0: 60000, 59210, 59210, 58441, 57692, 56962, 56250, 55555, michael@0: 54878, 54216, 53571, 52941, 52325, 51724, 51136, 50561, michael@0: 49450, 48387, 47368, 46875, 45918, 45000, 44554, 44117, michael@0: 43269, 42452, 41666, 40909, 40178, 39473, 38793, 38135, michael@0: 36885, 36290, 35714, 35156, 34615, 34090, 33582, 33088, michael@0: 32608, 32142, 31468, 31034, 30405, 29801, 29220, 28662, michael@0: }, michael@0: /* Inter case 285000/Qinter */ michael@0: { michael@0: 712500, 570000, 475000, 407142, 356250, 316666, 285000, 259090, michael@0: 237500, 219230, 203571, 190000, 178125, 167647, 158333, 150000, michael@0: 142500, 135714, 129545, 123913, 118750, 114000, 109615, 105555, michael@0: 101785, 98275, 95000, 91935, 89062, 86363, 83823, 81428, michael@0: 79166, 77027, 75000, 73076, 71250, 69512, 67857, 66279, michael@0: 64772, 63333, 61956, 60638, 59375, 58163, 57000, 55882, michael@0: 54807, 53773, 52777, 51818, 50892, 50000, 49137, 47500, michael@0: 45967, 44531, 43181, 41911, 40714, 39583, 38513, 37500, michael@0: 36538, 35625, 34756, 33928, 33139, 32386, 31666, 30978, michael@0: 30319, 29687, 29081, 28500, 27941, 27403, 26886, 26388, michael@0: 25909, 25446, 25000, 24568, 23949, 23360, 22800, 22265, michael@0: 21755, 21268, 20802, 20357, 19930, 19520, 19127, 18750, michael@0: 18387, 18037, 17701, 17378, 17065, 16764, 16473, 16101, michael@0: 15745, 15405, 15079, 14766, 14467, 14179, 13902, 13636, michael@0: 13380, 13133, 12895, 12666, 12445, 12179, 11924, 11632, michael@0: 11445, 11220, 11003, 10795, 10594, 10401, 10215, 10035, michael@0: } michael@0: }; michael@0: michael@0: static const int kf_boost_qadjustment[QINDEX_RANGE] = michael@0: { michael@0: 128, 129, 130, 131, 132, 133, 134, 135, michael@0: 136, 137, 138, 139, 140, 141, 142, 143, michael@0: 144, 145, 146, 147, 148, 149, 150, 151, michael@0: 152, 153, 154, 155, 156, 157, 158, 159, michael@0: 160, 161, 162, 163, 164, 165, 166, 167, michael@0: 168, 169, 170, 171, 172, 173, 174, 175, michael@0: 176, 177, 178, 179, 180, 181, 182, 183, michael@0: 184, 185, 186, 187, 188, 189, 190, 191, michael@0: 192, 193, 194, 195, 196, 197, 198, 199, michael@0: 200, 200, 201, 201, 202, 203, 203, 203, michael@0: 204, 204, 205, 205, 206, 206, 207, 207, michael@0: 208, 208, 209, 209, 210, 210, 211, 211, michael@0: 212, 212, 213, 213, 214, 214, 215, 215, michael@0: 216, 216, 217, 217, 218, 218, 219, 219, michael@0: 220, 220, 220, 220, 220, 220, 220, 220, michael@0: 220, 220, 220, 220, 220, 220, 220, 220, michael@0: }; michael@0: michael@0: /* #define GFQ_ADJUSTMENT (Q+100) */ michael@0: #define GFQ_ADJUSTMENT vp8_gf_boost_qadjustment[Q] michael@0: const int vp8_gf_boost_qadjustment[QINDEX_RANGE] = michael@0: { michael@0: 80, 82, 84, 86, 88, 90, 92, 94, michael@0: 96, 97, 98, 99, 100, 101, 102, 103, michael@0: 104, 105, 106, 107, 108, 109, 110, 111, michael@0: 112, 113, 114, 115, 116, 117, 118, 119, michael@0: 120, 121, 122, 123, 124, 125, 126, 127, michael@0: 128, 129, 130, 131, 132, 133, 134, 135, michael@0: 136, 137, 138, 139, 140, 141, 142, 143, michael@0: 144, 145, 146, 147, 148, 149, 150, 151, michael@0: 152, 153, 154, 155, 156, 157, 158, 159, michael@0: 160, 161, 162, 163, 164, 165, 166, 167, michael@0: 168, 169, 170, 171, 172, 173, 174, 175, michael@0: 176, 177, 178, 179, 180, 181, 182, 183, michael@0: 184, 184, 185, 185, 186, 186, 187, 187, michael@0: 188, 188, 189, 189, 190, 190, 191, 191, michael@0: 192, 192, 193, 193, 194, 194, 194, 194, michael@0: 195, 195, 196, 196, 197, 197, 198, 198 michael@0: }; michael@0: michael@0: /* michael@0: const int vp8_gf_boost_qadjustment[QINDEX_RANGE] = michael@0: { michael@0: 100,101,102,103,104,105,105,106, michael@0: 106,107,107,108,109,109,110,111, michael@0: 112,113,114,115,116,117,118,119, michael@0: 120,121,122,123,124,125,126,127, michael@0: 128,129,130,131,132,133,134,135, michael@0: 136,137,138,139,140,141,142,143, michael@0: 144,145,146,147,148,149,150,151, michael@0: 152,153,154,155,156,157,158,159, michael@0: 160,161,162,163,164,165,166,167, michael@0: 168,169,170,170,171,171,172,172, michael@0: 173,173,173,174,174,174,175,175, michael@0: 175,176,176,176,177,177,177,177, michael@0: 178,178,179,179,180,180,181,181, michael@0: 182,182,183,183,184,184,185,185, michael@0: 186,186,187,187,188,188,189,189, michael@0: 190,190,191,191,192,192,193,193, michael@0: }; michael@0: */ michael@0: michael@0: static const int kf_gf_boost_qlimits[QINDEX_RANGE] = michael@0: { michael@0: 150, 155, 160, 165, 170, 175, 180, 185, michael@0: 190, 195, 200, 205, 210, 215, 220, 225, michael@0: 230, 235, 240, 245, 250, 255, 260, 265, michael@0: 270, 275, 280, 285, 290, 295, 300, 305, michael@0: 310, 320, 330, 340, 350, 360, 370, 380, michael@0: 390, 400, 410, 420, 430, 440, 450, 460, michael@0: 470, 480, 490, 500, 510, 520, 530, 540, michael@0: 550, 560, 570, 580, 590, 600, 600, 600, michael@0: 600, 600, 600, 600, 600, 600, 600, 600, michael@0: 600, 600, 600, 600, 600, 600, 600, 600, michael@0: 600, 600, 600, 600, 600, 600, 600, 600, michael@0: 600, 600, 600, 600, 600, 600, 600, 600, michael@0: 600, 600, 600, 600, 600, 600, 600, 600, michael@0: 600, 600, 600, 600, 600, 600, 600, 600, michael@0: 600, 600, 600, 600, 600, 600, 600, 600, michael@0: 600, 600, 600, 600, 600, 600, 600, 600, michael@0: }; michael@0: michael@0: /* % adjustment to target kf size based on seperation from previous frame */ michael@0: static const int kf_boost_seperation_adjustment[16] = michael@0: { michael@0: 30, 40, 50, 55, 60, 65, 70, 75, michael@0: 80, 85, 90, 95, 100, 100, 100, 100, michael@0: }; michael@0: michael@0: michael@0: static const int gf_adjust_table[101] = michael@0: { michael@0: 100, michael@0: 115, 130, 145, 160, 175, 190, 200, 210, 220, 230, michael@0: 240, 260, 270, 280, 290, 300, 310, 320, 330, 340, michael@0: 350, 360, 370, 380, 390, 400, 400, 400, 400, 400, michael@0: 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, michael@0: 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, michael@0: 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, michael@0: 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, michael@0: 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, michael@0: 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, michael@0: 400, 400, 400, 400, 400, 400, 400, 400, 400, 400, michael@0: }; michael@0: michael@0: static const int gf_intra_usage_adjustment[20] = michael@0: { michael@0: 125, 120, 115, 110, 105, 100, 95, 85, 80, 75, michael@0: 70, 65, 60, 55, 50, 50, 50, 50, 50, 50, michael@0: }; michael@0: michael@0: static const int gf_interval_table[101] = michael@0: { michael@0: 7, michael@0: 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, michael@0: 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, michael@0: 7, 7, 7, 7, 7, 7, 7, 7, 7, 7, michael@0: 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, michael@0: 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, michael@0: 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, michael@0: 9, 9, 9, 9, 9, 9, 9, 9, 9, 9, michael@0: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, michael@0: 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, michael@0: 11, 11, 11, 11, 11, 11, 11, 11, 11, 11, michael@0: }; michael@0: michael@0: static const unsigned int prior_key_frame_weight[KEY_FRAME_CONTEXT] = { 1, 2, 3, 4, 5 }; michael@0: michael@0: michael@0: void vp8_save_coding_context(VP8_COMP *cpi) michael@0: { michael@0: CODING_CONTEXT *const cc = & cpi->coding_context; michael@0: michael@0: /* Stores a snapshot of key state variables which can subsequently be michael@0: * restored with a call to vp8_restore_coding_context. These functions are michael@0: * intended for use in a re-code loop in vp8_compress_frame where the michael@0: * quantizer value is adjusted between loop iterations. michael@0: */ michael@0: michael@0: cc->frames_since_key = cpi->frames_since_key; michael@0: cc->filter_level = cpi->common.filter_level; michael@0: cc->frames_till_gf_update_due = cpi->frames_till_gf_update_due; michael@0: cc->frames_since_golden = cpi->frames_since_golden; michael@0: michael@0: vp8_copy(cc->mvc, cpi->common.fc.mvc); michael@0: vp8_copy(cc->mvcosts, cpi->rd_costs.mvcosts); michael@0: michael@0: vp8_copy(cc->ymode_prob, cpi->common.fc.ymode_prob); michael@0: vp8_copy(cc->uv_mode_prob, cpi->common.fc.uv_mode_prob); michael@0: michael@0: vp8_copy(cc->ymode_count, cpi->mb.ymode_count); michael@0: vp8_copy(cc->uv_mode_count, cpi->mb.uv_mode_count); michael@0: michael@0: michael@0: /* Stats */ michael@0: #ifdef MODE_STATS michael@0: vp8_copy(cc->y_modes, y_modes); michael@0: vp8_copy(cc->uv_modes, uv_modes); michael@0: vp8_copy(cc->b_modes, b_modes); michael@0: vp8_copy(cc->inter_y_modes, inter_y_modes); michael@0: vp8_copy(cc->inter_uv_modes, inter_uv_modes); michael@0: vp8_copy(cc->inter_b_modes, inter_b_modes); michael@0: #endif michael@0: michael@0: cc->this_frame_percent_intra = cpi->this_frame_percent_intra; michael@0: } michael@0: michael@0: michael@0: void vp8_restore_coding_context(VP8_COMP *cpi) michael@0: { michael@0: CODING_CONTEXT *const cc = & cpi->coding_context; michael@0: michael@0: /* Restore key state variables to the snapshot state stored in the michael@0: * previous call to vp8_save_coding_context. michael@0: */ michael@0: michael@0: cpi->frames_since_key = cc->frames_since_key; michael@0: cpi->common.filter_level = cc->filter_level; michael@0: cpi->frames_till_gf_update_due = cc->frames_till_gf_update_due; michael@0: cpi->frames_since_golden = cc->frames_since_golden; michael@0: michael@0: vp8_copy(cpi->common.fc.mvc, cc->mvc); michael@0: michael@0: vp8_copy(cpi->rd_costs.mvcosts, cc->mvcosts); michael@0: michael@0: vp8_copy(cpi->common.fc.ymode_prob, cc->ymode_prob); michael@0: vp8_copy(cpi->common.fc.uv_mode_prob, cc->uv_mode_prob); michael@0: michael@0: vp8_copy(cpi->mb.ymode_count, cc->ymode_count); michael@0: vp8_copy(cpi->mb.uv_mode_count, cc->uv_mode_count); michael@0: michael@0: /* Stats */ michael@0: #ifdef MODE_STATS michael@0: vp8_copy(y_modes, cc->y_modes); michael@0: vp8_copy(uv_modes, cc->uv_modes); michael@0: vp8_copy(b_modes, cc->b_modes); michael@0: vp8_copy(inter_y_modes, cc->inter_y_modes); michael@0: vp8_copy(inter_uv_modes, cc->inter_uv_modes); michael@0: vp8_copy(inter_b_modes, cc->inter_b_modes); michael@0: #endif michael@0: michael@0: michael@0: cpi->this_frame_percent_intra = cc->this_frame_percent_intra; michael@0: } michael@0: michael@0: michael@0: void vp8_setup_key_frame(VP8_COMP *cpi) michael@0: { michael@0: /* Setup for Key frame: */ michael@0: michael@0: vp8_default_coef_probs(& cpi->common); michael@0: michael@0: vpx_memcpy(cpi->common.fc.mvc, vp8_default_mv_context, sizeof(vp8_default_mv_context)); michael@0: { michael@0: int flag[2] = {1, 1}; michael@0: vp8_build_component_cost_table(cpi->mb.mvcost, (const MV_CONTEXT *) cpi->common.fc.mvc, flag); michael@0: } michael@0: michael@0: /* Make sure we initialize separate contexts for altref,gold, and normal. michael@0: * TODO shouldn't need 3 different copies of structure to do this! michael@0: */ michael@0: vpx_memcpy(&cpi->lfc_a, &cpi->common.fc, sizeof(cpi->common.fc)); michael@0: vpx_memcpy(&cpi->lfc_g, &cpi->common.fc, sizeof(cpi->common.fc)); michael@0: vpx_memcpy(&cpi->lfc_n, &cpi->common.fc, sizeof(cpi->common.fc)); michael@0: michael@0: cpi->common.filter_level = cpi->common.base_qindex * 3 / 8 ; michael@0: michael@0: /* Provisional interval before next GF */ michael@0: if (cpi->auto_gold) michael@0: cpi->frames_till_gf_update_due = cpi->baseline_gf_interval; michael@0: else michael@0: cpi->frames_till_gf_update_due = DEFAULT_GF_INTERVAL; michael@0: michael@0: cpi->common.refresh_golden_frame = 1; michael@0: cpi->common.refresh_alt_ref_frame = 1; michael@0: } michael@0: michael@0: michael@0: static int estimate_bits_at_q(int frame_kind, int Q, int MBs, michael@0: double correction_factor) michael@0: { michael@0: int Bpm = (int)(.5 + correction_factor * vp8_bits_per_mb[frame_kind][Q]); michael@0: michael@0: /* Attempt to retain reasonable accuracy without overflow. The cutoff is michael@0: * chosen such that the maximum product of Bpm and MBs fits 31 bits. The michael@0: * largest Bpm takes 20 bits. michael@0: */ michael@0: if (MBs > (1 << 11)) michael@0: return (Bpm >> BPER_MB_NORMBITS) * MBs; michael@0: else michael@0: return (Bpm * MBs) >> BPER_MB_NORMBITS; michael@0: } michael@0: michael@0: michael@0: static void calc_iframe_target_size(VP8_COMP *cpi) michael@0: { michael@0: /* boost defaults to half second */ michael@0: int kf_boost; michael@0: uint64_t target; michael@0: michael@0: /* Clear down mmx registers to allow floating point in what follows */ michael@0: vp8_clear_system_state(); michael@0: michael@0: if (cpi->oxcf.fixed_q >= 0) michael@0: { michael@0: int Q = cpi->oxcf.key_q; michael@0: michael@0: target = estimate_bits_at_q(INTRA_FRAME, Q, cpi->common.MBs, michael@0: cpi->key_frame_rate_correction_factor); michael@0: } michael@0: else if (cpi->pass == 2) michael@0: { michael@0: /* New Two pass RC */ michael@0: target = cpi->per_frame_bandwidth; michael@0: } michael@0: /* First Frame is a special case */ michael@0: else if (cpi->common.current_video_frame == 0) michael@0: { michael@0: /* 1 Pass there is no information on which to base size so use michael@0: * bandwidth per second * fraction of the initial buffer michael@0: * level michael@0: */ michael@0: target = cpi->oxcf.starting_buffer_level / 2; michael@0: michael@0: if(target > cpi->oxcf.target_bandwidth * 3 / 2) michael@0: target = cpi->oxcf.target_bandwidth * 3 / 2; michael@0: } michael@0: else michael@0: { michael@0: /* if this keyframe was forced, use a more recent Q estimate */ michael@0: int Q = (cpi->common.frame_flags & FRAMEFLAGS_KEY) michael@0: ? cpi->avg_frame_qindex : cpi->ni_av_qi; michael@0: michael@0: int initial_boost = 32; /* |3.0 * per_frame_bandwidth| */ michael@0: /* Boost depends somewhat on frame rate: only used for 1 layer case. */ michael@0: if (cpi->oxcf.number_of_layers == 1) { michael@0: kf_boost = MAX(initial_boost, (int)(2 * cpi->output_framerate - 16)); michael@0: } michael@0: else { michael@0: /* Initial factor: set target size to: |3.0 * per_frame_bandwidth|. */ michael@0: kf_boost = initial_boost; michael@0: } michael@0: michael@0: /* adjustment up based on q: this factor ranges from ~1.2 to 2.2. */ michael@0: kf_boost = kf_boost * kf_boost_qadjustment[Q] / 100; michael@0: michael@0: /* frame separation adjustment ( down) */ michael@0: if (cpi->frames_since_key < cpi->output_framerate / 2) michael@0: kf_boost = (int)(kf_boost michael@0: * cpi->frames_since_key / (cpi->output_framerate / 2)); michael@0: michael@0: /* Minimal target size is |2* per_frame_bandwidth|. */ michael@0: if (kf_boost < 16) michael@0: kf_boost = 16; michael@0: michael@0: target = ((16 + kf_boost) * cpi->per_frame_bandwidth) >> 4; michael@0: } michael@0: michael@0: michael@0: if (cpi->oxcf.rc_max_intra_bitrate_pct) michael@0: { michael@0: unsigned int max_rate = cpi->per_frame_bandwidth michael@0: * cpi->oxcf.rc_max_intra_bitrate_pct / 100; michael@0: michael@0: if (target > max_rate) michael@0: target = max_rate; michael@0: } michael@0: michael@0: cpi->this_frame_target = (int)target; michael@0: michael@0: /* TODO: if we separate rate targeting from Q targetting, move this. michael@0: * Reset the active worst quality to the baseline value for key frames. michael@0: */ michael@0: if (cpi->pass != 2) michael@0: cpi->active_worst_quality = cpi->worst_quality; michael@0: michael@0: #if 0 michael@0: { michael@0: FILE *f; michael@0: michael@0: f = fopen("kf_boost.stt", "a"); michael@0: fprintf(f, " %8u %10d %10d %10d\n", michael@0: cpi->common.current_video_frame, cpi->gfu_boost, cpi->baseline_gf_interval, cpi->source_alt_ref_pending); michael@0: michael@0: fclose(f); michael@0: } michael@0: #endif michael@0: } michael@0: michael@0: michael@0: /* Do the best we can to define the parameters for the next GF based on what michael@0: * information we have available. michael@0: */ michael@0: static void calc_gf_params(VP8_COMP *cpi) michael@0: { michael@0: int Q = (cpi->oxcf.fixed_q < 0) ? cpi->last_q[INTER_FRAME] : cpi->oxcf.fixed_q; michael@0: int Boost = 0; michael@0: michael@0: int gf_frame_useage = 0; /* Golden frame useage since last GF */ michael@0: int tot_mbs = cpi->recent_ref_frame_usage[INTRA_FRAME] + michael@0: cpi->recent_ref_frame_usage[LAST_FRAME] + michael@0: cpi->recent_ref_frame_usage[GOLDEN_FRAME] + michael@0: cpi->recent_ref_frame_usage[ALTREF_FRAME]; michael@0: michael@0: int pct_gf_active = (100 * cpi->gf_active_count) / (cpi->common.mb_rows * cpi->common.mb_cols); michael@0: michael@0: if (tot_mbs) michael@0: gf_frame_useage = (cpi->recent_ref_frame_usage[GOLDEN_FRAME] + cpi->recent_ref_frame_usage[ALTREF_FRAME]) * 100 / tot_mbs; michael@0: michael@0: if (pct_gf_active > gf_frame_useage) michael@0: gf_frame_useage = pct_gf_active; michael@0: michael@0: /* Not two pass */ michael@0: if (cpi->pass != 2) michael@0: { michael@0: /* Single Pass lagged mode: TBD */ michael@0: if (0) michael@0: { michael@0: } michael@0: michael@0: /* Single Pass compression: Has to use current and historical data */ michael@0: else michael@0: { michael@0: #if 0 michael@0: /* Experimental code */ michael@0: int index = cpi->one_pass_frame_index; michael@0: int frames_to_scan = (cpi->max_gf_interval <= MAX_LAG_BUFFERS) ? cpi->max_gf_interval : MAX_LAG_BUFFERS; michael@0: michael@0: /* ************** Experimental code - incomplete */ michael@0: /* michael@0: double decay_val = 1.0; michael@0: double IIAccumulator = 0.0; michael@0: double last_iiaccumulator = 0.0; michael@0: double IIRatio; michael@0: michael@0: cpi->one_pass_frame_index = cpi->common.current_video_frame%MAX_LAG_BUFFERS; michael@0: michael@0: for ( i = 0; i < (frames_to_scan - 1); i++ ) michael@0: { michael@0: if ( index < 0 ) michael@0: index = MAX_LAG_BUFFERS; michael@0: index --; michael@0: michael@0: if ( cpi->one_pass_frame_stats[index].frame_coded_error > 0.0 ) michael@0: { michael@0: IIRatio = cpi->one_pass_frame_stats[index].frame_intra_error / cpi->one_pass_frame_stats[index].frame_coded_error; michael@0: michael@0: if ( IIRatio > 30.0 ) michael@0: IIRatio = 30.0; michael@0: } michael@0: else michael@0: IIRatio = 30.0; michael@0: michael@0: IIAccumulator += IIRatio * decay_val; michael@0: michael@0: decay_val = decay_val * cpi->one_pass_frame_stats[index].frame_pcnt_inter; michael@0: michael@0: if ( (i > MIN_GF_INTERVAL) && michael@0: ((IIAccumulator - last_iiaccumulator) < 2.0) ) michael@0: { michael@0: break; michael@0: } michael@0: last_iiaccumulator = IIAccumulator; michael@0: } michael@0: michael@0: Boost = IIAccumulator*100.0/16.0; michael@0: cpi->baseline_gf_interval = i; michael@0: michael@0: */ michael@0: #else michael@0: michael@0: /*************************************************************/ michael@0: /* OLD code */ michael@0: michael@0: /* Adjust boost based upon ambient Q */ michael@0: Boost = GFQ_ADJUSTMENT; michael@0: michael@0: /* Adjust based upon most recently measure intra useage */ michael@0: Boost = Boost * gf_intra_usage_adjustment[(cpi->this_frame_percent_intra < 15) ? cpi->this_frame_percent_intra : 14] / 100; michael@0: michael@0: /* Adjust gf boost based upon GF usage since last GF */ michael@0: Boost = Boost * gf_adjust_table[gf_frame_useage] / 100; michael@0: #endif michael@0: } michael@0: michael@0: /* golden frame boost without recode loop often goes awry. be michael@0: * safe by keeping numbers down. michael@0: */ michael@0: if (!cpi->sf.recode_loop) michael@0: { michael@0: if (cpi->compressor_speed == 2) michael@0: Boost = Boost / 2; michael@0: } michael@0: michael@0: /* Apply an upper limit based on Q for 1 pass encodes */ michael@0: if (Boost > kf_gf_boost_qlimits[Q] && (cpi->pass == 0)) michael@0: Boost = kf_gf_boost_qlimits[Q]; michael@0: michael@0: /* Apply lower limits to boost. */ michael@0: else if (Boost < 110) michael@0: Boost = 110; michael@0: michael@0: /* Note the boost used */ michael@0: cpi->last_boost = Boost; michael@0: michael@0: } michael@0: michael@0: /* Estimate next interval michael@0: * This is updated once the real frame size/boost is known. michael@0: */ michael@0: if (cpi->oxcf.fixed_q == -1) michael@0: { michael@0: if (cpi->pass == 2) /* 2 Pass */ michael@0: { michael@0: cpi->frames_till_gf_update_due = cpi->baseline_gf_interval; michael@0: } michael@0: else /* 1 Pass */ michael@0: { michael@0: cpi->frames_till_gf_update_due = cpi->baseline_gf_interval; michael@0: michael@0: if (cpi->last_boost > 750) michael@0: cpi->frames_till_gf_update_due++; michael@0: michael@0: if (cpi->last_boost > 1000) michael@0: cpi->frames_till_gf_update_due++; michael@0: michael@0: if (cpi->last_boost > 1250) michael@0: cpi->frames_till_gf_update_due++; michael@0: michael@0: if (cpi->last_boost >= 1500) michael@0: cpi->frames_till_gf_update_due ++; michael@0: michael@0: if (gf_interval_table[gf_frame_useage] > cpi->frames_till_gf_update_due) michael@0: cpi->frames_till_gf_update_due = gf_interval_table[gf_frame_useage]; michael@0: michael@0: if (cpi->frames_till_gf_update_due > cpi->max_gf_interval) michael@0: cpi->frames_till_gf_update_due = cpi->max_gf_interval; michael@0: } michael@0: } michael@0: else michael@0: cpi->frames_till_gf_update_due = cpi->baseline_gf_interval; michael@0: michael@0: /* ARF on or off */ michael@0: if (cpi->pass != 2) michael@0: { michael@0: /* For now Alt ref is not allowed except in 2 pass modes. */ michael@0: cpi->source_alt_ref_pending = 0; michael@0: michael@0: /*if ( cpi->oxcf.fixed_q == -1) michael@0: { michael@0: if ( cpi->oxcf.play_alternate && (cpi->last_boost > (100 + (AF_THRESH*cpi->frames_till_gf_update_due)) ) ) michael@0: cpi->source_alt_ref_pending = 1; michael@0: else michael@0: cpi->source_alt_ref_pending = 0; michael@0: }*/ michael@0: } michael@0: } michael@0: michael@0: michael@0: static void calc_pframe_target_size(VP8_COMP *cpi) michael@0: { michael@0: int min_frame_target; michael@0: int old_per_frame_bandwidth = cpi->per_frame_bandwidth; michael@0: michael@0: if ( cpi->current_layer > 0) michael@0: cpi->per_frame_bandwidth = michael@0: cpi->layer_context[cpi->current_layer].avg_frame_size_for_layer; michael@0: michael@0: min_frame_target = 0; michael@0: michael@0: if (cpi->pass == 2) michael@0: { michael@0: min_frame_target = cpi->min_frame_bandwidth; michael@0: michael@0: if (min_frame_target < (cpi->av_per_frame_bandwidth >> 5)) michael@0: min_frame_target = cpi->av_per_frame_bandwidth >> 5; michael@0: } michael@0: else if (min_frame_target < cpi->per_frame_bandwidth / 4) michael@0: min_frame_target = cpi->per_frame_bandwidth / 4; michael@0: michael@0: michael@0: /* Special alt reference frame case */ michael@0: if((cpi->common.refresh_alt_ref_frame) && (cpi->oxcf.number_of_layers == 1)) michael@0: { michael@0: if (cpi->pass == 2) michael@0: { michael@0: /* Per frame bit target for the alt ref frame */ michael@0: cpi->per_frame_bandwidth = cpi->twopass.gf_bits; michael@0: cpi->this_frame_target = cpi->per_frame_bandwidth; michael@0: } michael@0: michael@0: /* One Pass ??? TBD */ michael@0: } michael@0: michael@0: /* Normal frames (gf,and inter) */ michael@0: else michael@0: { michael@0: /* 2 pass */ michael@0: if (cpi->pass == 2) michael@0: { michael@0: cpi->this_frame_target = cpi->per_frame_bandwidth; michael@0: } michael@0: /* 1 pass */ michael@0: else michael@0: { michael@0: int Adjustment; michael@0: /* Make rate adjustment to recover bits spent in key frame michael@0: * Test to see if the key frame inter data rate correction michael@0: * should still be in force michael@0: */ michael@0: if (cpi->kf_overspend_bits > 0) michael@0: { michael@0: Adjustment = (cpi->kf_bitrate_adjustment <= cpi->kf_overspend_bits) ? cpi->kf_bitrate_adjustment : cpi->kf_overspend_bits; michael@0: michael@0: if (Adjustment > (cpi->per_frame_bandwidth - min_frame_target)) michael@0: Adjustment = (cpi->per_frame_bandwidth - min_frame_target); michael@0: michael@0: cpi->kf_overspend_bits -= Adjustment; michael@0: michael@0: /* Calculate an inter frame bandwidth target for the next michael@0: * few frames designed to recover any extra bits spent on michael@0: * the key frame. michael@0: */ michael@0: cpi->this_frame_target = cpi->per_frame_bandwidth - Adjustment; michael@0: michael@0: if (cpi->this_frame_target < min_frame_target) michael@0: cpi->this_frame_target = min_frame_target; michael@0: } michael@0: else michael@0: cpi->this_frame_target = cpi->per_frame_bandwidth; michael@0: michael@0: /* If appropriate make an adjustment to recover bits spent on a michael@0: * recent GF michael@0: */ michael@0: if ((cpi->gf_overspend_bits > 0) && (cpi->this_frame_target > min_frame_target)) michael@0: { michael@0: Adjustment = (cpi->non_gf_bitrate_adjustment <= cpi->gf_overspend_bits) ? cpi->non_gf_bitrate_adjustment : cpi->gf_overspend_bits; michael@0: michael@0: if (Adjustment > (cpi->this_frame_target - min_frame_target)) michael@0: Adjustment = (cpi->this_frame_target - min_frame_target); michael@0: michael@0: cpi->gf_overspend_bits -= Adjustment; michael@0: cpi->this_frame_target -= Adjustment; michael@0: } michael@0: michael@0: /* Apply small + and - boosts for non gf frames */ michael@0: if ((cpi->last_boost > 150) && (cpi->frames_till_gf_update_due > 0) && michael@0: (cpi->current_gf_interval >= (MIN_GF_INTERVAL << 1))) michael@0: { michael@0: /* % Adjustment limited to the range 1% to 10% */ michael@0: Adjustment = (cpi->last_boost - 100) >> 5; michael@0: michael@0: if (Adjustment < 1) michael@0: Adjustment = 1; michael@0: else if (Adjustment > 10) michael@0: Adjustment = 10; michael@0: michael@0: /* Convert to bits */ michael@0: Adjustment = (cpi->this_frame_target * Adjustment) / 100; michael@0: michael@0: if (Adjustment > (cpi->this_frame_target - min_frame_target)) michael@0: Adjustment = (cpi->this_frame_target - min_frame_target); michael@0: michael@0: if (cpi->frames_since_golden == (cpi->current_gf_interval >> 1)) michael@0: cpi->this_frame_target += ((cpi->current_gf_interval - 1) * Adjustment); michael@0: else michael@0: cpi->this_frame_target -= Adjustment; michael@0: } michael@0: } michael@0: } michael@0: michael@0: /* Sanity check that the total sum of adjustments is not above the michael@0: * maximum allowed That is that having allowed for KF and GF penalties michael@0: * we have not pushed the current interframe target to low. If the michael@0: * adjustment we apply here is not capable of recovering all the extra michael@0: * bits we have spent in the KF or GF then the remainder will have to michael@0: * be recovered over a longer time span via other buffer / rate control michael@0: * mechanisms. michael@0: */ michael@0: if (cpi->this_frame_target < min_frame_target) michael@0: cpi->this_frame_target = min_frame_target; michael@0: michael@0: if (!cpi->common.refresh_alt_ref_frame) michael@0: /* Note the baseline target data rate for this inter frame. */ michael@0: cpi->inter_frame_target = cpi->this_frame_target; michael@0: michael@0: /* One Pass specific code */ michael@0: if (cpi->pass == 0) michael@0: { michael@0: /* Adapt target frame size with respect to any buffering constraints: */ michael@0: if (cpi->buffered_mode) michael@0: { michael@0: int one_percent_bits = (int) michael@0: (1 + cpi->oxcf.optimal_buffer_level / 100); michael@0: michael@0: if ((cpi->buffer_level < cpi->oxcf.optimal_buffer_level) || michael@0: (cpi->bits_off_target < cpi->oxcf.optimal_buffer_level)) michael@0: { michael@0: int percent_low = 0; michael@0: michael@0: /* Decide whether or not we need to adjust the frame data michael@0: * rate target. michael@0: * michael@0: * If we are are below the optimal buffer fullness level michael@0: * and adherence to buffering constraints is important to michael@0: * the end usage then adjust the per frame target. michael@0: */ michael@0: if ((cpi->oxcf.end_usage == USAGE_STREAM_FROM_SERVER) && michael@0: (cpi->buffer_level < cpi->oxcf.optimal_buffer_level)) michael@0: { michael@0: percent_low = (int) michael@0: ((cpi->oxcf.optimal_buffer_level - cpi->buffer_level) / michael@0: one_percent_bits); michael@0: } michael@0: /* Are we overshooting the long term clip data rate... */ michael@0: else if (cpi->bits_off_target < 0) michael@0: { michael@0: /* Adjust per frame data target downwards to compensate. */ michael@0: percent_low = (int)(100 * -cpi->bits_off_target / michael@0: (cpi->total_byte_count * 8)); michael@0: } michael@0: michael@0: if (percent_low > cpi->oxcf.under_shoot_pct) michael@0: percent_low = cpi->oxcf.under_shoot_pct; michael@0: else if (percent_low < 0) michael@0: percent_low = 0; michael@0: michael@0: /* lower the target bandwidth for this frame. */ michael@0: cpi->this_frame_target -= michael@0: (cpi->this_frame_target * percent_low) / 200; michael@0: michael@0: /* Are we using allowing control of active_worst_allowed_q michael@0: * according to buffer level. michael@0: */ michael@0: if (cpi->auto_worst_q && cpi->ni_frames > 150) michael@0: { michael@0: int64_t critical_buffer_level; michael@0: michael@0: /* For streaming applications the most important factor is michael@0: * cpi->buffer_level as this takes into account the michael@0: * specified short term buffering constraints. However, michael@0: * hitting the long term clip data rate target is also michael@0: * important. michael@0: */ michael@0: if (cpi->oxcf.end_usage == USAGE_STREAM_FROM_SERVER) michael@0: { michael@0: /* Take the smaller of cpi->buffer_level and michael@0: * cpi->bits_off_target michael@0: */ michael@0: critical_buffer_level = michael@0: (cpi->buffer_level < cpi->bits_off_target) michael@0: ? cpi->buffer_level : cpi->bits_off_target; michael@0: } michael@0: /* For local file playback short term buffering constraints michael@0: * are less of an issue michael@0: */ michael@0: else michael@0: { michael@0: /* Consider only how we are doing for the clip as a michael@0: * whole michael@0: */ michael@0: critical_buffer_level = cpi->bits_off_target; michael@0: } michael@0: michael@0: /* Set the active worst quality based upon the selected michael@0: * buffer fullness number. michael@0: */ michael@0: if (critical_buffer_level < cpi->oxcf.optimal_buffer_level) michael@0: { michael@0: if ( critical_buffer_level > michael@0: (cpi->oxcf.optimal_buffer_level >> 2) ) michael@0: { michael@0: int64_t qadjustment_range = michael@0: cpi->worst_quality - cpi->ni_av_qi; michael@0: int64_t above_base = michael@0: (critical_buffer_level - michael@0: (cpi->oxcf.optimal_buffer_level >> 2)); michael@0: michael@0: /* Step active worst quality down from michael@0: * cpi->ni_av_qi when (critical_buffer_level == michael@0: * cpi->optimal_buffer_level) to michael@0: * cpi->worst_quality when michael@0: * (critical_buffer_level == michael@0: * cpi->optimal_buffer_level >> 2) michael@0: */ michael@0: cpi->active_worst_quality = michael@0: cpi->worst_quality - michael@0: (int)((qadjustment_range * above_base) / michael@0: (cpi->oxcf.optimal_buffer_level*3>>2)); michael@0: } michael@0: else michael@0: { michael@0: cpi->active_worst_quality = cpi->worst_quality; michael@0: } michael@0: } michael@0: else michael@0: { michael@0: cpi->active_worst_quality = cpi->ni_av_qi; michael@0: } michael@0: } michael@0: else michael@0: { michael@0: cpi->active_worst_quality = cpi->worst_quality; michael@0: } michael@0: } michael@0: else michael@0: { michael@0: int percent_high = 0; michael@0: michael@0: if ((cpi->oxcf.end_usage == USAGE_STREAM_FROM_SERVER) michael@0: && (cpi->buffer_level > cpi->oxcf.optimal_buffer_level)) michael@0: { michael@0: percent_high = (int)((cpi->buffer_level michael@0: - cpi->oxcf.optimal_buffer_level) michael@0: / one_percent_bits); michael@0: } michael@0: else if (cpi->bits_off_target > cpi->oxcf.optimal_buffer_level) michael@0: { michael@0: percent_high = (int)((100 * cpi->bits_off_target) michael@0: / (cpi->total_byte_count * 8)); michael@0: } michael@0: michael@0: if (percent_high > cpi->oxcf.over_shoot_pct) michael@0: percent_high = cpi->oxcf.over_shoot_pct; michael@0: else if (percent_high < 0) michael@0: percent_high = 0; michael@0: michael@0: cpi->this_frame_target += (cpi->this_frame_target * michael@0: percent_high) / 200; michael@0: michael@0: /* Are we allowing control of active_worst_allowed_q according michael@0: * to buffer level. michael@0: */ michael@0: if (cpi->auto_worst_q && cpi->ni_frames > 150) michael@0: { michael@0: /* When using the relaxed buffer model stick to the michael@0: * user specified value michael@0: */ michael@0: cpi->active_worst_quality = cpi->ni_av_qi; michael@0: } michael@0: else michael@0: { michael@0: cpi->active_worst_quality = cpi->worst_quality; michael@0: } michael@0: } michael@0: michael@0: /* Set active_best_quality to prevent quality rising too high */ michael@0: cpi->active_best_quality = cpi->best_quality; michael@0: michael@0: /* Worst quality obviously must not be better than best quality */ michael@0: if (cpi->active_worst_quality <= cpi->active_best_quality) michael@0: cpi->active_worst_quality = cpi->active_best_quality + 1; michael@0: michael@0: if(cpi->active_worst_quality > 127) michael@0: cpi->active_worst_quality = 127; michael@0: } michael@0: /* Unbuffered mode (eg. video conferencing) */ michael@0: else michael@0: { michael@0: /* Set the active worst quality */ michael@0: cpi->active_worst_quality = cpi->worst_quality; michael@0: } michael@0: michael@0: /* Special trap for constrained quality mode michael@0: * "active_worst_quality" may never drop below cq level michael@0: * for any frame type. michael@0: */ michael@0: if ( cpi->oxcf.end_usage == USAGE_CONSTRAINED_QUALITY && michael@0: cpi->active_worst_quality < cpi->cq_target_quality) michael@0: { michael@0: cpi->active_worst_quality = cpi->cq_target_quality; michael@0: } michael@0: } michael@0: michael@0: /* Test to see if we have to drop a frame michael@0: * The auto-drop frame code is only used in buffered mode. michael@0: * In unbufferd mode (eg vide conferencing) the descision to michael@0: * code or drop a frame is made outside the codec in response to real michael@0: * world comms or buffer considerations. michael@0: */ michael@0: if (cpi->drop_frames_allowed && michael@0: (cpi->oxcf.end_usage == USAGE_STREAM_FROM_SERVER) && michael@0: ((cpi->common.frame_type != KEY_FRAME))) michael@0: { michael@0: /* Check for a buffer underun-crisis in which case we have to drop michael@0: * a frame michael@0: */ michael@0: if ((cpi->buffer_level < 0)) michael@0: { michael@0: #if 0 michael@0: FILE *f = fopen("dec.stt", "a"); michael@0: fprintf(f, "%10d %10d %10d %10d ***** BUFFER EMPTY\n", michael@0: (int) cpi->common.current_video_frame, michael@0: cpi->decimation_factor, cpi->common.horiz_scale, michael@0: (cpi->buffer_level * 100) / cpi->oxcf.optimal_buffer_level); michael@0: fclose(f); michael@0: #endif michael@0: cpi->drop_frame = 1; michael@0: michael@0: /* Update the buffer level variable. */ michael@0: cpi->bits_off_target += cpi->av_per_frame_bandwidth; michael@0: if (cpi->bits_off_target > cpi->oxcf.maximum_buffer_size) michael@0: cpi->bits_off_target = (int)cpi->oxcf.maximum_buffer_size; michael@0: cpi->buffer_level = cpi->bits_off_target; michael@0: michael@0: if (cpi->oxcf.number_of_layers > 1) { michael@0: unsigned int i; michael@0: michael@0: // Propagate bits saved by dropping the frame to higher layers. michael@0: for (i = cpi->current_layer + 1; i < cpi->oxcf.number_of_layers; michael@0: i++) { michael@0: LAYER_CONTEXT *lc = &cpi->layer_context[i]; michael@0: lc->bits_off_target += (int)(lc->target_bandwidth / michael@0: lc->framerate); michael@0: if (lc->bits_off_target > lc->maximum_buffer_size) michael@0: lc->bits_off_target = lc->maximum_buffer_size; michael@0: lc->buffer_level = lc->bits_off_target; michael@0: } michael@0: } michael@0: } michael@0: } michael@0: michael@0: /* Adjust target frame size for Golden Frames: */ michael@0: if (cpi->oxcf.error_resilient_mode == 0 && michael@0: (cpi->frames_till_gf_update_due == 0) && !cpi->drop_frame) michael@0: { michael@0: int Q = (cpi->oxcf.fixed_q < 0) ? cpi->last_q[INTER_FRAME] : cpi->oxcf.fixed_q; michael@0: michael@0: int gf_frame_useage = 0; /* Golden frame useage since last GF */ michael@0: int tot_mbs = cpi->recent_ref_frame_usage[INTRA_FRAME] + michael@0: cpi->recent_ref_frame_usage[LAST_FRAME] + michael@0: cpi->recent_ref_frame_usage[GOLDEN_FRAME] + michael@0: cpi->recent_ref_frame_usage[ALTREF_FRAME]; michael@0: michael@0: int pct_gf_active = (100 * cpi->gf_active_count) / (cpi->common.mb_rows * cpi->common.mb_cols); michael@0: michael@0: if (tot_mbs) michael@0: gf_frame_useage = (cpi->recent_ref_frame_usage[GOLDEN_FRAME] + cpi->recent_ref_frame_usage[ALTREF_FRAME]) * 100 / tot_mbs; michael@0: michael@0: if (pct_gf_active > gf_frame_useage) michael@0: gf_frame_useage = pct_gf_active; michael@0: michael@0: /* Is a fixed manual GF frequency being used */ michael@0: if (cpi->auto_gold) michael@0: { michael@0: /* For one pass throw a GF if recent frame intra useage is michael@0: * low or the GF useage is high michael@0: */ michael@0: if ((cpi->pass == 0) && (cpi->this_frame_percent_intra < 15 || gf_frame_useage >= 5)) michael@0: cpi->common.refresh_golden_frame = 1; michael@0: michael@0: /* Two pass GF descision */ michael@0: else if (cpi->pass == 2) michael@0: cpi->common.refresh_golden_frame = 1; michael@0: } michael@0: michael@0: #if 0 michael@0: michael@0: /* Debug stats */ michael@0: if (0) michael@0: { michael@0: FILE *f; michael@0: michael@0: f = fopen("gf_useaget.stt", "a"); michael@0: fprintf(f, " %8ld %10ld %10ld %10ld %10ld\n", michael@0: cpi->common.current_video_frame, cpi->gfu_boost, GFQ_ADJUSTMENT, cpi->gfu_boost, gf_frame_useage); michael@0: fclose(f); michael@0: } michael@0: michael@0: #endif michael@0: michael@0: if (cpi->common.refresh_golden_frame == 1) michael@0: { michael@0: #if 0 michael@0: michael@0: if (0) michael@0: { michael@0: FILE *f; michael@0: michael@0: f = fopen("GFexit.stt", "a"); michael@0: fprintf(f, "%8ld GF coded\n", cpi->common.current_video_frame); michael@0: fclose(f); michael@0: } michael@0: michael@0: #endif michael@0: michael@0: if (cpi->auto_adjust_gold_quantizer) michael@0: { michael@0: calc_gf_params(cpi); michael@0: } michael@0: michael@0: /* If we are using alternate ref instead of gf then do not apply the michael@0: * boost It will instead be applied to the altref update Jims michael@0: * modified boost michael@0: */ michael@0: if (!cpi->source_alt_ref_active) michael@0: { michael@0: if (cpi->oxcf.fixed_q < 0) michael@0: { michael@0: if (cpi->pass == 2) michael@0: { michael@0: /* The spend on the GF is defined in the two pass michael@0: * code for two pass encodes michael@0: */ michael@0: cpi->this_frame_target = cpi->per_frame_bandwidth; michael@0: } michael@0: else michael@0: { michael@0: int Boost = cpi->last_boost; michael@0: int frames_in_section = cpi->frames_till_gf_update_due + 1; michael@0: int allocation_chunks = (frames_in_section * 100) + (Boost - 100); michael@0: int bits_in_section = cpi->inter_frame_target * frames_in_section; michael@0: michael@0: /* Normalize Altboost and allocations chunck down to michael@0: * prevent overflow michael@0: */ michael@0: while (Boost > 1000) michael@0: { michael@0: Boost /= 2; michael@0: allocation_chunks /= 2; michael@0: } michael@0: michael@0: /* Avoid loss of precision but avoid overflow */ michael@0: if ((bits_in_section >> 7) > allocation_chunks) michael@0: cpi->this_frame_target = Boost * (bits_in_section / allocation_chunks); michael@0: else michael@0: cpi->this_frame_target = (Boost * bits_in_section) / allocation_chunks; michael@0: } michael@0: } michael@0: else michael@0: cpi->this_frame_target = michael@0: (estimate_bits_at_q(1, Q, cpi->common.MBs, 1.0) michael@0: * cpi->last_boost) / 100; michael@0: michael@0: } michael@0: /* If there is an active ARF at this location use the minimum michael@0: * bits on this frame even if it is a contructed arf. michael@0: * The active maximum quantizer insures that an appropriate michael@0: * number of bits will be spent if needed for contstructed ARFs. michael@0: */ michael@0: else michael@0: { michael@0: cpi->this_frame_target = 0; michael@0: } michael@0: michael@0: cpi->current_gf_interval = cpi->frames_till_gf_update_due; michael@0: michael@0: } michael@0: } michael@0: michael@0: cpi->per_frame_bandwidth = old_per_frame_bandwidth; michael@0: } michael@0: michael@0: michael@0: void vp8_update_rate_correction_factors(VP8_COMP *cpi, int damp_var) michael@0: { michael@0: int Q = cpi->common.base_qindex; michael@0: int correction_factor = 100; michael@0: double rate_correction_factor; michael@0: double adjustment_limit; michael@0: michael@0: int projected_size_based_on_q = 0; michael@0: michael@0: /* Clear down mmx registers to allow floating point in what follows */ michael@0: vp8_clear_system_state(); michael@0: michael@0: if (cpi->common.frame_type == KEY_FRAME) michael@0: { michael@0: rate_correction_factor = cpi->key_frame_rate_correction_factor; michael@0: } michael@0: else michael@0: { michael@0: if (cpi->oxcf.number_of_layers == 1 && michael@0: (cpi->common.refresh_alt_ref_frame || michael@0: cpi->common.refresh_golden_frame)) michael@0: rate_correction_factor = cpi->gf_rate_correction_factor; michael@0: else michael@0: rate_correction_factor = cpi->rate_correction_factor; michael@0: } michael@0: michael@0: /* Work out how big we would have expected the frame to be at this Q michael@0: * given the current correction factor. Stay in double to avoid int michael@0: * overflow when values are large michael@0: */ michael@0: projected_size_based_on_q = (int)(((.5 + rate_correction_factor * vp8_bits_per_mb[cpi->common.frame_type][Q]) * cpi->common.MBs) / (1 << BPER_MB_NORMBITS)); michael@0: michael@0: /* Make some allowance for cpi->zbin_over_quant */ michael@0: if (cpi->mb.zbin_over_quant > 0) michael@0: { michael@0: int Z = cpi->mb.zbin_over_quant; michael@0: double Factor = 0.99; michael@0: double factor_adjustment = 0.01 / 256.0; michael@0: michael@0: while (Z > 0) michael@0: { michael@0: Z --; michael@0: projected_size_based_on_q = michael@0: (int)(Factor * projected_size_based_on_q); michael@0: Factor += factor_adjustment; michael@0: michael@0: if (Factor >= 0.999) michael@0: Factor = 0.999; michael@0: } michael@0: } michael@0: michael@0: /* Work out a size correction factor. */ michael@0: if (projected_size_based_on_q > 0) michael@0: correction_factor = (100 * cpi->projected_frame_size) / projected_size_based_on_q; michael@0: michael@0: /* More heavily damped adjustment used if we have been oscillating michael@0: * either side of target michael@0: */ michael@0: switch (damp_var) michael@0: { michael@0: case 0: michael@0: adjustment_limit = 0.75; michael@0: break; michael@0: case 1: michael@0: adjustment_limit = 0.375; michael@0: break; michael@0: case 2: michael@0: default: michael@0: adjustment_limit = 0.25; michael@0: break; michael@0: } michael@0: michael@0: if (correction_factor > 102) michael@0: { michael@0: /* We are not already at the worst allowable quality */ michael@0: correction_factor = (int)(100.5 + ((correction_factor - 100) * adjustment_limit)); michael@0: rate_correction_factor = ((rate_correction_factor * correction_factor) / 100); michael@0: michael@0: /* Keep rate_correction_factor within limits */ michael@0: if (rate_correction_factor > MAX_BPB_FACTOR) michael@0: rate_correction_factor = MAX_BPB_FACTOR; michael@0: } michael@0: else if (correction_factor < 99) michael@0: { michael@0: /* We are not already at the best allowable quality */ michael@0: correction_factor = (int)(100.5 - ((100 - correction_factor) * adjustment_limit)); michael@0: rate_correction_factor = ((rate_correction_factor * correction_factor) / 100); michael@0: michael@0: /* Keep rate_correction_factor within limits */ michael@0: if (rate_correction_factor < MIN_BPB_FACTOR) michael@0: rate_correction_factor = MIN_BPB_FACTOR; michael@0: } michael@0: michael@0: if (cpi->common.frame_type == KEY_FRAME) michael@0: cpi->key_frame_rate_correction_factor = rate_correction_factor; michael@0: else michael@0: { michael@0: if (cpi->oxcf.number_of_layers == 1 && michael@0: (cpi->common.refresh_alt_ref_frame || michael@0: cpi->common.refresh_golden_frame)) michael@0: cpi->gf_rate_correction_factor = rate_correction_factor; michael@0: else michael@0: cpi->rate_correction_factor = rate_correction_factor; michael@0: } michael@0: } michael@0: michael@0: michael@0: int vp8_regulate_q(VP8_COMP *cpi, int target_bits_per_frame) michael@0: { michael@0: int Q = cpi->active_worst_quality; michael@0: michael@0: /* Reset Zbin OQ value */ michael@0: cpi->mb.zbin_over_quant = 0; michael@0: michael@0: if (cpi->oxcf.fixed_q >= 0) michael@0: { michael@0: Q = cpi->oxcf.fixed_q; michael@0: michael@0: if (cpi->common.frame_type == KEY_FRAME) michael@0: { michael@0: Q = cpi->oxcf.key_q; michael@0: } michael@0: else if (cpi->oxcf.number_of_layers == 1 && michael@0: cpi->common.refresh_alt_ref_frame) michael@0: { michael@0: Q = cpi->oxcf.alt_q; michael@0: } michael@0: else if (cpi->oxcf.number_of_layers == 1 && michael@0: cpi->common.refresh_golden_frame) michael@0: { michael@0: Q = cpi->oxcf.gold_q; michael@0: } michael@0: michael@0: } michael@0: else michael@0: { michael@0: int i; michael@0: int last_error = INT_MAX; michael@0: int target_bits_per_mb; michael@0: int bits_per_mb_at_this_q; michael@0: double correction_factor; michael@0: michael@0: /* Select the appropriate correction factor based upon type of frame. */ michael@0: if (cpi->common.frame_type == KEY_FRAME) michael@0: correction_factor = cpi->key_frame_rate_correction_factor; michael@0: else michael@0: { michael@0: if (cpi->oxcf.number_of_layers == 1 && michael@0: (cpi->common.refresh_alt_ref_frame || michael@0: cpi->common.refresh_golden_frame)) michael@0: correction_factor = cpi->gf_rate_correction_factor; michael@0: else michael@0: correction_factor = cpi->rate_correction_factor; michael@0: } michael@0: michael@0: /* Calculate required scaling factor based on target frame size and michael@0: * size of frame produced using previous Q michael@0: */ michael@0: if (target_bits_per_frame >= (INT_MAX >> BPER_MB_NORMBITS)) michael@0: /* Case where we would overflow int */ michael@0: target_bits_per_mb = (target_bits_per_frame / cpi->common.MBs) << BPER_MB_NORMBITS; michael@0: else michael@0: target_bits_per_mb = (target_bits_per_frame << BPER_MB_NORMBITS) / cpi->common.MBs; michael@0: michael@0: i = cpi->active_best_quality; michael@0: michael@0: do michael@0: { michael@0: bits_per_mb_at_this_q = (int)(.5 + correction_factor * vp8_bits_per_mb[cpi->common.frame_type][i]); michael@0: michael@0: if (bits_per_mb_at_this_q <= target_bits_per_mb) michael@0: { michael@0: if ((target_bits_per_mb - bits_per_mb_at_this_q) <= last_error) michael@0: Q = i; michael@0: else michael@0: Q = i - 1; michael@0: michael@0: break; michael@0: } michael@0: else michael@0: last_error = bits_per_mb_at_this_q - target_bits_per_mb; michael@0: } michael@0: while (++i <= cpi->active_worst_quality); michael@0: michael@0: michael@0: /* If we are at MAXQ then enable Q over-run which seeks to claw michael@0: * back additional bits through things like the RD multiplier michael@0: * and zero bin size. michael@0: */ michael@0: if (Q >= MAXQ) michael@0: { michael@0: int zbin_oqmax; michael@0: michael@0: double Factor = 0.99; michael@0: double factor_adjustment = 0.01 / 256.0; michael@0: michael@0: if (cpi->common.frame_type == KEY_FRAME) michael@0: zbin_oqmax = 0; michael@0: else if (cpi->oxcf.number_of_layers == 1 && michael@0: (cpi->common.refresh_alt_ref_frame || michael@0: (cpi->common.refresh_golden_frame && michael@0: !cpi->source_alt_ref_active))) michael@0: zbin_oqmax = 16; michael@0: else michael@0: zbin_oqmax = ZBIN_OQ_MAX; michael@0: michael@0: /*{ michael@0: double Factor = (double)target_bits_per_mb/(double)bits_per_mb_at_this_q; michael@0: double Oq; michael@0: michael@0: Factor = Factor/1.2683; michael@0: michael@0: Oq = pow( Factor, (1.0/-0.165) ); michael@0: michael@0: if ( Oq > zbin_oqmax ) michael@0: Oq = zbin_oqmax; michael@0: michael@0: cpi->zbin_over_quant = (int)Oq; michael@0: }*/ michael@0: michael@0: /* Each incrment in the zbin is assumed to have a fixed effect michael@0: * on bitrate. This is not of course true. The effect will be michael@0: * highly clip dependent and may well have sudden steps. The michael@0: * idea here is to acheive higher effective quantizers than the michael@0: * normal maximum by expanding the zero bin and hence michael@0: * decreasing the number of low magnitude non zero coefficients. michael@0: */ michael@0: while (cpi->mb.zbin_over_quant < zbin_oqmax) michael@0: { michael@0: cpi->mb.zbin_over_quant ++; michael@0: michael@0: if (cpi->mb.zbin_over_quant > zbin_oqmax) michael@0: cpi->mb.zbin_over_quant = zbin_oqmax; michael@0: michael@0: /* Adjust bits_per_mb_at_this_q estimate */ michael@0: bits_per_mb_at_this_q = (int)(Factor * bits_per_mb_at_this_q); michael@0: Factor += factor_adjustment; michael@0: michael@0: if (Factor >= 0.999) michael@0: Factor = 0.999; michael@0: michael@0: /* Break out if we get down to the target rate */ michael@0: if (bits_per_mb_at_this_q <= target_bits_per_mb) michael@0: break; michael@0: } michael@0: michael@0: } michael@0: } michael@0: michael@0: return Q; michael@0: } michael@0: michael@0: michael@0: static int estimate_keyframe_frequency(VP8_COMP *cpi) michael@0: { michael@0: int i; michael@0: michael@0: /* Average key frame frequency */ michael@0: int av_key_frame_frequency = 0; michael@0: michael@0: /* First key frame at start of sequence is a special case. We have no michael@0: * frequency data. michael@0: */ michael@0: if (cpi->key_frame_count == 1) michael@0: { michael@0: /* Assume a default of 1 kf every 2 seconds, or the max kf interval, michael@0: * whichever is smaller. michael@0: */ michael@0: int key_freq = cpi->oxcf.key_freq>0 ? cpi->oxcf.key_freq : 1; michael@0: av_key_frame_frequency = 1 + (int)cpi->output_framerate * 2; michael@0: michael@0: if (cpi->oxcf.auto_key && av_key_frame_frequency > key_freq) michael@0: av_key_frame_frequency = key_freq; michael@0: michael@0: cpi->prior_key_frame_distance[KEY_FRAME_CONTEXT - 1] michael@0: = av_key_frame_frequency; michael@0: } michael@0: else michael@0: { michael@0: unsigned int total_weight = 0; michael@0: int last_kf_interval = michael@0: (cpi->frames_since_key > 0) ? cpi->frames_since_key : 1; michael@0: michael@0: /* reset keyframe context and calculate weighted average of last michael@0: * KEY_FRAME_CONTEXT keyframes michael@0: */ michael@0: for (i = 0; i < KEY_FRAME_CONTEXT; i++) michael@0: { michael@0: if (i < KEY_FRAME_CONTEXT - 1) michael@0: cpi->prior_key_frame_distance[i] michael@0: = cpi->prior_key_frame_distance[i+1]; michael@0: else michael@0: cpi->prior_key_frame_distance[i] = last_kf_interval; michael@0: michael@0: av_key_frame_frequency += prior_key_frame_weight[i] michael@0: * cpi->prior_key_frame_distance[i]; michael@0: total_weight += prior_key_frame_weight[i]; michael@0: } michael@0: michael@0: av_key_frame_frequency /= total_weight; michael@0: michael@0: } michael@0: // TODO (marpan): Given the checks above, |av_key_frame_frequency| michael@0: // should always be above 0. But for now we keep the sanity check in. michael@0: if (av_key_frame_frequency == 0) michael@0: av_key_frame_frequency = 1; michael@0: return av_key_frame_frequency; michael@0: } michael@0: michael@0: michael@0: void vp8_adjust_key_frame_context(VP8_COMP *cpi) michael@0: { michael@0: /* Clear down mmx registers to allow floating point in what follows */ michael@0: vp8_clear_system_state(); michael@0: michael@0: /* Do we have any key frame overspend to recover? */ michael@0: /* Two-pass overspend handled elsewhere. */ michael@0: if ((cpi->pass != 2) michael@0: && (cpi->projected_frame_size > cpi->per_frame_bandwidth)) michael@0: { michael@0: int overspend; michael@0: michael@0: /* Update the count of key frame overspend to be recovered in michael@0: * subsequent frames. A portion of the KF overspend is treated as gf michael@0: * overspend (and hence recovered more quickly) as the kf is also a michael@0: * gf. Otherwise the few frames following each kf tend to get more michael@0: * bits allocated than those following other gfs. michael@0: */ michael@0: overspend = (cpi->projected_frame_size - cpi->per_frame_bandwidth); michael@0: michael@0: if (cpi->oxcf.number_of_layers > 1) michael@0: cpi->kf_overspend_bits += overspend; michael@0: else michael@0: { michael@0: cpi->kf_overspend_bits += overspend * 7 / 8; michael@0: cpi->gf_overspend_bits += overspend * 1 / 8; michael@0: } michael@0: michael@0: /* Work out how much to try and recover per frame. */ michael@0: cpi->kf_bitrate_adjustment = cpi->kf_overspend_bits michael@0: / estimate_keyframe_frequency(cpi); michael@0: } michael@0: michael@0: cpi->frames_since_key = 0; michael@0: cpi->key_frame_count++; michael@0: } michael@0: michael@0: michael@0: void vp8_compute_frame_size_bounds(VP8_COMP *cpi, int *frame_under_shoot_limit, int *frame_over_shoot_limit) michael@0: { michael@0: /* Set-up bounds on acceptable frame size: */ michael@0: if (cpi->oxcf.fixed_q >= 0) michael@0: { michael@0: /* Fixed Q scenario: frame size never outranges target michael@0: * (there is no target!) michael@0: */ michael@0: *frame_under_shoot_limit = 0; michael@0: *frame_over_shoot_limit = INT_MAX; michael@0: } michael@0: else michael@0: { michael@0: if (cpi->common.frame_type == KEY_FRAME) michael@0: { michael@0: *frame_over_shoot_limit = cpi->this_frame_target * 9 / 8; michael@0: *frame_under_shoot_limit = cpi->this_frame_target * 7 / 8; michael@0: } michael@0: else michael@0: { michael@0: if (cpi->oxcf.number_of_layers > 1 || michael@0: cpi->common.refresh_alt_ref_frame || michael@0: cpi->common.refresh_golden_frame) michael@0: { michael@0: *frame_over_shoot_limit = cpi->this_frame_target * 9 / 8; michael@0: *frame_under_shoot_limit = cpi->this_frame_target * 7 / 8; michael@0: } michael@0: else michael@0: { michael@0: /* For CBR take buffer fullness into account */ michael@0: if (cpi->oxcf.end_usage == USAGE_STREAM_FROM_SERVER) michael@0: { michael@0: if (cpi->buffer_level >= ((cpi->oxcf.optimal_buffer_level + cpi->oxcf.maximum_buffer_size) >> 1)) michael@0: { michael@0: /* Buffer is too full so relax overshoot and tighten michael@0: * undershoot michael@0: */ michael@0: *frame_over_shoot_limit = cpi->this_frame_target * 12 / 8; michael@0: *frame_under_shoot_limit = cpi->this_frame_target * 6 / 8; michael@0: } michael@0: else if (cpi->buffer_level <= (cpi->oxcf.optimal_buffer_level >> 1)) michael@0: { michael@0: /* Buffer is too low so relax undershoot and tighten michael@0: * overshoot michael@0: */ michael@0: *frame_over_shoot_limit = cpi->this_frame_target * 10 / 8; michael@0: *frame_under_shoot_limit = cpi->this_frame_target * 4 / 8; michael@0: } michael@0: else michael@0: { michael@0: *frame_over_shoot_limit = cpi->this_frame_target * 11 / 8; michael@0: *frame_under_shoot_limit = cpi->this_frame_target * 5 / 8; michael@0: } michael@0: } michael@0: /* VBR and CQ mode */ michael@0: /* Note that tighter restrictions here can help quality michael@0: * but hurt encode speed michael@0: */ michael@0: else michael@0: { michael@0: /* Stron overshoot limit for constrained quality */ michael@0: if (cpi->oxcf.end_usage == USAGE_CONSTRAINED_QUALITY) michael@0: { michael@0: *frame_over_shoot_limit = cpi->this_frame_target * 11 / 8; michael@0: *frame_under_shoot_limit = cpi->this_frame_target * 2 / 8; michael@0: } michael@0: else michael@0: { michael@0: *frame_over_shoot_limit = cpi->this_frame_target * 11 / 8; michael@0: *frame_under_shoot_limit = cpi->this_frame_target * 5 / 8; michael@0: } michael@0: } michael@0: } michael@0: } michael@0: michael@0: /* For very small rate targets where the fractional adjustment michael@0: * (eg * 7/8) may be tiny make sure there is at least a minimum michael@0: * range. michael@0: */ michael@0: *frame_over_shoot_limit += 200; michael@0: *frame_under_shoot_limit -= 200; michael@0: if ( *frame_under_shoot_limit < 0 ) michael@0: *frame_under_shoot_limit = 0; michael@0: michael@0: } michael@0: } michael@0: michael@0: michael@0: /* return of 0 means drop frame */ michael@0: int vp8_pick_frame_size(VP8_COMP *cpi) michael@0: { michael@0: VP8_COMMON *cm = &cpi->common; michael@0: michael@0: if (cm->frame_type == KEY_FRAME) michael@0: calc_iframe_target_size(cpi); michael@0: else michael@0: { michael@0: calc_pframe_target_size(cpi); michael@0: michael@0: /* Check if we're dropping the frame: */ michael@0: if (cpi->drop_frame) michael@0: { michael@0: cpi->drop_frame = 0; michael@0: return 0; michael@0: } michael@0: } michael@0: return 1; michael@0: }