/* Eye Of MATE -- Affine Transformations
 *
 * Copyright (C) 2003-2009 The Free Software Foundation
 *
 * Portions based on code from libart_lgpl by Raph Levien.
 *
 * 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 2 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, write to the Free Software Foundation,
 * Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
 */

#ifdef HAVE_CONFIG_H
#include <config.h>
#endif

#include <time.h>
#include <stdlib.h>
#include <math.h>
#include <gtk/gtk.h>
#include <cairo/cairo.h>

#include "eom-transform.h"
#include "eom-jobs.h"

/* The number of progress updates per transformation */
#define EOM_TRANSFORM_N_PROG_UPDATES 20

struct _EomTransformPrivate {
	cairo_matrix_t affine;
};

typedef struct {
	gdouble x;
	gdouble y;
} EomPoint;

/* Convert degrees into radians */
#define EOM_DEG_TO_RAD(degree) ((degree) * (G_PI/180.0)) 

#define EOM_TRANSFORM_GET_PRIVATE(object) \
	(G_TYPE_INSTANCE_GET_PRIVATE ((object), EOM_TYPE_TRANSFORM, EomTransformPrivate))

G_DEFINE_TYPE (EomTransform, eom_transform, G_TYPE_OBJECT)

static void
eom_transform_init (EomTransform *trans)
{
	trans->priv = EOM_TRANSFORM_GET_PRIVATE (trans);
}

static void
eom_transform_class_init (EomTransformClass *klass)
{
	g_type_class_add_private (klass, sizeof (EomTransformPrivate));
}

GdkPixbuf*
eom_transform_apply (EomTransform *trans, GdkPixbuf *pixbuf, EomJob *job)
{
	EomPoint dest_top_left;
	EomPoint dest_bottom_right;
	EomPoint vertices[4] = { {0, 0}, {1, 0}, {1, 1}, {0, 1} };
	double r_det;
	int inverted [6];
	EomPoint dest;

	int src_width;
	int src_height;
	int src_rowstride;
	int src_n_channels;
	guchar *src_buffer;

	GdkPixbuf *dest_pixbuf;
	int dest_width;
	int dest_height;
	int dest_rowstride;
	int dest_n_channels;
	guchar *dest_buffer;

	guchar *src_pos;
	guchar *dest_pos;
	int dx, dy, sx, sy;
	int i, x, y;

	int progress_delta;

	g_return_val_if_fail (pixbuf != NULL, NULL);

	g_object_ref (pixbuf);

	src_width = gdk_pixbuf_get_width (pixbuf);
	src_height = gdk_pixbuf_get_height (pixbuf);
	src_rowstride = gdk_pixbuf_get_rowstride (pixbuf);
	src_n_channels = gdk_pixbuf_get_n_channels (pixbuf);
	src_buffer = gdk_pixbuf_get_pixels (pixbuf);

	/* find out the dimension of the destination pixbuf */
	dest_top_left.x = 100000;
	dest_top_left.y = 100000;
	dest_bottom_right.x = -100000;
	dest_bottom_right.y = -100000;

	for (i = 0; i < 4; i++) {
		dest.x = vertices[i].x * (src_width - 1);
		dest.y = vertices[i].y * (src_height -1);

		cairo_matrix_transform_point (&trans->priv->affine,
					      &dest.x, &dest.y);

		dest_top_left.x = MIN (dest_top_left.x, dest.x);
		dest_top_left.y = MIN (dest_top_left.y, dest.y);

		dest_bottom_right.x = MAX (dest_bottom_right.x, dest.x);
		dest_bottom_right.y = MAX (dest_bottom_right.y, dest.y);
	}

	/* create the resulting pixbuf */
	dest_width = abs (dest_bottom_right.x - dest_top_left.x + 1);
	dest_height = abs (dest_bottom_right.y - dest_top_left.y + 1);

	dest_pixbuf = gdk_pixbuf_new (GDK_COLORSPACE_RGB,
			       gdk_pixbuf_get_has_alpha (pixbuf),
			       gdk_pixbuf_get_bits_per_sample (pixbuf),
			       dest_width,
			       dest_height);
	dest_rowstride = gdk_pixbuf_get_rowstride (dest_pixbuf);
	dest_n_channels = gdk_pixbuf_get_n_channels (dest_pixbuf);
	dest_buffer = gdk_pixbuf_get_pixels (dest_pixbuf);

	/* invert the matrix so that we can compute the source pixel
	   from the target pixel and convert the values to integer
	   ones (faster!)  FIXME: Maybe we can do some more
	   improvements by using special mmx/3dnow features if
	   available.
	*/
	r_det = 1.0 / (trans->priv->affine.xx * trans->priv->affine.yy - trans->priv->affine.yx * trans->priv->affine.xy);
	inverted[0] =  trans->priv->affine.yy * r_det;
	inverted[1] = -trans->priv->affine.yx * r_det;
	inverted[2] = -trans->priv->affine.xy * r_det;
	inverted[3] =  trans->priv->affine.xx * r_det;
	inverted[4] = -trans->priv->affine.x0 * inverted[0] - trans->priv->affine.y0 * inverted[2];
	inverted[5] = -trans->priv->affine.x0 * inverted[1] - trans->priv->affine.y0 * inverted[3];

	progress_delta = MAX (1, dest_height / EOM_TRANSFORM_N_PROG_UPDATES);

	/*
	 * for every destination pixel (dx,dy) compute the source pixel (sx, sy)
	 * and copy the color values
	 */
	for (y = 0, dy = dest_top_left.y; y < dest_height; y++, dy++) {
		for (x = 0, dx = dest_top_left.x; x < dest_width; x++, dx++) {

			sx = dx * inverted[0] + dy * inverted[2] + inverted[4];
			sy = dx * inverted[1] + dy * inverted[3] + inverted[5];

			if (sx >= 0 && sx < src_width && sy >= 0 && sy < src_height) {
				src_pos  = src_buffer  + sy * src_rowstride  + sx * src_n_channels;
				dest_pos = dest_buffer +  y * dest_rowstride +  x * dest_n_channels;

				for (i = 0; i <  src_n_channels; i++) {
					dest_pos[i] = src_pos[i];
				}
			}
		}

		if (job != NULL && y % progress_delta == 0) {
			gfloat progress;

			progress = (gfloat) (y + 1.0) / (gfloat) dest_height;

			eom_job_set_progress (job, progress);
		}
	}

	g_object_unref (pixbuf);

	if (job != NULL) {
		eom_job_set_progress (job, 1.0);
	}

	return dest_pixbuf;
}

static void
_eom_cairo_matrix_copy (const cairo_matrix_t *src, cairo_matrix_t *dest)
{
	cairo_matrix_init (dest, src->xx, src->yx, src->xy, src->yy, src->x0, src->y0);
}

#define DOUBLE_EQUAL_MAX_DIFF 1e-6
#define DOUBLE_EQUAL(a,b) (fabs (a - b) < DOUBLE_EQUAL_MAX_DIFF)
/* art_affine_equal modified to work with cairo_matrix_t */
static gboolean
_eom_cairo_matrix_equal (const cairo_matrix_t *a, const cairo_matrix_t *b)
{
	return (DOUBLE_EQUAL (a->xx, b->xx) && DOUBLE_EQUAL (a->yx, b->yx) &&
		DOUBLE_EQUAL (a->xy, b->xy) && DOUBLE_EQUAL (a->yy, b->yy) &&
		DOUBLE_EQUAL (a->x0, b->x0) && DOUBLE_EQUAL (a->y0, b->y0) );
}

/* art_affine_flip modified to work with cairo_matrix_t */
static void
_eom_cairo_matrix_flip (cairo_matrix_t *dst, const cairo_matrix_t *src, gboolean horiz, gboolean vert)
{
	dst->xx = horiz ? -src->xx : src->xx;
	dst->yx = horiz ? -src->yx : src->yx;
	dst->xy = vert ? -src->xy : src->xy;
	dst->yy = vert ? -src->yy : src->yy;
	dst->x0 = horiz ? -src->x0 : src->x0;
	dst->y0 = vert ? -src->y0 : src->y0;
}

EomTransform*
eom_transform_reverse (EomTransform *trans)
{
	EomTransform *reverse;

	g_return_val_if_fail (EOM_IS_TRANSFORM (trans), NULL);

	reverse = EOM_TRANSFORM (g_object_new (EOM_TYPE_TRANSFORM, NULL));

	_eom_cairo_matrix_copy (&trans->priv->affine, &reverse->priv->affine);

	g_return_val_if_fail (cairo_matrix_invert (&reverse->priv->affine) == CAIRO_STATUS_SUCCESS, reverse);

	return reverse;
}

EomTransform*
eom_transform_compose (EomTransform *trans, EomTransform *compose)
{
	EomTransform *composition;

	g_return_val_if_fail (EOM_IS_TRANSFORM (trans), NULL);
	g_return_val_if_fail (EOM_IS_TRANSFORM (compose), NULL);

	composition = EOM_TRANSFORM (g_object_new (EOM_TYPE_TRANSFORM, NULL));

	cairo_matrix_multiply (&composition->priv->affine,
			       &trans->priv->affine,
			       &compose->priv->affine);

	return composition;
}

gboolean
eom_transform_is_identity (EomTransform *trans)
{
	static const cairo_matrix_t identity = { 1, 0, 0, 1, 0, 0 };

	g_return_val_if_fail (EOM_IS_TRANSFORM (trans), FALSE);

	return _eom_cairo_matrix_equal (&identity, &trans->priv->affine);
}

EomTransform*
eom_transform_identity_new (void)
{
	EomTransform *trans;

	trans = EOM_TRANSFORM (g_object_new (EOM_TYPE_TRANSFORM, NULL));

	cairo_matrix_init_identity (&trans->priv->affine);

	return trans;
}

EomTransform*
eom_transform_rotate_new (int degree)
{
	EomTransform *trans;

	trans = EOM_TRANSFORM (g_object_new (EOM_TYPE_TRANSFORM, NULL));

	cairo_matrix_init_rotate (&trans->priv->affine, EOM_DEG_TO_RAD(degree));

	return trans;
}

EomTransform*
eom_transform_flip_new   (EomTransformType type)
{
	EomTransform *trans;
	gboolean horiz, vert;

	trans = EOM_TRANSFORM (g_object_new (EOM_TYPE_TRANSFORM, NULL));

	cairo_matrix_init_identity (&trans->priv->affine);

	horiz = (type == EOM_TRANSFORM_FLIP_HORIZONTAL);
	vert = (type == EOM_TRANSFORM_FLIP_VERTICAL);

	_eom_cairo_matrix_flip (&trans->priv->affine,
				&trans->priv->affine,
				horiz, vert);

	return trans;
}

EomTransform*
eom_transform_new (EomTransformType type)
{
	EomTransform *trans = NULL;
	EomTransform *temp1 = NULL, *temp2 = NULL;

	switch (type) {
	case EOM_TRANSFORM_NONE:
		trans = eom_transform_identity_new ();
		break;
	case EOM_TRANSFORM_FLIP_HORIZONTAL:
		trans = eom_transform_flip_new (EOM_TRANSFORM_FLIP_HORIZONTAL);
		break;
	case EOM_TRANSFORM_ROT_180:
		trans = eom_transform_rotate_new (180);
		break;
	case EOM_TRANSFORM_FLIP_VERTICAL:
		trans = eom_transform_flip_new (EOM_TRANSFORM_FLIP_VERTICAL);
		break;
	case EOM_TRANSFORM_TRANSPOSE:
		temp1 = eom_transform_rotate_new (90);
		temp2 = eom_transform_flip_new (EOM_TRANSFORM_FLIP_HORIZONTAL);
		trans = eom_transform_compose (temp1, temp2);
		g_object_unref (temp1);
		g_object_unref (temp2);
		break;
	case EOM_TRANSFORM_ROT_90:
		trans = eom_transform_rotate_new (90);
		break;
	case EOM_TRANSFORM_TRANSVERSE:
		temp1 = eom_transform_rotate_new (90);
		temp2 = eom_transform_flip_new (EOM_TRANSFORM_FLIP_VERTICAL);
		trans = eom_transform_compose (temp1, temp2);
		g_object_unref (temp1);
		g_object_unref (temp2);
		break;
	case EOM_TRANSFORM_ROT_270:
		trans = eom_transform_rotate_new (270);
		break;
	default:
		trans = eom_transform_identity_new ();
		break;
	}

	return trans;
}

EomTransformType
eom_transform_get_transform_type (EomTransform *trans)
{
	cairo_matrix_t affine, a1, a2;
	EomTransformPrivate *priv;

	g_return_val_if_fail (EOM_IS_TRANSFORM (trans), EOM_TRANSFORM_NONE);

	priv = trans->priv;

	cairo_matrix_init_rotate (&affine, EOM_DEG_TO_RAD(90));
	if (_eom_cairo_matrix_equal (&affine, &priv->affine)) {
		return EOM_TRANSFORM_ROT_90;
	}

	cairo_matrix_init_rotate (&affine, EOM_DEG_TO_RAD(180));
	if (_eom_cairo_matrix_equal (&affine, &priv->affine)) {
		return EOM_TRANSFORM_ROT_180;
	}

	cairo_matrix_init_rotate (&affine, EOM_DEG_TO_RAD(270));
	if (_eom_cairo_matrix_equal (&affine, &priv->affine)) {
		return EOM_TRANSFORM_ROT_270;
	}

	cairo_matrix_init_identity (&affine);
	_eom_cairo_matrix_flip (&affine, &affine, TRUE, FALSE);
	if (_eom_cairo_matrix_equal (&affine, &priv->affine)) {
		return EOM_TRANSFORM_FLIP_HORIZONTAL;
	}

	cairo_matrix_init_identity (&affine);
	_eom_cairo_matrix_flip (&affine, &affine, FALSE, TRUE);
	if (_eom_cairo_matrix_equal (&affine, &priv->affine)) {
		return EOM_TRANSFORM_FLIP_VERTICAL;
	}

	cairo_matrix_init_rotate (&a1, EOM_DEG_TO_RAD(90));
	cairo_matrix_init_identity (&a2);
	_eom_cairo_matrix_flip (&a2, &a2, TRUE, FALSE);
	cairo_matrix_multiply(&affine, &a1, &a2);
	if (_eom_cairo_matrix_equal (&affine, &priv->affine)) {
		return EOM_TRANSFORM_TRANSPOSE;
	}

	/* A transversion is a 180° rotation followed by a transposition */
	/* Reuse the transposition from the previous step for this. */
	cairo_matrix_init_rotate (&a1, EOM_DEG_TO_RAD(180));
	cairo_matrix_multiply(&a2, &a1, &affine);
	if (_eom_cairo_matrix_equal (&a2, &priv->affine)) {
		return EOM_TRANSFORM_TRANSVERSE;
	}

	return EOM_TRANSFORM_NONE;
}

gboolean
eom_transform_get_affine (EomTransform *trans, cairo_matrix_t *affine)
{
	g_return_val_if_fail (EOM_IS_TRANSFORM (trans), FALSE);

	_eom_cairo_matrix_copy (&trans->priv->affine, affine);

	return TRUE;
}