+
+/*
+ *
+ *1 means failed to allocate point
+ *NOTE: label is duplicated if not NULL (using g_strdup)
+ */
+int graph_add_point(Graph* graph, double x, double y, const GString *label) {
+ struct graph_point *pt = (struct graph_point*)malloc(sizeof(struct graph_point));
+ if(pt == NULL) return 1;
+ pt->x = x;
+ pt->y = y;
+ pt->label = (label)?g_string_new(label->str):NULL;
+ graph_add_graph_point(graph, pt);
+ return 0;
+}
+
+int graph_add_graph_line(Graph* graph, struct graph_line *l) {
+ if(graph->lines == NULL) {
+ if((graph->lines = g_ptr_array_new()) == NULL) return 1;
+ }
+ g_ptr_array_add(graph->lines, l);
+ if((graph->points->len == 0) &&
+ (graph->lines->len == 1)) {
+ graph->maxx = graph->minx = l->p0_x;
+ graph->maxy = graph->miny = l->p0_y;
+ }else{
+ if(l->p0_x > graph->maxx) graph->maxx = l->p0_x;
+ if(l->p0_x < graph->minx) graph->minx = l->p0_x;
+ if(l->p0_y > graph->maxy) graph->maxy = l->p0_y;
+ if(l->p0_y < graph->miny) graph->miny = l->p0_y;
+ }
+ if(l->p3_x > graph->maxx) graph->maxx = l->p3_x;
+ if(l->p3_x < graph->minx) graph->minx = l->p3_x;
+ if(l->p3_y > graph->maxy) graph->maxy = l->p3_y;
+ if(l->p3_y < graph->miny) graph->miny = l->p3_y;
+ return 0;
+}
+
+
+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) {
+ #ifdef DEBUG
+ 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);
+ #endif
+ struct graph_line *l;
+ if((l=(struct graph_line*)malloc(sizeof(struct graph_line))) == NULL) return 2;
+ l->p0_x=p0_x;
+ l->p0_y=p0_y;
+ l->p1_x=p1_x;
+ l->p1_y=p1_y;
+ l->p2_x=p2_x;
+ l->p2_y=p2_y;
+ l->p3_x=p3_x;
+ l->p3_y=p3_y;
+ return graph_add_graph_line(graph, l);
+}
+
+struct cxt {
+ Graph *graph;
+ struct graph_point *p0;
+};
+
+static void linear_interpolate(gpointer data, gpointer user_data) {
+ struct cxt *cxt = user_data;
+ struct graph_point *p3 = data;
+ if(cxt->p0 != NULL) {
+ double p1_x = (p3->x + 2*cxt->p0->x)/3.0;
+ double p1_y = (p3->y + 2*cxt->p0->y)/3.0;
+ double p2_x = (2*p3->x + cxt->p0->x)/3.0;
+ double p2_y = (2*p3->y + cxt->p0->y)/3.0;
+ graph_add_line(cxt->graph, cxt->p0->x, cxt->p0->y, p1_x, p1_y, p2_x, p2_y, p3->x, p3->y);
+ }
+ cxt->p0 = p3;
+}
+
+void graph_add_linear_connectors(Graph* graph) {
+ struct cxt cxt;
+ cxt.graph = graph;
+ cxt.p0 = NULL;
+ g_ptr_array_foreach(graph->points, &linear_interpolate, &cxt);
+}
+
+void graph_bezier_linear_to_cubic(double x0, double x3, double *x1, double *x2) {
+ *x1 = (x3 + 2.0*x0)/3.0;
+ *x2 = (2.0*x3 + x0)/3.0;
+}
+
+void graph_bezier_linear_to_quadratic(double x0, double x3, double *x1) {
+ *x1 = (x3 + x0)/2.0;
+}
+
+void graph_bezier_quadratic_to_cubic(double x0, double x3, double *x1, double *x2) {
+ *x2 = (x3-2.0*(*x1))/3.0;
+ *x1 = (x0-2.0*(*x1))/3.0;
+}