src/math/chart-geometry.c

   1 /* PSPP - computes sample statistics.
   2    Copyright (C) 2004 Free Software Foundation, Inc.
   3    Written by John Darrington <john@darrington.wattle.id.au>
   4
   5    This program is free software; you can redistribute it and/or
   6    modify it under the terms of the GNU General Public License as
   7    published by the Free Software Foundation; either version 2 of the
   8    License, or (at your option) any later version.
   9
  10    This program is distributed in the hope that it will be useful, but
  11    WITHOUT ANY WARRANTY; without even the implied warranty of
  12    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  13    General Public License for more details.
  14
  15    You should have received a copy of the GNU General Public License
  16    along with this program; if not, write to the Free Software
  17    Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
  18    02110-1301, USA. */
  19
  20
  21 #include <math.h>
  22 #include <float.h>
  23
  24 #include "chart-geometry.h"
  25
  26 /* Adjust tick to be a sensible value
  27    ie:  ... 0.1,0.2,0.5,   1,2,5,  10,20,50 ... */
  28 double
  29 chart_rounded_tick(double tick)
  30 {
  31
  32   int i;
  33
  34   double diff = DBL_MAX;
  35   double t = tick;
  36
  37   static const double standard_ticks[] = {1, 2, 5, 10};
  38
  39   const double factor = pow(10,ceil(log10(standard_ticks[0] / tick))) ;
  40
  41   for (i = 3  ; i >= 0 ; --i)
  42     {
  43       const double d = fabs( tick - standard_ticks[i] / factor ) ;
  44
  45       if ( d < diff )
  46         {
  47           diff = d;
  48           t = standard_ticks[i] / factor ;
  49         }
  50     }
  51
  52   return t;
  53
  54 }
  55