src/output/charts/cartesian.c

   1 /* PSPP - computes sample statistics.
   2    Copyright (C) 2004 Free Software Foundation, Inc.
   3    Written by John Darrington <john@darrington.wattle.id.au>
   4
   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.
   9
  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.
  14
  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
  18    02110-1301, USA. */
  19
  20
  21 #include <math.h>
  22 #include <assert.h>
  23
  24 #include <output/chart.h>
  25
  26 #include <output/charts/plot-chart.h>
  27 #include <output/charts/cartesian.h>
  28 #include <libpspp/compiler.h>
  29
  30
  31
  32 struct dataset
  33 {
  34   int n_data;
  35   const char *label;
  36 };
  37
  38
  39 #define DATASETS 2
  40
  41 static const struct dataset dataset[DATASETS] =
  42   {
  43     { 13, "male"},
  44     { 11, "female"},
  45   };
  46
  47
  48 /* Plot a data point */
  49 void
  50 chart_datum(struct chart *ch, int dataset UNUSED, double x, double y)
  51 {
  52   if ( ! ch )
  53     return ;
  54
  55   {
  56     const double x_pos =
  57       (x - ch->x_min) * ch->abscissa_scale + ch->data_left ;
  58
  59     const double y_pos =
  60       (y - ch->y_min) * ch->ordinate_scale + ch->data_bottom ;
  61
  62     pl_savestate_r(ch->lp);
  63
  64     pl_fmarker_r(ch->lp, x_pos, y_pos, 6, 15);
  65
  66     pl_restorestate_r(ch->lp);
  67   }
  68 }
  69
  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
  74 */
  75 void
  76 chart_line(struct chart *ch, double slope, double intercept,
  77            double limit1, double limit2, enum CHART_DIM lim_dim)
  78 {
  79   double x1, y1;
  80   double x2, y2 ;
  81
  82   if ( ! ch )
  83     return ;
  84
  85
  86   if ( lim_dim == CHART_DIM_Y )
  87     {
  88       x1 = ( limit1 - intercept ) / slope ;
  89       x2 = ( limit2 - intercept ) / slope ;
  90       y1 = limit1;
  91       y2 = limit2;
  92     }
  93   else
  94     {
  95       x1 = limit1;
  96       x2 = limit2;
  97       y1 = slope * x1 + intercept;
  98       y2 = slope * x2 + intercept;
  99     }
 100
 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 ;
 105
 106   pl_savestate_r(ch->lp);
 107
 108   pl_fline_r(ch->lp, x1, y1, x2, y2);
 109
 110   pl_restorestate_r(ch->lp);
 111
 112 }