---

spectrum.cpp

Go to the documentation of this file.

#include <vector>

#include "spectrum.h"

void spectrum::spectralCalc( const double *a, 
                             const size_t N,
                             std::vector<double> &v )
{
  size_t end = 2;
  size_t start = 1;

  while (end <= N) {
    double power = 0;
    for (size_t i = start; i < end; i++) {
      power = power + (a[i] * a[i]);
    }
    v.push_back( power );
    start = end;
    end = end << 1;
  }
} // spectralCalc




void spectrum::copyBands( double *dest,
                          const double *src,
                          const size_t N,
                          size_t startBand,
                          const size_t endBand )
{
  size_t end = 1;
  size_t start = 0;
  size_t bandCnt = 0;

  while (end <= N) {
    double power = 0;
    for (size_t i = start; i < end; i++) {
      if (bandCnt >= startBand && bandCnt <= endBand)
        dest[i] = src[i];
      else
        dest[i] = 0;
    }
    start = end;
    end = end << 1;
    bandCnt++;
  }
} // copyBands

---