pintos-os.org Git - pspp/blob - src/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 "chart.h"
  23 #include <assert.h>
  24
  25
  26
  27 struct dataset
  28 {
  29   int n_data;
  30   const char *label;
  31 };
  32
  33
  34
  35
  36 #define DATASETS 2
  37
  38 static const struct dataset dataset[DATASETS] =
  39   {
  40     { 13, "male"},
  41     { 11, "female"},
  42   };
  43
  44
  45
  46
  47 static void
  48 write_legend(struct chart *chart, const char *heading, int n);
  49
  50
  51
  52 /* Write the abscissa label */
  53 void
  54 chart_write_xlabel(struct chart *ch, const char *label)
  55 {
  56   if ( ! ch )
  57     return ;
  58
  59   pl_savestate_r(ch->lp);
  60
  61   pl_move_r(ch->lp,ch->data_left, ch->abscissa_top);
  62   pl_alabel_r(ch->lp,0,'t',label);
  63
  64   pl_restorestate_r(ch->lp);
  65
  66 }
  67
  68
  69
  70 /* Write the ordinate label */
  71 void
  72 chart_write_ylabel(struct chart *ch, const char *label)
  73 {
  74   if ( ! ch )
  75     return ;
  76
  77   pl_savestate_r(ch->lp);
  78
  79   pl_move_r(ch->lp, ch->data_bottom, ch->ordinate_right);
  80   pl_textangle_r(ch->lp, 90);
  81   pl_alabel_r(ch->lp, 0, 0, label);
  82
  83   pl_restorestate_r(ch->lp);
  84 }
  85
  86
  87
  88 static void
  89 write_legend(struct chart *chart, const char *heading,
  90              int n)
  91 {
  92   int ds;
  93
  94   if ( ! chart )
  95     return ;
  96
  97
  98   pl_savestate_r(chart->lp);
  99
 100   pl_filltype_r(chart->lp,1);
 101
 102   pl_move_r(chart->lp, chart->legend_left,
 103             chart->data_bottom + chart->font_size * n * 1.5);
 104
 105   pl_alabel_r(chart->lp,0,'b',heading);
 106
 107   for (ds = 0 ; ds < n ; ++ds )
 108     {
 109       pl_fmove_r(chart->lp,
 110                  chart->legend_left,
 111                  chart->data_bottom + chart->font_size * ds  * 1.5);
 112
 113       pl_savestate_r(chart->lp);
 114       pl_fillcolorname_r(chart->lp,data_colour[ds]);
 115       pl_fboxrel_r (chart->lp,
 116                     0,0,
 117                     chart->font_size, chart->font_size);
 118       pl_restorestate_r(chart->lp);
 119
 120       pl_fmove_r(chart->lp,
 121                  chart->legend_left + chart->font_size * 1.5,
 122                  chart->data_bottom + chart->font_size * ds  * 1.5);
 123
 124       pl_alabel_r(chart->lp,'l','b',dataset[ds].label);
 125     }
 126
 127
 128   pl_restorestate_r(chart->lp);
 129 }
 130
 131
 132 /* Plot a data point */
 133 void
 134 chart_datum(struct chart *ch, int dataset UNUSED, double x, double y)
 135 {
 136   if ( ! ch )
 137     return ;
 138
 139   {
 140     const double x_pos =
 141       (x - ch->x_min) * ch->abscissa_scale + ch->data_left ;
 142
 143     const double y_pos =
 144       (y - ch->y_min) * ch->ordinate_scale + ch->data_bottom ;
 145
 146     pl_savestate_r(ch->lp);
 147
 148     pl_fmarker_r(ch->lp, x_pos, y_pos, 6, 15);
 149
 150     pl_restorestate_r(ch->lp);
 151   }
 152 }
 153
 154 /* Draw a line with slope SLOPE and intercept INTERCEPT.
 155    between the points limit1 and limit2.
 156    If lim_dim is CHART_DIM_Y then the limit{1,2} are on the
 157    y axis otherwise the x axis
 158 */
 159 void
 160 chart_line(struct chart *ch, double slope, double intercept,
 161            double limit1, double limit2, enum CHART_DIM lim_dim)
 162 {
 163   double x1, y1;
 164   double x2, y2 ;
 165
 166   if ( ! ch )
 167     return ;
 168
 169
 170   if ( lim_dim == CHART_DIM_Y )
 171     {
 172       x1 = ( limit1 - intercept ) / slope ;
 173       x2 = ( limit2 - intercept ) / slope ;
 174       y1 = limit1;
 175       y2 = limit2;
 176     }
 177   else
 178     {
 179       x1 = limit1;
 180       x2 = limit2;
 181       y1 = slope * x1 + intercept;
 182       y2 = slope * x2 + intercept;
 183     }
 184
 185   y1 = (y1 - ch->y_min) * ch->ordinate_scale + ch->data_bottom ;
 186   y2 = (y2 - ch->y_min) * ch->ordinate_scale + ch->data_bottom ;
 187   x1 = (x1 - ch->x_min) * ch->abscissa_scale + ch->data_left ;
 188   x2 = (x2 - ch->x_min) * ch->abscissa_scale + ch->data_left ;
 189
 190   pl_savestate_r(ch->lp);
 191
 192   pl_fline_r(ch->lp, x1, y1, x2, y2);
 193
 194   pl_restorestate_r(ch->lp);
 195
 196 }