sadwt_int.c

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

                                                         baseband_num_frames,
                                                         wavelet))
              {
                QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDyadicDWT3DInt): Error calling QccWAVWaveletShapeAdaptiveSynthesis1DInt()");
                goto Error;
              }
            
            for (frame = 0; frame < baseband_num_frames; frame++)
              {
                signal[frame][row][col] = output_signal_vector[frame];
                mask[frame][row][col] = output_mask_vector[frame];
              }
          }
      
      for (frame = 0; frame < baseband_num_frames; frame++)
        for (col = 0; col < baseband_num_cols; col++)
          {
            for (row = 0; row < baseband_num_rows; row++)
              {
                input_signal_vector[row] = signal[frame][row][col];
                input_mask_vector[row] = mask[frame][row][col];
              }
            
            if (QccWAVWaveletShapeAdaptiveSynthesis1DInt(input_signal_vector,
                                                         output_signal_vector,
                                                         subsequence,
                                                         input_mask_vector,
                                                         output_mask_vector,
                                                         baseband_num_rows,
                                                         wavelet))
              {
                QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDyadicDWT3DInt): Error calling QccWAVWaveletShapeAdaptiveSynthesis1DInt()");
                goto Error;
              }
            
            for (row = 0; row < baseband_num_rows; row++)
              {
                signal[frame][row][col] = output_signal_vector[row];
                mask[frame][row][col] = output_mask_vector[row];
              }
          }
      
      for (frame = 0; frame < baseband_num_frames; frame++)
        for (row = 0; row < baseband_num_rows; row++)
          {
            for (col = 0; col < baseband_num_cols; col++)
              {
                input_signal_vector[col] = signal[frame][row][col];
                input_mask_vector[col] = mask[frame][row][col];
              }
            
            if (QccWAVWaveletShapeAdaptiveSynthesis1DInt(input_signal_vector,
                                                         output_signal_vector,
                                                         subsequence,
                                                         input_mask_vector,
                                                         output_mask_vector,
                                                         baseband_num_cols,
                                                         wavelet))
              {
                QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptiveDyadicDWT3DInt): Error calling QccWAVWaveletShapeAdaptiveSynthesis1DInt()");
                goto Error;
              }
            
            for (col = 0; col < baseband_num_cols; col++)
              {
                signal[frame][row][col] = output_signal_vector[col];
                mask[frame][row][col] = output_mask_vector[col];
              }
          }
    }
  
  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 QccWAVWaveletShapeAdaptivePacketDWT3DInt(QccVolumeInt signal,
                                             QccVolumeInt mask,
                                             int num_frames,
                                             int num_rows,
                                             int num_cols,
                                             int temporal_num_scales,
                                             int spatial_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 frame, row, col;
  int baseband_num_frames;

  /* Allocate temporary working arrays */
  if ((input_signal_vector =
       QccVectorIntAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3DInt): Error calling VectorAlloc()");
      goto Error;
    }
  if ((output_signal_vector =
       QccVectorIntAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3DInt): Error calling VectorAlloc()");
      goto Error;
    }
  if ((subsequence =
       QccVectorIntAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3DInt): Error calling VectorAlloc()");
      goto Error;
    }
  if ((input_mask_vector =
       QccVectorIntAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3DInt): Error calling VectorAlloc()");
      goto Error;
    }
  if ((output_mask_vector =
       QccVectorIntAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWaveletShapeAdaptivePacketDWT3DInt): Error calling VectorAlloc()");
      goto Error;
    }
  
  for (scale = 0; scale < temporal_num_scales; scale++)
    {
      baseband_num_frames =
        QccWAVWaveletDWTSubbandLength(num_frames, scale, 0, 0, 0);

      for (row = 0; row < num_rows; row++)
        for (col = 0; col < num_cols; col++)
          {
            for (frame = 0; frame < baseband_num_frames; frame++)
              {
                input_signal_vector[frame] = signal[frame][row][col];
                input_mask_vector[frame] = mask[frame][row][col];
              }
            
            if (QccWAVWaveletShapeAdaptiveAnalysis1DInt(input_signal_vector,
                                                        output_signal_vector,
                                                        subsequence,
                                                        input_mask_vector,
                                                        output_mask_vector,
                                                        baseband_num_frames,
                                                        wavelet))
              {
                QccErrorAddMessage("(QccWAVWaveletShapeAdaptivePacketDWT3DInt): Error calling QccWAVWaveletShapeAdaptiveAnalysis1DInt()");
                goto Error;
              }
            
            for (frame = 0; frame < baseband_num_frames; frame++)
              {
                signal[frame][row][col] = output_signal_vector[frame];
                mask[frame][row][col] = output_mask_vector[frame];
              }
          }
    }
      
  for (frame = 0; frame < num_frames; frame++)
    if (QccWAVWaveletShapeAdaptiveDWT2DInt(signal[frame],
                                           mask[frame],
                                           num_rows,
                                           num_cols,
                                           spatial_num_scales,
                                           wavelet))
      {
        QccErrorAddMessage("(QccWAVWaveletShapeAdaptivePacketDWT3DInt): Error calling QccWAVWaveletShapeAdaptiveDWT2DInt()");
        goto Error;
      }
  
  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 QccWAVWaveletInverseShapeAdaptivePacketDWT3DInt(QccVolumeInt signal,
                                                    QccVolumeInt mask,
                                                    int num_frames,
                                                    int num_rows,
                                                    int num_cols,
                                                    int temporal_num_scales,
                                                    int spatial_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 frame, row, col;
  int baseband_num_frames;

  /* Allocate temporary working arrays */
  if ((input_signal_vector =
       QccVectorIntAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3DInt): Error calling QccVectorIntAlloc()");
      goto Error;
    }
  if ((output_signal_vector =
       QccVectorIntAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3DInt): Error calling QccVectorIntAlloc()");
      goto Error;
    }
  if ((subsequence =
       QccVectorIntAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3DInt): Error calling QccVectorIntAlloc()");
      goto Error;
    }
  if ((input_mask_vector =
       QccVectorIntAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3DInt): Error calling QccVectorIntAlloc()");
      goto Error;
    }
  if ((output_mask_vector =
       QccVectorIntAlloc(num_frames)) == NULL)
    {
      QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3DInt): Error calling QccVectorIntAlloc()");
      goto Error;
    }
  
  for (frame = 0; frame < num_frames; frame++)
    if (QccWAVWaveletInverseShapeAdaptiveDWT2DInt(signal[frame],
                                                  mask[frame],
                                                  num_rows,
                                                  num_cols,
                                                  spatial_num_scales,
                                                  wavelet))
      {
        QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3DInt): Error calling QccWAVWaveletInverseShapeAdaptiveDWT2DInt()");
        goto Error;
      }

  for (scale = temporal_num_scales - 1; scale >= 0; scale--)
    {
      baseband_num_frames =
        QccWAVWaveletDWTSubbandLength(num_frames, scale, 0, 0, 0);
      
      for (row = 0; row < num_rows; row++)
        for (col = 0; col < num_cols; col++)
          {
            for (frame = 0; frame < baseband_num_frames; frame++)
              {
                input_signal_vector[frame] = signal[frame][row][col];
                input_mask_vector[frame] = mask[frame][row][col];
              }
            
            if (QccWAVWaveletShapeAdaptiveSynthesis1DInt(input_signal_vector,
                                                         output_signal_vector,
                                                         subsequence,
                                                         input_mask_vector,
                                                         output_mask_vector,
                                                         baseband_num_frames,
                                                         wavelet))
              {
                QccErrorAddMessage("(QccWAVWaveletInverseShapeAdaptivePacketDWT3DInt): Error calling QccWAVWaveletShapeAdaptiveSynthesis1DInt()");
                goto Error;
              }
            
            for (frame = 0; frame < baseband_num_frames; frame++)
              {
                signal[frame][row][col] = output_signal_vector[frame];
                mask[frame][row][col] = output_mask_vector[frame];
              }
          }
    }
  
  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);
}