tce3d.c

spiht for linux this is used to decod and encode vedio i wich all enjoy

                  else
                    {
                      if (return_value)
                        {
                          QccErrorAddMessage("(QccWAVtceIPBand): QccWAVtceIPBand()");
                          goto Error;
                        }
                    } 
                  
                  if (col == 0)
                    p1[row][col] =
                      alpha * QCCWAVTCE3D_BOUNDARY_VALUE + filter_coef * v;
                  else
                    p1[row][col] = alpha * p1[row][col-1] + filter_coef * v;
                  p2[col] = p1[row][col] + alpha * p2[col];
                  
                }//a
              else
                { 
                  if (col == 0)
                    p1[row][col] =
                      alpha * QCCWAVTCE3D_BOUNDARY_VALUE  +
                      filter_coef * QCCWAVTCE3D_REFINE_HOLDER;
                  else
                    p1[row][col] =
                      alpha* p1[row][col-1] +
                      filter_coef * QCCWAVTCE3D_REFINE_HOLDER;
                  p2[col] = p1[row][col] + alpha * p2[col];
                  
                }
              p_estimation
                [subband_origin_frame+frame]
                [subband_origin_row + row]
                [subband_origin_col + col] =
                p_forward;
            }//aa
          
          
          for (col = subband_num_cols - 1; col >= 0; col--)
            {  
              if (col ==( subband_num_cols - 1)) 
                p2[col] = p2[col] + alpha * QCCWAVTCE3D_BOUNDARY_VALUE;
              else
                p2[col] = p2[col] + alpha * p4[row][col+1];
              p3[row][col]=p2[col]+alpha*p3[row][col];
              p_char = &(significance_map
                         [current_frame]
                         [current_row]
                         [subband_origin_col + col]);
              if (*p_char == QCCWAVTCE3D_S)
                {
                  if (col == (subband_num_cols - 1))
                    p4[row][col] =
                      alpha * QCCWAVTCE3D_BOUNDARY_VALUE +
                      filter_coef * QCCWAVTCE3D_REFINE_HOLDER; 
                  else 
                    p4[row][col] =
                      alpha * p4[row][col+1] +
                      filter_coef * QCCWAVTCE3D_REFINE_HOLDER;;
                } 
              else
                {
                  v = (*p_char == QCCWAVTCE3D_S_NEW);
                  if (col == subband_num_cols - 1)
                    p4[row][col] =
                      alpha * QCCWAVTCE3D_BOUNDARY_VALUE + filter_coef * v;
                  else
                    p4[row][col] = alpha * p4[row][col+1] + filter_coef * v;
                }
            }
        }//aaa
      //	printf("frame=%d\n",frame);
      for (col = 0; col < subband_num_cols; col++)
        p5[col]=0;
      for (row = subband_num_rows - 1; row >= 0; row--)
        {
          for (col = 0; col < subband_num_cols; col++)
            {
              p3[row][col]=p3[row][col]+alpha*p5[col];
              if (col ==( subband_num_cols - 1))
                p5[col] =
                  p1[row][col] + alpha*p5[col]+alpha*QCCWAVTCE3D_BOUNDARY_VALUE ;
              else 
                p5[col] = p1[row][col] + alpha*p5[col] + alpha*p4[row][col+1];
            }
        }
    }//aaaa
  
  for (frame = subband_num_frames-1;frame >= 0; frame--)
    {//bbb
      current_frame=subband_origin_frame+frame;

      for (row = subband_num_rows - 1; row >= 0; row--)
        {//bb
          
          current_row = subband_origin_row + row;
          for (col = subband_num_cols-1; col >= 0; col--)
            {   
              current_col = subband_origin_col + col; 
              p_char = &(significance_map
                         [current_frame][current_row][current_col]);
              if (col ==( subband_num_cols-1) )
                p_estimation[current_frame][current_row][current_col] +=
                  alpha * (QCCWAVTCE3D_BOUNDARY_VALUE + p2[col] + p3[row][col]);
              else
                p_estimation[current_frame][current_row][current_col] +=
                  alpha *(p1[row][col+1]+p2[col]+p3[row][col]);
              
              if (*p_char == QCCWAVTCE3D_S)
                {
                  p_estimation[current_frame][current_row][current_col]=
                    QCCWAVTCE3D_REFINE_HOLDER * weight[0] +
                    p_estimation[current_frame][current_row][current_col] *
                    weight[1];
                  if (col == (subband_num_cols - 1))
                    p1[row][col] = alpha * QCCWAVTCE3D_BOUNDARY_VALUE +
                      filter_coef * QCCWAVTCE3D_REFINE_HOLDER;
                  else
                    p1[row][col] = alpha * p1[row][col+1] +
                      filter_coef * QCCWAVTCE3D_REFINE_HOLDER;
                }
              else
                {
                  v = (*p_char == QCCWAVTCE3D_S_NEW);
                  p_estimation[current_frame][current_row][current_col] =
                    v * weight[0] +
                    p_estimation[current_frame][current_row][current_col] *
                    weight[1];
                  if (col == (subband_num_cols-1) )
                    p1[row][col] = alpha *QCCWAVTCE3D_BOUNDARY_VALUE +
                      filter_coef * v;
                  else
                    p1[row][col] = alpha * p1[row][col+1] + filter_coef * v;
                }    
              p2[col] = p1[row][col] + alpha * p2[col];   
            }
          
          for (col = 0; col < subband_num_cols; col++)
            {
	      if (col==0)
                p2[col] = p2[col] + alpha * QCCWAVTCE3D_BOUNDARY_VALUE;
	      else
                p2[col] = p2[col] + alpha * p4[row][col-1];
	      p3[row][col]=p2[col]+alpha*p3[row][col];
              current_col = subband_origin_col + col; 
              p_char = &(significance_map
                         [current_frame][current_row][current_col]);
              if (*p_char == QCCWAVTCE3D_S)
                {
	          if(col==0)
                    p4[row][col] =
                      alpha * QCCWAVTCE3D_BOUNDARY_VALUE +
                      filter_coef * QCCWAVTCE3D_REFINE_HOLDER;
                  else 
                    p4[row][col] =
                      alpha * p4[row][col-1] +
                      filter_coef * QCCWAVTCE3D_REFINE_HOLDER;
                }
              else
                {
                  v = (*p_char == QCCWAVTCE3D_S_NEW);
                  if(col==0)
                    p4[row][col] =
                      alpha * QCCWAVTCE3D_BOUNDARY_VALUE + filter_coef * v;
                  else 
                    p4[row][col] = alpha * p4[row][col - 1] + filter_coef *v;
                }        
            }
        }//bb
      for (row = 0; row < subband_num_rows; row++)
        for (col = 0; col < subband_num_cols; col++)
          {
            p3[row][col] = p3[row][col] + alpha*p5[col];
            if (col==0)
              p5[col] = p1[row][col] + alpha*p5[col] +
                alpha* QCCWAVTCE3D_BOUNDARY_VALUE; 
            else
              p5[col] = p1[row][col] + alpha*p5[col] + alpha*p4[row][col-1];
          }
    }//bbb
  
  return_value = 0;
  goto Return;
 Error:
  return_value = 1;
 Return: 
  QccVectorFree(p2); 
  QccVectorFree(p5);
  QccMatrixFree(p1, subband_num_rows); 
  QccMatrixFree(p3, subband_num_rows); 
  QccMatrixFree(p4, subband_num_rows); 
  return(return_value); 
}


static int QccWAVtce3DIPPass(QccWAVSubbandPyramid3D *coefficients,    
                             char ***significance_map,
                             char ***sign_array,
                             double threshold,
                             double ***p_estimation,
                             double *subband_significance,
                             QccENTArithmeticModel *model,
                             QccBitBuffer *buffer,
                             int *max_coefficient_bits,
                             double alpha)
{
  int subband,num_subbands; 
  int return_value;
  
  if (coefficients->transform_type ==
      QCCWAVSUBBANDPYRAMID3D_DYADIC)
    num_subbands =
      QccWAVSubbandPyramid3DNumLevelsToNumSubbandsDyadic(coefficients->spatial_num_levels);
  else
    num_subbands =
      QccWAVSubbandPyramid3DNumLevelsToNumSubbandsPacket(coefficients->temporal_num_levels,
                                                         coefficients->spatial_num_levels);
  
  for (subband = 0; subband < num_subbands; subband++)
    {
      if ((threshold - 0.000001) <
          pow((double)2, (double)max_coefficient_bits[subband]))
        {
	  // printf("IP=%d\n",subband);
          return_value =
            QccWAVtce3DRevEst(coefficients,
                              significance_map,
                              subband,
                              subband_significance,
                              p_estimation,
                              alpha);
          if (return_value == 1)
            {
              QccErrorAddMessage("(QccWAVtce3DIPPass): Error calling QccWAVtce3DRevEst()");
              return(1);
            }
          else
            {
              if (return_value == 2)
                return(2);
            }
          
          return_value = QccWAVtce3DIPBand(coefficients,
                                           significance_map,
                                           sign_array,
                                           p_estimation,
                                           subband_significance,
                                           subband,
                                           threshold,
                                           model,
                                           buffer,
                                           alpha);
          
          if (return_value == 1)
            {
              QccErrorAddMessage("(QccWAVtce3DIPPass): Error calling QccWAVtce3DIPBand()");
              return(1);
            }
          else
            {
              if (return_value == 2)
                return(2);
            }
        }
    } 
  
  return(0);
}


int static QccWAVtce3DSPass(QccWAVSubbandPyramid3D *coefficients,    
                            char ***significance_map,
                            double threshold,
                            QccENTArithmeticModel *model,
                            QccBitBuffer *buffer,
                            int *max_coefficient_bits)
{
  int subband,num_subbands; 
  int return_value;
  int subband_origin_frame;
  int subband_origin_row;
  int subband_origin_col;
  int subband_num_frames;
  int subband_num_rows;
  int subband_num_cols;
  int frame,row, col;
  int current_frame, current_row, current_col;
  char *p_char;
  int symbol;
  double p = 0.4;
  
  if (coefficients->transform_type ==
      QCCWAVSUBBANDPYRAMID3D_DYADIC)
    num_subbands =
      QccWAVSubbandPyramid3DNumLevelsToNumSubbandsDyadic(coefficients->spatial_num_levels);
  else
    num_subbands =
      QccWAVSubbandPyramid3DNumLevelsToNumSubbandsPacket(coefficients->temporal_num_levels,
                                                         coefficients->spatial_num_levels);
  
  for (subband = 0; subband < num_subbands; subband++)
    {
      if ((threshold - 0.000001) <
          pow((double)2, (double)(max_coefficient_bits[subband])))
        { 
          if (QccWAVSubbandPyramid3DSubbandSize(coefficients,
                                                subband,
                                                &subband_num_frames,
                                                &subband_num_rows,
                                                &subband_num_cols))
            {
              QccErrorAddMessage("(QccWAVtce3DSPass): Error calling QccWAVSubbandPyramid3DSubbandSize()");
              return(1);
            }
          if (QccWAVSubbandPyramid3DSubbandOffsets(coefficients,
                                                   subband,
                                                   &subband_origin_frame,
                                                   &subband_origin_row,
                                                   &subband_origin_col))
            {
              QccErrorAddMessage("(QccWAVtce3DSPass): Error calling QccWAVSubbandPyramid3DSubbandOffsets()");
              return(1);
            }  
          
          for (frame = 0; frame < subband_num_frames; frame++) 
            {
              current_frame = subband_origin_frame + frame;
              for (row = 0; row < subband_num_rows; row++)
                {
                  current_row = subband_origin_row + row;
                  for (col = 0; col < subband_num_cols; col++)
                    { 
                      current_col = subband_origin_col + col;
                      p_char = &(significance_map
                                 [current_frame][current_row][current_col]);
                      
                      if (*p_char == QCCWAVTCE3D_S)
                        { 
                          
                          if (QccWAVtce3DUpdateModel(model, p))
                            {
                              QccErrorAddMessage("(QccWAVtce3DSPass): Error calling QccWAVTce3DUpdateModel()");
                              return(1);
                            }
                          
                          if (buffer->type == QCCBITBUFFER_OUTPUT)
                            {