16 #include "TMatrixDSym.h"
17 #include "TMatrixDSymEigen.h"
82 void solve(
double th,
int nE=0,
char c=
's');
84 void apply(
double threshold=0.,
char c=
'a');
135 int tsize=
vrank[
n].size();
143 chList.clear(); std::vector< WSeries<double> >().swap(
chList);
145 chMask.clear(); std::vector< wavearray<int> >().swap(
chMask);
147 matrix.clear(); std::vector<TMatrixDSym>().swap(
matrix);
148 vCROSS.clear(); std::vector< wavearray<double> >().swap(
vCROSS);
149 vEIGEN.clear(); std::vector< wavearray<double> >().swap(
vEIGEN);
150 vrank.clear(); std::vector< wavearray<double> >().swap(
vrank);
170 std::vector< wavearray<double> >
vrank;
184 #endif // REGRESSION_HH
void setMatrix(double edge=0., double f=1.)
void _apply_(int n, std::vector< wavearray< double > > &w, std::vector< wavearray< double > > &W)
WSeries< double > * getTFmap(int n=0)
WSeries< double > getWNoise()
wavearray< double > target
wavearray< double > getNoise()
TMatrixTSym< double > TMatrixDSym
cout<< endl;cout<< "ts size = "<< ts.size()<< " ts rate = "<< ts.rate()<< endl;tf.Forward(ts, wdm);int levels=tf.getLevel();cout<< "tf size = "<< tf.size()<< endl;double dF=tf.resolution();double dT=1./(2 *dF);cout<< "rate(hz) : "<< RATE<< "\t layers : "<< nLAYERS<< "\t dF(hz) : "<< dF<< "\t dT(ms) : "<< dT *1000.<< endl;int itime=TIME_PIXEL_INDEX;int ifreq=FREQ_PIXEL_INDEX;int index=(levels+1)*itime+ifreq;double time=itime *dT;double freq=(ifreq >0)?ifreq *dF:dF/4;cout<< endl;cout<< "PIXEL TIME = "<< time<< " sec "<< endl;cout<< "PIXEL FREQ = "<< freq<< " Hz "<< endl;cout<< endl;wavearray< double > x
TMatrixDSym getMatrix(size_t n=0)
std::vector< wavearray< double > > vCROSS
size_t add(WSeries< double > &target, char *name, double fL=0., double fH=0.)
std::vector< wavearray< double > > vEIGEN
std::vector< Wiener > FILTER
std::vector< WSeries< double > > chList
wavearray< double > rank(int nbins=0, double fL=0., double fH=0.)
std::vector< double > vectorD
std::vector< vectorD > filter00
void apply(double threshold=0., char c='a')
wavearray< double > getRank(int n)
std::vector< TMatrixDSym > matrix
std::vector< wavearray< int > > chMask
void solve(double th, int nE=0, char c='s')
std::vector< int > channel
wavearray< double > vfreq
std::vector< vectorD > filter90
wavearray< double > getVEIGEN(int n=-1)
wavearray< double > getClean()
std::vector< char * > chName
std::vector< double > norm
wavearray< double > getFILTER(char c='a', int nT=-1, int nW=-1)
void unmask(int n, double flow=0., double fhigh=0.)
void mask(int n, double flow=0., double fhigh=0.)
wavearray< double > getVCROSS(size_t n=0)
std::vector< wavearray< double > > vrank
regression & operator=(const regression &)
wavearray< double > rnoise
virtual void resize(unsigned int)
void Inverse(int n=-1)
param: n - number of steps (-1 means full reconstruction)
wavearray< double > channel(size_t n)