website/archive/binaries/mac/src/disembl/Tisean_3.0.1/source_c/routines/rescale_data.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: Sep 5, 2004*/
  21 /* Changes:
  22  * Sep 5 2004: + include <stdlib.h>
  23  */
  24
  25 #include <stdio.h>
  26 #include "tisean_cec.h"
  27 #include <stdlib.h>
  28
  29 void rescale_data(double *x,unsigned long l,double *min,double *interval)
  30 {
  31   int i;
  32
  33   *min=*interval=x[0];
  34
  35   for (i=1;i<l;i++) {
  36     if (x[i] < *min) *min=x[i];
  37     if (x[i] > *interval) *interval=x[i];
  38   }
  39   *interval -= *min;
  40
  41   if (*interval != 0.0) {
  42     for (i=0;i<l;i++)
  43       x[i]=(x[i]- *min)/ *interval;
  44   }
  45   else {
  46     fprintf(stderr,"rescale_data: data ranges from %e to %e. It makes\n"
  47             "\t\tno sense to continue. Exiting!\n\n",*min,*min+(*interval));
  48     exit(RESCALE_DATA_ZERO_INTERVAL);
  49   }
  50 }