IUPred wrapper, tester and new binary for X64 linux systems
[jabaws.git] / binaries / src / disembl / Tisean_3.0.1 / source_c / rescale.c
1 /*
2  *   This file is part of TISEAN
3  *
4  *   Copyright (c) 1998-2007 Rainer Hegger, Holger Kantz, Thomas Schreiber
5  *
6  *   TISEAN is free software; you can redistribute it and/or modify
7  *   it under the terms of the GNU General Public License as published by
8  *   the Free Software Foundation; either version 2 of the License, or
9  *   (at your option) any later version.
10  *
11  *   TISEAN is distributed in the hope that it will be useful,
12  *   but WITHOUT ANY WARRANTY; without even the implied warranty of
13  *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
14  *   GNU General Public License for more details.
15  *
16  *   You should have received a copy of the GNU General Public License
17  *   along with TISEAN; if not, write to the Free Software
18  *   Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
19  */
20 /*Author: Rainer Hegger. Last modified: Nov 23, 2000 */
21 #include <stdio.h>
22 #include <stdlib.h>
23 #include <string.h>
24 #include <limits.h>
25 #include <math.h>
26 #include "routines/tsa.h"
27
28 #define WID_STR "Rescales the data"
29
30 unsigned long length=ULONG_MAX,exclude=0;
31 unsigned int dim=1;
32 unsigned int verbosity=0xff;
33 char *column=NULL;
34 char *outfile=NULL,stdo=1,set_av=0,set_var=0,dimset=0;
35 char *infile=NULL;
36 double **series;
37 double xmin=0.0,xmax=1.0;
38
39 void show_options(char *progname)
40 {
41   what_i_do(progname,WID_STR);
42   fprintf(stderr," Usage: %s [options]\n",progname);
43   fprintf(stderr," Options:\n");
44   fprintf(stderr,"Everything not being a valid option will be interpreted"
45           " as a possible"
46           " datafile.\nIf no datafile is given stdin is read. Just - also"
47           " means stdin\n");
48   fprintf(stderr,"\t-l # of data to use [default: whole file]\n");
49   fprintf(stderr,"\t-x # of lines to ignore [default: 0]\n");
50   fprintf(stderr,"\t-m # of components to be read [default: %u]\n",dim);
51   fprintf(stderr,"\t-c columns to read [default: 1,...,# of components]\n");
52   fprintf(stderr,"\t-z minimum of the new series [default: 0.0]\n");
53   fprintf(stderr,"\t-Z maximum of the new series [default: 1.0]\n");
54   fprintf(stderr,"\t-a create a series with average value equals 0\n");
55   fprintf(stderr,"\t-v create a series with variance 1\n");
56   fprintf(stderr,"\t-o output file name [default: 'datafile'.res']\n");
57   fprintf(stderr,"\t-V verbosity level [default: 1]\n\t\t"
58           "0='only panic messages'\n\t\t"
59           "1='+ input/output messages'\n");
60   fprintf(stderr,"\t-h show these options\n");
61   exit(0);
62 }
63
64 void scan_options(int n,char **in)
65 {
66   char *out;
67
68   if ((out=check_option(in,n,'l','u')) != NULL)
69     sscanf(out,"%lu",&length);
70   if ((out=check_option(in,n,'x','u')) != NULL)
71     sscanf(out,"%lu",&exclude);
72   if ((out=check_option(in,n,'m','u')) != NULL) {
73     sscanf(out,"%u",&dim);
74     dimset=1;
75   }
76   if ((out=check_option(in,n,'c','s')) != NULL)
77     column=out;
78   if ((out=check_option(in,n,'V','u')) != NULL)
79     sscanf(out,"%u",&verbosity);
80   if ((out=check_option(in,n,'z','f')) != NULL)
81     sscanf(out,"%lf",&xmin);
82   if ((out=check_option(in,n,'Z','f')) != NULL)
83     sscanf(out,"%lf",&xmax);
84   if ((out=check_option(in,n,'a','n')) != NULL)
85     set_av=1;
86   if ((out=check_option(in,n,'v','n')) != NULL)
87     set_var=1;
88   if ((out=check_option(in,n,'o','o')) != NULL) {
89     stdo=0;
90     if (strlen(out) > 0)
91       outfile=out;
92   }
93 }
94
95 int main(int argc,char **argv)
96 {
97   char stdi=0;
98   FILE *file;
99   double min,max;
100   double av,varianz;
101   long i,n;
102     
103   if (scan_help(argc,argv))
104     show_options(argv[0]);
105
106   scan_options(argc,argv);
107 #ifndef OMIT_WHAT_I_DO
108   if (verbosity&VER_INPUT)
109     what_i_do(argv[0],WID_STR);
110 #endif
111
112   infile=search_datafile(argc,argv,NULL,verbosity);
113   if (infile == NULL)
114     stdi=1;
115
116   if (outfile == NULL) {
117     if (!stdi) {
118       check_alloc(outfile=(char*)calloc(strlen(infile)+5,(size_t)1));
119       strcpy(outfile,infile);
120       strcat(outfile,".res");
121     }
122     else {
123       check_alloc(outfile=(char*)calloc((size_t)10,(size_t)1));
124       strcpy(outfile,"stdin.res");
125     }
126   }
127   if (!stdo)
128     test_outfile(outfile);
129
130   if (xmin >= xmax) {
131     fprintf(stderr,"Choosing the minimum larger or equal the maximum\n"
132             "makes no sense. Exiting!\n");
133     exit(RESCALE__WRONG_INTERVAL);
134   }
135
136   if (column == NULL)
137     series=(double**)get_multi_series(infile,&length,exclude,&dim,"",dimset,
138                                       verbosity);
139   else
140     series=(double**)get_multi_series(infile,&length,exclude,&dim,column,
141                                       dimset,verbosity);
142
143   for (n=0;n<dim;n++) {
144     variance(series[n],length,&av,&varianz);
145     
146     if (set_av)
147       for (i=0;i<length;i++)
148         series[n][i] -= av;
149
150     if (set_var)
151       for (i=0;i<length;i++)
152         series[n][i] /= varianz;
153   
154     if (!set_var && !set_av) {
155       rescale_data(series[n],length,&min,&max);
156       for (i=0;i<length;i++)
157         series[n][i]=series[n][i]*(xmax-xmin)+xmin;
158     }
159   }
160
161   if (stdo) {
162     if (verbosity&VER_INPUT)
163       fprintf(stderr,"Writing to stdout\n");
164     for (i=0;i<length;i++) {
165       fprintf(stdout,"%e",series[0][i]);
166       for (n=1;n<dim;n++)
167         fprintf(stdout," %e",series[n][i]);
168       fprintf(stdout,"\n");
169     }
170   }
171   else {
172     file=fopen(outfile,"w");
173     if (verbosity&VER_INPUT)
174       fprintf(stderr,"Opened %s for writing\n",outfile);
175     for (i=0;i<length;i++) {
176       fprintf(file,"%e",series[0][i]);
177       for (n=1;n<dim;n++)
178         fprintf(file," %e",series[n][i]);
179       fprintf(file,"\n");
180     }
181     fclose(file);
182   }
183
184   return 0;
185 }