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