lifting_int_cdf5_3.c

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

              
            default:
              QccErrorAddMessage("(QccWAVLiftingAnalysisIntCohenDaubechiesFeauveau5_3): Undefined boundary method");
              return(1);
            }
          break;
        }
      break;
    }
  
  return(0);
}


static void QccWAVCDF53SynthesisIntSymmetricEvenEven(QccVectorInt signal,
                                                     int signal_length)
{
  int index;
  
  for (index = 0; index < signal_length; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index + 1], 75, 6);

  for (index = 1; index < signal_length; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index - 1], 181, 7);

  for (index = 0; index < signal_length; index += 2)
    signal[index] -=
      QccWAVLiftingIntRound(signal[index + 1], 53, 6);

  for (index = 1; index < signal_length; index += 2)
    signal[index] = signal[index - 1] - signal[index];

  signal[0] -= QccWAVLiftingIntRound(signal[1], 1, 0);
  for (index = 2; index < signal_length; index += 2)
    signal[index] -=
      QccWAVLiftingIntRound(signal[index + 1] + signal[index - 1], 1, 1);

  for (index = 1; index < signal_length - 2; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index - 1] + signal[index + 1], 1, 0);
  signal[signal_length - 1] += signal[signal_length - 2];
  
}


static void QccWAVCDF53SynthesisIntSymmetricEvenOdd(QccVectorInt signal,
                                                    int signal_length)
{
  int index;
  
  for (index = 1; index < signal_length; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index - 1], 75, 6);

  for (index = 0; index < signal_length; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index + 1], 181, 7);

  for (index = 1; index < signal_length; index += 2)
    signal[index] -=
      QccWAVLiftingIntRound(signal[index - 1], 53, 6);

  for (index = 0; index < signal_length; index += 2)
    signal[index] = signal[index + 1] - signal[index];

  for (index = 1; index < signal_length - 2; index += 2)
    signal[index] -=
      QccWAVLiftingIntRound(signal[index + 1] + signal[index - 1], 1, 1);
  signal[signal_length - 1] -=
    QccWAVLiftingIntRound(signal[signal_length - 2], 1, 0);
  
  signal[0] += signal[1];
  for (index = 2; index < signal_length; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index - 1] + signal[index + 1], 1, 0);
  
}


static void QccWAVCDF53SynthesisIntSymmetricOddEven(QccVectorInt signal,
                                                    int signal_length)
{
  int index;
  
  for (index = 0; index < signal_length - 2; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index + 1], 75, 6);
  signal[signal_length - 1] +=
      QccWAVLiftingIntRound(signal[signal_length - 2], 75, 6);

  for (index = 1; index < signal_length; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index - 1], 181, 7);

  for (index = 0; index < signal_length - 2; index += 2)
    signal[index] -=
      QccWAVLiftingIntRound(signal[index + 1], 53, 6);
  signal[signal_length - 1] -=
      QccWAVLiftingIntRound(signal[signal_length - 2], 53, 6);

  for (index = 1; index < signal_length; index += 2)
    signal[index] = signal[index - 1] - signal[index];

  signal[0] -= QccWAVLiftingIntRound(signal[1], 1, 0);
  for (index = 2; index < signal_length - 2; index += 2)
    signal[index] -=
      QccWAVLiftingIntRound(signal[index + 1] + signal[index - 1], 1, 1);
  signal[signal_length - 1] -=
    QccWAVLiftingIntRound(signal[signal_length - 2], 1, 0);
  
  for (index = 1; index < signal_length; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index - 1] + signal[index + 1], 1, 0);
  
}


static void QccWAVCDF53SynthesisIntSymmetricOddOdd(QccVectorInt signal,
                                                   int signal_length)
{
  int index;
  
  for (index = 1; index < signal_length; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index - 1], 75, 6);

  for (index = 0; index < signal_length - 2; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index + 1], 181, 7);
  signal[signal_length - 1] +=
    QccWAVLiftingIntRound(signal[signal_length - 2], 181, 7);

  for (index = 1; index < signal_length; index += 2)
    signal[index] -=
      QccWAVLiftingIntRound(signal[index - 1], 53, 6);

  for (index = 0; index < signal_length - 2; index += 2)
    signal[index] = signal[index + 1] - signal[index];
  signal[signal_length - 1] =
    signal[signal_length - 2] - signal[signal_length - 1];

  for (index = 1; index < signal_length; index += 2)
    signal[index] -=
      QccWAVLiftingIntRound(signal[index + 1] + signal[index - 1], 1, 1);

  signal[0] += signal[1];
  for (index = 2; index < signal_length - 2; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index - 1] + signal[index + 1], 1, 0);
  signal[signal_length - 1] += signal[signal_length - 2];

}


static void QccWAVCDF53SynthesisIntPeriodicEvenEven(QccVectorInt signal,
                                                    int signal_length)
{
  int index;
  
  for (index = 0; index < signal_length; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index + 1], 75, 6);

  for (index = 1; index < signal_length; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index - 1], 181, 7);

  for (index = 0; index < signal_length; index += 2)
    signal[index] -=
      QccWAVLiftingIntRound(signal[index + 1], 53, 6);

  for (index = 1; index < signal_length; index += 2)
    signal[index] = signal[index - 1] - signal[index];

  signal[0] -=
    QccWAVLiftingIntRound(signal[1] + signal[signal_length - 1], 1, 1);
  for (index = 2; index < signal_length; index += 2)
    signal[index] -=
      QccWAVLiftingIntRound(signal[index + 1] + signal[index - 1], 1, 1);
  
  for (index = 1; index < signal_length - 2; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index - 1] + signal[index + 1], 1, 0);
  signal[signal_length - 1] +=
    QccWAVLiftingIntRound(signal[signal_length - 2] + signal[0], 1, 0);
  
}


static void QccWAVCDF53SynthesisIntPeriodicEvenOdd(QccVectorInt signal,
                                                   int signal_length)
{
  int index;

  for (index = 1; index < signal_length; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index - 1], 75, 6);

  for (index = 0; index < signal_length; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index + 1], 181, 7);

  for (index = 1; index < signal_length; index += 2)
    signal[index] -=
      QccWAVLiftingIntRound(signal[index - 1], 53, 6);

  for (index = 0; index < signal_length; index += 2)
    signal[index] = signal[index + 1] - signal[index];

  for (index = 1; index < signal_length - 1; index += 2)
    signal[index] -=
      QccWAVLiftingIntRound(signal[index + 1] + signal[index - 1], 1, 1);
  signal[signal_length - 1] -=
    QccWAVLiftingIntRound(signal[0] + signal[signal_length - 2], 1, 1);

  signal[0] +=
    QccWAVLiftingIntRound(signal[signal_length - 1] + signal[1], 1, 0);
  for (index = 2; index < signal_length; index += 2)
    signal[index] +=
      QccWAVLiftingIntRound(signal[index - 1] + signal[index + 1], 1, 0);
  
}


int QccWAVLiftingSynthesisIntCohenDaubechiesFeauveau5_3(QccVectorInt signal,
                                                        int signal_length,
                                                        int phase,
                                                        int boundary)
{
  if (signal == NULL)
    return(0);
  
  if (!signal_length)
    return(0);
  
  if (signal_length == 1)
    {
      if (phase == QCCWAVWAVELET_PHASE_EVEN)
        signal[0] /= M_SQRT2;
      else
        signal[0] *= M_SQRT2;
      return(0);
    }
  
  if ((boundary == QCCWAVWAVELET_BOUNDARY_BOUNDARY_WAVELET) &&
      ((signal_length == 2) || (signal_length == 3)))
    boundary = QCCWAVWAVELET_BOUNDARY_SYMMETRIC_EXTENSION;
  
  switch (phase)
    {
    case QCCWAVWAVELET_PHASE_EVEN:
      switch (signal_length % 2)
        {
        case 0:
          switch (boundary)
            {
            case QCCWAVWAVELET_BOUNDARY_SYMMETRIC_EXTENSION:
              QccWAVCDF53SynthesisIntSymmetricEvenEven(signal, signal_length);
              break;
              
            case QCCWAVWAVELET_BOUNDARY_PERIODIC_EXTENSION:
              QccWAVCDF53SynthesisIntPeriodicEvenEven(signal, signal_length);
              break;
              
            default:
              QccErrorAddMessage("(QccWAVLiftingSynthesisIntCohenDaubechiesFeauveau5_3): Undefined boundary method");
              return(1);
            }
          break;
          
        case 1:
          switch (boundary)
            {
            case QCCWAVWAVELET_BOUNDARY_SYMMETRIC_EXTENSION:
              QccWAVCDF53SynthesisIntSymmetricOddEven(signal, signal_length);
              break;
              
            case QCCWAVWAVELET_BOUNDARY_PERIODIC_EXTENSION:
              QccErrorAddMessage("(QccWAVLiftingSynthesisIntCohenDaubechiesFeauveau5_3): Signal length must be even for periodic extension");
              return(1);
              break;
              
            default:
              QccErrorAddMessage("(QccWAVLiftingSynthesisIntCohenDaubechiesFeauveau5_3): Undefined boundary method");
              return(1);
            }
          break;
        }
      break;

    case QCCWAVWAVELET_PHASE_ODD:
      switch (signal_length % 2)
        {
        case 0:
          switch (boundary)
            {
            case QCCWAVWAVELET_BOUNDARY_SYMMETRIC_EXTENSION:
              QccWAVCDF53SynthesisIntSymmetricEvenOdd(signal, signal_length);
              break;
              
            case QCCWAVWAVELET_BOUNDARY_PERIODIC_EXTENSION:
              QccWAVCDF53SynthesisIntPeriodicEvenOdd(signal, signal_length);
              break;
              
            default:
              QccErrorAddMessage("(QccWAVLiftingSynthesisIntCohenDaubechiesFeauveau5_3): Undefined boundary method");
              return(1);
            }
          break;
          
        case 1:
          switch (boundary)
            {
            case QCCWAVWAVELET_BOUNDARY_SYMMETRIC_EXTENSION:
              QccWAVCDF53SynthesisIntSymmetricOddOdd(signal, signal_length);
              break;
              
            case QCCWAVWAVELET_BOUNDARY_PERIODIC_EXTENSION:
              QccErrorAddMessage("(QccWAVLiftingSynthesisIntCohenDaubechiesFeauveau5_3): Signal length must be even for periodic extension");
              return(1);
              break;
              
            default:
              QccErrorAddMessage("(QccWAVLiftingSynthesisIntCohenDaubechiesFeauveau5_3): Undefined boundary method");
              return(1);
            }
          break;
        }
      break;
    }

  return(0);
}