4c75464caaffc2334e6c54e21d4dbd30f5d3e4b0
[scdataviz] / graph.c
1 /*
2 ** Graph.c
3 ** 
4 ** Made by (Johnny Q. Hacker)
5 ** Login   <solarion@johnathan>
6 ** 
7
8
9 Copyright (C) 2008 Joseph Pingenot
10
11 This program is free software: you can redistribute it and/or modify
12 it under the terms of the GNU Affero General Public License as published by
13 the Free Software Foundation, either version 3 of the License, or
14 (at your option) any later version.
15
16 This program is distributed in the hope that it will be useful,
17 but WITHOUT ANY WARRANTY; without even the implied warranty of
18 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
19 GNU Affero General Public License for more details.
20
21 You should have received a copy of the GNU Affero General Public License
22 along with this program.  If not, see <http://www.gnu.org/licenses/>.
23
24
25 **Many thanks to Davyd Madeley for the excellent Cairo  Tutorial
26 **  at http://gnomejournal.org/article/34/writing-a--using-cairo-and-gtk28
27
28 ** Started on  Thu Jul 17 10:51:32 2008 Johnny Q. Hacker
29 ** Last update Sun May 12 01:17:25 2002 Speed Blue
30 */
31
32 #include <graph.h>
33 #include <stdlib.h>
34
35 G_DEFINE_TYPE(Graph, graph, G_TYPE_OBJECT);
36
37 static void graph_class_init(GraphClass *klass) {
38   GObjectClass *gobject_klass;
39   gobject_klass = G_OBJECT_CLASS(klass);
40 }
41
42 static void graph_init(Graph *graph) {
43   graph->symbol = CIRCLE;
44   graph->points = g_ptr_array_new();
45   graph->lines = NULL;
46 }
47
48 Graph *graph_new(void) {
49   return g_object_new(GRAPH_TYPE, NULL);
50 }
51
52 /*1 means you're stupid and passed NULL*/
53 int graph_add_graph_point(Graph* graph, struct graph_point *pt) {
54   if(pt == NULL) return 1;
55   g_ptr_array_add(graph->points, pt);
56   if((graph->points->len == 1) && ((graph->lines == NULL) || (graph->lines->len == 0))) {
57     graph->maxx = graph->minx = pt->x;
58     graph->maxy = graph->miny = pt->y;
59   }else{
60     if(pt->x > graph->maxx) graph->maxx = pt->x;
61     if(pt->x < graph->minx) graph->minx = pt->x;
62     if(pt->y > graph->maxy) graph->maxy = pt->y;
63     if(pt->y < graph->miny) graph->miny = pt->y;
64   }
65   return 0;
66 }
67
68 /*
69  *
70  *1 means failed to allocate point
71  *NOTE: label is duplicated if not NULL (using g_strdup)
72  */
73 int graph_add_point(Graph* graph, double x, double y, const GString *label) {
74   struct graph_point *pt = (struct graph_point*)malloc(sizeof(struct graph_point));
75   if(pt == NULL) return 1;
76   pt->x = x;
77   pt->y = y;
78   pt->label = (label)?g_string_new(label->str):NULL;
79   graph_add_graph_point(graph, pt);
80   return 0;
81 }
82
83 int graph_add_line(Graph* graph, double p0_x, double p0_y, double p1_x, double p1_y, double p2_x, double p2_y, double p3_x, double p3_y) {
84   if(graph->lines == NULL) {
85     if((graph->points = g_ptr_array_new()) == NULL) return 1;
86   }
87   struct graph_line *l;
88   if((l=(struct graph_line*)malloc(sizeof(struct graph_line))) == NULL) return 2;
89   l->p0_x=p0_x;
90   l->p0_y=p0_y;
91   l->p1_x=p1_x;
92   l->p1_y=p1_y;
93   l->p2_x=p2_x;
94   l->p2_y=p2_y;
95   l->p3_x=p3_x;
96   l->p3_y=p3_y;
97   g_ptr_array_add(graph->lines, l);
98   if((graph->points->len == 0) && (graph->lines->len == 1)) {
99     graph->maxx = graph->minx = p0_x;
100     graph->maxy = graph->miny = p0_y;
101   }else{
102     if(p0_x > graph->maxx) graph->maxx = p0_x;
103     if(p0_x < graph->minx) graph->minx = p0_x;
104     if(p0_y > graph->maxy) graph->maxy = p0_y;
105     if(p0_y < graph->miny) graph->miny = p0_y;
106   }
107   if(p3_x > graph->maxx) graph->maxx = p3_x;
108   if(p3_x < graph->minx) graph->minx = p3_x;
109   if(p3_y > graph->maxy) graph->maxy = p3_y;
110   if(p3_y < graph->miny) graph->miny = p3_y;
111   return 0;
112 }
113     
114 struct cxt {
115   Graph *graph;
116   struct graph_point *p0;
117 };
118     
119 static void linear_interpolate(gpointer data, gpointer user_data) {
120   struct cxt *cxt = user_data;
121   struct graph_point *p3 = data;
122   if(cxt->p0 != NULL) {
123     double p1_x = (p3->x - cxt->p0->x)/3.0;
124     double p1_y = (p3->y - cxt->p0->y)/3.0;
125     double p2_x = (p3->x + 2*cxt->p0->x)/3.0;
126     double p2_y = (p3->y + 2*cxt->p0->y)/3.0;
127     graph_add_line(cxt->graph, cxt->p0->x, cxt->p0->y, p1_x, p1_y, p2_x, p2_y, p3->x, p3->y);
128   }
129   cxt->p0 = p3;
130 }
131     
132 void graph_add_linear_connectors(Graph* graph) {
133   struct cxt cxt;
134   cxt.graph = graph;
135   cxt.p0 = NULL;
136   g_ptr_array_foreach(graph->points, &linear_interpolate, &cxt);
137 }