git.maemo.org Git - scdataviz/blob - 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 }
  48
  49 Graph *graph_new(void) {
  50   return g_object_new(GRAPH_TYPE, NULL);
  51 }
  52
  53 /*1 means you're stupid and passed NULL*/
  54 int graph_add_graph_point(Graph* graph, struct graph_point *pt) {
  55   if(pt == NULL) return 1;
  56   g_ptr_array_add(graph->points, pt);
  57   if((graph->points->len == 1) && ((graph->lines == NULL) || (graph->lines->len == 0))) {
  58     graph->maxx = graph->minx = pt->x;
  59     graph->maxy = graph->miny = pt->y;
  60   }else{
  61     if(pt->x > graph->maxx) graph->maxx = pt->x;
  62     if(pt->x < graph->minx) graph->minx = pt->x;
  63     if(pt->y > graph->maxy) graph->maxy = pt->y;
  64     if(pt->y < graph->miny) graph->miny = pt->y;
  65   }
  66   return 0;
  67 }
  68
  69 /*
  70  *
  71  *1 means failed to allocate point
  72  *NOTE: label is duplicated if not NULL (using g_strdup)
  73  */
  74 int graph_add_point(Graph* graph, double x, double y, const GString *label) {
  75   struct graph_point *pt = (struct graph_point*)malloc(sizeof(struct graph_point));
  76   if(pt == NULL) return 1;
  77   pt->x = x;
  78   pt->y = y;
  79   pt->label = (label)?g_string_new(label->str):NULL;
  80   graph_add_graph_point(graph, pt);
  81   return 0;
  82 }
  83
  84 int graph_add_graph_line(Graph* graph, struct graph_line *l) {
  85   if(graph->lines == NULL) {
  86     if((graph->lines = g_ptr_array_new()) == NULL) return 1;
  87   }
  88   g_ptr_array_add(graph->lines, l);
  89   if((graph->points->len == 0) &&
  90      (graph->lines->len == 1)) {
  91     graph->maxx = graph->minx = l->p0_x;
  92     graph->maxy = graph->miny = l->p0_y;
  93   }else{
  94     if(l->p0_x > graph->maxx) graph->maxx = l->p0_x;
  95     if(l->p0_x < graph->minx) graph->minx = l->p0_x;
  96     if(l->p0_y > graph->maxy) graph->maxy = l->p0_y;
  97     if(l->p0_y < graph->miny) graph->miny = l->p0_y;
  98   }
  99   if(l->p3_x > graph->maxx) graph->maxx = l->p3_x;
 100   if(l->p3_x < graph->minx) graph->minx = l->p3_x;
 101   if(l->p3_y > graph->maxy) graph->maxy = l->p3_y;
 102   if(l->p3_y < graph->miny) graph->miny = l->p3_y;
 103   return 0;
 104 }
 105
 106
 107 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) {
 108   #ifdef DEBUG
 109   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);
 110   #endif
 111   struct graph_line *l;
 112   if((l=(struct graph_line*)malloc(sizeof(struct graph_line))) == NULL) return 2;
 113   l->p0_x=p0_x;
 114   l->p0_y=p0_y;
 115   l->p1_x=p1_x;
 116   l->p1_y=p1_y;
 117   l->p2_x=p2_x;
 118   l->p2_y=p2_y;
 119   l->p3_x=p3_x;
 120   l->p3_y=p3_y;
 121   return graph_add_graph_line(graph, l);
 122 }
 123
 124 struct cxt {
 125   Graph *graph;
 126   struct graph_point *p0;
 127 };
 128
 129 static void linear_interpolate(gpointer data, gpointer user_data) {
 130   struct cxt *cxt = user_data;
 131   struct graph_point *p3 = data;
 132   if(cxt->p0 != NULL) {
 133     double p1_x = (p3->x + 2*cxt->p0->x)/3.0;
 134     double p1_y = (p3->y + 2*cxt->p0->y)/3.0;
 135     double p2_x = (2*p3->x + cxt->p0->x)/3.0;
 136     double p2_y = (2*p3->y + cxt->p0->y)/3.0;
 137     graph_add_line(cxt->graph, cxt->p0->x, cxt->p0->y, p1_x, p1_y, p2_x, p2_y, p3->x, p3->y);
 138   }
 139   cxt->p0 = p3;
 140 }
 141
 142 void graph_add_linear_connectors(Graph* graph) {
 143   struct cxt cxt;
 144   cxt.graph = graph;
 145   cxt.p0 = NULL;
 146   g_ptr_array_foreach(graph->points, &linear_interpolate, &cxt);
 147 }
 148
 149 void graph_bezier_linear_to_cubic(double x0, double x3, double *x1, double *x2) {
 150   *x1 = (x3 + 2.0*x0)/3.0;
 151   *x2 = (2.0*x3 + x0)/3.0;
 152 }
 153
 154 void graph_bezier_linear_to_quadratic(double x0, double x3, double *x1) {
 155   *x1 = (x3 + x0)/2.0;
 156 }
 157
 158 void graph_bezier_quadratic_to_cubic(double x0, double x3, double *x1, double *x2) {
 159   *x2 = (x3-2.0*(*x1))/3.0;
 160   *x1 = (x0-2.0*(*x1))/3.0;
 161 }