PIKApp/app/core/pikacurve.c

1297 lines
36 KiB
C
Raw Permalink Normal View History

2023-09-26 00:35:21 +02:00
/* PIKA - Photo and Image Kooker Application
* a rebranding of The GNU Image Manipulation Program (created with heckimp)
* A derived work which may be trivial. However, any changes may be (C)2023 by Aldercone Studio
*
* Original copyright, applying to most contents (license remains unchanged):
* Copyright (C) 1995 Spencer Kimball and Peter Mattis
*
* This program is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 3 of the License, or
* (at your option) any later version.
*
* This program 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 General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program. If not, see <https://www.gnu.org/licenses/>.
*/
#include "config.h"
#include <stdlib.h>
#include <string.h> /* memcmp */
#include <gdk-pixbuf/gdk-pixbuf.h>
#include <gegl.h>
#include "libpikabase/pikabase.h"
#include "libpikamath/pikamath.h"
#include "libpikaconfig/pikaconfig.h"
#include "core-types.h"
#include "pikacurve.h"
#include "pikacurve-load.h"
#include "pikacurve-save.h"
#include "pikaparamspecs.h"
#include "pika-intl.h"
#define EPSILON 1e-6
enum
{
PROP_0,
PROP_CURVE_TYPE,
PROP_N_POINTS,
PROP_POINTS,
PROP_POINT_TYPES,
PROP_N_SAMPLES,
PROP_SAMPLES,
N_PROPS
};
static GParamSpec *obj_props[N_PROPS] = { NULL, };
/* local function prototypes */
static void pika_curve_config_iface_init (PikaConfigInterface *iface);
static void pika_curve_finalize (GObject *object);
static void pika_curve_set_property (GObject *object,
guint property_id,
const GValue *value,
GParamSpec *pspec);
static void pika_curve_get_property (GObject *object,
guint property_id,
GValue *value,
GParamSpec *pspec);
static gint64 pika_curve_get_memsize (PikaObject *object,
gint64 *gui_size);
static void pika_curve_get_preview_size (PikaViewable *viewable,
gint size,
gboolean popup,
gboolean dot_for_dot,
gint *width,
gint *height);
static gboolean pika_curve_get_popup_size (PikaViewable *viewable,
gint width,
gint height,
gboolean dot_for_dot,
gint *popup_width,
gint *popup_height);
static PikaTempBuf * pika_curve_get_new_preview (PikaViewable *viewable,
PikaContext *context,
gint width,
gint height);
static gchar * pika_curve_get_description (PikaViewable *viewable,
gchar **tooltip);
static void pika_curve_dirty (PikaData *data);
static const gchar * pika_curve_get_extension (PikaData *data);
static void pika_curve_data_copy (PikaData *data,
PikaData *src_data);
static gboolean pika_curve_serialize (PikaConfig *config,
PikaConfigWriter *writer,
gpointer data);
static gboolean pika_curve_deserialize (PikaConfig *config,
GScanner *scanner,
gint nest_level,
gpointer data);
static gboolean pika_curve_equal (PikaConfig *a,
PikaConfig *b);
static void _pika_curve_reset (PikaConfig *config);
static gboolean pika_curve_config_copy (PikaConfig *src,
PikaConfig *dest,
GParamFlags flags);
static void pika_curve_calculate (PikaCurve *curve);
static void pika_curve_plot (PikaCurve *curve,
gint p1,
gint p2,
gint p3,
gint p4);
G_DEFINE_TYPE_WITH_CODE (PikaCurve, pika_curve, PIKA_TYPE_DATA,
G_IMPLEMENT_INTERFACE (PIKA_TYPE_CONFIG,
pika_curve_config_iface_init))
#define parent_class pika_curve_parent_class
static void
pika_curve_class_init (PikaCurveClass *klass)
{
GObjectClass *object_class = G_OBJECT_CLASS (klass);
PikaObjectClass *pika_object_class = PIKA_OBJECT_CLASS (klass);
PikaViewableClass *viewable_class = PIKA_VIEWABLE_CLASS (klass);
PikaDataClass *data_class = PIKA_DATA_CLASS (klass);
GParamSpec *array_spec;
object_class->finalize = pika_curve_finalize;
object_class->set_property = pika_curve_set_property;
object_class->get_property = pika_curve_get_property;
pika_object_class->get_memsize = pika_curve_get_memsize;
viewable_class->default_icon_name = "FIXME icon name";
viewable_class->get_preview_size = pika_curve_get_preview_size;
viewable_class->get_popup_size = pika_curve_get_popup_size;
viewable_class->get_new_preview = pika_curve_get_new_preview;
viewable_class->get_description = pika_curve_get_description;
data_class->dirty = pika_curve_dirty;
data_class->save = pika_curve_save;
data_class->get_extension = pika_curve_get_extension;
data_class->copy = pika_curve_data_copy;
obj_props[PROP_CURVE_TYPE] =
g_param_spec_enum ("curve-type",
"Curve Type",
"The curve type",
PIKA_TYPE_CURVE_TYPE,
PIKA_CURVE_SMOOTH,
PIKA_CONFIG_PARAM_FLAGS);
obj_props[PROP_N_POINTS] =
g_param_spec_int ("n-points",
"Number of Points",
"The number of points",
0, G_MAXINT, 0,
/* for backward compatibility */
PIKA_CONFIG_PARAM_IGNORE | PIKA_CONFIG_PARAM_FLAGS);
array_spec = g_param_spec_double ("point", NULL, NULL,
-1.0, 1.0, 0.0, PIKA_PARAM_READWRITE);
obj_props[PROP_POINTS] =
pika_param_spec_value_array ("points",
NULL, NULL,
array_spec,
PIKA_CONFIG_PARAM_FLAGS);
array_spec = g_param_spec_enum ("point-type", NULL, NULL,
PIKA_TYPE_CURVE_POINT_TYPE,
PIKA_CURVE_POINT_SMOOTH,
PIKA_PARAM_READWRITE);
obj_props[PROP_POINT_TYPES] =
pika_param_spec_value_array ("point-types",
NULL, NULL,
array_spec,
PIKA_CONFIG_PARAM_FLAGS);
obj_props[PROP_N_SAMPLES] =
g_param_spec_int ("n-samples",
"Number of Samples",
"The number of samples",
256, 256, 256,
PIKA_CONFIG_PARAM_FLAGS);
array_spec = g_param_spec_double ("sample", NULL, NULL,
0.0, 1.0, 0.0, PIKA_PARAM_READWRITE);
obj_props[PROP_SAMPLES] =
pika_param_spec_value_array ("samples",
NULL, NULL,
array_spec,
PIKA_CONFIG_PARAM_FLAGS);
g_object_class_install_properties (object_class, N_PROPS, obj_props);
}
static void
pika_curve_config_iface_init (PikaConfigInterface *iface)
{
iface->serialize = pika_curve_serialize;
iface->deserialize = pika_curve_deserialize;
iface->equal = pika_curve_equal;
iface->reset = _pika_curve_reset;
iface->copy = pika_curve_config_copy;
}
static void
pika_curve_init (PikaCurve *curve)
{
curve->n_points = 0;
curve->points = NULL;
curve->n_samples = 0;
curve->samples = NULL;
curve->identity = FALSE;
}
static void
pika_curve_finalize (GObject *object)
{
PikaCurve *curve = PIKA_CURVE (object);
g_clear_pointer (&curve->points, g_free);
curve->n_points = 0;
g_clear_pointer (&curve->samples, g_free);
curve->n_samples = 0;
G_OBJECT_CLASS (parent_class)->finalize (object);
}
static void
pika_curve_set_property (GObject *object,
guint property_id,
const GValue *value,
GParamSpec *pspec)
{
PikaCurve *curve = PIKA_CURVE (object);
switch (property_id)
{
case PROP_CURVE_TYPE:
pika_curve_set_curve_type (curve, g_value_get_enum (value));
break;
case PROP_N_POINTS:
/* ignored */
break;
case PROP_POINTS:
{
PikaValueArray *array = g_value_get_boxed (value);
PikaCurvePoint *points;
gint length;
gint n_points;
gint i;
if (! array)
{
pika_curve_clear_points (curve);
break;
}
length = pika_value_array_length (array) / 2;
n_points = 0;
points = g_new0 (PikaCurvePoint, length);
for (i = 0; i < length; i++)
{
GValue *x = pika_value_array_index (array, i * 2);
GValue *y = pika_value_array_index (array, i * 2 + 1);
/* for backward compatibility */
if (g_value_get_double (x) < 0.0)
continue;
points[n_points].x = CLAMP (g_value_get_double (x), 0.0, 1.0);
points[n_points].y = CLAMP (g_value_get_double (y), 0.0, 1.0);
if (n_points > 0)
{
points[n_points].x = MAX (points[n_points].x,
points[n_points - 1].x);
}
if (n_points < curve->n_points)
points[n_points].type = curve->points[n_points].type;
else
points[n_points].type = PIKA_CURVE_POINT_SMOOTH;
n_points++;
}
g_free (curve->points);
curve->n_points = n_points;
curve->points = points;
g_object_notify_by_pspec (object, obj_props[PROP_N_POINTS]);
g_object_notify_by_pspec (object, obj_props[PROP_POINT_TYPES]);
}
break;
case PROP_POINT_TYPES:
{
PikaValueArray *array = g_value_get_boxed (value);
PikaCurvePoint *points;
gint length;
gdouble x = 0.0;
gdouble y = 0.0;
gint i;
if (! array)
{
pika_curve_clear_points (curve);
break;
}
length = pika_value_array_length (array);
points = g_new0 (PikaCurvePoint, length);
for (i = 0; i < length; i++)
{
GValue *type = pika_value_array_index (array, i);
points[i].type = g_value_get_enum (type);
if (i < curve->n_points)
{
x = curve->points[i].x;
y = curve->points[i].y;
}
points[i].x = x;
points[i].y = y;
}
g_free (curve->points);
curve->n_points = length;
curve->points = points;
g_object_notify_by_pspec (object, obj_props[PROP_N_POINTS]);
g_object_notify_by_pspec (object, obj_props[PROP_POINTS]);
}
break;
case PROP_N_SAMPLES:
pika_curve_set_n_samples (curve, g_value_get_int (value));
break;
case PROP_SAMPLES:
{
PikaValueArray *array = g_value_get_boxed (value);
gint length;
gint i;
if (! array)
break;
length = pika_value_array_length (array);
for (i = 0; i < curve->n_samples && i < length; i++)
{
GValue *v = pika_value_array_index (array, i);
curve->samples[i] = CLAMP (g_value_get_double (v), 0.0, 1.0);
}
}
break;
default:
G_OBJECT_WARN_INVALID_PROPERTY_ID (object, property_id, pspec);
break;
}
}
static void
pika_curve_get_property (GObject *object,
guint property_id,
GValue *value,
GParamSpec *pspec)
{
PikaCurve *curve = PIKA_CURVE (object);
switch (property_id)
{
case PROP_CURVE_TYPE:
g_value_set_enum (value, curve->curve_type);
break;
case PROP_N_POINTS:
g_value_set_int (value, curve->n_points);
break;
case PROP_POINTS:
{
PikaValueArray *array = pika_value_array_new (curve->n_points * 2);
GValue v = G_VALUE_INIT;
gint i;
g_value_init (&v, G_TYPE_DOUBLE);
for (i = 0; i < curve->n_points; i++)
{
g_value_set_double (&v, curve->points[i].x);
pika_value_array_append (array, &v);
g_value_set_double (&v, curve->points[i].y);
pika_value_array_append (array, &v);
}
g_value_unset (&v);
g_value_take_boxed (value, array);
}
break;
case PROP_POINT_TYPES:
{
PikaValueArray *array = pika_value_array_new (curve->n_points);
GValue v = G_VALUE_INIT;
gint i;
g_value_init (&v, PIKA_TYPE_CURVE_POINT_TYPE);
for (i = 0; i < curve->n_points; i++)
{
g_value_set_enum (&v, curve->points[i].type);
pika_value_array_append (array, &v);
}
g_value_unset (&v);
g_value_take_boxed (value, array);
}
break;
case PROP_N_SAMPLES:
g_value_set_int (value, curve->n_samples);
break;
case PROP_SAMPLES:
{
PikaValueArray *array = pika_value_array_new (curve->n_samples);
GValue v = G_VALUE_INIT;
gint i;
g_value_init (&v, G_TYPE_DOUBLE);
for (i = 0; i < curve->n_samples; i++)
{
g_value_set_double (&v, curve->samples[i]);
pika_value_array_append (array, &v);
}
g_value_unset (&v);
g_value_take_boxed (value, array);
}
break;
default:
G_OBJECT_WARN_INVALID_PROPERTY_ID (object, property_id, pspec);
break;
}
}
static gint64
pika_curve_get_memsize (PikaObject *object,
gint64 *gui_size)
{
PikaCurve *curve = PIKA_CURVE (object);
gint64 memsize = 0;
memsize += curve->n_points * sizeof (PikaCurvePoint);
memsize += curve->n_samples * sizeof (gdouble);
return memsize + PIKA_OBJECT_CLASS (parent_class)->get_memsize (object,
gui_size);
}
static void
pika_curve_get_preview_size (PikaViewable *viewable,
gint size,
gboolean popup,
gboolean dot_for_dot,
gint *width,
gint *height)
{
*width = size;
*height = size;
}
static gboolean
pika_curve_get_popup_size (PikaViewable *viewable,
gint width,
gint height,
gboolean dot_for_dot,
gint *popup_width,
gint *popup_height)
{
*popup_width = width * 2;
*popup_height = height * 2;
return TRUE;
}
static PikaTempBuf *
pika_curve_get_new_preview (PikaViewable *viewable,
PikaContext *context,
gint width,
gint height)
{
return NULL;
}
static gchar *
pika_curve_get_description (PikaViewable *viewable,
gchar **tooltip)
{
PikaCurve *curve = PIKA_CURVE (viewable);
return g_strdup_printf ("%s", pika_object_get_name (curve));
}
static void
pika_curve_dirty (PikaData *data)
{
PikaCurve *curve = PIKA_CURVE (data);
curve->identity = FALSE;
pika_curve_calculate (curve);
PIKA_DATA_CLASS (parent_class)->dirty (data);
}
static const gchar *
pika_curve_get_extension (PikaData *data)
{
return PIKA_CURVE_FILE_EXTENSION;
}
static void
pika_curve_data_copy (PikaData *data,
PikaData *src_data)
{
pika_data_freeze (data);
pika_config_copy (PIKA_CONFIG (src_data),
PIKA_CONFIG (data), 0);
pika_data_thaw (data);
}
static gboolean
pika_curve_serialize (PikaConfig *config,
PikaConfigWriter *writer,
gpointer data)
{
return pika_config_serialize_properties (config, writer);
}
static gboolean
pika_curve_deserialize (PikaConfig *config,
GScanner *scanner,
gint nest_level,
gpointer data)
{
gboolean success;
success = pika_config_deserialize_properties (config, scanner, nest_level);
PIKA_CURVE (config)->identity = FALSE;
return success;
}
static gboolean
pika_curve_equal (PikaConfig *a,
PikaConfig *b)
{
PikaCurve *a_curve = PIKA_CURVE (a);
PikaCurve *b_curve = PIKA_CURVE (b);
if (a_curve->curve_type != b_curve->curve_type)
return FALSE;
if (a_curve->n_points != b_curve->n_points ||
memcmp (a_curve->points, b_curve->points,
sizeof (PikaCurvePoint) * a_curve->n_points))
{
return FALSE;
}
if (a_curve->n_samples != b_curve->n_samples ||
memcmp (a_curve->samples, b_curve->samples,
sizeof (gdouble) * a_curve->n_samples))
{
return FALSE;
}
return TRUE;
}
static void
_pika_curve_reset (PikaConfig *config)
{
pika_curve_reset (PIKA_CURVE (config), TRUE);
}
static gboolean
pika_curve_config_copy (PikaConfig *src,
PikaConfig *dest,
GParamFlags flags)
{
PikaCurve *src_curve = PIKA_CURVE (src);
PikaCurve *dest_curve = PIKA_CURVE (dest);
/* make sure the curve type is copied *before* the points, so that we don't
* overwrite the copied points when changing the type
*/
dest_curve->curve_type = src_curve->curve_type;
pika_config_sync (G_OBJECT (src), G_OBJECT (dest), flags);
dest_curve->identity = src_curve->identity;
pika_data_dirty (PIKA_DATA (dest));
return TRUE;
}
/* public functions */
PikaData *
pika_curve_new (const gchar *name)
{
g_return_val_if_fail (name != NULL, NULL);
g_return_val_if_fail (*name != '\0', NULL);
return g_object_new (PIKA_TYPE_CURVE,
"name", name,
NULL);
}
PikaData *
pika_curve_get_standard (void)
{
static PikaData *standard_curve = NULL;
if (! standard_curve)
{
standard_curve = pika_curve_new ("Standard");
pika_data_clean (standard_curve);
pika_data_make_internal (standard_curve,
"pika-curve-standard");
g_object_ref (standard_curve);
}
return standard_curve;
}
void
pika_curve_reset (PikaCurve *curve,
gboolean reset_type)
{
gint i;
g_return_if_fail (PIKA_IS_CURVE (curve));
g_object_freeze_notify (G_OBJECT (curve));
for (i = 0; i < curve->n_samples; i++)
curve->samples[i] = (gdouble) i / (gdouble) (curve->n_samples - 1);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_SAMPLES]);
g_free (curve->points);
curve->n_points = 2;
curve->points = g_new0 (PikaCurvePoint, 2);
curve->points[0].x = 0.0;
curve->points[0].y = 0.0;
curve->points[0].type = PIKA_CURVE_POINT_SMOOTH;
curve->points[1].x = 1.0;
curve->points[1].y = 1.0;
curve->points[1].type = PIKA_CURVE_POINT_SMOOTH;
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_N_POINTS]);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_POINTS]);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_POINT_TYPES]);
if (reset_type)
{
curve->curve_type = PIKA_CURVE_SMOOTH;
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_CURVE_TYPE]);
}
curve->identity = TRUE;
g_object_thaw_notify (G_OBJECT (curve));
pika_data_dirty (PIKA_DATA (curve));
}
void
pika_curve_set_curve_type (PikaCurve *curve,
PikaCurveType curve_type)
{
g_return_if_fail (PIKA_IS_CURVE (curve));
if (curve->curve_type != curve_type)
{
pika_data_freeze (PIKA_DATA (curve));
g_object_freeze_notify (G_OBJECT (curve));
curve->curve_type = curve_type;
if (curve_type == PIKA_CURVE_SMOOTH)
{
gint i;
g_free (curve->points);
/* pick some points from the curve and make them control
* points
*/
curve->n_points = 9;
curve->points = g_new0 (PikaCurvePoint, 9);
for (i = 0; i < curve->n_points; i++)
{
gint sample = i * (curve->n_samples - 1) / (curve->n_points - 1);
curve->points[i].x = (gdouble) sample /
(gdouble) (curve->n_samples - 1);
curve->points[i].y = curve->samples[sample];
curve->points[i].type = PIKA_CURVE_POINT_SMOOTH;
}
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_N_POINTS]);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_POINTS]);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_POINT_TYPES]);
}
else
{
pika_curve_clear_points (curve);
}
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_CURVE_TYPE]);
g_object_thaw_notify (G_OBJECT (curve));
pika_data_thaw (PIKA_DATA (curve));
}
}
PikaCurveType
pika_curve_get_curve_type (PikaCurve *curve)
{
g_return_val_if_fail (PIKA_IS_CURVE (curve), PIKA_CURVE_SMOOTH);
return curve->curve_type;
}
gint
pika_curve_get_n_points (PikaCurve *curve)
{
g_return_val_if_fail (PIKA_IS_CURVE (curve), 0);
return curve->n_points;
}
void
pika_curve_set_n_samples (PikaCurve *curve,
gint n_samples)
{
g_return_if_fail (PIKA_IS_CURVE (curve));
g_return_if_fail (n_samples >= 256);
g_return_if_fail (n_samples <= 4096);
if (n_samples != curve->n_samples)
{
gint i;
g_object_freeze_notify (G_OBJECT (curve));
curve->n_samples = n_samples;
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_N_SAMPLES]);
curve->samples = g_renew (gdouble, curve->samples, curve->n_samples);
for (i = 0; i < curve->n_samples; i++)
curve->samples[i] = (gdouble) i / (gdouble) (curve->n_samples - 1);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_SAMPLES]);
if (curve->curve_type == PIKA_CURVE_FREE)
curve->identity = TRUE;
g_object_thaw_notify (G_OBJECT (curve));
}
}
gint
pika_curve_get_n_samples (PikaCurve *curve)
{
g_return_val_if_fail (PIKA_IS_CURVE (curve), 0);
return curve->n_samples;
}
gint
pika_curve_get_point_at (PikaCurve *curve,
gdouble x)
{
gint closest_point = -1;
gdouble distance = EPSILON;
gint i;
g_return_val_if_fail (PIKA_IS_CURVE (curve), -1);
for (i = 0; i < curve->n_points; i++)
{
gdouble point_distance;
point_distance = fabs (x - curve->points[i].x);
if (point_distance <= distance)
{
closest_point = i;
distance = point_distance;
}
}
return closest_point;
}
gint
pika_curve_get_closest_point (PikaCurve *curve,
gdouble x,
gdouble y,
gdouble max_distance)
{
gint closest_point = -1;
gdouble distance2 = G_MAXDOUBLE;
gint i;
g_return_val_if_fail (PIKA_IS_CURVE (curve), -1);
if (max_distance >= 0.0)
distance2 = SQR (max_distance);
for (i = curve->n_points - 1; i >= 0; i--)
{
gdouble point_distance2;
point_distance2 = SQR (x - curve->points[i].x) +
SQR (y - curve->points[i].y);
if (point_distance2 <= distance2)
{
closest_point = i;
distance2 = point_distance2;
}
}
return closest_point;
}
gint
pika_curve_add_point (PikaCurve *curve,
gdouble x,
gdouble y)
{
PikaCurvePoint *points;
gint point;
g_return_val_if_fail (PIKA_IS_CURVE (curve), -1);
if (curve->curve_type == PIKA_CURVE_FREE)
return -1;
x = CLAMP (x, 0.0, 1.0);
y = CLAMP (y, 0.0, 1.0);
for (point = 0; point < curve->n_points; point++)
{
if (curve->points[point].x > x)
break;
}
points = g_new0 (PikaCurvePoint, curve->n_points + 1);
memcpy (points, curve->points,
point * sizeof (PikaCurvePoint));
memcpy (points + point + 1, curve->points + point,
(curve->n_points - point) * sizeof (PikaCurvePoint));
points[point].x = x;
points[point].y = y;
points[point].type = PIKA_CURVE_POINT_SMOOTH;
g_free (curve->points);
curve->n_points++;
curve->points = points;
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_N_POINTS]);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_POINTS]);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_POINT_TYPES]);
pika_data_dirty (PIKA_DATA (curve));
return point;
}
void
pika_curve_delete_point (PikaCurve *curve,
gint point)
{
PikaCurvePoint *points;
g_return_if_fail (PIKA_IS_CURVE (curve));
g_return_if_fail (point >= 0 && point < curve->n_points);
points = g_new0 (PikaCurvePoint, curve->n_points - 1);
memcpy (points, curve->points,
point * sizeof (PikaCurvePoint));
memcpy (points + point, curve->points + point + 1,
(curve->n_points - point - 1) * sizeof (PikaCurvePoint));
g_free (curve->points);
curve->n_points--;
curve->points = points;
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_N_POINTS]);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_POINTS]);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_POINT_TYPES]);
pika_data_dirty (PIKA_DATA (curve));
}
void
pika_curve_set_point (PikaCurve *curve,
gint point,
gdouble x,
gdouble y)
{
g_return_if_fail (PIKA_IS_CURVE (curve));
g_return_if_fail (point >= 0 && point < curve->n_points);
curve->points[point].x = CLAMP (x, 0.0, 1.0);
curve->points[point].y = CLAMP (y, 0.0, 1.0);
if (point > 0)
curve->points[point].x = MAX (x, curve->points[point - 1].x);
if (point < curve->n_points - 1)
curve->points[point].x = MIN (x, curve->points[point + 1].x);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_POINTS]);
pika_data_dirty (PIKA_DATA (curve));
}
void
pika_curve_move_point (PikaCurve *curve,
gint point,
gdouble y)
{
g_return_if_fail (PIKA_IS_CURVE (curve));
g_return_if_fail (point >= 0 && point < curve->n_points);
curve->points[point].y = CLAMP (y, 0.0, 1.0);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_POINTS]);
pika_data_dirty (PIKA_DATA (curve));
}
void
pika_curve_get_point (PikaCurve *curve,
gint point,
gdouble *x,
gdouble *y)
{
g_return_if_fail (PIKA_IS_CURVE (curve));
g_return_if_fail (point >= 0 && point < curve->n_points);
if (x) *x = curve->points[point].x;
if (y) *y = curve->points[point].y;
}
void
pika_curve_set_point_type (PikaCurve *curve,
gint point,
PikaCurvePointType type)
{
g_return_if_fail (PIKA_IS_CURVE (curve));
g_return_if_fail (point >= 0 && point < curve->n_points);
curve->points[point].type = type;
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_POINT_TYPES]);
pika_data_dirty (PIKA_DATA (curve));
}
PikaCurvePointType
pika_curve_get_point_type (PikaCurve *curve,
gint point)
{
g_return_val_if_fail (PIKA_IS_CURVE (curve), PIKA_CURVE_POINT_SMOOTH);
g_return_val_if_fail (point >= 0 && point < curve->n_points, PIKA_CURVE_POINT_SMOOTH);
return curve->points[point].type;
}
void
pika_curve_clear_points (PikaCurve *curve)
{
g_return_if_fail (PIKA_IS_CURVE (curve));
if (curve->points)
{
g_clear_pointer (&curve->points, g_free);
curve->n_points = 0;
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_N_POINTS]);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_POINTS]);
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_POINT_TYPES]);
pika_data_dirty (PIKA_DATA (curve));
}
}
void
pika_curve_set_curve (PikaCurve *curve,
gdouble x,
gdouble y)
{
g_return_if_fail (PIKA_IS_CURVE (curve));
g_return_if_fail (x >= 0 && x <= 1.0);
g_return_if_fail (y >= 0 && y <= 1.0);
if (curve->curve_type == PIKA_CURVE_SMOOTH)
return;
curve->samples[ROUND (x * (gdouble) (curve->n_samples - 1))] = y;
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_SAMPLES]);
pika_data_dirty (PIKA_DATA (curve));
}
/**
* pika_curve_is_identity:
* @curve: a #PikaCurve object
*
* If this function returns %TRUE, then the curve maps each value to
* itself. If it returns %FALSE, then this assumption can not be made.
*
* Returns: %TRUE if the curve is an identity mapping, %FALSE otherwise.
**/
gboolean
pika_curve_is_identity (PikaCurve *curve)
{
g_return_val_if_fail (PIKA_IS_CURVE (curve), FALSE);
return curve->identity;
}
void
pika_curve_get_uchar (PikaCurve *curve,
gint n_samples,
guchar *samples)
{
gint i;
g_return_if_fail (PIKA_IS_CURVE (curve));
/* FIXME: support n_samples != curve->n_samples */
g_return_if_fail (n_samples == curve->n_samples);
g_return_if_fail (samples != NULL);
for (i = 0; i < curve->n_samples; i++)
samples[i] = curve->samples[i] * 255.999;
}
/* private functions */
static void
pika_curve_calculate (PikaCurve *curve)
{
gint i;
gint p1, p2, p3, p4;
if (pika_data_is_frozen (PIKA_DATA (curve)))
return;
switch (curve->curve_type)
{
case PIKA_CURVE_SMOOTH:
/* Initialize boundary curve points */
if (curve->n_points > 0)
{
PikaCurvePoint point;
gint boundary;
point = curve->points[0];
boundary = ROUND (point.x * (gdouble) (curve->n_samples - 1));
for (i = 0; i < boundary; i++)
curve->samples[i] = point.y;
point = curve->points[curve->n_points - 1];
boundary = ROUND (point.x * (gdouble) (curve->n_samples - 1));
for (i = boundary; i < curve->n_samples; i++)
curve->samples[i] = point.y;
}
for (i = 0; i < curve->n_points - 1; i++)
{
p1 = MAX (i - 1, 0);
p2 = i;
p3 = i + 1;
p4 = MIN (i + 2, curve->n_points - 1);
if (curve->points[p2].type == PIKA_CURVE_POINT_CORNER)
p1 = p2;
if (curve->points[p3].type == PIKA_CURVE_POINT_CORNER)
p4 = p3;
pika_curve_plot (curve, p1, p2, p3, p4);
}
/* ensure that the control points are used exactly */
for (i = 0; i < curve->n_points; i++)
{
gdouble x = curve->points[i].x;
gdouble y = curve->points[i].y;
curve->samples[ROUND (x * (gdouble) (curve->n_samples - 1))] = y;
}
g_object_notify_by_pspec (G_OBJECT (curve), obj_props[PROP_SAMPLES]);
break;
case PIKA_CURVE_FREE:
break;
}
}
/*
* This function calculates the curve values between the control points
* p2 and p3, taking the potentially existing neighbors p1 and p4 into
* account.
*
* This function uses a cubic bezier curve for the individual segments and
* calculates the necessary intermediate control points depending on the
* neighbor curve control points.
*/
static void
pika_curve_plot (PikaCurve *curve,
gint p1,
gint p2,
gint p3,
gint p4)
{
gint i;
gdouble x0, x3;
gdouble y0, y1, y2, y3;
gdouble dx, dy;
gdouble slope;
/* the outer control points for the bezier curve. */
x0 = curve->points[p2].x;
y0 = curve->points[p2].y;
x3 = curve->points[p3].x;
y3 = curve->points[p3].y;
/*
* the x values of the inner control points are fixed at
* x1 = 2/3*x0 + 1/3*x3 and x2 = 1/3*x0 + 2/3*x3
* this ensures that the x values increase linearly with the
* parameter t and enables us to skip the calculation of the x
* values altogether - just calculate y(t) evenly spaced.
*/
dx = x3 - x0;
dy = y3 - y0;
if (dx <= EPSILON)
{
gint index;
index = ROUND (x0 * (gdouble) (curve->n_samples - 1));
curve->samples[index] = y3;
return;
}
if (p1 == p2 && p3 == p4)
{
/* No information about the neighbors,
* calculate y1 and y2 to get a straight line
*/
y1 = y0 + dy / 3.0;
y2 = y0 + dy * 2.0 / 3.0;
}
else if (p1 == p2 && p3 != p4)
{
/* only the right neighbor is available. Make the tangent at the
* right endpoint parallel to the line between the left endpoint
* and the right neighbor. Then point the tangent at the left towards
* the control handle of the right tangent, to ensure that the curve
* does not have an inflection point.
*/
slope = (curve->points[p4].y - y0) / (curve->points[p4].x - x0);
y2 = y3 - slope * dx / 3.0;
y1 = y0 + (y2 - y0) / 2.0;
}
else if (p1 != p2 && p3 == p4)
{
/* see previous case */
slope = (y3 - curve->points[p1].y) / (x3 - curve->points[p1].x);
y1 = y0 + slope * dx / 3.0;
y2 = y3 + (y1 - y3) / 2.0;
}
else /* (p1 != p2 && p3 != p4) */
{
/* Both neighbors are available. Make the tangents at the endpoints
* parallel to the line between the opposite endpoint and the adjacent
* neighbor.
*/
slope = (y3 - curve->points[p1].y) / (x3 - curve->points[p1].x);
y1 = y0 + slope * dx / 3.0;
slope = (curve->points[p4].y - y0) / (curve->points[p4].x - x0);
y2 = y3 - slope * dx / 3.0;
}
/*
* finally calculate the y(t) values for the given bezier values. We can
* use homogeneously distributed values for t, since x(t) increases linearly.
*/
for (i = 0; i <= ROUND (dx * (gdouble) (curve->n_samples - 1)); i++)
{
gdouble y, t;
gint index;
t = i / dx / (gdouble) (curve->n_samples - 1);
y = y0 * (1-t) * (1-t) * (1-t) +
3 * y1 * (1-t) * (1-t) * t +
3 * y2 * (1-t) * t * t +
y3 * t * t * t;
index = i + ROUND (x0 * (gdouble) (curve->n_samples - 1));
if (index < curve->n_samples)
curve->samples[index] = CLAMP (y, 0.0, 1.0);
}
}