# HG changeset patch # User Matti Hamalainen # Date 1131356944 0 # Node ID 4f247b19c9ea4f2f4c2a34efedc94c5e7510c178 # Parent 70ab1c28d4b2573a0c99579b70662fe1b9d21a3f Add new files diff -r 70ab1c28d4b2 -r 4f247b19c9ea src/xs_curve.c --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/src/xs_curve.c Mon Nov 07 09:49:04 2005 +0000 @@ -0,0 +1,963 @@ +/* GTK - The GIMP Toolkit + * Copyright (C) 1997 David Mosberger + * + * This library is free software; you can redistribute it and/or + * modify it under the terms of the GNU Library General Public + * License as published by the Free Software Foundation; either + * version 2 of the License, or (at your option) any later version. + * + * This library is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Library General Public License for more details. + * + * You should have received a copy of the GNU Library General Public + * License along with this library; if not, write to the + * Free Software Foundation, Inc., 59 Temple Place - Suite 330, + * Boston, MA 02111-1307, USA. + */ + +/* + * Modified by the GTK+ Team and others 1997-1999. See the AUTHORS + * file for a list of people on the GTK+ Team. See the ChangeLog + * files for a list of changes. These files are distributed with + * GTK+ at ftp://ftp.gtk.org/pub/gtk/. + */ + +#include +#include +#include + +#include "xs_curve.h" +#include +#include +#include +/* +#include +#include +*/ + +#define RADIUS 3 /* radius of the control points */ +#define MIN_DISTANCE 8 /* min distance between control points */ + +#define GRAPH_MASK (GDK_EXPOSURE_MASK | \ + GDK_POINTER_MOTION_MASK | \ + GDK_POINTER_MOTION_HINT_MASK | \ + GDK_ENTER_NOTIFY_MASK | \ + GDK_BUTTON_PRESS_MASK | \ + GDK_BUTTON_RELEASE_MASK | \ + GDK_BUTTON1_MOTION_MASK) + +enum { + ARG_0, + ARG_CURVE_TYPE, + ARG_MIN_X, + ARG_MAX_X, + ARG_MIN_Y, + ARG_MAX_Y +}; + +static GtkDrawingAreaClass *parent_class = NULL; +static guint curve_type_changed_signal = 0; + + +/* forward declarations: */ +static void gtk_curve_class_init (GtkCurveClass *class); +static void gtk_curve_init (GtkCurve *curve); +static void gtk_curve_set_arg (GtkObject *object, + GtkArg *arg, + guint arg_id); +static void gtk_curve_get_arg (GtkObject *object, + GtkArg *arg, + guint arg_id); +static void gtk_curve_finalize (GtkObject *object); +static gint gtk_curve_graph_events (GtkWidget *widget, + GdkEvent *event, + GtkCurve *c); +static void gtk_curve_size_graph (GtkCurve *curve); + +GtkType +gtk_curve_get_type (void) +{ + static GtkType curve_type = 0; + + if (!curve_type) + { + static const GtkTypeInfo curve_info = + { + "GtkCurve", + sizeof (GtkCurve), + sizeof (GtkCurveClass), + (GtkClassInitFunc) gtk_curve_class_init, + (GtkObjectInitFunc) gtk_curve_init, + /* reserved_1 */ NULL, + /* reserved_2 */ NULL, + (GtkClassInitFunc) NULL, + }; + + curve_type = gtk_type_unique (GTK_TYPE_DRAWING_AREA, &curve_info); + } + return curve_type; +} + +static void +gtk_curve_class_init (GtkCurveClass *class) +{ + GtkObjectClass *object_class; + + parent_class = gtk_type_class (GTK_TYPE_DRAWING_AREA); + + object_class = (GtkObjectClass *) class; + + object_class->set_arg = gtk_curve_set_arg; + object_class->get_arg = gtk_curve_get_arg; + object_class->finalize = gtk_curve_finalize; + + curve_type_changed_signal = + gtk_signal_new ("curve_type_changed", GTK_RUN_FIRST, object_class->type, + GTK_SIGNAL_OFFSET (GtkCurveClass, curve_type_changed), + gtk_marshal_NONE__NONE, GTK_TYPE_NONE, 0); + gtk_object_class_add_signals (object_class, &curve_type_changed_signal, 1); + + gtk_object_add_arg_type ("GtkCurve::curve_type", GTK_TYPE_CURVE_TYPE, + GTK_ARG_READWRITE, ARG_CURVE_TYPE); + gtk_object_add_arg_type ("GtkCurve::min_x", GTK_TYPE_FLOAT, + GTK_ARG_READWRITE, ARG_MIN_X); + gtk_object_add_arg_type ("GtkCurve::max_x", GTK_TYPE_FLOAT, + GTK_ARG_READWRITE, ARG_MAX_X); + gtk_object_add_arg_type ("GtkCurve::min_y", GTK_TYPE_FLOAT, + GTK_ARG_READWRITE, ARG_MIN_Y); + gtk_object_add_arg_type ("GtkCurve::max_y", GTK_TYPE_FLOAT, + GTK_ARG_READWRITE, ARG_MAX_Y); +} + +static void +gtk_curve_init (GtkCurve *curve) +{ + gint old_mask; + + curve->cursor_type = GDK_TOP_LEFT_ARROW; + curve->pixmap = NULL; + curve->curve_type = GTK_CURVE_TYPE_SPLINE; + curve->height = 0; + curve->grab_point = -1; + + curve->num_points = 0; + curve->point = 0; + + curve->num_ctlpoints = 0; + curve->ctlpoint = NULL; + + curve->min_x = 0.0; + curve->max_x = 1.0; + curve->min_y = 0.0; + curve->max_y = 1.0; + + old_mask = gtk_widget_get_events (GTK_WIDGET (curve)); + gtk_widget_set_events (GTK_WIDGET (curve), old_mask | GRAPH_MASK); + gtk_signal_connect (GTK_OBJECT (curve), "event", + (GtkSignalFunc) gtk_curve_graph_events, curve); + gtk_curve_size_graph (curve); +} + +static void +gtk_curve_set_arg (GtkObject *object, + GtkArg *arg, + guint arg_id) +{ + GtkCurve *curve = GTK_CURVE (object); + + switch (arg_id) + { + case ARG_CURVE_TYPE: + gtk_curve_set_curve_type (curve, GTK_VALUE_ENUM (*arg)); + break; + case ARG_MIN_X: + gtk_curve_set_range (curve, GTK_VALUE_FLOAT (*arg), curve->max_x, + curve->min_y, curve->max_y); + break; + case ARG_MAX_X: + gtk_curve_set_range (curve, curve->min_x, GTK_VALUE_FLOAT (*arg), + curve->min_y, curve->max_y); + break; + case ARG_MIN_Y: + gtk_curve_set_range (curve, curve->min_x, curve->max_x, + GTK_VALUE_FLOAT (*arg), curve->max_y); + break; + case ARG_MAX_Y: + gtk_curve_set_range (curve, curve->min_x, curve->max_x, + curve->min_y, GTK_VALUE_FLOAT (*arg)); + break; + } +} + +static void +gtk_curve_get_arg (GtkObject *object, + GtkArg *arg, + guint arg_id) +{ + GtkCurve *curve = GTK_CURVE (object); + + switch (arg_id) + { + case ARG_CURVE_TYPE: + GTK_VALUE_ENUM (*arg) = curve->curve_type; + break; + case ARG_MIN_X: + GTK_VALUE_FLOAT (*arg) = curve->min_x; + break; + case ARG_MAX_X: + GTK_VALUE_FLOAT (*arg) = curve->max_x; + break; + case ARG_MIN_Y: + GTK_VALUE_FLOAT (*arg) = curve->min_y; + break; + case ARG_MAX_Y: + GTK_VALUE_FLOAT (*arg) = curve->max_y; + break; + default: + arg->type = GTK_TYPE_INVALID; + break; + } +} + +static int +project (gfloat value, gfloat min, gfloat max, int norm) +{ + return (norm - 1) * ((value - min) / (max - min)) + 0.5; +} + +static gfloat +unproject (gint value, gfloat min, gfloat max, int norm) +{ + return value / (gfloat) (norm - 1) * (max - min) + min; +} + +/* Solve the tridiagonal equation system that determines the second + derivatives for the interpolation points. (Based on Numerical + Recipies 2nd Edition.) */ +static void +spline_solve (int n, gfloat x[], gfloat y[], gfloat y2[]) +{ + gfloat p, sig, *u; + gint i, k; + + u = g_malloc ((n - 1) * sizeof (u[0])); + + y2[0] = u[0] = 0.0; /* set lower boundary condition to "natural" */ + + for (i = 1; i < n - 1; ++i) + { + sig = (x[i] - x[i - 1]) / (x[i + 1] - x[i - 1]); + p = sig * y2[i - 1] + 2.0; + y2[i] = (sig - 1.0) / p; + u[i] = ((y[i + 1] - y[i]) + / (x[i + 1] - x[i]) - (y[i] - y[i - 1]) / (x[i] - x[i - 1])); + u[i] = (6.0 * u[i] / (x[i + 1] - x[i - 1]) - sig * u[i - 1]) / p; + } + + y2[n - 1] = 0.0; + for (k = n - 2; k >= 0; --k) + y2[k] = y2[k] * y2[k + 1] + u[k]; + + g_free (u); +} + +static gfloat +spline_eval (int n, gfloat x[], gfloat y[], gfloat y2[], gfloat val) +{ + gint k_lo, k_hi, k; + gfloat h, b, a; + + /* do a binary search for the right interval: */ + k_lo = 0; k_hi = n - 1; + while (k_hi - k_lo > 1) + { + k = (k_hi + k_lo) / 2; + if (x[k] > val) + k_hi = k; + else + k_lo = k; + } + + h = x[k_hi] - x[k_lo]; + g_assert (h > 0.0); + + a = (x[k_hi] - val) / h; + b = (val - x[k_lo]) / h; + return a*y[k_lo] + b*y[k_hi] + + ((a*a*a - a)*y2[k_lo] + (b*b*b - b)*y2[k_hi]) * (h*h)/6.0; +} + +static void +gtk_curve_interpolate (GtkCurve *c, gint width, gint height) +{ + gfloat *vector; + int i; + + vector = g_malloc (width * sizeof (vector[0])); + + gtk_curve_get_vector (c, width, vector); + + c->height = height; + if (c->num_points != width) + { + c->num_points = width; + if (c->point) + g_free (c->point); + c->point = g_malloc (c->num_points * sizeof (c->point[0])); + } + + for (i = 0; i < width; ++i) + { + c->point[i].x = RADIUS + i; + c->point[i].y = RADIUS + height + - project (vector[i], c->min_y, c->max_y, height); + } + + g_free (vector); +} + +static void +gtk_curve_draw (GtkCurve *c, gint width, gint height) +{ + GtkStateType state; + GtkStyle *style; + gint i; + + if (!c->pixmap) + return; + + if (c->height != height || c->num_points != width) + gtk_curve_interpolate (c, width, height); + + state = GTK_STATE_NORMAL; + if (!GTK_WIDGET_IS_SENSITIVE (GTK_WIDGET (c))) + state = GTK_STATE_INSENSITIVE; + + style = GTK_WIDGET (c)->style; + + /* clear the pixmap: */ + gtk_paint_flat_box (style, c->pixmap, GTK_STATE_NORMAL, GTK_SHADOW_NONE, + NULL, GTK_WIDGET(c), "curve_bg", + 0, 0, width + RADIUS * 2, height + RADIUS * 2); + /* draw the grid lines: (XXX make more meaningful) */ + for (i = 0; i < 5; i++) + { + gdk_draw_line (c->pixmap, style->dark_gc[state], + RADIUS, i * (height / 4.0) + RADIUS, + width + RADIUS, i * (height / 4.0) + RADIUS); + gdk_draw_line (c->pixmap, style->dark_gc[state], + i * (width / 4.0) + RADIUS, RADIUS, + i * (width / 4.0) + RADIUS, height + RADIUS); + } + + gdk_draw_points (c->pixmap, style->fg_gc[state], c->point, c->num_points); + if (c->curve_type != GTK_CURVE_TYPE_FREE) + for (i = 0; i < c->num_ctlpoints; ++i) + { + gint x, y; + + if (c->ctlpoint[i][0] < c->min_x) + continue; + + x = project (c->ctlpoint[i][0], c->min_x, c->max_x, + width); + y = height - + project (c->ctlpoint[i][1], c->min_y, c->max_y, + height); + + /* draw a bullet: */ + gdk_draw_arc (c->pixmap, style->fg_gc[state], TRUE, x, y, + RADIUS * 2, RADIUS*2, 0, 360*64); + } + gdk_draw_pixmap (GTK_WIDGET (c)->window, style->fg_gc[state], c->pixmap, + 0, 0, 0, 0, width + RADIUS * 2, height + RADIUS * 2); +} + +static gint +gtk_curve_graph_events (GtkWidget *widget, GdkEvent *event, GtkCurve *c) +{ + GdkCursorType new_type = c->cursor_type; + gint i, src, dst, leftbound, rightbound; + GdkEventButton *bevent; + GdkEventMotion *mevent; + GtkWidget *w; + gint tx, ty; + gint cx, x, y, width, height; + gint closest_point = 0; + gfloat rx, ry, min_x; + guint distance; + gint x1, x2, y1, y2; + + w = GTK_WIDGET (c); + width = w->allocation.width - RADIUS * 2; + height = w->allocation.height - RADIUS * 2; + + if ((width < 0) || (height < 0)) + return FALSE; + + /* get the pointer position */ + gdk_window_get_pointer (w->window, &tx, &ty, NULL); + x = CLAMP ((tx - RADIUS), 0, width-1); + y = CLAMP ((ty - RADIUS), 0, height-1); + + min_x = c->min_x; + + distance = ~0U; + for (i = 0; i < c->num_ctlpoints; ++i) + { + cx = project (c->ctlpoint[i][0], min_x, c->max_x, width); + if ((guint) abs (x - cx) < distance) + { + distance = abs (x - cx); + closest_point = i; + } + } + + switch (event->type) + { + case GDK_CONFIGURE: + if (c->pixmap) + gdk_pixmap_unref (c->pixmap); + c->pixmap = 0; + /* fall through */ + case GDK_EXPOSE: + if (!c->pixmap) + c->pixmap = gdk_pixmap_new (w->window, + w->allocation.width, + w->allocation.height, -1); + gtk_curve_draw (c, width, height); + break; + + case GDK_BUTTON_PRESS: + gtk_grab_add (widget); + + bevent = (GdkEventButton *) event; + new_type = GDK_TCROSS; + + switch (c->curve_type) + { + case GTK_CURVE_TYPE_LINEAR: + case GTK_CURVE_TYPE_SPLINE: + if (distance > MIN_DISTANCE) + { + /* insert a new control point */ + if (c->num_ctlpoints > 0) + { + cx = project (c->ctlpoint[closest_point][0], min_x, + c->max_x, width); + if (x > cx) + ++closest_point; + } + ++c->num_ctlpoints; + c->ctlpoint = + g_realloc (c->ctlpoint, + c->num_ctlpoints * sizeof (*c->ctlpoint)); + for (i = c->num_ctlpoints - 1; i > closest_point; --i) + memcpy (c->ctlpoint + i, c->ctlpoint + i - 1, + sizeof (*c->ctlpoint)); + } + c->grab_point = closest_point; + c->ctlpoint[c->grab_point][0] = + unproject (x, min_x, c->max_x, width); + c->ctlpoint[c->grab_point][1] = + unproject (height - y, c->min_y, c->max_y, height); + + gtk_curve_interpolate (c, width, height); + break; + + case GTK_CURVE_TYPE_FREE: + c->point[x].x = RADIUS + x; + c->point[x].y = RADIUS + y; + c->grab_point = x; + c->last = y; + break; + } + gtk_curve_draw (c, width, height); + break; + + case GDK_BUTTON_RELEASE: + gtk_grab_remove (widget); + + /* delete inactive points: */ + if (c->curve_type != GTK_CURVE_TYPE_FREE) + { + for (src = dst = 0; src < c->num_ctlpoints; ++src) + { + if (c->ctlpoint[src][0] >= min_x) + { + memcpy (c->ctlpoint + dst, c->ctlpoint + src, + sizeof (*c->ctlpoint)); + ++dst; + } + } + if (dst < src) + { + c->num_ctlpoints -= (src - dst); + if (c->num_ctlpoints <= 0) + { + c->num_ctlpoints = 1; + c->ctlpoint[0][0] = min_x; + c->ctlpoint[0][1] = c->min_y; + gtk_curve_interpolate (c, width, height); + gtk_curve_draw (c, width, height); + } + c->ctlpoint = + g_realloc (c->ctlpoint, + c->num_ctlpoints * sizeof (*c->ctlpoint)); + } + } + new_type = GDK_FLEUR; + c->grab_point = -1; + break; + + case GDK_MOTION_NOTIFY: + mevent = (GdkEventMotion *) event; + + switch (c->curve_type) + { + case GTK_CURVE_TYPE_LINEAR: + case GTK_CURVE_TYPE_SPLINE: + if (c->grab_point == -1) + { + /* if no point is grabbed... */ + if (distance <= MIN_DISTANCE) + new_type = GDK_FLEUR; + else + new_type = GDK_TCROSS; + } + else + { + /* drag the grabbed point */ + new_type = GDK_TCROSS; + + leftbound = -MIN_DISTANCE; + if (c->grab_point > 0) + leftbound = project (c->ctlpoint[c->grab_point - 1][0], + min_x, c->max_x, width); + + rightbound = width + RADIUS * 2 + MIN_DISTANCE; + if (c->grab_point + 1 < c->num_ctlpoints) + rightbound = project (c->ctlpoint[c->grab_point + 1][0], + min_x, c->max_x, width); + + if (tx <= leftbound || tx >= rightbound + || ty > height + RADIUS * 2 + MIN_DISTANCE + || ty < -MIN_DISTANCE) + c->ctlpoint[c->grab_point][0] = min_x - 1.0; + else + { + rx = unproject (x, min_x, c->max_x, width); + ry = unproject (height - y, c->min_y, c->max_y, height); + c->ctlpoint[c->grab_point][0] = rx; + c->ctlpoint[c->grab_point][1] = ry; + } + gtk_curve_interpolate (c, width, height); + gtk_curve_draw (c, width, height); + } + break; + + case GTK_CURVE_TYPE_FREE: + if (c->grab_point != -1) + { + if (c->grab_point > x) + { + x1 = x; + x2 = c->grab_point; + y1 = y; + y2 = c->last; + } + else + { + x1 = c->grab_point; + x2 = x; + y1 = c->last; + y2 = y; + } + + if (x2 != x1) + for (i = x1; i <= x2; i++) + { + c->point[i].x = RADIUS + i; + c->point[i].y = RADIUS + + (y1 + ((y2 - y1) * (i - x1)) / (x2 - x1)); + } + else + { + c->point[x].x = RADIUS + x; + c->point[x].y = RADIUS + y; + } + c->grab_point = x; + c->last = y; + gtk_curve_draw (c, width, height); + } + if (mevent->state & GDK_BUTTON1_MASK) + new_type = GDK_TCROSS; + else + new_type = GDK_PENCIL; + break; + } + if (new_type != (GdkCursorType) c->cursor_type) + { + GdkCursor *cursor; + + c->cursor_type = new_type; + + cursor = gdk_cursor_new (c->cursor_type); + gdk_window_set_cursor (w->window, cursor); + gdk_cursor_destroy (cursor); + } + break; + + default: + break; + } + return FALSE; +} + +void +gtk_curve_set_curve_type (GtkCurve *c, GtkCurveType new_type) +{ + gfloat rx, dx; + gint x, i; + + if (new_type != c->curve_type) + { + gint width, height; + + width = GTK_WIDGET(c)->allocation.width - RADIUS * 2; + height = GTK_WIDGET(c)->allocation.height - RADIUS * 2; + + if (new_type == GTK_CURVE_TYPE_FREE) + { + gtk_curve_interpolate (c, width, height); + c->curve_type = new_type; + } + else if (c->curve_type == GTK_CURVE_TYPE_FREE) + { + if (c->ctlpoint) + g_free (c->ctlpoint); + c->num_ctlpoints = 9; + c->ctlpoint = g_malloc (c->num_ctlpoints * sizeof (*c->ctlpoint)); + + rx = 0.0; + dx = (width - 1) / (gfloat) (c->num_ctlpoints - 1); + + for (i = 0; i < c->num_ctlpoints; ++i, rx += dx) + { + x = (int) (rx + 0.5); + c->ctlpoint[i][0] = + unproject (x, c->min_x, c->max_x, width); + c->ctlpoint[i][1] = + unproject (RADIUS + height - c->point[x].y, + c->min_y, c->max_y, height); + } + c->curve_type = new_type; + gtk_curve_interpolate (c, width, height); + } + else + { + c->curve_type = new_type; + gtk_curve_interpolate (c, width, height); + } + gtk_signal_emit (GTK_OBJECT (c), curve_type_changed_signal); + gtk_curve_draw (c, width, height); + } +} + +static void +gtk_curve_size_graph (GtkCurve *curve) +{ + gint width, height; + gfloat aspect; + + width = (curve->max_x - curve->min_x) + 1; + height = (curve->max_y - curve->min_y) + 1; + aspect = width / (gfloat) height; + if (width > gdk_screen_width () / 4) + width = gdk_screen_width () / 4; + if (height > gdk_screen_height () / 4) + height = gdk_screen_height () / 4; + + if (aspect < 1.0) + width = height * aspect; + else + height = width / aspect; + + gtk_drawing_area_size (GTK_DRAWING_AREA (curve), + width + RADIUS * 2, height + RADIUS * 2); +} + +static void +gtk_curve_reset_vector (GtkCurve *curve) +{ + if (curve->ctlpoint) + g_free (curve->ctlpoint); + + curve->num_ctlpoints = 2; + curve->ctlpoint = g_malloc (2 * sizeof (curve->ctlpoint[0])); + curve->ctlpoint[0][0] = curve->min_x; + curve->ctlpoint[0][1] = curve->min_y; + curve->ctlpoint[1][0] = curve->max_x; + curve->ctlpoint[1][1] = curve->max_y; + + if (curve->pixmap) + { + gint width, height; + + width = GTK_WIDGET (curve)->allocation.width - RADIUS * 2; + height = GTK_WIDGET (curve)->allocation.height - RADIUS * 2; + + if (curve->curve_type == GTK_CURVE_TYPE_FREE) + { + curve->curve_type = GTK_CURVE_TYPE_LINEAR; + gtk_curve_interpolate (curve, width, height); + curve->curve_type = GTK_CURVE_TYPE_FREE; + } + else + gtk_curve_interpolate (curve, width, height); + gtk_curve_draw (curve, width, height); + } +} + +void +gtk_curve_reset (GtkCurve *c) +{ + GtkCurveType old_type; + + old_type = c->curve_type; + c->curve_type = GTK_CURVE_TYPE_SPLINE; + gtk_curve_reset_vector (c); + + if (old_type != GTK_CURVE_TYPE_SPLINE) + gtk_signal_emit (GTK_OBJECT (c), curve_type_changed_signal); +} + +void +gtk_curve_set_gamma (GtkCurve *c, gfloat gamma) +{ + gfloat x, one_over_gamma, height, one_over_width; + GtkCurveType old_type; + gint i; + + if (c->num_points < 2) + return; + + old_type = c->curve_type; + c->curve_type = GTK_CURVE_TYPE_FREE; + + if (gamma <= 0) + one_over_gamma = 1.0; + else + one_over_gamma = 1.0 / gamma; + one_over_width = 1.0 / (c->num_points - 1); + height = c->height; + for (i = 0; i < c->num_points; ++i) + { + x = (gfloat) i / (c->num_points - 1); + c->point[i].x = RADIUS + i; + c->point[i].y = + RADIUS + (height * (1.0 - pow (x, one_over_gamma)) + 0.5); + } + + if (old_type != GTK_CURVE_TYPE_FREE) + gtk_signal_emit (GTK_OBJECT (c), curve_type_changed_signal); + + gtk_curve_draw (c, c->num_points, c->height); +} + +void +gtk_curve_set_range (GtkCurve *curve, + gfloat min_x, gfloat max_x, gfloat min_y, gfloat max_y) +{ + curve->min_x = min_x; + curve->max_x = max_x; + curve->min_y = min_y; + curve->max_y = max_y; + + gtk_curve_size_graph (curve); + gtk_curve_reset_vector (curve); +} + +void +gtk_curve_set_vector (GtkCurve *c, int veclen, gfloat vector[]) +{ + GtkCurveType old_type; + gfloat rx, dx, ry; + gint i, height; + + old_type = c->curve_type; + c->curve_type = GTK_CURVE_TYPE_FREE; + + if (c->point) + height = GTK_WIDGET (c)->allocation.height - RADIUS * 2; + else + { + height = (c->max_y - c->min_y); + if (height > gdk_screen_height () / 4) + height = gdk_screen_height () / 4; + + c->height = height; + c->num_points = veclen; + c->point = g_malloc (c->num_points * sizeof (c->point[0])); + } + rx = 0; + dx = (veclen - 1.0) / (c->num_points - 1.0); + + for (i = 0; i < c->num_points; ++i, rx += dx) + { + ry = vector[(int) (rx + 0.5)]; + if (ry > c->max_y) ry = c->max_y; + if (ry < c->min_y) ry = c->min_y; + c->point[i].x = RADIUS + i; + c->point[i].y = + RADIUS + height - project (ry, c->min_y, c->max_y, height); + } + if (old_type != GTK_CURVE_TYPE_FREE) + gtk_signal_emit (GTK_OBJECT (c), curve_type_changed_signal); + + gtk_curve_draw (c, c->num_points, height); +} + +void +gtk_curve_get_vector (GtkCurve *c, int veclen, gfloat vector[]) +{ + gfloat rx, ry, dx, dy, min_x, delta_x, *mem, *xv, *yv, *y2v, prev; + gint dst, i, x, next, num_active_ctlpoints = 0, first_active = -1; + + min_x = c->min_x; + + if (c->curve_type != GTK_CURVE_TYPE_FREE) + { + /* count active points: */ + prev = min_x - 1.0; + for (i = num_active_ctlpoints = 0; i < c->num_ctlpoints; ++i) + if (c->ctlpoint[i][0] > prev) + { + if (first_active < 0) + first_active = i; + prev = c->ctlpoint[i][0]; + ++num_active_ctlpoints; + } + + /* handle degenerate case: */ + if (num_active_ctlpoints < 2) + { + if (num_active_ctlpoints > 0) + ry = c->ctlpoint[first_active][1]; + else + ry = c->min_y; + if (ry < c->min_y) ry = c->min_y; + if (ry > c->max_y) ry = c->max_y; + for (x = 0; x < veclen; ++x) + vector[x] = ry; + return; + } + } + + switch (c->curve_type) + { + case GTK_CURVE_TYPE_SPLINE: + mem = g_malloc (3 * num_active_ctlpoints * sizeof (gfloat)); + xv = mem; + yv = mem + num_active_ctlpoints; + y2v = mem + 2*num_active_ctlpoints; + + prev = min_x - 1.0; + for (i = dst = 0; i < c->num_ctlpoints; ++i) + if (c->ctlpoint[i][0] > prev) + { + prev = c->ctlpoint[i][0]; + xv[dst] = c->ctlpoint[i][0]; + yv[dst] = c->ctlpoint[i][1]; + ++dst; + } + + spline_solve (num_active_ctlpoints, xv, yv, y2v); + + rx = min_x; + dx = (c->max_x - min_x) / (veclen - 1); + for (x = 0; x < veclen; ++x, rx += dx) + { + ry = spline_eval (num_active_ctlpoints, xv, yv, y2v, rx); + if (ry < c->min_y) ry = c->min_y; + if (ry > c->max_y) ry = c->max_y; + vector[x] = ry; + } + + g_free (mem); + break; + + case GTK_CURVE_TYPE_LINEAR: + dx = (c->max_x - min_x) / (veclen - 1); + rx = min_x; + ry = c->min_y; + dy = 0.0; + i = first_active; + for (x = 0; x < veclen; ++x, rx += dx) + { + if (rx >= c->ctlpoint[i][0]) + { + if (rx > c->ctlpoint[i][0]) + ry = c->min_y; + dy = 0.0; + next = i + 1; + while (next < c->num_ctlpoints + && c->ctlpoint[next][0] <= c->ctlpoint[i][0]) + ++next; + if (next < c->num_ctlpoints) + { + delta_x = c->ctlpoint[next][0] - c->ctlpoint[i][0]; + dy = ((c->ctlpoint[next][1] - c->ctlpoint[i][1]) + / delta_x); + dy *= dx; + ry = c->ctlpoint[i][1]; + i = next; + } + } + vector[x] = ry; + ry += dy; + } + break; + + case GTK_CURVE_TYPE_FREE: + if (c->point) + { + rx = 0.0; + dx = c->num_points / (double) veclen; + for (x = 0; x < veclen; ++x, rx += dx) + vector[x] = unproject (RADIUS + c->height - c->point[(int) rx].y, + c->min_y, c->max_y, + c->height); + } + else + memset (vector, 0, veclen * sizeof (vector[0])); + break; + } +} + +GtkWidget* +gtk_curve_new (void) +{ + return gtk_type_new (gtk_curve_get_type ()); +} + +static void +gtk_curve_finalize (GtkObject *object) +{ + GtkCurve *curve; + + g_return_if_fail (object != NULL); + g_return_if_fail (GTK_IS_CURVE (object)); + + curve = GTK_CURVE (object); + if (curve->pixmap) + gdk_pixmap_unref (curve->pixmap); + if (curve->point) + g_free (curve->point); + if (curve->ctlpoint) + g_free (curve->ctlpoint); + + (*GTK_OBJECT_CLASS (parent_class)->finalize) (object); +} diff -r 70ab1c28d4b2 -r 4f247b19c9ea src/xs_curve.h --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/src/xs_curve.h Mon Nov 07 09:49:04 2005 +0000 @@ -0,0 +1,101 @@ +/* GTK - The GIMP Toolkit + * Copyright (C) 1995-1997 Peter Mattis, Spencer Kimball and Josh MacDonald + * + * This library is free software; you can redistribute it and/or + * modify it under the terms of the GNU Library General Public + * License as published by the Free Software Foundation; either + * version 2 of the License, or (at your option) any later version. + * + * This library is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Library General Public License for more details. + * + * You should have received a copy of the GNU Library General Public + * License along with this library; if not, write to the + * Free Software Foundation, Inc., 59 Temple Place - Suite 330, + * Boston, MA 02111-1307, USA. + */ + +/* + * Modified by the GTK+ Team and others 1997-1999. See the AUTHORS + * file for a list of people on the GTK+ Team. See the ChangeLog + * files for a list of changes. These files are distributed with + * GTK+ at ftp://ftp.gtk.org/pub/gtk/. + */ + +#ifndef __GTK_CURVE_H__ +#define __GTK_CURVE_H__ + + +#include +#include + +#ifdef __cplusplus +extern "C" { +#endif /* __cplusplus */ + + +#define GTK_TYPE_CURVE (gtk_curve_get_type ()) +#define GTK_CURVE(obj) (GTK_CHECK_CAST ((obj), GTK_TYPE_CURVE, GtkCurve)) +#define GTK_CURVE_CLASS(klass) (GTK_CHECK_CLASS_CAST ((klass), GTK_TYPE_CURVE, GtkCurveClass)) +#define GTK_IS_CURVE(obj) (GTK_CHECK_TYPE ((obj), GTK_TYPE_CURVE)) +#define GTK_IS_CURVE_CLASS(klass) (GTK_CHECK_CLASS_TYPE ((klass), GTK_TYPE_CURVE)) + + +typedef struct _GtkCurve GtkCurve; +typedef struct _GtkCurveClass GtkCurveClass; + + +struct _GtkCurve +{ + GtkDrawingArea graph; + + gint cursor_type; + gfloat min_x; + gfloat max_x; + gfloat min_y; + gfloat max_y; + GdkPixmap *pixmap; + GtkCurveType curve_type; + gint height; /* (cached) graph height in pixels */ + gint grab_point; /* point currently grabbed */ + gint last; + + /* (cached) curve points: */ + gint num_points; + GdkPoint *point; + + /* control points: */ + gint num_ctlpoints; /* number of control points */ + gfloat (*ctlpoint)[2]; /* array of control points */ +}; + +struct _GtkCurveClass +{ + GtkDrawingAreaClass parent_class; + + void (* curve_type_changed) (GtkCurve *curve); +}; + + +GtkType gtk_curve_get_type (void); +GtkWidget* gtk_curve_new (void); +void gtk_curve_reset (GtkCurve *curve); +void gtk_curve_set_gamma (GtkCurve *curve, gfloat gamma); +void gtk_curve_set_range (GtkCurve *curve, + gfloat min_x, gfloat max_x, + gfloat min_y, gfloat max_y); +void gtk_curve_get_vector (GtkCurve *curve, + int veclen, gfloat vector[]); +void gtk_curve_set_vector (GtkCurve *curve, + int veclen, gfloat vector[]); +void gtk_curve_set_curve_type (GtkCurve *curve, GtkCurveType type); + + +#ifdef __cplusplus +} +#endif /* __cplusplus */ + + +#endif /* __GTK_CURVE_H__ */