1 /* PSPP - computes sample statistics.
2 Copyright (C) 2004 Free Software Foundation, Inc.
3 Written by John Darrington <john@darrington.wattle.id.au>
5 This program is free software; you can redistribute it and/or
6 modify it under the terms of the GNU General Public License as
7 published by the Free Software Foundation; either version 2 of the
8 License, or (at your option) any later version.
10 This program is distributed in the hope that it will be useful, but
11 WITHOUT ANY WARRANTY; without even the implied warranty of
12 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
13 General Public License for more details.
15 You should have received a copy of the GNU General Public License
16 along with this program; if not, write to the Free Software
17 Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
24 #include <output/chart.h>
26 #include <output/charts/plot-chart.h>
27 #include <output/charts/cartesian.h>
28 #include <libpspp/compiler.h>
41 static const struct dataset dataset[DATASETS] =
48 /* Plot a data point */
50 chart_datum(struct chart *ch, int dataset UNUSED, double x, double y)
57 (x - ch->x_min) * ch->abscissa_scale + ch->data_left ;
60 (y - ch->y_min) * ch->ordinate_scale + ch->data_bottom ;
62 pl_savestate_r(ch->lp);
64 pl_fmarker_r(ch->lp, x_pos, y_pos, 6, 15);
66 pl_restorestate_r(ch->lp);
70 /* Draw a line with slope SLOPE and intercept INTERCEPT.
71 between the points limit1 and limit2.
72 If lim_dim is CHART_DIM_Y then the limit{1,2} are on the
73 y axis otherwise the x axis
76 chart_line(struct chart *ch, double slope, double intercept,
77 double limit1, double limit2, enum CHART_DIM lim_dim)
86 if ( lim_dim == CHART_DIM_Y )
88 x1 = ( limit1 - intercept ) / slope ;
89 x2 = ( limit2 - intercept ) / slope ;
97 y1 = slope * x1 + intercept;
98 y2 = slope * x2 + intercept;
101 y1 = (y1 - ch->y_min) * ch->ordinate_scale + ch->data_bottom ;
102 y2 = (y2 - ch->y_min) * ch->ordinate_scale + ch->data_bottom ;
103 x1 = (x1 - ch->x_min) * ch->abscissa_scale + ch->data_left ;
104 x2 = (x2 - ch->x_min) * ch->abscissa_scale + ch->data_left ;
106 pl_savestate_r(ch->lp);
108 pl_fline_r(ch->lp, x1, y1, x2, y2);
110 pl_restorestate_r(ch->lp);