eb0ad9282a854885768472319a899dbfba6c14ea
[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 static const double standard_ticks[] = {1, 2, 5, 10};
24
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   int i;
32
33   double diff = DBL_MAX;
34   double t = tick;
35
36   double factor;
37
38   /* Avoid arithmetic problems with very small values */
39   if (fabs (tick) < DBL_EPSILON)
40      return 0;
41
42   factor = pow (10,ceil (log10 (standard_ticks[0] / tick)));
43
44   for (i = 3  ; i >= 0 ; --i)
45     {
46       const double d = fabs (tick - standard_ticks[i] / factor);
47
48       if ( d < diff )
49         {
50           diff = d;
51           t = standard_ticks[i] / factor ;
52         }
53     }
54
55   return t;
56 }
57