pintos-os.org Git - pspp/blob - src/output/charts/cartesian.c

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