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 <cairo/cairo.h>
  23 #include <math.h>
  24 #include <assert.h>
  25
  26 #include <output/chart.h>
  27 #include <output/chart-provider.h>
  28 #include <output/charts/plot-chart.h>
  29 #include <libpspp/compiler.h>
  30
  31 #include "xalloc.h"
  32
  33 /* Start a new vector called NAME */
  34 void
  35 chart_vector_start (cairo_t *cr, struct chart_geometry *geom, const char *name)
  36 {
  37   const struct chart_colour *colour;
  38
  39   cairo_save (cr);
  40
  41   colour = &data_colour[geom->n_datasets % N_CHART_COLOURS];
  42   cairo_set_source_rgb (cr,
  43                         colour->red / 255.0,
  44                         colour->green / 255.0,
  45                         colour->blue / 255.0);
  46
  47   geom->n_datasets++;
  48   geom->dataset = xrealloc (geom->dataset,
  49                             geom->n_datasets * sizeof (*geom->dataset));
  50
  51   geom->dataset[geom->n_datasets - 1] = strdup (name);
  52 }
  53
  54 /* Plot a data point */
  55 void
  56 chart_datum (cairo_t *cr, const struct chart_geometry *geom,
  57              int dataset UNUSED, double x, double y)
  58 {
  59   double x_pos = (x - geom->x_min) * geom->abscissa_scale + geom->data_left;
  60   double y_pos = (y - geom->y_min) * geom->ordinate_scale + geom->data_bottom;
  61
  62   chart_draw_marker (cr, x_pos, y_pos, MARKER_SQUARE, 15);
  63 }
  64
  65 void
  66 chart_vector_end (cairo_t *cr, struct chart_geometry *geom)
  67 {
  68   cairo_stroke (cr);
  69   cairo_restore (cr);
  70   geom->in_path = false;
  71 }
  72
  73 /* Plot a data point */
  74 void
  75 chart_vector (cairo_t *cr, struct chart_geometry *geom, double x, double y)
  76 {
  77   const double x_pos =
  78     (x - geom->x_min) * geom->abscissa_scale + geom->data_left ;
  79
  80   const double y_pos =
  81     (y - geom->y_min) * geom->ordinate_scale + geom->data_bottom ;
  82
  83   if (geom->in_path)
  84     cairo_line_to (cr, x_pos, y_pos);
  85   else
  86     {
  87       cairo_move_to (cr, x_pos, y_pos);
  88       geom->in_path = true;
  89     }
  90 }
  91
  92
  93
  94 /* Draw a line with slope SLOPE and intercept INTERCEPT.
  95    between the points limit1 and limit2.
  96    If lim_dim is CHART_DIM_Y then the limit{1,2} are on the
  97    y axis otherwise the x axis
  98 */
  99 void
 100 chart_line(cairo_t *cr, const struct chart_geometry *geom,
 101            double slope, double intercept,
 102            double limit1, double limit2, enum CHART_DIM lim_dim)
 103 {
 104   double x1, y1;
 105   double x2, y2;
 106
 107   if ( lim_dim == CHART_DIM_Y )
 108     {
 109       x1 = ( limit1 - intercept ) / slope;
 110       x2 = ( limit2 - intercept ) / slope;
 111       y1 = limit1;
 112       y2 = limit2;
 113     }
 114   else
 115     {
 116       x1 = limit1;
 117       x2 = limit2;
 118       y1 = slope * x1 + intercept;
 119       y2 = slope * x2 + intercept;
 120     }
 121
 122   y1 = (y1 - geom->y_min) * geom->ordinate_scale + geom->data_bottom;
 123   y2 = (y2 - geom->y_min) * geom->ordinate_scale + geom->data_bottom;
 124   x1 = (x1 - geom->x_min) * geom->abscissa_scale + geom->data_left;
 125   x2 = (x2 - geom->x_min) * geom->abscissa_scale + geom->data_left;
 126
 127   cairo_move_to (cr, x1, y1);
 128   cairo_line_to (cr, x2, y2);
 129   cairo_stroke (cr);
 130 }