output: Add initial support for PSPP output via Cairo.
[pspp] / src / output / charts / cartesian.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
18 #include <config.h>
19
20 #include <math.h>
21 #include <assert.h>
22
23 #include <output/chart.h>
24
25 #include <output/charts/plot-chart.h>
26 #include <output/charts/cartesian.h>
27 #include <libpspp/compiler.h>
28
29
30
31 struct dataset
32 {
33   int n_data;
34   const char *label;
35 };
36
37
38 #define DATASETS 2
39
40 static const struct dataset dataset[DATASETS] =
41   {
42     { 13, "male"},
43     { 11, "female"},
44   };
45
46
47 /* Plot a data point */
48 void
49 chart_datum(struct chart *ch, int dataset UNUSED, double x, double y)
50 {
51   if ( ! ch )
52     return ;
53
54   {
55     const double x_pos =
56       (x - ch->x_min) * ch->abscissa_scale + ch->data_left ;
57
58     const double y_pos =
59       (y - ch->y_min) * ch->ordinate_scale + ch->data_bottom ;
60
61     pl_savestate_r(ch->lp);
62
63     pl_fmarker_r(ch->lp, x_pos, y_pos, 6, 15);
64
65     pl_restorestate_r(ch->lp);
66   }
67 }
68
69 /* Draw a line with slope SLOPE and intercept INTERCEPT.
70    between the points limit1 and limit2.
71    If lim_dim is CHART_DIM_Y then the limit{1,2} are on the
72    y axis otherwise the x axis
73 */
74 void
75 chart_line(struct chart *ch, double slope, double intercept,
76            double limit1, double limit2, enum CHART_DIM lim_dim)
77 {
78   double x1, y1;
79   double x2, y2 ;
80
81   if ( ! ch )
82     return ;
83
84
85   if ( lim_dim == CHART_DIM_Y )
86     {
87       x1 = ( limit1 - intercept ) / slope ;
88       x2 = ( limit2 - intercept ) / slope ;
89       y1 = limit1;
90       y2 = limit2;
91     }
92   else
93     {
94       x1 = limit1;
95       x2 = limit2;
96       y1 = slope * x1 + intercept;
97       y2 = slope * x2 + intercept;
98     }
99
100   y1 = (y1 - ch->y_min) * ch->ordinate_scale + ch->data_bottom ;
101   y2 = (y2 - ch->y_min) * ch->ordinate_scale + ch->data_bottom ;
102   x1 = (x1 - ch->x_min) * ch->abscissa_scale + ch->data_left ;
103   x2 = (x2 - ch->x_min) * ch->abscissa_scale + ch->data_left ;
104
105   pl_savestate_r(ch->lp);
106
107   pl_fline_r(ch->lp, x1, y1, x2, y2);
108
109   pl_restorestate_r(ch->lp);
110
111 }