开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 多媒体 > Audio > 查看源码
report.c
资源名称：chapter15.rar　[点击查看]
上传用户：hjq518
上传日期：2021-12-09
资源大小：5084k
文件大小：58k
源码类别：

Audio

开发平台：

Visual C++

report.c：源码内容
  1. /*!
  2.  *  file
  3.  *     report.c
  4.  *  brief
  5.  *     Report related files()
  6.  *  author
  7.  *   Main contributors (see contributors.h for copyright, address and affiliation details)
  8.  *   - Karsten Suehring                <suehring@hhi.de>
  9.  *   - Alexis Michael Tourapis         <alexismt@ieee.org>
  10.  ***********************************************************************
  11.  */
  13. #include "contributors.h"
  15. #include <time.h>
  16. #include <math.h>
  17. #include <sys/timeb.h>
  18. #include "global.h"
  20. #include "context_ini.h"
  21. #include "explicit_gop.h"
  22. #include "filehandle.h"
  23. #include "fmo.h"
  24. #include "image.h"
  25. #include "intrarefresh.h"
  26. #include "leaky_bucket.h"
  27. #include "me_epzs.h"
  28. #include "output.h"
  29. #include "parset.h"
  30. #include "report.h"
  33. StatParameters   statistics,  *stats = &statistics;
  34. DistortionParams distortions, *dist  = &distortions;
  35. static char DistortionType[3][20] = {"SAD", "SSE", "Hadamard SAD"};
  37. extern int frame_statistic_start;
  38. extern ColocatedParams *Co_located;
  39. extern ColocatedParams *Co_located_JV[MAX_PLANE];  //!< Co_located to be used during 4:4:4 independent mode encoding
  40. extern void Clear_Motion_Search_Module (void);
  41. void   report_log_mode(InputParameters *params, StatParameters *stats, int64 bit_use[NUM_SLICE_TYPES][2]);
  42. /*!
  43.  ************************************************************************
  44.  * brief
  45.  *    Reports frame statistical data to a stats file
  46.  ************************************************************************
  47.  */
  48. void report_frame_statistic()
  49. {
  50.   FILE *p_stat_frm = NULL;
  51.   static int64 last_bit_ctr_n = 0;
  52.   int i;
  53.   char name[30];
  54.   int bitcounter;
  55.   StatParameters *cur_stats = &enc_picture->stats;
  57. #ifndef WIN32
  58.   time_t now;
  59.   struct tm *l_time;
  60.   char string[1000];
  61. #else
  62.   char timebuf[128];
  63. #endif
  65.   // write to log file
  66.   if ((p_stat_frm = fopen("stat_frame.dat", "r")) == 0)            // check if file exists
  67.   {
  68.     if ((p_stat_frm = fopen("stat_frame.dat", "a")) == NULL)       // append new statistic at the end
  69.     {
  70.       snprintf(errortext, ET_SIZE, "Error open file %s  n", "stat_frame.dat.dat");
  71.       error(errortext, 500);
  72.     }
  73.     else                                            // Create header for new log file
  74.     {
  75.       fprintf(p_stat_frm, " --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- n");
  76.       fprintf(p_stat_frm, "|            Encoder statistics. This file is generated during first encoding session, new sessions will be appended                                                                                                                                                                                                                                                                                                                                                              |n");
  77.       fprintf(p_stat_frm, " --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- n");
  78.     }
  79.   }
  80.   else
  81.   {
  82.     fclose (p_stat_frm);
  83.     if ((p_stat_frm = fopen("stat_frame.dat", "a")) == NULL)       // File exists, just open for appending
  84.     {
  85.       snprintf(errortext, ET_SIZE, "Error open file %s  n", "stat_frame.dat.dat");
  86.       error(errortext, 500);
  87.     }
  88.   }
  90.   if (frame_statistic_start)
  91.   {
  92.     fprintf(p_stat_frm, "|     ver     | Date  | Time  |    Sequence                  |Frm | QP |P/MbInt|   Bits   |  SNRY  |  SNRU  |  SNRV  |  I4  |  I8  | I16  | IC0  | IC1  | IC2  | IC3  | PI4  | PI8  | PI16 |  P0  |  P1  |  P2  |  P3  | P1*4*| P1*8*| P2*4*| P2*8*| P3*4*| P3*8*|  P8  | P8:4 | P4*4*| P4*8*| P8:5 | P8:6 | P8:7 | BI4  | BI8  | BI16 |  B0  |  B1  |  B2  |  B3  | B0*4*| B0*8*| B1*4*| B1*8*| B2*4*| B2*8*| B3*4*| B3*8*|  B8  | B8:0 |B80*4*|B80*8*| B8:4 | B4*4*| B4*8*| B8:5 | B8:6 | B8:7 |n");
  93.     fprintf(p_stat_frm, " ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ n");
  94.   }
  96.   //report
  97.   fprintf(p_stat_frm, "|%4s/%s", VERSION, EXT_VERSION);
  99. #ifdef WIN32
  100.   _strdate( timebuf );
  101.   fprintf(p_stat_frm, "| %1.5s |", timebuf);
  103.   _strtime( timebuf);
  104.   fprintf(p_stat_frm, " % 1.5s |", timebuf);
  105. #else
  106.   now = time ((time_t *) NULL); // Get the system time and put it into 'now' as 'calender time'
  107.   time (&now);
  108.   l_time = localtime (&now);
  109.   strftime (string, sizeof string, "%d-%b-%Y", l_time);
  110.   fprintf(p_stat_frm, "| %1.5s |", string );
  112.   strftime (string, sizeof string, "%H:%M:%S", l_time);
  113.   fprintf(p_stat_frm, " %1.5s |", string);
  114. #endif
  116.   for (i=0;i<30;i++)
  117.     name[i]=params->infile[i + imax(0,(int) (strlen(params->infile)- 30))]; // write last part of path, max 30 chars
  119.   fprintf(p_stat_frm, "%30.30s|", name);
  120.   fprintf(p_stat_frm, "%3d |", frame_no);
  121.   fprintf(p_stat_frm, "%3d |", img->qp);
  122.   fprintf(p_stat_frm, "  %d/%d  |", params->PicInterlace, params->MbInterlace);
  125.   if (img->frm_number == 0 && img->frame_num == 0)
  126.   {
  127.     bitcounter = (int) stats->bit_counter[I_SLICE];
  128.   }
  129.   else
  130.   {
  131.     bitcounter = (int) (stats->bit_ctr_n - last_bit_ctr_n);
  132.     last_bit_ctr_n = stats->bit_ctr_n;
  133.   }
  135.   //report bitrate
  136.   fprintf(p_stat_frm, " %9d|", bitcounter);
  138.   //report snr's  
  139.   fprintf(p_stat_frm, " %2.4f| %2.4f| %2.4f|", dist->metric[PSNR].value[0], dist->metric[PSNR].value[1], dist->metric[PSNR].value[2]);
  141.   //report modes
  142.   //I-Modes
  143.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[I_SLICE][I4MB ]);
  144.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[I_SLICE][I8MB ]);
  145.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[I_SLICE][I16MB]);
  147.   //chroma intra mode
  148.   fprintf(p_stat_frm, " %5d|", cur_stats->intra_chroma_mode[0]);
  149.   fprintf(p_stat_frm, " %5d|", cur_stats->intra_chroma_mode[1]);
  150.   fprintf(p_stat_frm, " %5d|", cur_stats->intra_chroma_mode[2]);
  151.   fprintf(p_stat_frm, " %5d|", cur_stats->intra_chroma_mode[3]);
  153.   //P-Modes
  154.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[P_SLICE][I4MB ]);
  155.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[P_SLICE][I8MB ]);
  156.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[P_SLICE][I16MB]);
  157.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[P_SLICE][0    ]);
  159.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[P_SLICE][1    ]);
  160.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[P_SLICE][2    ]);
  161.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[P_SLICE][3    ]);
  162.   
  163.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][1][0]);
  164.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][1][1]);
  165.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][2][0]);
  166.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][2][1]);
  167.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][3][0]);
  168.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][3][1]);
  170.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[P_SLICE][P8x8 ]);
  171.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[P_SLICE][4    ]);
  172.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use_transform[P_SLICE][4][0]);
  173.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use_transform[P_SLICE][4][1]);
  174.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[P_SLICE][5    ]);
  175.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[P_SLICE][6    ]);
  176.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[P_SLICE][7    ]);
  178.   //B-Modes
  179.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[B_SLICE][I4MB ]);
  180.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[B_SLICE][I8MB ]);
  181.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[B_SLICE][I16MB]);
  182.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[B_SLICE][0    ]);
  183.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[B_SLICE][1    ]);
  184.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[B_SLICE][2    ]);
  185.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[B_SLICE][3    ]);
  186.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use_transform[B_SLICE][0][0]);
  187.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use_transform[B_SLICE][0][1]);
  188.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use_transform[B_SLICE][1][0]);
  189.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use_transform[B_SLICE][1][1]);
  190.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use_transform[B_SLICE][2][0]);
  191.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use_transform[B_SLICE][2][1]);
  192.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use_transform[B_SLICE][3][0]);
  193.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use_transform[B_SLICE][3][1]);
  195.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[B_SLICE][P8x8]);
  196.   fprintf(p_stat_frm, " %d|", (cur_stats->b8_mode_0_use [B_SLICE][0] + cur_stats->b8_mode_0_use [B_SLICE][1]));
  197.   fprintf(p_stat_frm, " %5d|", cur_stats->b8_mode_0_use [B_SLICE][0]);
  198.   fprintf(p_stat_frm, " %5d|", cur_stats->b8_mode_0_use [B_SLICE][1]);
  199.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[B_SLICE][4   ]);
  200.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use_transform[B_SLICE][4][0]);
  201.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use_transform[B_SLICE][4][1]);
  202.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[B_SLICE][5   ]);
  203.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[B_SLICE][6   ]);
  204.   fprintf(p_stat_frm, " %5" FORMAT_OFF_T  "|", cur_stats->mode_use[B_SLICE][7   ]);
  206.   fprintf(p_stat_frm, "n");
  208.   //save the last results
  209.   frame_statistic_start = 0;
  210.   fclose(p_stat_frm);
  211. }
  214. double report_slice_pred_stats(FILE *p_stat, StatParameters *stats, int slice_type, double bit_use, char *slice_name)
  215. {
  216.   fprintf(p_stat,"n ---------------------|----------------|-----------------|");
  217.   fprintf(p_stat,"n   %8s           |   Mode used    | MotionInfo bits |", slice_name);
  218.   fprintf(p_stat,"n ---------------------|----------------|-----------------|");
  219.   fprintf(p_stat,"n Mode  0  (copy)      |  %5" FORMAT_OFF_T  "         |    %8.2f     |", stats->mode_use[slice_type][0   ], (double)stats->bit_use_mode[slice_type][0   ] / bit_use);
  220.   fprintf(p_stat,"n Mode  1  (16x16)     |  %5" FORMAT_OFF_T  "         |    %8.2f     |", stats->mode_use[slice_type][1   ], (double)stats->bit_use_mode[slice_type][1   ] / bit_use);
  221.   fprintf(p_stat,"n Mode  2  (16x8)      |  %5" FORMAT_OFF_T  "         |    %8.2f     |", stats->mode_use[slice_type][2   ], (double)stats->bit_use_mode[slice_type][2   ] / bit_use);
  222.   fprintf(p_stat,"n Mode  3  (8x16)      |  %5" FORMAT_OFF_T  "         |    %8.2f     |", stats->mode_use[slice_type][3   ], (double)stats->bit_use_mode[slice_type][3   ] / bit_use);
  223.   fprintf(p_stat,"n Mode  4  (8x8)       |  %5" FORMAT_OFF_T  "         |    %8.2f     |", stats->mode_use[slice_type][P8x8], (double)stats->bit_use_mode[slice_type][P8x8] / bit_use);
  224.   fprintf(p_stat,"n Mode  5  intra 4x4   |  %5" FORMAT_OFF_T  "         |-----------------|", stats->mode_use[slice_type][I4MB]);
  225.   fprintf(p_stat,"n Mode  6  intra 8x8   |  %5" FORMAT_OFF_T  "         |", stats->mode_use[slice_type][I8MB]);
  226.   fprintf(p_stat,"n Mode  7+ intra 16x16 |  %5" FORMAT_OFF_T  "         |", stats->mode_use[slice_type][I16MB]);
  227.   fprintf(p_stat,"n Mode     intra IPCM  |  %5" FORMAT_OFF_T  "         |", stats->mode_use[slice_type][IPCM ]);
  229.   return (double)(stats->bit_use_mode[slice_type][0] + stats->bit_use_mode[slice_type][1] + stats->bit_use_mode[slice_type][2]
  230.   + stats->bit_use_mode[slice_type][3] + stats->bit_use_mode[slice_type][P8x8]) / bit_use;
  231. }
  233. /*!
  234.  ***********************************************************************
  235.  * brief
  236.  *    Terminates and reports statistics on error.
  237.  *
  238.  ***********************************************************************
  239.  */
  240. void report_stats_on_error(void)
  241. {
  242.   params->no_frames = img->frm_number;
  243.   free_encoder_memory(img);
  244.   exit (-1);
  245. }
  248. void report_stats(InputParameters *params, StatParameters *stats, int64 bit_use[NUM_SLICE_TYPES][2], float frame_rate)
  249. {
  250.   FILE *p_stat;                    //!< status file for the last encoding session
  251.   double mean_motion_info_bit_use[NUM_SLICE_TYPES] = {0.0};
  252.   int i;
  254.   if (strlen(params->StatsFile) == 0)
  255.     strcpy (params->StatsFile,"stats.dat");
  257.   if ((p_stat = fopen(params->StatsFile, "wt")) == 0)
  258.   {
  259.     snprintf(errortext, ET_SIZE, "Error open file %s", params->StatsFile);
  260.     error(errortext, 500);
  261.   }
  263.   fprintf(p_stat," -------------------------------------------------------------- n");
  264.   fprintf(p_stat,"  This file contains statistics for the last encoded sequence   n");
  265.   fprintf(p_stat," -------------------------------------------------------------- n");
  266.   fprintf(p_stat,   " Sequence                     : %sn", params->infile);
  267.   fprintf(p_stat,   " No.of coded pictures         : %4dn", stats->frame_counter);
  268.   fprintf(p_stat,   " Freq. for encoded bitstream  : %4.0fn", frame_rate);
  270.   fprintf(p_stat,   " I Slice Bitrate(kb/s)        : %6.2fn", stats->bitrate_st[I_SLICE] / 1000);
  271.   fprintf(p_stat,   " P Slice Bitrate(kb/s)        : %6.2fn", stats->bitrate_st[P_SLICE] / 1000);
  272.   fprintf(p_stat,   " B Slice Bitrate(kb/s)        : %6.2fn", stats->bitrate_st[B_SLICE] / 1000);
  273.   fprintf(p_stat,   " Total Bitrate(kb/s)          : %6.2fn", stats->bitrate / 1000);
  275.   for (i = 0; i < 3; i++)
  276.   {
  277.     fprintf(p_stat," ME Level %1d Metric            : %sn", i, DistortionType[params->MEErrorMetric[i]]);
  278.   }
  279.   fprintf(p_stat," Mode Decision Metric         : %sn", DistortionType[params->ModeDecisionMetric]);
  281.   switch ( params->ChromaMEEnable )
  282.   {
  283.   case 1:
  284.     fprintf(p_stat," ME for components            : YCbCrn");
  285.     break;
  286.   default:
  287.     fprintf(p_stat," ME for components            : Yn");
  288.     break;
  289.   }
  291.   fprintf(p_stat,  " Image format                 : %dx%dn", params->output.width, params->output.height);
  293.   if (params->intra_upd)
  294.     fprintf(p_stat," Error robustness             : Onn");
  295.   else
  296.     fprintf(p_stat," Error robustness             : Offn");
  298.   fprintf(p_stat,  " Search range                 : %dn", params->search_range);
  300.   fprintf(p_stat,   " Total number of references   : %dn", params->num_ref_frames);
  301.   fprintf(p_stat,   " References for P slices      : %dn", params->P_List0_refs ? params->P_List0_refs : params->num_ref_frames);
  303.   if (stats->frame_ctr[B_SLICE]!=0)
  304.   {
  305.     fprintf(p_stat, " List0 refs for B slices      : %dn", params->B_List0_refs ? params->B_List0_refs : params->num_ref_frames);
  306.     fprintf(p_stat, " List1 refs for B slices      : %dn", params->B_List1_refs ? params->B_List1_refs : params->num_ref_frames);
  307.   }
  309.   fprintf(p_stat,   " Profile/Level IDC            : (%d,%d)n", params->ProfileIDC, params->LevelIDC);
  310.   if (params->symbol_mode == CAVLC)
  311.     fprintf(p_stat,   " Entropy coding method        : CAVLCn");
  312.   else
  313.     fprintf(p_stat,   " Entropy coding method        : CABACn");
  315.   if (params->MbInterlace)
  316.     fprintf(p_stat, " MB Field Coding : On n");
  318.   if (params->SearchMode == EPZS)
  319.     EPZSOutputStats(params, p_stat, 1);
  321.   if (params->full_search == 2)
  322.     fprintf(p_stat," Search range restrictions    : nonen");
  323.   else if (params->full_search == 1)
  324.     fprintf(p_stat," Search range restrictions    : older reference framesn");
  325.   else
  326.     fprintf(p_stat," Search range restrictions    : smaller blocks and older reference framesn");
  328.   if (params->rdopt)
  329.     fprintf(p_stat," RD-optimized mode decision   : usedn");
  330.   else
  331.     fprintf(p_stat," RD-optimized mode decision   : not usedn");
  333.   fprintf(p_stat,"n ---------------------|----------------|---------------|");
  334.   fprintf(p_stat,"n     Item             |     Intra      |   All frames  |");
  335.   fprintf(p_stat,"n ---------------------|----------------|---------------|");
  336.   fprintf(p_stat,"n SNR Y(dB)            |");
  337.   fprintf(p_stat," %5.2f          |", dist->metric[PSNR].avslice[I_SLICE][0]);
  338.   fprintf(p_stat," %5.2f         |", dist->metric[PSNR].average[0]);
  339.   fprintf(p_stat,"n SNR U/V (dB)         |");
  340.   fprintf(p_stat," %5.2f/%5.2f    |", dist->metric[PSNR].avslice[I_SLICE][1], dist->metric[PSNR].avslice[I_SLICE][2]);
  341.   fprintf(p_stat," %5.2f/%5.2f   |", dist->metric[PSNR].average[1], dist->metric[PSNR].average[2]);
  342.   fprintf(p_stat,"n ---------------------|----------------|---------------|");
  343.   fprintf(p_stat,"n");
  345.   fprintf(p_stat,"n ---------------------|----------------|---------------|---------------|");
  346.   fprintf(p_stat,"n     SNR              |        I       |       P       |       B       |");
  347.   fprintf(p_stat,"n ---------------------|----------------|---------------|---------------|");
  348.   fprintf(p_stat,"n SNR Y(dB)            |      %5.3f    |     %5.3f    |     %5.3f    |",
  349.     dist->metric[PSNR].avslice[I_SLICE][0], dist->metric[PSNR].avslice[P_SLICE][0], dist->metric[PSNR].avslice[B_SLICE][0]);
  350.   fprintf(p_stat,"n SNR U(dB)            |      %5.3f    |     %5.3f    |     %5.3f    |",
  351.     dist->metric[PSNR].avslice[I_SLICE][1], dist->metric[PSNR].avslice[P_SLICE][1], dist->metric[PSNR].avslice[B_SLICE][1]);
  352.   fprintf(p_stat,"n SNR V(dB)            |      %5.3f    |     %5.3f    |     %5.3f    |",
  353.     dist->metric[PSNR].avslice[I_SLICE][2], dist->metric[PSNR].avslice[P_SLICE][2], dist->metric[PSNR].avslice[B_SLICE][2]);
  354.   fprintf(p_stat,"n ---------------------|----------------|---------------|---------------|");  
  355.   fprintf(p_stat,"n");
  357.   // QUANT.
  358.   fprintf(p_stat,"n ---------------------|----------------|---------------|---------------|");
  359.   fprintf(p_stat,"n     Ave Quant        |        I       |       P       |       B       |");
  360.   fprintf(p_stat,"n ---------------------|----------------|---------------|---------------|");
  361.   fprintf(p_stat,"n        QP            |      %5.3f    |     %5.3f    |     %5.3f    |",
  362.     (float)stats->quant[I_SLICE]/dmax(1.0,(float)stats->num_macroblocks[I_SLICE]),
  363.     (float)stats->quant[P_SLICE]/dmax(1.0,(float)stats->num_macroblocks[P_SLICE]),
  364.     (float)stats->quant[B_SLICE]/dmax(1.0,(float)stats->num_macroblocks[B_SLICE]));
  365.   fprintf(p_stat,"n ---------------------|----------------|---------------|---------------|");  
  366.   fprintf(p_stat,"n");
  368.   // MODE
  369.   fprintf(p_stat,"n ---------------------|----------------|");
  370.   fprintf(p_stat,"n   Intra              |   Mode used    |");
  371.   fprintf(p_stat,"n ---------------------|----------------|");
  372.   fprintf(p_stat,"n Mode 0  intra 4x4    |  %5" FORMAT_OFF_T  "         |", stats->mode_use[I_SLICE][I4MB ]);
  373.   fprintf(p_stat,"n Mode 1  intra 8x8    |  %5" FORMAT_OFF_T  "         |", stats->mode_use[I_SLICE][I8MB ]);
  374.   fprintf(p_stat,"n Mode 2+ intra 16x16  |  %5" FORMAT_OFF_T  "         |", stats->mode_use[I_SLICE][I16MB]);
  375.   fprintf(p_stat,"n Mode    intra IPCM   |  %5" FORMAT_OFF_T  "         |", stats->mode_use[I_SLICE][IPCM ]);
  377.   // P slices
  378.   if (stats->frame_ctr[P_SLICE]!=0)
  379.   {    
  380.     mean_motion_info_bit_use[P_SLICE] = report_slice_pred_stats(p_stat, stats, P_SLICE,(double) bit_use[P_SLICE][0], "P Slice ");
  381.   }
  382.   // B slices
  383.   if (stats->frame_ctr[B_SLICE]!=0)
  384.   {
  385.     mean_motion_info_bit_use[B_SLICE] = report_slice_pred_stats(p_stat, stats, B_SLICE,(double) bit_use[B_SLICE][0], "B Slice ");
  386.   }
  387.   // SP slices
  388.   if (stats->frame_ctr[SP_SLICE]!=0)
  389.   {    
  390.     mean_motion_info_bit_use[SP_SLICE] = report_slice_pred_stats(p_stat, stats, SP_SLICE,(double) bit_use[SP_SLICE][0], "SP Slice");
  391.   }
  394.   fprintf(p_stat,"n ---------------------|----------------|");
  395.   fprintf(p_stat,"n");
  397.   fprintf(p_stat,"n ---------------------|----------------|----------------|----------------|");
  398.   fprintf(p_stat,"n  Bit usage:          |      Intra     |      Inter     |    B frame     |");
  399.   fprintf(p_stat,"n ---------------------|----------------|----------------|----------------|");
  401.   fprintf(p_stat,"n Header               |");
  402.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_header[I_SLICE] / bit_use[I_SLICE][0]);
  403.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_header[P_SLICE] / bit_use[P_SLICE][0]);
  404.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_header[B_SLICE] / bit_use[B_SLICE][0]);
  406.   fprintf(p_stat,"n Mode                 |");
  407.   fprintf(p_stat," %10.2f     |", (float)stats->bit_use_mb_type[I_SLICE] / bit_use[I_SLICE][0]);
  408.   fprintf(p_stat," %10.2f     |", (float)stats->bit_use_mb_type[P_SLICE] / bit_use[P_SLICE][0]);
  409.   fprintf(p_stat," %10.2f     |", (float)stats->bit_use_mb_type[B_SLICE] / bit_use[B_SLICE][0]);
  411.   fprintf(p_stat,"n Motion Info          |");
  412.   fprintf(p_stat,"        ./.     |");
  413.   fprintf(p_stat," %10.2f     |", mean_motion_info_bit_use[P_SLICE]);
  414.   fprintf(p_stat," %10.2f     |", mean_motion_info_bit_use[B_SLICE]);
  416.   fprintf(p_stat,"n CBP Y/C              |");
  417.   fprintf(p_stat," %10.2f     |", (float) stats->tmp_bit_use_cbp[I_SLICE] / bit_use[I_SLICE][0]);
  418.   fprintf(p_stat," %10.2f     |", (float) stats->tmp_bit_use_cbp[P_SLICE] / bit_use[P_SLICE][0]);   
  419.   fprintf(p_stat," %10.2f     |", (float) stats->tmp_bit_use_cbp[B_SLICE] / bit_use[B_SLICE][0]);
  421.   fprintf(p_stat,"n Coeffs. Y            |");
  422.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_coeff[0][I_SLICE] / bit_use[I_SLICE][0]);
  423.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_coeff[0][P_SLICE] / bit_use[P_SLICE][0]);   
  424.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_coeff[0][B_SLICE] / bit_use[B_SLICE][0]);
  426.   fprintf(p_stat,"n Coeffs. C            |");
  427.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_coeffC[I_SLICE] / bit_use[I_SLICE][0]);
  428.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_coeffC[P_SLICE] / bit_use[P_SLICE][0]);   
  429.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_coeffC[B_SLICE] / bit_use[B_SLICE][0]);
  431.   fprintf(p_stat,"n Coeffs. CB           |");
  432.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_coeff[1][I_SLICE] / bit_use[I_SLICE][0]);
  433.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_coeff[1][P_SLICE] / bit_use[P_SLICE][0]);   
  434.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_coeff[1][B_SLICE] / bit_use[B_SLICE][0]);
  436.   fprintf(p_stat,"n Coeffs. CB           |");
  437.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_coeff[2][I_SLICE] / bit_use[I_SLICE][0]);
  438.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_coeff[2][P_SLICE] / bit_use[P_SLICE][0]);   
  439.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_coeff[2][B_SLICE] / bit_use[B_SLICE][0]);
  441.   fprintf(p_stat,"n Delta quant          |");
  442.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_delta_quant[I_SLICE] / bit_use[I_SLICE][0]);
  443.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_delta_quant[P_SLICE] / bit_use[P_SLICE][0]);   
  444.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_delta_quant[B_SLICE] / bit_use[B_SLICE][0]);
  446.   fprintf(p_stat,"n Stuffing Bits        |");
  447.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_stuffingBits[I_SLICE] / bit_use[I_SLICE][0]);
  448.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_stuffingBits[P_SLICE] / bit_use[P_SLICE][0]);   
  449.   fprintf(p_stat," %10.2f     |", (float) stats->bit_use_stuffingBits[B_SLICE] / bit_use[B_SLICE][0]);
  451.   fprintf(p_stat,"n ---------------------|----------------|----------------|----------------|");
  452.   fprintf(p_stat,"n average bits/frame   |");
  453.   fprintf(p_stat," %10.2f     |", (float) bit_use[I_SLICE][1] / (float) bit_use[I_SLICE][0] );
  454.   fprintf(p_stat," %10.2f     |", (float) bit_use[P_SLICE][1] / (float) bit_use[P_SLICE][0] );
  455.   fprintf(p_stat," %10.2f     |", (float) bit_use[B_SLICE][1] / (float) bit_use[B_SLICE][0] );
  456.   fprintf(p_stat,"n ---------------------|----------------|----------------|----------------|");
  457.   fprintf(p_stat,"n");
  459.   fclose(p_stat);
  460. }
  463. void report_log(InputParameters *params, StatParameters *stats, float frame_rate)
  464. {
  465.   char name[40];
  466.   int i;
  467. #ifndef WIN32
  468.   time_t now;
  469.   struct tm *l_time;
  470.   char string[1000];
  471. #else
  472.   char timebuf[128];
  473. #endif
  475.   if ((p_log = fopen("log.dat", "r")) == 0)         // check if file exists
  476.   {
  477.     if ((p_log = fopen("log.dat", "a")) == NULL)    // append new statistic at the end
  478.     {
  479.       snprintf(errortext, ET_SIZE, "Error open file %s  n", "log.dat");
  480.       error(errortext, 500);
  481.     }
  482.     else                                            // Create header for new log file
  483.     {
  484.       fprintf(p_log," ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ n");
  485.       fprintf(p_log,"|                          Encoder statistics. This file is generated during first encoding session, new sessions will be appended                                                                                                                                                                                 |n");
  486.       fprintf(p_log," ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ n");
  487.       fprintf(p_log,"|     ver     | Date  | Time  |               Sequence                 | #Img |P/MbInt| QPI| QPP| QPB| Format  |Iperiod| #B | FMES | Hdmd | S.R |#Ref | Freq |Coding|RD-opt|Intra upd|8x8Tr| SNRY 1| SNRU 1| SNRV 1| SNRY N| SNRU N| SNRV N|#Bitr I|#Bitr P|#Bitr B|#Bitr IPB|     Total Time   |      Me Time     |n");
  488.       fprintf(p_log," ------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------ n");
  490.     }
  491.   }
  492.   else
  493.   {
  494.     fclose (p_log);
  495.     if ((p_log = fopen("log.dat", "a")) == NULL)         // File exists, just open for appending
  496.     {
  497.       snprintf(errortext, ET_SIZE, "Error open file %s  n", "log.dat");
  498.       error(errortext, 500);
  499.     }
  500.   }
  501.   fprintf(p_log,"|%5s/%-5s", VERSION, EXT_VERSION);
  503. #ifdef WIN32
  504.   _strdate( timebuf );
  505.   fprintf(p_log,"| %1.5s |", timebuf );
  507.   _strtime( timebuf);
  508.   fprintf(p_log," % 1.5s |", timebuf);
  509. #else
  510.   now = time ((time_t *) NULL); // Get the system time and put it into 'now' as 'calender time'
  511.   time (&now);
  512.   l_time = localtime (&now);
  513.   strftime (string, sizeof string, "%d-%b-%Y", l_time);
  514.   fprintf(p_log,"| %1.5s |", string );
  516.   strftime (string, sizeof string, "%H:%M:%S", l_time);
  517.   fprintf(p_log," %1.5s |", string );
  518. #endif
  520.   for (i=0; i < 40; i++)
  521.     name[i] = params->infile[i + imax(0, ((int) strlen(params->infile)) - 40)]; // write last part of path, max 40 chars
  522.   fprintf(p_log,"%40.40s|",name);
  524.   fprintf(p_log,"%5d |  %d/%d  |", params->no_frames, params->PicInterlace, params->MbInterlace);
  525.   fprintf(p_log," %-3d| %-3d| %-3d|", params->qp0, params->qpN, params->qpB);
  527.   fprintf(p_log,"%4dx%-4d|", params->output.width, params->output.height);
  528.   fprintf(p_log,"  %3d  |%3d |", params->intra_period, stats->successive_Bframe);
  531.   switch( params->SearchMode ) 
  532.   {
  533.   case UM_HEX:
  534.     fprintf(p_log,"  HEX |");
  535.     break;
  536.   case UM_HEX_SIMPLE:
  537.     fprintf(p_log," SHEX |");
  538.     break;
  539.   case EPZS:
  540.     fprintf(p_log," EPZS |");
  541.     break;
  542.   case FAST_FULL_SEARCH:
  543.     fprintf(p_log,"  FFS |");
  544.     break;
  545.   default:
  546.     fprintf(p_log,"  FS  |");
  547.     break;
  548.   }
  550.   fprintf(p_log,"  %1d%1d%1d |", params->MEErrorMetric[F_PEL], params->MEErrorMetric[H_PEL], params->MEErrorMetric[Q_PEL]);
  552.   fprintf(p_log," %3d | %2d  |", params->search_range, params->num_ref_frames );
  554.   fprintf(p_log," %5.2f|", (img->framerate *(float) (stats->successive_Bframe + 1)) / (float)(params->jumpd + 1));
  556.   if (params->symbol_mode == CAVLC)
  557.     fprintf(p_log," CAVLC|");
  558.   else
  559.     fprintf(p_log," CABAC|");
  561.   fprintf(p_log,"   %d  |", params->rdopt);
  563.   if (params->intra_upd == 1)
  564.     fprintf(p_log,"   ON    |");
  565.   else
  566.     fprintf(p_log,"   OFF   |");
  568.   fprintf(p_log,"  %d  |", params->Transform8x8Mode);
  570.   fprintf(p_log,"%7.3f|%7.3f|%7.3f|", 
  571.     dist->metric[PSNR].avslice[I_SLICE][0],
  572.     dist->metric[PSNR].avslice[I_SLICE][1],
  573.     dist->metric[PSNR].avslice[I_SLICE][2]);
  574.   fprintf(p_log,"%7.3f|%7.3f|%7.3f|", dist->metric[PSNR].average[0],dist->metric[PSNR].average[1],dist->metric[PSNR].average[2]);
  575.   /*
  576.   fprintf(p_log,"%-5.3f|%-5.3f|%-5.3f|", dist->metric[PSNR].avslice[I_SLICE][0], dist->metric[PSNR].avslice[I_SLICE][1], dist->metric[PSNR].avslice[I_SLICE][2]);
  577.   fprintf(p_log,"%-5.3f|%-5.3f|%-5.3f|", dist->metric[PSNR].avslice[P_SLICE][0], dist->metric[PSNR].avslice[P_SLICE][1], dist->metric[PSNR].avslice[P_SLICE][2]);
  578.   fprintf(p_log,"%-5.3f|%-5.3f|%-5.3f|", dist->metric[PSNR].avslice[B_SLICE][0], dist->metric[PSNR].avslice[B_SLICE][1], dist->metric[PSNR].avslice[B_SLICE][2]);
  579.   */
  580.   fprintf(p_log,"%7.0f|%7.0f|%7.0f|%9.0f|", stats->bitrate_st[I_SLICE],stats->bitrate_st[P_SLICE],stats->bitrate_st[B_SLICE], stats->bitrate);
  582.   fprintf(p_log,"   %12d   |   %12d   |", (int)tot_time,(int)me_tot_time);
  585.   fprintf(p_log,"n");
  587.   fclose(p_log);
  589.   p_log = fopen("data.txt", "a");
  591.   if ((stats->successive_Bframe != 0) && (stats->frame_ctr[B_SLICE] != 0)) // B picture used
  592.   {
  593.     fprintf(p_log, "%3d %2d %2d %2.2f %2.2f %2.2f %5" FORMAT_OFF_T  " "
  594.       "%2.2f %2.2f %2.2f %5d "
  595.       "%2.2f %2.2f %2.2f %5" FORMAT_OFF_T  " %5" FORMAT_OFF_T  " %.3fn",
  596.       params->no_frames, params->qp0, params->qpN,
  597.       dist->metric[PSNR].avslice[I_SLICE][0],
  598.       dist->metric[PSNR].avslice[I_SLICE][1],
  599.       dist->metric[PSNR].avslice[I_SLICE][2],
  600.       stats->bit_counter[I_SLICE],
  601.       0.0,
  602.       0.0,
  603.       0.0,
  604.       0,
  605.       dist->metric[PSNR].average[0],
  606.       dist->metric[PSNR].average[1],
  607.       dist->metric[PSNR].average[2],
  608.       (stats->bit_counter[I_SLICE] + stats->bit_ctr) / stats->frame_counter,
  609.       stats->bit_counter[B_SLICE] / stats->frame_ctr[B_SLICE],
  610.       (double) 0.001 * tot_time / (stats->frame_counter));
  611.   }
  612.   else
  613.   {
  614.     if (params->no_frames != 0)
  615.       fprintf(p_log, "%3d %2d %2d %2.2f %2.2f %2.2f %5" FORMAT_OFF_T  " "
  616.       "%2.2f %2.2f %2.2f %5d "
  617.       "%2.2f %2.2f %2.2f %5" FORMAT_OFF_T  " %5d %.3fn",
  618.       params->no_frames, params->qp0, params->qpN,
  619.       dist->metric[PSNR].avslice[I_SLICE][0],
  620.       dist->metric[PSNR].avslice[I_SLICE][1],
  621.       dist->metric[PSNR].avslice[I_SLICE][2],
  622.       stats->bit_counter[I_SLICE],
  623.       0.0,
  624.       0.0,
  625.       0.0,
  626.       0,
  627.       dist->metric[PSNR].average[0],
  628.       dist->metric[PSNR].average[1],
  629.       dist->metric[PSNR].average[2],
  630.       (stats->bit_counter[I_SLICE] + stats->bit_ctr)/ stats->frame_counter,
  631.       0,
  632.       (double)0.001*tot_time/params->no_frames);
  633.   }
  635.   fclose(p_log);
  636. }
  638. /*!
  639.  ************************************************************************
  640.  * brief
  641.  *    Reports the gathered information to appropriate outputs
  642.  * par Input:
  643.  *    struct inp_par *inp,                                            n
  644.  *    ImageParameters *img,                                            n
  645.  *    struct stat_par *stats,                                          n
  646.  *
  647.  * par Output:
  648.  *    None
  649.  ************************************************************************
  650.  */
  651. void report( ImageParameters *img, InputParameters *params, StatParameters *stats)
  652. {
  653.   int64 bit_use[NUM_SLICE_TYPES][2];
  654.   int i,j;
  655.   int64 total_bits;
  656.   float frame_rate;
  658.   bit_use[ I_SLICE][0] = stats->frame_ctr[I_SLICE];
  659.   bit_use[ P_SLICE][0] = imax(stats->frame_ctr[P_SLICE ] + stats->frame_ctr[SP_SLICE], 1);
  660.   bit_use[ B_SLICE][0] = imax(stats->frame_ctr[B_SLICE ], 1);
  661.   bit_use[SP_SLICE][0] = imax(stats->frame_ctr[SP_SLICE], 1);
  663.   //  Accumulate bit usage for inter and intra frames
  664.   for (j=0; j < NUM_SLICE_TYPES; j++)
  665.   {
  666.     bit_use[j][1] = 0;
  667.   }
  669.   for (j=0; j < NUM_SLICE_TYPES; j++)
  670.   {
  671.     for(i=0; i < MAXMODE; i++)
  672.       bit_use[j][1] += stats->bit_use_mode[j][i];
  674.     bit_use[j][1] += stats->bit_use_mb_type[j];
  675.     bit_use[j][1] += stats->bit_use_header[j];    
  676.     bit_use[j][1] += stats->tmp_bit_use_cbp[j];
  677.     bit_use[j][1] += stats->bit_use_coeffC[j];
  678.     bit_use[j][1] += stats->bit_use_coeff[0][j];   
  679.     bit_use[j][1] += stats->bit_use_coeff[1][j]; 
  680.     bit_use[j][1] += stats->bit_use_coeff[2][j]; 
  681.     bit_use[j][1] += stats->bit_use_delta_quant[j];
  682.     bit_use[j][1] += stats->bit_use_stuffingBits[j];
  683.   }
  685.   frame_rate = (img->framerate *(float)(stats->successive_Bframe + 1)) / (float) (params->jumpd + 1);
  687.   //! Currently adding NVB bits on P rate. Maybe additional stats info should be created instead and added in log file  
  688.   stats->bitrate_st[ I_SLICE] = (stats->bit_counter[ I_SLICE]) * (frame_rate) / (float) (stats->frame_counter);
  689.   stats->bitrate_st[ P_SLICE] = (stats->bit_counter[ P_SLICE]) * (frame_rate) / (float) (stats->frame_counter);
  690.   stats->bitrate_st[ B_SLICE] = (stats->bit_counter[ B_SLICE]) * (frame_rate) / (float) (stats->frame_counter);
  691.   stats->bitrate_st[SP_SLICE] = (stats->bit_counter[SP_SLICE]) * (frame_rate) / (float) (stats->frame_counter);
  693.   switch (params->Verbose)
  694.   {
  695.   case 0:
  696.   case 1:
  697.   default:
  698.     fprintf(stdout,"------------------ Average data all frames  -----------------------------------nn");
  699.     break;
  700.   case 2:
  701.     fprintf(stdout,"------------------------------------  Average data all frames  ---------------------------------nn");
  702.     break;
  703.   }
  705.   if (params->Verbose != 0)
  706.   {
  707.     DistMetric *snr     = &dist->metric[PSNR    ];
  708.     DistMetric *sse     = &dist->metric[SSE     ];
  709.     DistMetric *snr_rgb = &dist->metric[PSNR_RGB];
  710.     DistMetric *sse_rgb = &dist->metric[SSE_RGB ];
  712.     int  impix    = params->output.size_cmp[0];
  713.     int  impix_cr = params->output.size_cmp[1];
  715.     float csnr_y = psnr(img->max_imgpel_value_comp_sq[0], impix   , sse->average[0]);
  716.     float csnr_u = psnr(img->max_imgpel_value_comp_sq[1], impix_cr, sse->average[1]);
  717.     float csnr_v = psnr(img->max_imgpel_value_comp_sq[2], impix_cr, sse->average[2]);
  719.     fprintf(stdout,  " Total encoding time for the seq.  : %.3f sec (%.2f fps)n", tot_time*0.001, 1000.0 * (stats->frame_counter) / tot_time);
  720.     fprintf(stdout,  " Total ME time for sequence        : %.3f sec nn", me_tot_time*0.001);
  722.     fprintf(stdout," Y { PSNR (dB), cSNR (dB), MSE }   : { %5.2f, %5.2f, %5.2f }n", 
  723.       snr->average[0], csnr_y, sse->average[0]/impix);
  724.     fprintf(stdout," U { PSNR (dB), cSNR (dB), MSE }   : { %5.2f, %5.2f, %5.2f }n",
  725.       snr->average[1], csnr_u, sse->average[1]/impix_cr);
  726.     fprintf(stdout," V { PSNR (dB), cSNR (dB), MSE }   : { %5.2f, %5.2f, %5.2f }n",
  727.       snr->average[2], csnr_v, sse->average[2]/impix_cr);
  729.     if(params->DistortionYUVtoRGB == 1)
  730.     {
  731.       float csnr_r = psnr(img->max_imgpel_value_comp_sq[0], impix, sse_rgb->average[0]);
  732.       float csnr_g = psnr(img->max_imgpel_value_comp_sq[1], impix, sse_rgb->average[1]);
  733.       float csnr_b = psnr(img->max_imgpel_value_comp_sq[2], impix, sse_rgb->average[2]);
  735.       fprintf(stdout," R { PSNR (dB), cSNR (dB), MSE }   : { %5.2f, %5.2f, %5.2f }n", 
  736.         snr_rgb->average[0], csnr_r, sse_rgb->average[0]/impix);
  737.       fprintf(stdout," G { PSNR (dB), cSNR (dB), MSE }   : { %5.2f, %5.2f, %5.2f }n",
  738.         snr_rgb->average[1], csnr_g, sse_rgb->average[1]/impix);
  739.       fprintf(stdout," B { PSNR (dB), cSNR (dB), MSE }   : { %5.2f, %5.2f, %5.2f }n",
  740.         snr_rgb->average[2], csnr_b, sse_rgb->average[2]/impix);
  741.     }
  743.     if (params->Distortion[SSIM] == 1)
  744.     {
  745.       if(params->DistortionYUVtoRGB == 1)
  746.       {
  747.         fprintf(stdout," SSIM {Y, R}                       : { %5.4f, %5.4f }n", dist->metric[SSIM].average[0], dist->metric[SSIM_RGB].average[0]);
  748.         fprintf(stdout," SSIM {U, G}                       : { %5.4f, %5.4f }n", dist->metric[SSIM].average[1], dist->metric[SSIM_RGB].average[1]);
  749.         fprintf(stdout," SSIM {V, B}                       : { %5.4f, %5.4f }n", dist->metric[SSIM].average[2], dist->metric[SSIM_RGB].average[2]);
  750.       }
  751.       else
  752.       {
  753.         fprintf(stdout," Y SSIM                            : %5.4fn", dist->metric[SSIM].average[0]);
  754.         fprintf(stdout," U SSIM                            : %5.4fn", dist->metric[SSIM].average[1]);
  755.         fprintf(stdout," V SSIM                            : %5.4fn", dist->metric[SSIM].average[2]);
  756.       }
  757.     }
  758.     if (params->Distortion[MS_SSIM] == 1)
  759.     {
  760.       if(params->DistortionYUVtoRGB == 1)
  761.       {
  762.         fprintf(stdout," MS-SSIM {Y, R}                    : { %5.4f, %5.4f }n", dist->metric[MS_SSIM].average[0], dist->metric[MS_SSIM_RGB].average[0]);
  763.         fprintf(stdout," MS-SSIM {U, G}                    : { %5.4f, %5.4f }n", dist->metric[MS_SSIM].average[1], dist->metric[MS_SSIM_RGB].average[1]);
  764.         fprintf(stdout," MS-SSIM {V, B}                    : { %5.4f, %5.4f }n", dist->metric[MS_SSIM].average[2], dist->metric[MS_SSIM_RGB].average[2]);
  765.       }
  766.       else
  767.       {
  768.         fprintf(stdout," Y MS-SSIM                         : %5.4fn", dist->metric[MS_SSIM].average[0]);
  769.         fprintf(stdout," U MS-SSIM                         : %5.4fn", dist->metric[MS_SSIM].average[1]);
  770.         fprintf(stdout," V MS-SSIM                         : %5.4fn", dist->metric[MS_SSIM].average[2]);
  771.       }
  772.     }
  773.     fprintf(stdout,"n");
  774.   }
  775.   else
  776.     fprintf(stdout,  " Total encoding time for the seq.  : %.3f sec (%.2f fps)nn", tot_time*0.001, 1000.0 * (stats->frame_counter) / tot_time);
  778.   total_bits = stats->bit_ctr_parametersets;
  779.   for (i = 0; i < NUM_SLICE_TYPES; i++)
  780.     total_bits += stats->bit_counter[i];
  782.   if (stats->frame_ctr[B_SLICE] != 0)
  783.   {
  784.     fprintf(stdout, " Total bits                        : %" FORMAT_OFF_T  " (I %" FORMAT_OFF_T  ", P %" FORMAT_OFF_T  ", B %" FORMAT_OFF_T  " NVB %d) n",
  785.       total_bits,  stats->bit_counter[I_SLICE], stats->bit_counter[P_SLICE], stats->bit_counter[B_SLICE], stats->bit_ctr_parametersets);
  786.   }
  787.   else if (params->sp_periodicity == 0)
  788.   {
  789.     fprintf(stdout, " Total bits                        : %" FORMAT_OFF_T  " (I %" FORMAT_OFF_T  ", P %" FORMAT_OFF_T  ", NVB %d) n",
  790.       total_bits, stats->bit_counter[I_SLICE], stats->bit_counter[P_SLICE], stats->bit_ctr_parametersets);
  791.   }
  792.   else
  793.   {
  794.     fprintf(stdout, " Total bits                        : %" FORMAT_OFF_T  " (I %" FORMAT_OFF_T  ", P %" FORMAT_OFF_T  ", NVB %d) n",
  795.       total_bits, stats->bit_counter[I_SLICE], stats->bit_counter[P_SLICE], stats->bit_ctr_parametersets);
  796.   }
  798.   frame_rate = (img->framerate *(float)(stats->successive_Bframe + 1)) / (float) (params->jumpd + 1);
  800.   stats->bitrate= ((float) total_bits * frame_rate) / ((float) (stats->frame_counter));
  801.   fprintf(stdout, " Bit rate (kbit/s)  @ %2.2f Hz     : %5.2fn", frame_rate, stats->bitrate / 1000);
  802.   
  803.   for (i = 0; i < 5; i++)
  804.     stats->bit_ctr_emulationprevention += stats->bit_use_stuffingBits[i];
  806.   fprintf(stdout, " Bits to avoid Startcode Emulation : %" FORMAT_OFF_T  " n", stats->bit_ctr_emulationprevention);
  807.   fprintf(stdout, " Bits for parameter sets           : %d nn", stats->bit_ctr_parametersets);
  809.   switch (params->Verbose)
  810.   {
  811.   case 0:
  812.   case 1:
  813.   default:
  814.     fprintf(stdout,"-------------------------------------------------------------------------------n");
  815.     break;
  816.   case 2:
  817.     fprintf(stdout,"------------------------------------------------------------------------------------------------n");
  818.     break;
  819.   }  
  820.   fprintf(stdout,"Exit JM %s encoder ver %s ", JM, VERSION);
  821.   fprintf(stdout,"n");
  823.   // status file
  824.   report_stats(params, stats, bit_use, frame_rate);
  826.   // write to log file
  827.   report_log(params, stats, frame_rate);
  829.   if (params->ReportFrameStats)
  830.   {
  831.     if ((p_log = fopen("stat_frame.dat", "a")) == NULL)       // append new statistic at the end
  832.     {
  833.       snprintf(errortext, ET_SIZE, "Error open file %s  n", "stat_frame.dat.dat");
  834.       //    error(errortext, 500);
  835.     }
  836.     else
  837.     {
  838.       fprintf(p_log," --------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- n");
  839.       fclose(p_log);
  840.     }
  841.     report_log_mode(params, stats, bit_use);
  842.   }
  843. }
  846. /*!
  847.  ************************************************************************
  848.  * brief
  849.  *    Prints the header of the protocol.
  850.  * par Input:
  851.  *    struct inp_par *inp
  852.  * par Output:
  853.  *    none
  854.  ************************************************************************
  855.  */
  856. void information_init ( ImageParameters *img, InputParameters *params, StatParameters *stats)
  857. {
  858.   int i;
  859.   static char yuv_types[4][10] = {"YUV 4:0:0", "YUV 4:2:0", "YUV 4:2:2", "YUV 4:4:4"};
  860.   switch (params->Verbose)
  861.   {
  862.   case 0:
  863.   case 1:
  864.   default:
  865.     printf("------------------------------- JM %4.4s %7.7s -------------------------------n", VERSION, EXT_VERSION);
  866.     break;
  867.   case 2:
  868.     printf("--------------------------------------- JM %4.4s %7.7s ----------------------------------------n", VERSION, EXT_VERSION);
  869.     break;
  870.   }
  872.   fprintf(stdout,  " Input YUV file                    : %s n", params->infile);
  873.   fprintf(stdout,  " Output H.264 bitstream            : %s n", params->outfile);
  874.   if (p_dec != -1)
  875.     fprintf(stdout,  " Output YUV file                   : %s n", params->ReconFile);
  876.   fprintf(stdout,  " YUV Format                        : %s n", &yuv_types[img->yuv_format][0]);//img->yuv_format==YUV422?"YUV 4:2:2":(img->yuv_format==YUV444)?"YUV 4:4:4":"YUV 4:2:0");
  877.   fprintf(stdout,  " Frames to be encoded I-P/B        : %d/%dn", params->no_frames, (params->successive_Bframe*(params->no_frames-1)));
  878.   if (params->Verbose != 0)
  879.   {
  880.     fprintf(stdout,  " Freq. for encoded bitstream       : %1.0fn", img->framerate/(float)(params->jumpd+1));
  881.     fprintf(stdout,  " PicInterlace / MbInterlace        : %d/%dn", params->PicInterlace, params->MbInterlace);
  882.     fprintf(stdout,  " Transform8x8Mode                  : %dn", params->Transform8x8Mode);
  884.     for (i=0; i<3; i++)
  885.     {
  886.       fprintf(stdout," ME Metric for Refinement Level %1d  : %sn", i, DistortionType[params->MEErrorMetric[i]]);
  887.     }
  888.     fprintf(stdout,  " Mode Decision Metric              : %sn", DistortionType[params->ModeDecisionMetric]);
  890.     switch ( params->ChromaMEEnable )
  891.     {
  892.     case 1:
  893.       fprintf(stdout," Motion Estimation for components  : YCbCrn");
  894.       break;
  895.     default:
  896.       fprintf(stdout," Motion Estimation for components  : Yn");
  897.       break;
  898.     }
  900.     fprintf(stdout,  " Image format                      : %dx%d (%dx%d)n", params->output.width, params->output.height, img->width,img->height);
  902.     if (params->intra_upd)
  903.       fprintf(stdout," Error robustness                  : Onn");
  904.     else
  905.       fprintf(stdout," Error robustness                  : Offn");
  906.     fprintf(stdout,  " Search range                      : %dn", params->search_range);
  908.     fprintf(stdout,  " Total number of references        : %dn", params->num_ref_frames);
  909.     fprintf(stdout,  " References for P slices           : %dn", params->P_List0_refs ? params->P_List0_refs : params->num_ref_frames);
  910.     fprintf(stdout,  " List0 references for B slices     : %dn", params->B_List0_refs ? params->B_List0_refs : params->num_ref_frames);
  911.     fprintf(stdout,  " List1 references for B slices     : %dn", params->B_List1_refs ? params->B_List1_refs : params->num_ref_frames);
  913.     // Sequence Type
  914.     fprintf(stdout,  " Sequence type                     :");
  915.     if (stats->successive_Bframe > 0 && params->HierarchicalCoding)
  916.     {
  917.       fprintf(stdout, " Hierarchy (QP: I %d, P %d, B %d) n",
  918.         params->qp0, params->qpN, params->qpB);
  919.     }
  920.     else if (stats->successive_Bframe > 0)
  921.     {
  922.       char seqtype[80];
  923.       int i,j;
  925.       strcpy (seqtype,"I");
  927.       for (j=0; j < 2; j++)
  928.       {
  929.         for (i=0; i < stats->successive_Bframe; i++)
  930.         {
  931.           if (params->BRefPictures)
  932.             strncat(seqtype,"-RB", imax(0, (int) (79 - strlen(seqtype))));
  933.           else
  934.             strncat(seqtype,"-B", imax(0, (int) (79 - strlen(seqtype))));
  935.         }
  936.         strncat(seqtype,"-P", imax(0, (int) (79 - strlen(seqtype))));
  937.       }
  938.       if (params->BRefPictures)
  939.         fprintf(stdout, " %s (QP: I %d, P %d, RB %d) n", seqtype, params->qp0, params->qpN, iClip3(0, 51, params->qpB + params->qpBRSOffset));
  940.       else
  941.         fprintf(stdout, " %s (QP: I %d, P %d, B %d) n", seqtype, params->qp0, params->qpN, params->qpB);
  942.     }
  943.     else if (stats->successive_Bframe == 0 && params->sp_periodicity == 0) 
  944.       fprintf(stdout, " IPPP (QP: I %d, P %d) n", params->qp0, params->qpN);
  945.     else 
  946.       fprintf(stdout, " I-P-P-SP-P (QP: I %d, P %d, SP (%d, %d)) n",  params->qp0, params->qpN, params->qpsp, params->qpsp_pred);
  948.     // report on entropy coding  method
  949.     if (params->symbol_mode == CAVLC)
  950.       fprintf(stdout," Entropy coding method             : CAVLCn");
  951.     else
  952.       fprintf(stdout," Entropy coding method             : CABACn");
  954.     fprintf(stdout,  " Profile/Level IDC                 : (%d,%d)n", params->ProfileIDC, params->LevelIDC);
  956.     if (params->SearchMode == UM_HEX)
  957.       fprintf(stdout,  " Motion Estimation Scheme          : HEXn");
  958.     else if (params->SearchMode == UM_HEX_SIMPLE)
  959.       fprintf(stdout,  " Motion Estimation Scheme          : SHEXn");
  960.     else if (params->SearchMode == EPZS)
  961.     {
  962.       fprintf(stdout,  " Motion Estimation Scheme          : EPZSn");
  963.       EPZSOutputStats(params, stdout, 0);
  964.     }
  965.     else if (params->SearchMode == FAST_FULL_SEARCH)
  966.       fprintf(stdout,  " Motion Estimation Scheme          : Fast Full Searchn");
  967.     else
  968.       fprintf(stdout,  " Motion Estimation Scheme          : Full Searchn");
  970.     if (params->full_search == 2)
  971.       fprintf(stdout," Search range restrictions         : nonen");
  972.     else if (params->full_search == 1)
  973.       fprintf(stdout," Search range restrictions         : older reference framesn");
  974.     else
  975.       fprintf(stdout," Search range restrictions         : smaller blocks and older reference framesn");
  977.     if (params->rdopt)
  978.       fprintf(stdout," RD-optimized mode decision        : usedn");
  979.     else
  980.       fprintf(stdout," RD-optimized mode decision        : not usedn");
  982.     switch(params->partition_mode)
  983.     {
  984.     case PAR_DP_1:
  985.       fprintf(stdout," Data Partitioning Mode            : 1 partition n");
  986.       break;
  987.     case PAR_DP_3:
  988.       fprintf(stdout," Data Partitioning Mode            : 3 partitions n");
  989.       break;
  990.     default:
  991.       fprintf(stdout," Data Partitioning Mode            : not supportedn");
  992.       break;
  993.     }
  995.     switch(params->of_mode)
  996.     {
  997.     case PAR_OF_ANNEXB:
  998.       fprintf(stdout," Output File Format                : H.264/AVC Annex B Byte Stream Format n");
  999.       break;
  1000.     case PAR_OF_RTP:
  1001.       fprintf(stdout," Output File Format                : RTP Packet File Format n");
  1002.       break;
  1003.     default:
  1004.       fprintf(stdout," Output File Format                : not supportedn");
  1005.       break;
  1006.     }
  1007.   }
  1010.   switch (params->Verbose)
  1011.   {
  1012.   case 0:
  1013.   default:
  1014.     printf("-------------------------------------------------------------------------------n");
  1015.     printf("nEncoding. Please Wait.nn");
  1016.     break;    
  1017.   case 1:
  1018.     printf("-------------------------------------------------------------------------------n");
  1019.     printf("  Frame  Bit/pic    QP   SnrY    SnrU    SnrV    Time(ms) MET(ms) Frm/Fld Ref  n");
  1020.     printf("-------------------------------------------------------------------------------n");
  1021.     break;
  1022.   case 2:
  1023.     if (params->Distortion[SSIM] == 1)
  1024.     {
  1025.       printf("------------------------------------------------------------------------------------------------------------------------n");
  1026.       printf("  Frame  Bit/pic WP QP QL   SnrY    SnrU    SnrV   SsimY   SsimU   SsimV    Time(ms) MET(ms) Frm/Fld   I D L0 L1 RDP Refn");
  1027.       printf("------------------------------------------------------------------------------------------------------------------------n");
  1028.     }
  1029.     else
  1030.     {
  1031.       printf("------------------------------------------------------------------------------------------------n");
  1032.       printf("  Frame  Bit/pic WP QP QL   SnrY    SnrU    SnrV    Time(ms) MET(ms) Frm/Fld   I D L0 L1 RDP Refn");
  1033.       printf("------------------------------------------------------------------------------------------------n");
  1034.     }
  1035.     break;
  1036.    
  1037.   }
  1038. }
  1040. /*!
  1041.  ************************************************************************
  1042.  * brief
  1043.  *    Report mode distribution of the sequence to log_mode.dat
  1044.  ************************************************************************
  1045.  */
  1046. void report_log_mode(InputParameters *params, StatParameters *stats, int64 bit_use[NUM_SLICE_TYPES][2])
  1047. {
  1048.   FILE *p_stat;                    //!< status file for the last encoding session
  1049.   int i;
  1050.   char name[40];
  1051. #ifndef WIN32
  1052.   time_t now;
  1053.   struct tm *l_time;
  1054.   char string[1000];
  1055. #else
  1056.   char timebuf[128];
  1057. #endif
  1059.   if ((p_stat = fopen("log_mode.dat", "r")) == 0)         // check if file exists
  1060.   {
  1061.     if ((p_stat = fopen("log_mode.dat", "a")) == NULL)    // append new statistic at the end
  1062.     {
  1063.       snprintf(errortext, ET_SIZE, "Error open file %s  n", "log_mode.dat");
  1064.       error(errortext, 500);
  1065.     }
  1066.     else                                            // Create header for new log file
  1067.     {
  1068.       fprintf(p_stat, " ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- n");
  1069.       fprintf(p_stat,"|                          Encoder statistics. This file is generated during first encoding session, new sessions will be appended                                                                                                                                                                                 |n");
  1070.       fprintf(p_stat, " ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- n");
  1071.       fprintf(p_stat, "|     ver     | Date  | Time  |    Sequence                  | QP |  I4  |  I8  | I16  | IC0  | IC1  | IC2  | IC3  | PI4  | PI8  | PI16 |  P0  |  P1  |  P2  |  P3  | P1*4*| P1*8*| P2*4*| P2*8*| P3*4*| P3*8*|  P8  | P8:4 | P4*4*| P4*8*| P8:5 | P8:6 | P8:7 | BI4  | BI8  | BI16 |  B0  |  B1  |  B2  |  B3  | B0*4*| B0*8*| B1*4*| B1*8*| B2*4*| B2*8*| B3*4*| B3*8*|  B8  | B8:0 |B80*4*|B80*8*| B8:4 | B4*4*| B4*8*| B8:5 | B8:6 | B8:7 |n");
  1072.       fprintf(p_stat, " ---------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------- n");
  1073.     }
  1074.   }
  1075.   else
  1076.   {
  1077.     fclose (p_stat);
  1078.     if ((p_stat = fopen("log_mode.dat", "a")) == NULL)         // File exists, just open for appending
  1079.     {
  1080.       snprintf(errortext, ET_SIZE, "Error open file %s  n", "log_mode.dat");
  1081.       error(errortext, 500);
  1082.     }
  1083.   }
  1085.   //report
  1086.   fprintf(p_stat, "|%4s/%s", VERSION, EXT_VERSION);
  1088. #ifdef WIN32
  1089.   _strdate( timebuf );
  1090.   fprintf(p_stat, "| %1.5s |", timebuf);
  1092.   _strtime( timebuf);
  1093.   fprintf(p_stat, " % 1.5s |", timebuf);
  1094. #else
  1095.   now = time ((time_t *) NULL); // Get the system time and put it into 'now' as 'calender time'
  1096.   time (&now);
  1097.   l_time = localtime (&now);
  1098.   strftime (string, sizeof string, "%d-%b-%Y", l_time);
  1099.   fprintf(p_stat, "| %1.5s |", string );
  1101.   strftime (string, sizeof string, "%H:%M:%S", l_time);
  1102.   fprintf(p_stat, " %1.5s |", string);
  1103. #endif
  1105.   for (i=0;i<30;i++)
  1106.     name[i]=params->infile[i + imax(0,(int) (strlen(params->infile)- 30))]; // write last part of path, max 30 chars
  1108.   fprintf(p_stat, "%30.30s|", name);
  1109.   fprintf(p_stat, "%3d |", img->qp);
  1111.   //report modes
  1112.   //I-Modes
  1113.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[I_SLICE][I4MB ]);
  1114.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[I_SLICE][I8MB ]);
  1115.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[I_SLICE][I16MB]);
  1117.   //chroma intra mode
  1118.   fprintf(p_stat, " %5d|", stats->intra_chroma_mode[0]);
  1119.   fprintf(p_stat, " %5d|", stats->intra_chroma_mode[1]);
  1120.   fprintf(p_stat, " %5d|", stats->intra_chroma_mode[2]);
  1121.   fprintf(p_stat, " %5d|", stats->intra_chroma_mode[3]);
  1123.   //P-Modes
  1124.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[P_SLICE][I4MB ]);
  1125.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[P_SLICE][I8MB ]);
  1126.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[P_SLICE][I16MB]);
  1127.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[P_SLICE][0    ]);
  1129.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[P_SLICE][1    ]);
  1130.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[P_SLICE][2    ]);
  1131.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[P_SLICE][3    ]);
  1132.   
  1133.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][1][0]);
  1134.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][1][1]);
  1135.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][2][0]);
  1136.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][2][1]);
  1137.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][3][0]);
  1138.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][3][1]);
  1140.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[P_SLICE][P8x8 ]);
  1141.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[P_SLICE][4    ]);
  1142.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][4][0]);
  1143.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[P_SLICE][4][1]);
  1144.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[P_SLICE][5    ]);
  1145.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[P_SLICE][6    ]);
  1146.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[P_SLICE][7    ]);
  1148.   //B-Modes
  1149.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[B_SLICE][I4MB ]);
  1150.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[B_SLICE][I8MB ]);
  1151.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[B_SLICE][I16MB]);
  1152.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[B_SLICE][0    ]);
  1153.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[B_SLICE][1    ]);
  1154.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[B_SLICE][2    ]);
  1155.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[B_SLICE][3    ]);
  1156.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[B_SLICE][0][0]);
  1157.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[B_SLICE][0][1]);
  1158.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[B_SLICE][1][0]);
  1159.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[B_SLICE][1][1]);
  1160.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[B_SLICE][2][0]);
  1161.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[B_SLICE][2][1]);
  1162.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[B_SLICE][3][0]);
  1163.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[B_SLICE][3][1]);
  1165.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[B_SLICE][P8x8]);
  1166.   fprintf(p_stat, " %d|", (stats->b8_mode_0_use [B_SLICE][0]+stats->b8_mode_0_use [B_SLICE][1]));
  1167.   fprintf(p_stat, " %5d|", stats->b8_mode_0_use [B_SLICE][0]);
  1168.   fprintf(p_stat, " %5d|", stats->b8_mode_0_use [B_SLICE][1]);
  1169.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[B_SLICE][4   ]);
  1170.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[B_SLICE][4][0]);
  1171.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use_transform[B_SLICE][4][1]);
  1172.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[B_SLICE][5   ]);
  1173.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[B_SLICE][6   ]);
  1174.   fprintf(p_stat, " %5" FORMAT_OFF_T  "|", stats->mode_use[B_SLICE][7   ]);
  1176.   fprintf(p_stat, "n");
  1177.   fclose(p_stat);
  1178. }