Recent spread-sheet-widget to version 0.2
[pspp] / src / output / charts / scree-cairo.c
1 /* PSPP - a program for statistical analysis.
2    Copyright (C) 2009, 2011 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 #include <config.h>
18
19 #include "output/charts/scree.h"
20
21 #include <math.h>
22
23 #include "output/cairo-chart.h"
24
25 #include "gettext.h"
26 #define _(msgid) gettext (msgid)
27
28 void
29 xrchart_draw_scree (const struct chart_item *chart_item, cairo_t *cr,
30                     struct xrchart_geometry *geom)
31 {
32   const struct scree *rc = to_scree (chart_item);
33   size_t i;
34   double min, max;
35
36   xrchart_write_title (cr, geom, _("Scree Plot"));
37   xrchart_write_xlabel (cr, geom, rc->xlabel);
38   xrchart_write_ylabel (cr, geom, _("Eigenvalue"));
39
40   gsl_vector_minmax (rc->eval, &min, &max);
41
42   if ( fabs (max) > fabs (min))
43     max = fabs (max);
44   else
45     max = fabs (min);
46
47   xrchart_write_yscale (cr, geom, 0, max);
48   xrchart_write_xscale (cr, geom, 0, rc->eval->size + 1);
49
50   xrchart_vector_start (cr, geom, "");
51   for (i = 0 ; i < rc->eval->size; ++i)
52     {
53       const double x = 1 + i;
54       const double y = gsl_vector_get (rc->eval, i);
55       xrchart_vector (cr, geom, x, y);
56     }
57   xrchart_vector_end (cr, geom);
58
59   for (i = 0 ; i < rc->eval->size; ++i)
60     {
61       const double x = 1 + i;
62       const double y = gsl_vector_get (rc->eval, i);
63       xrchart_datum (cr, geom, 0, x, y);
64     }
65 }