// g++ `root-config --cflags` `root-config --libs` -L/usr/lib/x86_64-linux-gnu/root5.34 -lFFTW -lTreeViewer -I/home/ecm/Notes/13chirp_cwb/watrepo/tools/install/inc  test_writewdm.cc /home/ecm/Notes/13chirp_cwb/watrepo/tools/install/lib/wavelet.so -o test_writewdm

#include <math.h>
#include <iostream>
#include <stdio.h>
#include "TFile.h"
#include "wseries.hh"
#include "WDM.hh"

using namespace std;
 
void test_wdm()
{

  int N=16364;
  wavearray<double> x(N); 
  for(int i=0; i<N; i++) x[i] = sin(2*3.1415*i/100.0);

  WDM<double> wdm(512, 512, 4, 10);
  WSeries<double> data;
  data.Forward(x, wdm);
  
  printf("metadata=%d\n",data.sizeZero());
  printf("data[0]=%g\n",data.pWavelet->pWWS[0]);

  remove("test_wdm.root");
  TFile *file = new TFile("test_wdm.root", "CREATE");
  data.Write("data");
  file->Close();
}
 
# ifndef __CINT__
int main()
{
  test_wdm();
  return 0;
}
# endif