int main(int argc,char **argv)
{
long i;
char stdi=0;
double min,max,maxmax;
if (scan_help(argc,argv))
show_options(argv[0]);
scan_options(argc,argv);
#ifndef OMIT_WHAT_I_DO
if (verbosity&VER_INPUT)
what_i_do(argv[0],WID_STR);
#endif
infile=search_datafile(argc,argv,NULL,verbosity);
if (infile == NULL)
stdi=1;
if (outfile == NULL) {
if (!stdi) {
check_alloc(outfile=(char*)calloc(strlen(infile)+5,(size_t)1));
strcpy(outfile,infile);
strcat(outfile,".rec");
}
else {
check_alloc(outfile=(char*)calloc((size_t)10,(size_t)1));
strcpy(outfile,"stdin.rec");
}
}
if (!stdo)
test_outfile(outfile);
if (columns == NULL)
series=(double**)get_multi_series(infile,&length,exclude,&dim,"",dimset,
verbosity);
else
series=(double**)get_multi_series(infile,&length,exclude,&dim,columns,
dimset,verbosity);
maxmax=0.0;
for (i=0;i<dim;i++) {
rescale_data(series[i],length,&min,&max);
if (max > maxmax)
maxmax=max;
}
if (epsset)
eps /= maxmax;
check_alloc(list=(long*)malloc(sizeof(long)*length));
check_alloc(box=(long**)malloc(sizeof(long*)*BOX));
for (i=0;i<BOX;i++)
check_alloc(box[i]=(long*)malloc(sizeof(long)*BOX));
make_multi_box(series,box,list,length,BOX,dim,embed,delay,eps);
lfind_neighbors();
return 0;
}
int main(int argc,char **argv)
{
char stdi=0,*done,alldone;
int i;
long n;
long maxlength;
double min,max;
FILE *file;
if (scan_help(argc,argv))
show_options(argv[0]);
scan_options(argc,argv);
#ifndef OMIT_WHAT_I_DO
if (verbosity&VER_INPUT)
what_i_do(argv[0],WID_STR);
#endif
infile=search_datafile(argc,argv,&column,verbosity);
if (infile == NULL)
stdi=1;
if (outfile == NULL) {
if (!stdi) {
check_alloc(outfile=(char*)calloc(strlen(infile)+5,(size_t)1));
strcpy(outfile,infile);
strcat(outfile,".ros");
}
else {
check_alloc(outfile=(char*)calloc((size_t)10,(size_t)1));
strcpy(outfile,"stdin.ros");
}
}
test_outfile(outfile);
series=(double*)get_series(infile,&length,exclude,column,verbosity);
rescale_data(series,length,&min,&max);
if (epsset)
eps0 /= max;
check_alloc(list=(long*)malloc(length*sizeof(long)));
check_alloc(lyap=(double*)malloc((steps+1)*sizeof(double)));
check_alloc(found=(long*)malloc((steps+1)*sizeof(long)));
check_alloc(done=(char*)malloc(length));
for (i=0;i<=steps;i++) {
lyap[i]=0.0;
found[i]=0;
}
for (i=0;i<length;i++)
done[i]=0;
maxlength=length-delay*(dim-1)-steps-1-mindist;
alldone=0;
file=fopen(outfile,"w");
if (verbosity&VER_INPUT)
fprintf(stderr,"Opened %s for writing\n",outfile);
for (eps=eps0;!alldone;eps*=1.1) {
epsinv=1.0/eps;
put_in_boxes();
alldone=1;
for (n=0;n<=maxlength;n++) {
if (!done[n])
done[n]=make_iterate(n);
alldone &= done[n];
}
if (verbosity&VER_USR1)
fprintf(stderr,"epsilon: %e already found: %ld\n",eps*max,found[0]);
}
for (i=0;i<=steps;i++)
if (found[i])
fprintf(file,"%d %e\n",i,lyap[i]/found[i]/2.0);
fclose(file);
return 0;
}
int main(int argc,char **argv)
{
char stdi=0;
FILE *file;
double min,max;
double av,varianz;
long i,n;
if (scan_help(argc,argv))
show_options(argv[0]);
scan_options(argc,argv);
#ifndef OMIT_WHAT_I_DO
if (verbosity&VER_INPUT)
what_i_do(argv[0],WID_STR);
#endif
infile=search_datafile(argc,argv,NULL,verbosity);
if (infile == NULL)
stdi=1;
if (outfile == NULL) {
if (!stdi) {
check_alloc(outfile=(char*)calloc(strlen(infile)+5,(size_t)1));
strcpy(outfile,infile);
strcat(outfile,".res");
}
else {
check_alloc(outfile=(char*)calloc((size_t)10,(size_t)1));
strcpy(outfile,"stdin.res");
}
}
if (!stdo)
test_outfile(outfile);
if (xmin >= xmax) {
fprintf(stderr,"Choosing the minimum larger or equal the maximum\n"
"makes no sense. Exiting!\n");
exit(RESCALE__WRONG_INTERVAL);
}
if (column == NULL)
series=(double**)get_multi_series(infile,&length,exclude,&dim,"",dimset,
verbosity);
else
series=(double**)get_multi_series(infile,&length,exclude,&dim,column,
dimset,verbosity);
for (n=0;n<dim;n++) {
variance(series[n],length,&av,&varianz);
if (ddof!=0)
varianz *= sqrt((double) length/(length-ddof));
if (set_av)
for (i=0;i<length;i++)
series[n][i] -= av;
if (set_var)
for (i=0;i<length;i++)
series[n][i] /= varianz;
if (!set_var && !set_av) {
rescale_data(series[n],length,&min,&max);
for (i=0;i<length;i++)
series[n][i]=series[n][i]*(xmax-xmin)+xmin;
}
}
if (stdo) {
if (verbosity&VER_INPUT)
fprintf(stderr,"Writing to stdout\n");
for (i=0;i<length;i++) {
fprintf(stdout,"%e",series[0][i]);
for (n=1;n<dim;n++)
fprintf(stdout," %e",series[n][i]);
fprintf(stdout,"\n");
}
}
else {
file=fopen(outfile,"w");
if (verbosity&VER_INPUT)
fprintf(stderr,"Opened %s for writing\n",outfile);
for (i=0;i<length;i++) {
fprintf(file,"%e",series[0][i]);
for (n=1;n<dim;n++)
fprintf(file," %e",series[n][i]);
fprintf(file,"\n");
}
fclose(file);
}
return 0;
}
int main(int argc,char **argv)
{
char stdi=0,done;
long i,j,hdim,actfound;
unsigned long count=1;
double *swap,*newcast,maxinterval,*min,*interval,dummy,epsilon0;
FILE *file=NULL;
if (scan_help(argc,argv))
show_options(argv[0]);
scan_options(argc,argv);
#ifndef OMIT_WHAT_I_DO
if (verbosity&VER_INPUT)
what_i_do(argv[0],WID_STR);
#endif
infile=search_datafile(argc,argv,NULL,verbosity);
if (infile == NULL)
stdi=1;
if (outfile == NULL) {
if (!stdi) {
check_alloc(outfile=(char*)calloc(strlen(infile)+5,(size_t)1));
strcpy(outfile,infile);
strcat(outfile,".lzr");
}
else {
check_alloc(outfile=(char*)calloc((size_t)10,(size_t)1));
strcpy(outfile,"stdin.lzr");
}
}
if (!onscreen)
test_outfile(outfile);
hdim=(embed-1)*DELAY+1;
if (column == NULL)
series=(double**)get_multi_series(infile,&LENGTH,exclude,&dim,"",dimset,
verbosity);
else
series=(double**)get_multi_series(infile,&LENGTH,exclude,&dim,column,
dimset,verbosity);
dim1=dim-1;
check_alloc(min=(double*)malloc(sizeof(double)*dim));
check_alloc(interval=(double*)malloc(sizeof(double)*dim));
check_alloc(var=(double*)malloc(sizeof(double)*dim));
maxinterval=0.0;
for (i=0;i<dim;i++) {
rescale_data(series[i],LENGTH,&min[i],&interval[i]);
variance(series[i],LENGTH,&dummy,&var[i]);
if (interval[i] > maxinterval)
maxinterval=interval[i];
}
if (epsset)
EPS0 /= maxinterval;
check_alloc(cast=(double**)malloc(sizeof(double*)*hdim));
for (i=0;i<hdim;i++)
check_alloc(cast[i]=(double*)malloc(sizeof(double)*dim));
check_alloc(newcast=(double*)malloc(sizeof(double)*dim));
check_alloc(newav=(double*)malloc(sizeof(double)*dim));
check_alloc(list=(long*)malloc(sizeof(long)*LENGTH));
check_alloc(found=(long*)malloc(sizeof(long)*LENGTH));
check_alloc(abstand=(double*)malloc(sizeof(double)*LENGTH));
check_alloc(box=(long**)malloc(sizeof(long*)*NMAX));
for (i=0;i<NMAX;i++)
check_alloc(box[i]=(long*)malloc(sizeof(long)*NMAX));
check_alloc(vec=(double*)malloc(sizeof(double)*dim));
check_alloc(hsum=(double*)malloc(sizeof(double)*dim));
check_alloc(mat=(double**)malloc(sizeof(double*)*dim));
for (i=0;i<dim;i++) {
check_alloc(mat[i]=(double*)malloc(sizeof(double)*dim));
}
for (j=0;j<dim;j++)
for (i=0;i<hdim;i++)
cast[i][j]=series[j][LENGTH-hdim+i];
indexes=make_multi_index(dim,embed,DELAY);
if (!onscreen) {
file=fopen(outfile,"w");
if (verbosity&VER_INPUT)
fprintf(stderr,"Opened %s for writing\n",outfile);
}
else {
if (verbosity&VER_INPUT)
fprintf(stderr,"Writing to stdout\n");
}
rnd_init(seed);
epsilon0=EPS0/EPSF;
if (setnoise)
Q /= 100.0;
for (i=0;i<FLENGTH;i++) {
done=0;
if (setsort)
epsilon= epsilon0/((double)count*EPSF);
else
epsilon=epsilon0;
while (!done) {
epsilon*=EPSF;
put_in_boxes();
actfound=hfind_neighbors();
if (actfound >= MINN) {
if (setsort) {
epsilon0 += epsilon;
count++;
sort(actfound);
actfound=MINN;
}
make_zeroth(actfound,newcast);
if (onscreen) {
for (j=0;j<dim-1;j++)
printf("%e ",newcast[j]*interval[j]+min[j]);
printf("%e\n",newcast[dim-1]*interval[dim-1]+min[dim-1]);
fflush(stdout);
}
else {
for (j=0;j<dim-1;j++)
fprintf(file,"%e ",newcast[j]*interval[j]+min[j]);
fprintf(file,"%e\n",newcast[dim-1]*interval[dim-1]+min[dim-1]);
fflush(file);
}
done=1;
swap=cast[0];
for (j=0;j<hdim-1;j++)
cast[j]=cast[j+1];
cast[hdim-1]=swap;
for (j=0;j<dim;j++)
cast[hdim-1][j]=newcast[j];
}
}
}
if (!onscreen)
fclose(file);
if (outfile != NULL)
free(outfile);
for (i=0;i<dim;i++)
free(mat[i]);
free(mat);
for (i=0;i<hdim;i++)
free(cast[i]);
free(cast);
free(newcast);
free(found);
free(list);
for (i=0;i<NMAX;i++)
free(box[i]);
free(box);
free(vec);
free(newav);
for (i=0;i<dim;i++)
free(series[i]);
free(series);
return 0;
}
int main(int argc,char **argv)
{
char stdi=0;
FILE *file=NULL;
double min,inter=0.0,ind_inter,epsilon,av,ind_var;
char *nearest,alldone;
long i;
unsigned int dim,emb;
unsigned long donesofar;
if (scan_help(argc,argv))
show_options(argv[0]);
scan_options(argc,argv);
#ifndef OMIT_WHAT_I_DO
if (verbosity&VER_INPUT)
what_i_do(argv[0],WID_STR);
#endif
infile=search_datafile(argc,argv,NULL,verbosity);
if (infile == NULL)
stdi=1;
if (outfile == NULL) {
if (!stdi) {
check_alloc(outfile=(char*)calloc(strlen(infile)+5,(size_t)1));
strcpy(outfile,infile);
strcat(outfile,".fnn");
}
else {
check_alloc(outfile=(char*)calloc((size_t)10,(size_t)1));
strcpy(outfile,"stdin.fnn");
}
}
if (!stdo)
test_outfile(outfile);
if (column == NULL)
series=(double**)get_multi_series(infile,&length,exclude,&comp,"",dimset,
verbosity);
else
series=(double**)get_multi_series(infile,&length,exclude,&comp,column,
dimset,verbosity);
for (i=0;i<comp;i++) {
rescale_data(series[i],length,&min,&ind_inter);
variance(series[i],length,&av,&ind_var);
if (i == 0) {
varianz=ind_var;
inter=ind_inter;
}
else {
varianz=(varianz>ind_var)?ind_var:varianz;
inter=(inter<ind_inter)?ind_inter:inter;
}
}
check_alloc(list=(long*)malloc(sizeof(long)*length));
check_alloc(nearest=(char*)malloc(length));
check_alloc(box=(long**)malloc(sizeof(long*)*BOX));
for (i=0;i<BOX;i++)
check_alloc(box[i]=(long*)malloc(sizeof(long)*BOX));
if (!stdo) {
file=fopen(outfile,"w");
if (verbosity&VER_INPUT)
fprintf(stderr,"Opened %s for writing\n",outfile);
}
else {
if (verbosity&VER_INPUT)
fprintf(stderr,"Writing to stdout\n");
}
check_alloc(vcomp=(unsigned int*)malloc(sizeof(int)*(maxdim)));
check_alloc(vemb=(unsigned int*)malloc(sizeof(int)*(maxdim)));
for (i=0;i<maxdim;i++) {
if (comp == 1) {
vcomp[i]=0;
vemb[i]=i;
}
else {
vcomp[i]=i%comp;
vemb[i]=(i/comp)*delay;
}
}
for (emb=minemb;emb<=maxemb;emb++) {
dim=emb*comp-1;
epsilon=eps0;
toolarge=0;
alldone=0;
donesofar=0;
aveps=0.0;
vareps=0.0;
for (i=0;i<length;i++)
nearest[i]=0;
if (verbosity&VER_USR1)
fprintf(stderr,"Start for dimension=%u\n",dim+1);
while (!alldone && (epsilon < 2.*varianz/rt)) {
alldone=1;
mmb(vcomp[dim],vemb[dim],epsilon);
for (i=0;i<length-maxemb*delay;i++)
if (!nearest[i]) {
nearest[i]=find_nearest(i,dim,epsilon);
alldone &= nearest[i];
donesofar += (unsigned long)nearest[i];
}
if (verbosity&VER_USR1)
fprintf(stderr,"Found %lu up to epsilon=%e\n",donesofar,epsilon*inter);
epsilon*=sqrt(2.0);
if (!donesofar)
eps0=epsilon;
}
if (donesofar == 0) {
fprintf(stderr,"Not enough points found!\n");
exit(FALSE_NEAREST_NOT_ENOUGH_POINTS);
}
aveps *= (1./(double)donesofar);
vareps *= (1./(double)donesofar);
if (stdo) {
fprintf(stdout,"%u %e %e %e\n",dim+1,(double)toolarge/(double)donesofar,
aveps*inter,sqrt(vareps)*inter);
fflush(stdout);
}
else {
fprintf(file,"%u %e %e %e\n",dim+1,(double)toolarge/(double)donesofar,
aveps*inter,sqrt(vareps)*inter);
fflush(file);
}
}
if (!stdo)
fclose(file);
if (infile != NULL)
free(infile);
if (outfile != NULL)
free(outfile);
free(series);
free(list);
free(nearest);
for (i=0;i<BOX;i++)
free(box[i]);
free(box);
return 0;
}
int main(int argc,char **argv)
{
unsigned int outstrlen=0;
int i,j,k,n;
unsigned int offset,wdim=1,alldim;
unsigned int has_future,*future,**useries,h;
unsigned long count;
long *segment;
char *outstring,*infile,*wcol,*use_array;
double **series,*sermin,*serinter,**dfuture;
ptree *root;
FILE *fout;
if (scan_help(argc,argv))
show_options(argv[0]);
for (i=0;i<argc;i++) {
outstrlen += strlen(argv[i]);
outstrlen++;
}
check_alloc(outstring=(char*)malloc(sizeof(char)*(outstrlen+8)));
sprintf(outstring,"#Prog: ");
offset=7;
for (i=0;i<argc;i++) {
sprintf(outstring+offset,"%s ",argv[i]);
offset += (strlen(argv[i])+1);
}
scan_options(argc,argv);
if (!compdimset)
compdim=pars.DIM;
else
if (compdim > pars.DIM) {
fprintf(stderr,"-C value larger -m Value.\n");
exit(LANGEVIN_MAIN_C_TOO_LARGE);
}
infile=search_datafile(argc,argv,NULL,VERBOSITY);
if (infile == NULL) {
fprintf(stderr,"No input datafile found.\n");
exit(LANGEVIN_MAIN_NO_INPUTFILE);
}
if (OUTFILE == NULL) {
check_alloc(OUTFILE=(char*)calloc(strlen(infile)+5,(size_t)1));
sprintf(OUTFILE,"%s.pru",infile);
}
OUTFILE=test_outfile(OUTFILE);
if (COLUMN == NULL)
series=(double**)get_multi_series(infile,&(pars.LENGTH),EXCLUDE,
&(pars.DIM),"",dimset,VERBOSITY);
else
series=(double**)get_multi_series(infile,&(pars.LENGTH),EXCLUDE,
&(pars.DIM),COLUMN,dimset,VERBOSITY);
check_alloc(sermin=malloc(sizeof(double)*pars.DIM));
check_alloc(serinter=malloc(sizeof(double)*pars.DIM));
rescale_data(series,pars,sermin,serinter);
alldim=pars.DIM;
pars.DIM=compdim;
check_alloc(useries=(unsigned int**)malloc(sizeof(int*)*pars.DIM));
check_alloc(PART=(unsigned int*)malloc(sizeof(int)*pars.DIM));
if (SPART == NULL)
for (i=0;i<pars.DIM;i++)
PART[i]=2;
else
set_part(PART,SPART,pars);
for (i=0;i<pars.DIM;i++) {
check_alloc(useries[i]=(unsigned int*)malloc(sizeof(int)*pars.LENGTH));
for (j=0;j<pars.LENGTH;j++) {
h=(unsigned int)((series[i][j]*(double)PART[i])/serinter[i]);
useries[i][j]=((h<PART[i])?h:(PART[i]-1));
}
}
check_alloc(future=(unsigned int*)malloc(sizeof(int)*pars.LENGTH));
if (WHICHFUTURE>0) {
check_alloc(wcol=(char*)calloc((size_t)10,(size_t)1));
sprintf(wcol,"%u",WHICHFUTURE);
dfuture=(double**)get_multi_series(infile,&(pars.LENGTH),EXCLUDE,
&wdim,wcol,1,0L);
for (i=0;i<pars.LENGTH-1;i++) {
future[i] = (dfuture[0][i]>0.0);
}
future[pars.LENGTH-1]=1;
free(dfuture);
}
else {
for (i=0;i<pars.LENGTH;i++)
future[i]=1;
}
root=make_ptree(PART[0]);
count=0;
check_alloc(segment=(long*)malloc(sizeof(long)*pars.LENGTH));
for (i=0;i<pars.LENGTH;i += ODEPTH) {
has_future=1;
for (j=0;j<ODEPTH;j++)
has_future &= future[i+j];
if (has_future) {
fill_tree(root,useries,pars,i,0,DEPTH,PART);
}
segment[count++]=i;
}
if (scrambleset)
lscramble(segment,count);
check_alloc(use_array=(char*)malloc(sizeof(char)*pars.LENGTH));
for (i=0;i<count;i++) {
n=segment[i];
has_future=1;
for (j=0;j<ODEPTH;j++)
has_future &= future[n+j];
if (has_future) {
h=read_tree(root,useries,pars,n,0,DEPTH);
if (h <= MAXOUT)
use_array[i]=1;
else
use_array[i]=0;
}
else
use_array[i]=2;
}
pars.DIM=alldim;
fout=fopen(OUTFILE,"w");
fprintf(fout,"%s\n",outstring);
fprintf(fout,"#Content: ");
for (i=0;i<pars.DIM;i++)
fprintf(fout,"x%d ",i+1);
fprintf(fout,"future time_index\n");
fflush(fout);
for (i=0;i<count;i++) {
n=segment[i];
if (use_array[i] == 1) {
for (j=0;j<(ODEPTH-1);j++) {
for (k=0;k<pars.DIM;k++)
fprintf(fout,"%f ",series[k][n+j]+sermin[k]);
fprintf(fout,"1 %u\n",n+j);
}
for (k=0;k<pars.DIM;k++)
fprintf(fout,"%f ",series[k][n+ODEPTH-1]+sermin[k]);
if (i<(count-1)) {
if (segment[i+1] == (n+ODEPTH) && use_array[i+1])
fprintf(fout,"1 %u\n",n+ODEPTH-1);
else
fprintf(fout,"0 %u\n",n+ODEPTH-1);
}
else
fprintf(fout,"0 %u\n",n+ODEPTH-1);
}
}
fclose(fout);
}
int main(int argc,char **argv)
{
char stdi=0,done;
long i,j,hdim,actfound;
double maxinterval,*swap,*newcast;
FILE *file=NULL;
if (scan_help(argc,argv))
show_options(argv[0]);
scan_options(argc,argv);
#ifndef OMIT_WHAT_I_DO
if (verbosity&VER_INPUT)
what_i_do(argv[0],WID_STR);
#endif
infile=search_datafile(argc,argv,NULL,verbosity);
if (infile == NULL)
stdi=1;
if (outfile == NULL) {
if (!stdi) {
check_alloc(outfile=(char*)calloc(strlen(infile)+6,(size_t)1));
strcpy(outfile,infile);
strcat(outfile,".cast");
}
else {
check_alloc(outfile=(char*)calloc((size_t)11,(size_t)1));
strcpy(outfile,"stdin.cast");
}
}
if (!onscreen)
test_outfile(outfile);
hdim=(embed-1)*DELAY+1;
if (column == NULL)
series=(double**)get_multi_series(infile,&LENGTH,exclude,&dim,"",dimset,
verbosity);
else
series=(double**)get_multi_series(infile,&LENGTH,exclude,&dim,column,
dimset,verbosity);
check_alloc(min=(double*)malloc(sizeof(double)*dim));
check_alloc(interval=(double*)malloc(sizeof(double)*dim));
dim1=dim-1;
maxinterval=0.0;
for (i=0;i<dim;i++) {
rescale_data(series[i],LENGTH,&min[i],&interval[i]);
if (interval[i] > maxinterval)
maxinterval=interval[i];
}
check_alloc(cast=(double**)malloc(sizeof(double*)*hdim));
for (i=0;i<hdim;i++)
check_alloc(cast[i]=(double*)malloc(sizeof(double)*dim));
check_alloc(newcast=(double*)malloc(sizeof(double)*dim));
check_alloc(list=(long*)malloc(sizeof(long)*LENGTH));
check_alloc(found=(long*)malloc(sizeof(long)*LENGTH));
check_alloc(box=(long**)malloc(sizeof(long*)*NMAX));
for (i=0;i<NMAX;i++)
check_alloc(box[i]=(long*)malloc(sizeof(long)*NMAX));
check_alloc(localav=(double*)malloc(sizeof(double)*dim*embed));
check_alloc(foreav=(double*)malloc(sizeof(double)*dim));
check_alloc(vec=(double*)malloc(sizeof(double)*dim*embed));
check_alloc(mat=(double**)malloc(sizeof(double*)*dim*embed));
for (i=0;i<dim*embed;i++)
check_alloc(mat[i]=(double*)malloc(sizeof(double)*dim*embed));
if (epsset)
EPS0 /= maxinterval;
for (j=0;j<dim;j++)
for (i=0;i<hdim;i++)
cast[i][j]=series[j][LENGTH-hdim+i];
if (!onscreen) {
file=fopen(outfile,"w");
if (verbosity&VER_INPUT)
fprintf(stderr,"Opened %s for writing\n",outfile);
}
else {
if (verbosity&VER_INPUT)
fprintf(stderr,"Writing to stdout\n");
}
for (i=0;i<FLENGTH;i++) {
done=0;
epsilon=EPS0/EPSF;
while (!done) {
epsilon*=EPSF;
put_in_boxes();
actfound=hfind_neighbors();
if (actfound >= MINN) {
if (!do_zeroth)
make_fit(actfound,newcast);
else
make_zeroth(actfound,newcast);
if (onscreen) {
for (j=0;j<dim-1;j++)
printf("%e ",newcast[j]*interval[j]+min[j]);
printf("%e\n",newcast[dim-1]*interval[dim-1]+min[dim-1]);
fflush(stdout);
}
else {
for (j=0;j<dim-1;j++)
fprintf(file,"%e ",newcast[j]*interval[j]+min[j]);
fprintf(file,"%e\n",newcast[dim-1]*interval[dim-1]+min[dim-1]);
fflush(file);
}
done=1;
for (j=0;j<dim;j++) {
if ((newcast[j] > 2.0) || (newcast[j] < -1.0)) {
fprintf(stderr,"Forecast failed. Escaping data region!\n");
exit(NSTEP__ESCAPE_REGION);
}
}
swap=cast[0];
for (j=0;j<hdim-1;j++)
cast[j]=cast[j+1];
cast[hdim-1]=swap;
for (j=0;j<dim;j++)
cast[hdim-1][j]=newcast[j];
}
}
}
if (!onscreen)
fclose(file);
if (outfile != NULL)
free(outfile);
for (i=0;i<embed*dim;i++)
free(mat[i]);
free(mat);
for (i=0;i<hdim;i++)
free(cast[i]);
free(cast);
free(newcast);
free(found);
free(list);
for (i=0;i<NMAX;i++)
free(box[i]);
free(box);
free(vec);
free(localav);
free(foreav);
free(min);
free(interval);
for (i=0;i<dim;i++)
free(series[i]);
free(series);
return 0;
}
int main(int argc,char **argv)
{
unsigned int outstrlen=0,offset,hdim,wdim=1;
int i,j,minnadd;
unsigned int *future,has_future;
long li,lj,n,refn;
char *outstring,*infile,*wcol;
char ntoolarge;
double **series,**rseries,*sermin,*serinter,**dfuture,max,**cast,sweights;
struct sfound sf;
FILE *fout;
if (scan_help(argc,argv))
show_options(argv[0]);
for (i=0;i<argc;i++) {
outstrlen += strlen(argv[i]);
outstrlen++;
}
check_alloc(outstring=(char*)malloc(sizeof(char)*(outstrlen+8)));
sprintf(outstring,"#Prog: ");
offset=7;
for (i=0;i<argc;i++) {
sprintf(outstring+offset,"%s ",argv[i]);
offset += (strlen(argv[i])+1);
}
scan_options(argc,argv);
infile=search_datafile(argc,argv,NULL,VERBOSITY);
if (infile == NULL) {
fprintf(stderr,"No input datafile found.\n");
exit(LANGEVIN_MAIN_NO_INPUTFILE);
}
if (reference) {
fout=fopen(refname,"r");
if (fout == NULL) {
fprintf(stderr,"No test point datafile found.\n");
exit(LANGEVIN_MAIN_NO_INPUTFILE);
}
else
fclose(fout);
}
if (OUTFILE == NULL) {
check_alloc(OUTFILE=(char*)calloc(strlen(infile)+6,(size_t)1));
sprintf(OUTFILE,"%s.osn",infile);
}
OUTFILE=test_outfile(OUTFILE);
if (!reference) {
refname=infile;
ignore += EXCLUDE;
}
if (COLUMN == NULL)
series=(double**)get_multi_series(infile,&(pars.LENGTH),EXCLUDE,
&(pars.DIM),"",dimset,VERBOSITY);
else
series=(double**)get_multi_series(infile,&(pars.LENGTH),EXCLUDE,
&(pars.DIM),COLUMN,dimset,VERBOSITY);
switch (neighborkind) {
case 0: {hdim=0;minnadd=1;}break;
case 1: {hdim=0;minnadd=2;}break;
case 2: {hdim=1;minnadd=3;}
}
if (!reference && (ignore < hdim) && (neighborkind == 2)) {
ignore=hdim;
}
pars.hdim=hdim;
if (COLUMN == NULL)
rseries=(double**)get_multi_series(refname,&FLENGTH,ignore,
&(pars.DIM),"",dimset,VERBOSITY);
else
rseries=(double**)get_multi_series(refname,&FLENGTH,ignore,
&(pars.DIM),COLUMN,dimset,VERBOSITY);
if (VERBOSITY)
fprintf(stderr,"%lf\n",rseries[0][0]);
check_alloc(sermin=malloc(sizeof(double)*pars.DIM));
check_alloc(serinter=malloc(sizeof(double)*pars.DIM));
rescale_data(series,pars,sermin,serinter);
max=serinter[0];
for (i=1;i<pars.DIM;i++)
if (serinter[i] > max)
max=serinter[i];
for (i=0;i<pars.DIM;i++) {
for (j=0;j<pars.LENGTH;j++)
series[i][j] /= max;
for (j=0;j<FLENGTH;j++)
rseries[i][j]=(rseries[i][j]-sermin[i])/max;
}
if (pars.maxr > 0.0)
pars.maxr /= max;
check_alloc(future=(unsigned int*)malloc(sizeof(int)*pars.LENGTH));
if (neighborkind == 0) {
for (li=0;li<pars.LENGTH;li++)
future[li]=1;
}
else {
if (WHICHFUTURE>0) {
check_alloc(wcol=(char*)calloc((size_t)10,(size_t)1));
sprintf(wcol,"%u",WHICHFUTURE);
dfuture=(double**)get_multi_series(infile,&(pars.LENGTH),EXCLUDE,
&wdim,wcol,1,0L);
if (neighborkind == 2) {
for (li=0;li<hdim;li++)
future[li]=0;
}
for (li=hdim;li<pars.LENGTH;li++) {
has_future= (dfuture[0][li]>0.0);
for (lj=1;lj<=hdim;lj++)
has_future &= (dfuture[0][li-lj]>0.0);
future[li]=has_future;
}
free(dfuture);
}
else {
for (li=0;li<pars.LENGTH-1;li++)
future[li]=1;
future[pars.LENGTH-1]=0;
if (neighborkind == 2) {
for (li=0;li<hdim;li++)
future[li]=0;
}
}
}
check_alloc(sf.found=(unsigned long*)
malloc(sizeof(unsigned long)*pars.LENGTH));
check_alloc(sf.distance=(double*)malloc(sizeof(double)*pars.LENGTH));
check_alloc(sf.weight=(double*)malloc(sizeof(double)*pars.LENGTH));
check_alloc(sf.count=(unsigned long*)malloc(sizeof(unsigned long)));
check_alloc(sf.aveps=(double*)malloc(sizeof(double)));
sf.count[0]=0;
sf.aveps[0]=0.0;
check_alloc(cast=(double**)malloc(sizeof(double*)*pars.DIM));
for (i=0;i<pars.DIM;i++)
check_alloc(cast[i]=(double*)malloc(sizeof(double)*(hdim+1)));
pars.minminn=1;
init_neighbor_search(series,pars,future);
/* print program call and column labels */
fout=fopen(OUTFILE,"w");
fprintf(fout,"%s\n",outstring);
#if defined(MAXNORM)
fprintf(fout,"#Norm: Maxnorm\n");
#else
fprintf(fout,"#Norm: L2 Norm\n");
#endif
#if defined(WEIGHTS)
if (pars.maxr > 0.0)
fprintf(fout,"#Weights: on WFACT = %lf maxr = %lf\n",WFACT,pars.maxr*max);
else
fprintf(fout,"#Weights: on WFACT = %lf maxr not set\n",WFACT);
#else
fprintf(fout,"#Weights: off\n");
#endif
fprintf(fout,"#Content: ");
for (i=0;i<pars.DIM;i++) {
fprintf(fout,"x%d ",i+1);
}
for (i=0;i<pars.MINN;i++)
fprintf(fout,"n%d ",i+1);
fprintf(fout,"\n");
fflush(fout);
n=0;
while (n < FLENGTH) {
refn=n+ignore-EXCLUDE;
if (!reference) {
while (!future[refn]) {
n++;
refn++;
if (n>=(pars.LENGTH-1))
goto ntoolarge;
}
}
for (j=0;j<pars.DIM;j++)
cast[j][hdim]=rseries[j][n];
/* search neighbors */
if (!reference)
pars.MINN += minnadd;
search_neighbors(series,cast,pars,sf);
if (!reference) {
j=0;
switch (neighborkind) {
case 0: {
for (i=0;i<pars.MINN;i++) {
if (sf.found[i] == refn) {
j++;
}
sf.found[i]=sf.found[j];
sf.distance[i]=sf.distance[j];
j++;
}
} break;
case 1: {
for (i=0;i<pars.MINN;i++) {
if ((sf.found[i] == refn) || (sf.found[i] == (refn-1))) {
j++;
}
sf.found[i]=sf.found[j];
sf.distance[i]=sf.distance[j];
j++;
}
} break;
case 2: {
for (i=0;i<pars.MINN;i++) {
if ((sf.found[i] == refn) || (sf.found[i] == (refn-1)) ||
(sf.found[i] == (refn+1))) {
j++;
}
sf.found[i]=sf.found[j];
sf.distance[i]=sf.distance[j];
j++;
}
} break;
}
pars.MINN -= minnadd;
}
/* print coordinate */
for (i=0;i<pars.DIM;i++)
fprintf(fout,"%e ",cast[i][hdim]*max+sermin[i]);
/* print neighbors */
for (i=0;i<pars.MINN;i++)
fprintf(fout,"%lu ",sf.found[i]+1);
fprintf(fout,"\n");
fflush(fout);
n++;
}
ntoolarge: fclose(fout);
return 0;
}
int main(int argc,char **argv)
{
char stdi=0;
unsigned long i,j;
double x,norm,size=1.0,size2=1.0;
FILE *fout;
if (scan_help(argc,argv))
show_options(argv[0]);
scan_options(argc,argv);
#ifndef OMIT_WHAT_I_DO
if (verbosity&VER_INPUT)
what_i_do(argv[0],WID_STR);
#endif
infile=search_datafile(argc,argv,&column,verbosity);
if (infile == NULL)
stdi=1;
if (outfile == NULL) {
if (!stdi) {
check_alloc(outfile=(char*)calloc(strlen(infile)+5,1));
strcpy(outfile,infile);
strcat(outfile,".his");
}
else {
check_alloc(outfile=(char*)calloc((size_t)10,1));
strcpy(outfile,"stdin.his");
}
}
if (!my_stdout)
test_outfile(outfile);
series=(double*)get_series(infile,&length,exclude,column,verbosity);
variance(series,length,&average,&var);
rescale_data(series,length,&min,&max);
if (base > 0) {
check_alloc(box=(long*)malloc(sizeof(long)*base));
for (i=0;i<base;i++)
box[i]=0;
size=1./base;
size2=size/2.0;
for (i=0;i<length;i++) {
if (series[i] > (1.0-size2))
series[i]=1.0-size2;
j=(long)(series[i]*base);
box[j]++;
}
}
norm=1.0/(double)length;
if (!my_stdout) {
fout=fopen(outfile,"w");
if (verbosity&VER_INPUT)
fprintf(stderr,"Opened %s for writing\n",outfile);
fprintf(fout,"#interval of data: [%e:%e]\n",min,max+min);
fprintf(fout,"#average= %e\n",average);
fprintf(fout,"#standard deviation= %e\n",var);
for (i=0;i<base;i++) {
x=(double)(i*size);
fprintf(fout,"%e %e\n",(x+size2)*max+min,(double)box[i]*norm);
}
fclose(fout);
}
else {
if (verbosity&VER_INPUT)
fprintf(stderr,"Writing to stdout\n");
fprintf(stdout,"#interval of data: [%e:%e]\n",min,max+min);
fprintf(stdout,"#average= %e\n",average);
fprintf(stdout,"#standard deviation= %e\n",var);
for (i=0;i<base;i++) {
x=(double)(i*size);
fprintf(stdout,"%e %e\n",(x+size2)*max+min,(double)box[i]*norm);
fflush(stdout);
}
}
return 0;
}
