35 #undef NDEBUG // Always check asserts, the speed effect is far too small to disable them.
39 #define M_E 2.718281828
46 snprintf(s->
avctx->
stats_out, 256,
"in:%d out:%d type:%d q:%d itex:%d ptex:%d mv:%d misc:%d fcode:%d bcode:%d mc-var:%d var:%d icount:%d skipcount:%d hbits:%d;\n",
97 static double (*
const func1[])(
void *, double)={
102 static const char *
const func1_names[]={
161 next= strchr(p,
';');
166 e= sscanf(p,
" in:%d ", &picture_number);
168 assert(picture_number >= 0);
169 assert(picture_number < rcc->num_entries);
170 rce= &rcc->
entry[picture_number];
172 e+=sscanf(p,
" in:%*d out:%*d type:%d q:%f itex:%d ptex:%d mv:%d misc:%d fcode:%d bcode:%d mc-var:%d var:%d icount:%d skipcount:%d hbits:%d",
210 for(i=0; i<60*30; i++){
274 buffer_size, rcc->
buffer_index, frame_size, min_rate, max_rate);
288 int stuffing= ceil((rcc->
buffer_index - buffer_size)/8);
311 const double mb_num= s->
mb_num;
350 if(bits<0.0) bits=0.0;
357 if(rco[i].
end_frame < frame_num)
continue;
360 bits=
qp2bits(rce, rco[i].qscale);
395 if (q > last_q + maxdiff) q= last_q + maxdiff;
396 else if(q < last_q - maxdiff) q= last_q - maxdiff;
414 assert(qmin <= qmax);
427 if(qmax<qmin) qmax= qmin;
454 double d= 2*(buffer_size - expected_size)/buffer_size;
456 else if(d<0.0001) d=0.0001;
469 double d= 2*expected_size/buffer_size;
471 else if(d<0.0001) d=0.0001;
483 av_dlog(s,
"q:%f max:%f min:%f size:%f index:%f agr:%f\n",
488 else if(q>qmax) q=qmax;
490 double min2= log(qmin);
491 double max2= log(qmax);
494 q= (q - min2)/(max2-min2) - 0.5;
496 q= 1.0/(1.0 + exp(q));
497 q= q*(max2-min2) + min2;
515 double new_coeff= size*q / (var + 1);
521 p->
coeff+= new_coeff;
542 for(i=0; i<s->
mb_num; i++){
544 float temp_cplx= sqrt(pic->
mc_mb_var[mb_xy]);
545 float spat_cplx= sqrt(pic->
mb_var[mb_xy]);
546 const int lumi= pic->
mb_mean[mb_xy];
547 float bits, cplx, factor;
551 float mb_factor = 0.0;
552 if(spat_cplx < 4) spat_cplx= 4;
553 if(temp_cplx < 4) temp_cplx= 4;
557 factor= 1.0 + p_masking;
560 factor= pow(temp_cplx, - temp_cplx_masking);
562 factor*=pow(spat_cplx, - spatial_cplx_masking);
565 factor*= (1.0 - (lumi-128)*(lumi-128)*lumi_masking);
567 factor*= (1.0 - (lumi-128)*(lumi-128)*dark_masking);
569 if(mb_x < mb_width/5){
570 mb_distance = mb_width/5 - mb_x;
571 mb_factor = (float)mb_distance / (
float)(mb_width/5);
572 }
else if(mb_x > 4*mb_width/5){
573 mb_distance = mb_x - 4*mb_width/5;
574 mb_factor = (float)mb_distance / (
float)(mb_width/5);
576 if(mb_y < mb_height/5){
577 mb_distance = mb_height/5 - mb_y;
578 mb_factor =
FFMAX(mb_factor, (
float)mb_distance / (
float)(mb_height/5));
579 }
else if(mb_y > 4*mb_height/5){
580 mb_distance = mb_y - 4*mb_height/5;
581 mb_factor =
FFMAX(mb_factor, (
float)mb_distance / (
float)(mb_height/5));
584 factor*= 1.0 - border_masking*mb_factor;
586 if(factor<0.00001) factor= 0.00001;
597 float factor= bits_sum/cplx_sum;
598 for(i=0; i<s->
mb_num; i++){
599 float newq= q*cplx_tab[i]/bits_tab[i];
603 bits_sum -= bits_tab[i];
604 cplx_sum -= cplx_tab[i]*q/qmax;
606 else if(newq < qmin){
607 bits_sum -= bits_tab[i];
608 cplx_sum -= cplx_tab[i]*q/qmin;
611 if(bits_sum < 0.001) bits_sum= 0.001;
612 if(cplx_sum < 0.001) cplx_sum= 0.001;
615 for(i=0; i<s->
mb_num; i++){
617 float newq= q*cplx_tab[i]/bits_tab[i];
621 newq*= bits_sum/cplx_sum;
624 intq= (int)(newq + 0.5);
626 if (intq > qmax) intq= qmax;
627 else if(intq < qmin) intq= qmin;
637 rce= &rcc->
entry[picture_number];
648 float br_compensation;
673 if(picture_number>2 && !dry_run){
679 assert(picture_number>=0);
680 assert(picture_number<rcc->num_entries);
681 rce= &rcc->
entry[picture_number];
695 wanted_bits= (uint64_t)(s->
bit_rate*(
double)picture_number/fps);
697 wanted_bits = (uint64_t)(s->
bit_rate*(
double)dts_pic->
f.
pts / fps);
702 if(br_compensation<=0.0) br_compensation=0.001;
713 br_compensation, s->
frame_bits, var, pict_type);
745 q=
get_qscale(s, rce, rate_factor, picture_number);
771 av_log(s->
avctx,
AV_LOG_DEBUG,
"%c qp:%d<%2.1f<%d %d want:%d total:%d comp:%f st_q:%2.2f size:%d var:%d/%d br:%d fps:%d\n",
778 else if(q>qmax) q=qmax;
802 double complexity[5]={0,0,0,0,0};
803 uint64_t const_bits[5]={0,0,0,0,0};
804 uint64_t all_const_bits;
806 double rate_factor=0;
809 const int filter_size= (int)(a->
qblur*4) | 1;
810 double expected_bits;
811 double *qscale, *blurred_qscale, qscale_sum;
828 if(all_available_bits < all_const_bits){
837 for(step=256*256; step>0.0000001; step*=0.5){
849 assert(filter_size%2==1);
863 double q=0.0, sum=0.0;
865 for(j=0; j<filter_size; j++){
866 int index= i+j-filter_size/2;
870 if(index < 0 || index >= rcc->
num_entries)
continue;
872 q+= qscale[
index] * coeff;
875 blurred_qscale[i]= q/sum;
887 expected_bits +=
bits;
891 "expected_bits: %f all_available_bits: %d rate_factor: %f\n",
892 expected_bits, (
int)all_available_bits, rate_factor);
893 if(expected_bits > all_available_bits) {
904 av_dlog(s,
"[lavc rc] entry[%d].new_qscale = %.3f qp = %.3f\n",
910 assert(toobig <= 40);
912 "[lavc rc] requested bitrate: %d bps expected bitrate: %d bps\n",
914 (
int)(expected_bits / ((
double)all_available_bits/s->
bit_rate)));
916 "[lavc rc] estimated target average qp: %.3f\n",
920 "[lavc rc] Using all of requested bitrate is not "
921 "necessary for this video with these parameters.\n");
922 }
else if (toobig == 40) {
924 "[lavc rc] Error: bitrate too low for this video "
925 "with these parameters.\n");
927 }
else if (fabs(expected_bits/all_available_bits - 1.0) > 0.01) {
929 "[lavc rc] Error: 2pass curve failed to converge\n");