* Bowing params are working better.
[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 #include <stdio.h>
35
36 G_DEFINE_TYPE(Graph, graph, G_TYPE_OBJECT);
37
38 static void graph_class_init(GraphClass *klass) {
39   GObjectClass *gobject_klass;
40   gobject_klass = G_OBJECT_CLASS(klass);
41 }
42
43 static void graph_init(Graph *graph) {
44   graph->symbol = CIRCLE;
45   graph->points = g_ptr_array_new();
46   graph->lines = NULL;
47   graph->xaxis = graph->yaxis = NULL;
48 }
49
50 Graph *graph_new(void) {
51   return g_object_new(GRAPH_TYPE, NULL);
52 }
53
54 /*1 means you're stupid and passed NULL*/
55 int graph_add_graph_point(Graph* graph, struct graph_point *pt) {
56   if(pt == NULL) return 1;
57   g_ptr_array_add(graph->points, pt);
58   if((graph->points->len == 1) && ((graph->lines == NULL) || (graph->lines->len == 0))) {
59     graph->maxx = graph->minx = pt->x;
60     graph->maxy = graph->miny = pt->y;
61   }else{
62     if(pt->x > graph->maxx) graph->maxx = pt->x;
63     if(pt->x < graph->minx) graph->minx = pt->x;
64     if(pt->y > graph->maxy) graph->maxy = pt->y;
65     if(pt->y < graph->miny) graph->miny = pt->y;
66   }
67   return 0;
68 }
69
70 /*
71  *
72  *1 means failed to allocate point
73  *NOTE: label is duplicated if not NULL (using g_strdup)
74  */
75 int graph_add_point(Graph* graph, double x, double y, const GString *label) {
76   struct graph_point *pt = (struct graph_point*)malloc(sizeof(struct graph_point));
77   if(pt == NULL) return 1;
78   pt->x = x;
79   pt->y = y;
80   pt->label = (label)?g_string_new(label->str):NULL;
81   graph_add_graph_point(graph, pt);
82   return 0;
83 }
84
85 int graph_add_graph_line(Graph* graph, struct graph_line *l) {
86   if(l == NULL) return 1;
87   if(graph->lines == NULL) {
88     if((graph->lines = g_ptr_array_new()) == NULL) return 1;
89   }
90   g_ptr_array_add(graph->lines, l);
91   if((graph->points->len == 0) && 
92      (graph->lines->len == 1)) {
93     graph->maxx = graph->minx = l->p0_x;
94     graph->maxy = graph->miny = l->p0_y;
95   }else{
96     if(l->p0_x > graph->maxx) graph->maxx = l->p0_x;
97     if(l->p0_x < graph->minx) graph->minx = l->p0_x;
98     if(l->p0_y > graph->maxy) graph->maxy = l->p0_y;
99     if(l->p0_y < graph->miny) graph->miny = l->p0_y;
100   }
101   if(l->p3_x > graph->maxx) graph->maxx = l->p3_x;
102   if(l->p3_x < graph->minx) graph->minx = l->p3_x;
103   if(l->p3_y > graph->maxy) graph->maxy = l->p3_y;
104   if(l->p3_y < graph->miny) graph->miny = l->p3_y;
105   return 0;
106 }
107
108
109 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) {
110   #ifdef DEBUG
111   fprintf(stderr, "graph_add_line(%d, p0_x=%g, p0_y=%g, p1_x=%g, p1_y=%g, p2_x=%g, p2_y=%g, p3_x=%g, p3_y=%g\n", (unsigned int) graph, p0_x, p0_y, p1_x, p1_y, p2_x, p2_y, p3_x, p3_y);
112   #endif
113   struct graph_line *l;
114   if((l=(struct graph_line*)malloc(sizeof(struct graph_line))) == NULL) return 2;
115   l->p0_x=p0_x;
116   l->p0_y=p0_y;
117   l->p1_x=p1_x;
118   l->p1_y=p1_y;
119   l->p2_x=p2_x;
120   l->p2_y=p2_y;
121   l->p3_x=p3_x;
122   l->p3_y=p3_y;
123   return graph_add_graph_line(graph, l);
124 }
125     
126 struct cxt {
127   Graph *graph;
128   struct graph_point *p0;
129 };
130     
131 static void linear_interpolate(gpointer data, gpointer user_data) {
132   struct cxt *cxt = user_data;
133   struct graph_point *p3 = data;
134   if(cxt->p0 != NULL) {
135     double p1_x = (p3->x + 2*cxt->p0->x)/3.0;
136     double p1_y = (p3->y + 2*cxt->p0->y)/3.0;
137     double p2_x = (2*p3->x + cxt->p0->x)/3.0;
138     double p2_y = (2*p3->y + cxt->p0->y)/3.0;
139     graph_add_line(cxt->graph, cxt->p0->x, cxt->p0->y, p1_x, p1_y, p2_x, p2_y, p3->x, p3->y);
140   }
141   cxt->p0 = p3;
142 }
143     
144 void graph_add_linear_connectors(Graph* graph) {
145   struct cxt cxt;
146   cxt.graph = graph;
147   cxt.p0 = NULL;
148   g_ptr_array_foreach(graph->points, &linear_interpolate, &cxt);
149 }
150
151 void graph_bezier_linear_to_cubic(double x0, double x3, double *x1, double *x2) {
152   *x1 = (x3 + 2.0*x0)/3.0;
153   *x2 = (2.0*x3 + x0)/3.0;
154 }
155
156 void graph_bezier_linear_to_quadratic(double x0, double x3, double *x1) {
157   *x1 = (x3 + x0)/2.0;
158 }
159
160 void graph_bezier_quadratic_to_cubic(double x0, double x3, double *x1, double *x2) {
161   *x2 = (x3-2.0*(*x1))/3.0;
162   *x1 = (x0-2.0*(*x1))/3.0;
163 }
164
165 void graph_set_xaxis(Graph *g, struct graph_axis *axis) {
166   if(g->xaxis != NULL) {
167     free(g->xaxis);
168   }
169   g->xaxis=axis;
170 }
171 void graph_set_yaxis(Graph *g, struct graph_axis *axis) {
172   if(g->yaxis != NULL) {
173     free(g->yaxis);
174   }
175   g->yaxis=axis;
176 }
177
178 void graph_autoset_xaxis(Graph *g) {
179   struct graph_axis *axis;
180   if((axis=(struct graph_axis *)malloc(sizeof(struct graph_axis))) == NULL) {
181     return 1;
182   }
183   graph_set_xaxis(g, axis);
184 }
185
186 void graph_autoset_yaxis(Graph *g) {
187   struct graph_axis *axis;
188   if((axis=(struct graph_axis *)malloc(sizeof(struct graph_axis))) == NULL) {
189     return 1;
190   }
191   graph_set_xaxis(g, axis);
192 }