---

line_norm.h

Go to the documentation of this file.

#ifndef _LINE_NORM_H_
#define _LINE_NORM_H_

#include "line.h"

template <class T>
class line_norm : public line<T> {

protected:
  void normalize( T& vec, int N, transDirection direction )
  {
    const double sqrt2 = sqrt( 2.0 );
    int half = N >> 1;

    for (int i = 0; i < half; i++) {
      int j = i + half;

      if (direction == forward) {
        vec[i] = sqrt2 * vec[i];
        vec[j] = vec[j]/sqrt2;
      }
      else if (direction == inverse) {
        vec[i] = vec[i]/sqrt2;
        vec[j] = sqrt2 * vec[j];
      }
      else {
        printf("normalize: bad direction value\n");
      }
    } // for
  } // normalize

  void inverseStep( T& vec, const int n )
  {
    normalize( vec, n, inverse );
    update( vec, n, inverse );
    predict( vec, n, inverse );
    merge( vec, n );
  }  // inverseStep

  void forwardStep( T& vec, const int n )
  {
    split( vec, n );
    predict( vec, n, forward );
    update( vec, n, forward );
    normalize( vec, n, forward );
  } // forwardStep

};

#endif

---