Recent spread-sheet-widget to version 0.2
[pspp] / src / output / charts / spreadlevel-plot.c
1 /* PSPP - a program for statistical analysis.
2    Copyright (C) 2012 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
19 #include "output/charts/spreadlevel-plot.h"
20
21 #include "libpspp/cast.h"
22 #include "output/chart-item-provider.h"
23
24 #include "gl/xalloc.h"
25 #include "gl/minmax.h"
26
27 #include <math.h>
28 #include <float.h>
29 #include <stdlib.h>
30
31 struct chart_item *
32 spreadlevel_plot_create (const char *label, double tx_pwr)
33 {
34   struct spreadlevel_plot_chart *sl = xzalloc (sizeof *sl);
35   chart_item_init (&sl->chart_item, &spreadlevel_plot_chart_class, label);
36
37   sl->x_lower = DBL_MAX;
38   sl->x_upper = -DBL_MAX;
39
40   sl->y_lower = DBL_MAX;
41   sl->y_upper = -DBL_MAX;
42
43   sl->tx_pwr = tx_pwr;
44
45   sl->n_data = 0;
46   sl->data = NULL;
47
48   return &sl->chart_item;
49 }
50
51 void
52 spreadlevel_plot_add (struct chart_item *ci, double spread, double level)
53 {
54   struct spreadlevel_plot_chart *sl = to_spreadlevel_plot_chart (ci);
55
56   if ( sl->tx_pwr == 0)
57     {
58       spread = log (spread);
59       level = log (level);
60     }
61   else
62     {
63       spread = pow (spread, sl->tx_pwr);
64       level = pow (level, sl->tx_pwr);
65     }
66
67   sl->x_lower = MIN (sl->x_lower, level);
68   sl->x_upper = MAX (sl->x_upper, level);
69
70   sl->y_lower = MIN (sl->y_lower, spread);
71   sl->y_upper = MAX (sl->y_upper, spread);
72
73   sl->n_data++;
74   sl->data = xrealloc (sl->data, sizeof (*sl->data) * sl->n_data);
75   sl->data[sl->n_data - 1].x = level;
76   sl->data[sl->n_data - 1].y = spread;
77 }
78
79
80 static void
81 spreadlevel_plot_chart_destroy (struct chart_item *chart_item)
82 {
83   struct spreadlevel_plot_chart *sl = to_spreadlevel_plot_chart (chart_item);
84
85   free (sl->data);
86   free (sl);
87 }
88
89 const struct chart_item_class spreadlevel_plot_chart_class =
90   {
91     spreadlevel_plot_chart_destroy
92   };