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