sadwt_int.c

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

      
      current_length =
        QccWAVWaveletDWTSubbandLength(signal_length, scale + 1, 0, 0, 0);
      
      if (QccVectorIntCopy(temp_signal, signal, current_length))
        {
          QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveDWT1DInt): Error calling QccVectorIntCopy()");
          goto Error;
        }
      if (QccVectorIntCopy(temp_mask, mask, current_length))
        {
          QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveDWT1DInt): Error calling QccVectorIntCopy()");
          goto Error;
        }
    }
  
  return_value = 0;
  goto Return;
 Error:
  return_value = 1;
 Return:
  QccVectorIntFree(temp_signal);
  QccVectorIntFree(temp_mask);
  QccVectorIntFree(subsequence);
  return(return_value);
}


int QccWAVWaveletInverseShapeAdaptiveDWT1DInt(QccVectorInt signal,
                                              QccVectorInt mask,
                                              int signal_length,
                                              int num_scales,
                                              const QccWAVWavelet *wavelet)
{
  int return_value;
  int scale;
  QccVectorInt temp_signal = NULL;
  QccVectorInt temp_mask = NULL;
  QccVectorInt subsequence = NULL;
  int current_length;
  
  if (signal == NULL)
    return(0);
  if (mask == NULL)
    return(0);
  if (wavelet == NULL)
    return(0);
  
  if (num_scales <= 0)
    return(0);
  
  if ((temp_signal = QccVectorIntAlloc(signal_length)) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveInverseDWT1DInt): Error calling QccVectorIntAlloc()");
      goto Error;
    }
  if ((temp_mask = QccVectorIntAlloc(signal_length)) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveInverseDWT1DInt): Error calling QccVectorIntAlloc()");
      goto Error;
    }
  if ((subsequence = QccVectorIntAlloc(signal_length)) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveInverseDWT1DInt): Error calling QccVectorIntAlloc()");
      goto Error;
    }
  
  if (QccVectorIntCopy(temp_signal, signal, signal_length))
    {
      QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveInverseDWT1DInt): Error calling QccVectorIntCopy()");
      goto Error;
    }
  if (QccVectorIntCopy(temp_mask, mask, signal_length))
    {
      QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveInverseDWT1DInt): Error calling QccVectorIntCopy()");
      goto Error;
    }
  
  for (scale = num_scales - 1; scale >= 0 ; scale--)
    {
      current_length = 
        QccWAVWaveletDWTSubbandLength(signal_length, scale, 0, 0, 0);
      
      if (QccWAVWaveletShapeAdaptiveSynthesis1DInt(temp_signal,
                                                   signal,
                                                   subsequence,
                                                   temp_mask,
                                                   mask,
                                                   current_length,
                                                   wavelet))
        {
          QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveInverseDWT1DInt): Error calling QccWAVWaveletShapeAdaptiveSynthesis1DInt()");
          goto Error;
        }
      
      if (QccVectorIntCopy(temp_signal, signal, current_length))
        {
          QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveInverseDWT1DInt): Error calling QccVectorIntCopy()");
          goto Error;
        }
      if (QccVectorIntCopy(temp_mask, mask, current_length))
        {
          QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveInverseDWT1DInt): Error calling QccVectorIntCopy()");
          goto Error;
        }
    }
  
  return_value = 0;
  goto Return;
 Error:
  return_value = 1;
 Return:
  QccVectorIntFree(temp_signal);
  QccVectorIntFree(temp_mask);
  QccVectorIntFree(subsequence);
  return(return_value);
  
}


int QccWAVWaveletShapeAdaptiveDWT2DInt(QccMatrixInt signal, 
                                       QccMatrixInt mask, 
                                       int num_rows,
                                       int num_cols,
                                       int num_scales,
                                       const QccWAVWavelet *wavelet)
{
  int return_value;
  QccVectorInt input_signal_vector = NULL;
  QccVectorInt output_signal_vector = NULL;
  QccVectorInt subsequence = NULL;
  QccVectorInt input_mask_vector = NULL;
  QccVectorInt output_mask_vector = NULL;
  int scale;
  int row, col;
  int baseband_num_rows;
  int baseband_num_cols;
  
  /* Allocate temporary working arrays */
  if ((input_signal_vector =
       QccVectorIntAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDWT2DInt): Error calling VectorAlloc()");
      goto Error;
    }
  if ((output_signal_vector =
       QccVectorIntAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDWT2DInt): Error calling VectorAlloc()");
      goto Error;
    }
  if ((subsequence =
       QccVectorIntAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDWT2DInt): Error calling VectorAlloc()");
      goto Error;
    }
  if ((input_mask_vector =
       QccVectorIntAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDWT2DInt): Error calling VectorAlloc()");
      goto Error;
    }
  if ((output_mask_vector =
       QccVectorIntAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptiveDWT2DInt): Error calling VectorAlloc()");
      goto Error;
    }
  
  for (scale = 0; scale < num_scales; scale++)
    {
      baseband_num_rows =
        QccWAVWaveletDWTSubbandLength(num_rows, scale, 0, 0, 0);
      baseband_num_cols =
        QccWAVWaveletDWTSubbandLength(num_cols, scale, 0, 0, 0);
      
      for (row = 0; row < baseband_num_rows; row++)
        {
          for (col = 0; col < baseband_num_cols; col++)
            {
              input_signal_vector[col] = signal[row][col];
              input_mask_vector[col] = mask[row][col];
            }
          
          if (QccWAVWaveletShapeAdaptiveAnalysis1DInt(input_signal_vector,
                                                      output_signal_vector,
                                                      subsequence,
                                                      input_mask_vector,
                                                      output_mask_vector,
                                                      baseband_num_cols,
                                                      wavelet))
            {
              QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveDWT2DInt): Error calling QccWAVWaveletShapeAdaptiveAnalysis1DInt()");
              goto Error;
            }
          
          for (col = 0; col < baseband_num_cols; col++)
            {
              signal[row][col] = output_signal_vector[col];
              mask[row][col] = output_mask_vector[col];
            }
        }
      
      for (col = 0; col < baseband_num_cols; col++)
        {
          for (row = 0; row < baseband_num_rows; row++)
            {
              input_signal_vector[row] = signal[row][col];
              input_mask_vector[row] = mask[row][col];
            }
          
          if (QccWAVWaveletShapeAdaptiveAnalysis1DInt(input_signal_vector,
                                                      output_signal_vector,
                                                      subsequence,
                                                      input_mask_vector,
                                                      output_mask_vector,
                                                      baseband_num_rows,
                                                      wavelet))
            {
              QccErrorAddMessage("(QccWAVWaveletShapeAdaptiveDWT2DInt): Error calling QccWAVWaveletShapeAdaptiveAnalysis1DInt()");
              goto Error;
            }
          
          for (row = 0; row < baseband_num_rows; row++)
            {
              signal[row][col] = output_signal_vector[row];
              mask[row][col] = output_mask_vector[row];
            }
        }
    }
  
  return_value = 0;
  goto Return;
 Error:
  return_value = 1;
 Return:
  QccVectorIntFree(input_signal_vector);
  QccVectorIntFree(output_signal_vector);
  QccVectorIntFree(subsequence);
  QccVectorIntFree(input_mask_vector);
  QccVectorIntFree(output_mask_vector);
  return(return_value);
} 


int QccWAVWaveletInverseShapeAdaptiveDWT2DInt(QccMatrixInt signal, 
                                              QccMatrixInt mask, 
                                              int num_rows,
                                              int num_cols,
                                              int num_scales,
                                              const QccWAVWavelet *wavelet)
{
  int return_value;
  QccVectorInt input_signal_vector = NULL;
  QccVectorInt output_signal_vector = NULL;
  QccVectorInt subsequence = NULL;
  QccVectorInt input_mask_vector = NULL;
  QccVectorInt output_mask_vector = NULL;
  int scale;
  int row, col;
  int baseband_num_rows;
  int baseband_num_cols;
  
  /* Allocate temporary working arrays */
  if ((input_signal_vector =
       QccVectorIntAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDWT2DInt): Error calling QccVectorIntAlloc()");
      goto Error;
    }
  if ((output_signal_vector =
       QccVectorIntAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDWT2DInt): Error calling QccVectorIntAlloc()");
      goto Error;
    }
  if ((subsequence =
       QccVectorIntAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDWT2DInt): Error calling QccVectorIntAlloc()");
      goto Error;
    }
  if ((input_mask_vector =
       QccVectorIntAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDWT2DInt): Error calling QccVectorIntAlloc()");
      goto Error;
    }
  if ((output_mask_vector =
       QccVectorIntAlloc(QccMathMax(num_rows, num_cols))) == NULL)
    
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDWT2DInt): Error calling QccVectorIntAlloc()");
      goto Error;
    }
  
  for (scale = num_scales - 1; scale >= 0; scale--)
    {
      baseband_num_rows =
        QccWAVWaveletDWTSubbandLength(num_rows, scale, 0, 0, 0);
      baseband_num_cols =
        QccWAVWaveletDWTSubbandLength(num_cols, scale, 0, 0, 0);
      
      for (col = 0; col < baseband_num_cols; col++)
        {
          for (row = 0; row < baseband_num_rows; row++)