Adopt use of gnulib for portability.
[pspp] / src / chart.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.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