32 if(NET==NULL) {cout <<
"ced::ced Error : NET is NULL" << endl;
exit(1);}
33 if(EVT==NULL) {cout <<
"ced::CED Error : EVT is NULL" << endl;
exit(1);}
34 if(sbasedirCED==NULL) {cout <<
"ced::ced Error : sbasedirCED is NULL" << endl;
exit(1);}
37 this->sbasedirCED = sbasedirCED;
38 this->rbasedirCED = NULL;
45 this->useSparse=
false;
50 if(NET==NULL) {cout <<
"ced::ced Error : NET is NULL" << endl;
exit(1);}
51 if(EVT==NULL) {cout <<
"ced::ced Error : EVT is NULL" << endl;
exit(1);}
52 if(rbasedirCED==NULL) {cout <<
"ced::ced Error : rbasedirCED is NULL" << endl;
exit(1);}
55 this->sbasedirCED = NULL;
56 this->rbasedirCED = rbasedirCED;
63 this->useSparse=
false;
73 if(inRate<=0) {cout <<
"ced::SetOptions : Error - inrate must be > 0" << endl;
exit(1);}
75 this->paletteId=paletteId;
76 this->useSparse=useSparse;
84 void Write(
double factor,
size_t iID,
int LAG,
char* dirCED);
cout<< "skymap size : "<< L<< endl;for(int l=0;l< L;l++) sm.set(l, l);sm > const_cast< char * >("skymap.dat")
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
void SetOptions(int simulation, double rho, double inRate, bool useSparse=false, char *gtype=const_cast< char * >("png"), int paletteId=0)
void SetChannelName(char *chName)
ced(network *NET, netevent *EVT, char *sbasedirCED)
ced(network *NET, netevent *EVT, TDirectory *rbasedirCED)
netevent EVT(itree, nifo)
strcpy(RunLabel, RUN_LABEL)