/* -*- Mode: c; c-basic-offset: 4; indent-tabs-mode: t; tab-width: 8; -*- */

/* cairo - a vector graphics library with display and print output

 *

 * Copyright © 2002 University of Southern California

 *

 * This library is free software; you can redistribute it and/or

 * modify it either under the terms of the GNU Lesser General Public

 * License version 2.1 as published by the Free Software Foundation

 * (the "LGPL") or, at your option, under the terms of the Mozilla

 * Public License Version 1.1 (the "MPL"). If you do not alter this

 * notice, a recipient may use your version of this file under either

 * the MPL or the LGPL.

 *

 * You should have received a copy of the LGPL along with this library

 * in the file COPYING-LGPL-2.1; if not, write to the Free Software

 * Foundation, Inc., 51 Franklin Street, Suite 500, Boston, MA 02110-1335, USA

 * You should have received a copy of the MPL along with this library

 * in the file COPYING-MPL-1.1

 *

 * The contents of this file are subject to the Mozilla Public License

 * Version 1.1 (the "License"); you may not use this file except in

 * compliance with the License. You may obtain a copy of the License at

 * http://www.mozilla.org/MPL/

 *

 * This software is distributed on an "AS IS" basis, WITHOUT WARRANTY

 * OF ANY KIND, either express or implied. See the LGPL or the MPL for

 * the specific language governing rights and limitations.

 *

 * The Original Code is the cairo graphics library.

 *

 * The Initial Developer of the Original Code is University of Southern

 * California.

 *

 * Contributor(s):

 *  Carl D. Worth <cworth@cworth.org>

 */

#include "cairoint.h"

#include "cairo-boxes-private.h"

#include "cairo-contour-private.h"

#include "cairo-error-private.h"

#define DEBUG_POLYGON 0

#if DEBUG_POLYGON && !NDEBUG

static void

assert_last_edge_is_valid(cairo_polygon_t *polygon,

        const cairo_box_t *limit)

{

    cairo_edge_t *edge;

    cairo_fixed_t x;

    edge = &polygon->edges[polygon->num_edges-1];

    assert (edge->bottom > edge->top);

    assert (edge->top >= limit->p1.y);

    assert (edge->bottom <= limit->p2.y);

    x = _cairo_edge_compute_intersection_x_for_y (&edge->line.p1,

              &edge->line.p2,

              edge->top);

    assert (x >= limit->p1.x);

    assert (x <= limit->p2.x);

    x = _cairo_edge_compute_intersection_x_for_y (&edge->line.p1,

              &edge->line.p2,

              edge->bottom);

    assert (x >= limit->p1.x);

    assert (x <= limit->p2.x);

}

#else

#define assert_last_edge_is_valid(p, l)

#endif

static void

_cairo_polygon_add_edge (cairo_polygon_t *polygon,

       const cairo_point_t *p1,

       const cairo_point_t *p2,

       int dir);

void

_cairo_polygon_limit (cairo_polygon_t *polygon,

         const cairo_box_t *limits,

         int num_limits)

{

    int n;

    polygon->limits = limits;

    polygon->num_limits = num_limits;

    if (polygon->num_limits) {

  polygon->limit = limits[0];

  for (n = 1; n < num_limits; n++) {

      if (limits[n].p1.x < polygon->limit.p1.x)

    polygon->limit.p1.x = limits[n].p1.x;

      if (limits[n].p1.y < polygon->limit.p1.y)

    polygon->limit.p1.y = limits[n].p1.y;

      if (limits[n].p2.x > polygon->limit.p2.x)

    polygon->limit.p2.x = limits[n].p2.x;

      if (limits[n].p2.y > polygon->limit.p2.y)

    polygon->limit.p2.y = limits[n].p2.y;

  }

    }

}

void

_cairo_polygon_limit_to_clip (cairo_polygon_t *polygon,

            const cairo_clip_t *clip)

{

    if (clip)

  _cairo_polygon_limit (polygon, clip->boxes, clip->num_boxes);

    else

  _cairo_polygon_limit (polygon, 0, 0);

}

void

_cairo_polygon_init (cairo_polygon_t *polygon,

         const cairo_box_t *limits,

         int num_limits)

{

    VG (VALGRIND_MAKE_MEM_UNDEFINED (polygon, sizeof (cairo_polygon_t)));

    polygon->status = CAIRO_STATUS_SUCCESS;

    polygon->num_edges = 0;

    polygon->edges = polygon->edges_embedded;

    polygon->edges_size = ARRAY_LENGTH (polygon->edges_embedded);

    polygon->extents.p1.x = polygon->extents.p1.y = INT32_MAX;

    polygon->extents.p2.x = polygon->extents.p2.y = INT32_MIN;

    _cairo_polygon_limit (polygon, limits, num_limits);

}

void

_cairo_polygon_init_with_clip (cairo_polygon_t *polygon,

             const cairo_clip_t *clip)

{

    if (clip)

  _cairo_polygon_init (polygon, clip->boxes, clip->num_boxes);

    else

  _cairo_polygon_init (polygon, 0, 0);

}

cairo_status_t

_cairo_polygon_init_boxes (cairo_polygon_t *polygon,

         const cairo_boxes_t *boxes)

{

    const struct _cairo_boxes_chunk *chunk;

    int i;

    VG (VALGRIND_MAKE_MEM_UNDEFINED (polygon, sizeof (cairo_polygon_t)));

    polygon->status = CAIRO_STATUS_SUCCESS;

    polygon->num_edges = 0;

    polygon->edges = polygon->edges_embedded;

    polygon->edges_size = ARRAY_LENGTH (polygon->edges_embedded);

    if (boxes->num_boxes > ARRAY_LENGTH (polygon->edges_embedded)/2) {

  polygon->edges_size = 2 * boxes->num_boxes;

  polygon->edges = _cairo_malloc_ab (polygon->edges_size,

             2*sizeof(cairo_edge_t));

  if (unlikely (polygon->edges == NULL))

      return polygon->status = _cairo_error (CAIRO_STATUS_NO_MEMORY);

    }

    polygon->extents.p1.x = polygon->extents.p1.y = INT32_MAX;

    polygon->extents.p2.x = polygon->extents.p2.y = INT32_MIN;

    polygon->limits = NULL;

    polygon->num_limits = 0;

    for (chunk = &boxes->chunks; chunk != NULL; chunk = chunk->next) {

  for (i = 0; i < chunk->count; i++) {

      cairo_point_t p1, p2;

      p1 = chunk->base[i].p1;

      p2.x = p1.x;

      p2.y = chunk->base[i].p2.y;

      _cairo_polygon_add_edge (polygon, &p1, &p2, 1);

      p1 = chunk->base[i].p2;

      p2.x = p1.x;

      p2.y = chunk->base[i].p1.y;

      _cairo_polygon_add_edge (polygon, &p1, &p2, 1);

  }

    }

    return polygon->status;

}

cairo_status_t

_cairo_polygon_init_box_array (cairo_polygon_t *polygon,

             cairo_box_t *boxes,

             int num_boxes)

{

    int i;

    VG (VALGRIND_MAKE_MEM_UNDEFINED (polygon, sizeof (cairo_polygon_t)));

    polygon->status = CAIRO_STATUS_SUCCESS;

    polygon->num_edges = 0;

    polygon->edges = polygon->edges_embedded;

    polygon->edges_size = ARRAY_LENGTH (polygon->edges_embedded);

    if (num_boxes > ARRAY_LENGTH (polygon->edges_embedded)/2) {

  polygon->edges_size = 2 * num_boxes;

  polygon->edges = _cairo_malloc_ab (polygon->edges_size,

             2*sizeof(cairo_edge_t));

  if (unlikely (polygon->edges == NULL))

      return polygon->status = _cairo_error (CAIRO_STATUS_NO_MEMORY);

    }

    polygon->extents.p1.x = polygon->extents.p1.y = INT32_MAX;

    polygon->extents.p2.x = polygon->extents.p2.y = INT32_MIN;

    polygon->limits = NULL;

    polygon->num_limits = 0;

    for (i = 0; i < num_boxes; i++) {

  cairo_point_t p1, p2;

  p1 = boxes[i].p1;

  p2.x = p1.x;

  p2.y = boxes[i].p2.y;

  _cairo_polygon_add_edge (polygon, &p1, &p2, 1);

  p1 = boxes[i].p2;

  p2.x = p1.x;

  p2.y = boxes[i].p1.y;

  _cairo_polygon_add_edge (polygon, &p1, &p2, 1);

    }

    return polygon->status;

}

void

_cairo_polygon_fini (cairo_polygon_t *polygon)

{

    if (polygon->edges != polygon->edges_embedded)

  free (polygon->edges);

    VG (VALGRIND_MAKE_MEM_UNDEFINED (polygon, sizeof (cairo_polygon_t)));

}

/* make room for at least one more edge */

static cairo_bool_t

_cairo_polygon_grow (cairo_polygon_t *polygon)

{

    cairo_edge_t *new_edges;

    int old_size = polygon->edges_size;

    int new_size = 4 * old_size;

    if (CAIRO_INJECT_FAULT ()) {

  polygon->status = _cairo_error (CAIRO_STATUS_NO_MEMORY);

  return FALSE;

    }

    if (polygon->edges == polygon->edges_embedded) {

  new_edges = _cairo_malloc_ab (new_size, sizeof (cairo_edge_t));

  if (new_edges != NULL)

      memcpy (new_edges, polygon->edges, old_size * sizeof (cairo_edge_t));

    } else {

  new_edges = _cairo_realloc_ab (polygon->edges,

                           new_size, sizeof (cairo_edge_t));

    }

    if (unlikely (new_edges == NULL)) {

  polygon->status = _cairo_error (CAIRO_STATUS_NO_MEMORY);

  return FALSE;

    }

    polygon->edges = new_edges;

    polygon->edges_size = new_size;

    return TRUE;

}

static void

_add_edge (cairo_polygon_t *polygon,

     const cairo_point_t *p1,

     const cairo_point_t *p2,

     int top, int bottom,

     int dir)

{

    cairo_edge_t *edge;

    assert (top < bottom);

    if (unlikely (polygon->num_edges == polygon->edges_size)) {

  if (! _cairo_polygon_grow (polygon))

      return;

    }

    edge = &polygon->edges[polygon->num_edges++];

    edge->line.p1 = *p1;

    edge->line.p2 = *p2;

    edge->top = top;

    edge->bottom = bottom;

    edge->dir = dir;

    if (top < polygon->extents.p1.y)

  polygon->extents.p1.y = top;

    if (bottom > polygon->extents.p2.y)

  polygon->extents.p2.y = bottom;

    if (p1->x < polygon->extents.p1.x || p1->x > polygon->extents.p2.x) {

  cairo_fixed_t x = p1->x;

  if (top != p1->y)

      x = _cairo_edge_compute_intersection_x_for_y (p1, p2, top);

  if (x < polygon->extents.p1.x)

      polygon->extents.p1.x = x;

  if (x > polygon->extents.p2.x)

      polygon->extents.p2.x = x;

    }

    if (p2->x < polygon->extents.p1.x || p2->x > polygon->extents.p2.x) {

  cairo_fixed_t x = p2->x;

  if (bottom != p2->y)

      x = _cairo_edge_compute_intersection_x_for_y (p1, p2, bottom);

  if (x < polygon->extents.p1.x)

      polygon->extents.p1.x = x;

  if (x > polygon->extents.p2.x)

      polygon->extents.p2.x = x;

    }

}

static void

_add_clipped_edge (cairo_polygon_t *polygon,

       const cairo_point_t *p1,

       const cairo_point_t *p2,

       const int top, const int bottom,

       const int dir)

{

    cairo_point_t bot_left, top_right;

    cairo_fixed_t top_y, bot_y;

    int n;

    for (n = 0; n < polygon->num_limits; n++) {

  const cairo_box_t *limits = &polygon->limits[n];

  cairo_fixed_t pleft, pright;

  if (top >= limits->p2.y)

      continue;

  if (bottom <= limits->p1.y)

      continue;

  bot_left.x = limits->p1.x;

  bot_left.y = limits->p2.y;

  top_right.x = limits->p2.x;

  top_right.y = limits->p1.y;

  /* The useful region */

  top_y = MAX (top, limits->p1.y);

  bot_y = MIN (bottom, limits->p2.y);

  /* The projection of the edge on the horizontal axis */

  pleft = MIN (p1->x, p2->x);

  pright = MAX (p1->x, p2->x);

  if (limits->p1.x <= pleft && pright <= limits->p2.x) {

      /* Projection of the edge completely contained in the box:

       * clip vertically by restricting top and bottom */

      _add_edge (polygon, p1, p2, top_y, bot_y, dir);

      assert_last_edge_is_valid (polygon, limits);

  } else if (pright <= limits->p1.x) {

      /* Projection of the edge to the left of the box:

       * replace with the left side of the box (clipped top/bottom) */

      _add_edge (polygon, &limits->p1, &bot_left, top_y, bot_y, dir);

      assert_last_edge_is_valid (polygon, limits);

  } else if (limits->p2.x <= pleft) {

      /* Projection of the edge to the right of the box:

       * replace with the right side of the box (clipped top/bottom) */

      _add_edge (polygon, &top_right, &limits->p2, top_y, bot_y, dir);

      assert_last_edge_is_valid (polygon, limits);

  } else {

      /* The edge and the box intersect in a generic way */

      cairo_fixed_t left_y, right_y;

      cairo_bool_t top_left_to_bottom_right;

      /*

       * The edge intersects the lines corresponding to the left

       * and right sides of the limit box at left_y and right_y,

       * but we need to add edges for the range from top_y to

       * bot_y.

       *

       * For both intersections, there are three cases:

       *

       *  1) It is outside the vertical range of the limit

       *     box. In this case we can simply further clip the

       *     edge we will be emitting (i.e. restrict its

       *     top/bottom limits to those of the limit box).

       *

       *  2) It is inside the vertical range of the limit

       *     box. In this case, we need to add the vertical edge

       *     connecting the correct vertex to the intersection,

       *     in order to preserve the winding count.

       *

       *  3) It is exactly on the box. In this case, do nothing.

       *

       * These operations restrict the active range (stored in

       * top_y/bot_y) so that the p1-p2 edge is completely

       * inside the box if it is clipped to this vertical range.

       */

      top_left_to_bottom_right = (p1->x <= p2->x) == (p1->y <= p2->y);

      if (top_left_to_bottom_right) {

    if (pleft >= limits->p1.x) {

        left_y = top_y;

    } else {

        left_y = _cairo_edge_compute_intersection_y_for_x (p1, p2,

                       limits->p1.x);

        if (_cairo_edge_compute_intersection_x_for_y (p1, p2, left_y) < limits->p1.x)

      left_y++;

    }

    left_y = MIN (left_y, bot_y);

    if (top_y < left_y) {

        _add_edge (polygon, &limits->p1, &bot_left,

             top_y, left_y, dir);

        assert_last_edge_is_valid (polygon, limits);

        top_y = left_y;

    }

    if (pright <= limits->p2.x) {

        right_y = bot_y;

    } else {

        right_y = _cairo_edge_compute_intersection_y_for_x (p1, p2,

                  limits->p2.x);

        if (_cairo_edge_compute_intersection_x_for_y (p1, p2, right_y) > limits->p2.x)

      right_y--;

    }

    right_y = MAX (right_y, top_y);

    if (bot_y > right_y) {

        _add_edge (polygon, &top_right, &limits->p2,

             right_y, bot_y, dir);

        assert_last_edge_is_valid (polygon, limits);

        bot_y = right_y;

    }

      } else {

    if (pright <= limits->p2.x) {

        right_y = top_y;

    } else {

        right_y = _cairo_edge_compute_intersection_y_for_x (p1, p2,

                  limits->p2.x);

        if (_cairo_edge_compute_intersection_x_for_y (p1, p2, right_y) > limits->p2.x)

      right_y++;

    }

    right_y = MIN (right_y, bot_y);

    if (top_y < right_y) {

        _add_edge (polygon, &top_right, &limits->p2,

             top_y, right_y, dir);

        assert_last_edge_is_valid (polygon, limits);

        top_y = right_y;

    }

    if (pleft >= limits->p1.x) {

        left_y = bot_y;

    } else {

        left_y = _cairo_edge_compute_intersection_y_for_x (p1, p2,

                       limits->p1.x);

        if (_cairo_edge_compute_intersection_x_for_y (p1, p2, left_y) < limits->p1.x)

      left_y--;

    }

    left_y = MAX (left_y, top_y);

    if (bot_y > left_y) {

        _add_edge (polygon, &limits->p1, &bot_left,

             left_y, bot_y, dir);

        assert_last_edge_is_valid (polygon, limits);

        bot_y = left_y;

    }

      }

      if (top_y != bot_y) {

    _add_edge (polygon, p1, p2, top_y, bot_y, dir);

    assert_last_edge_is_valid (polygon, limits);

      }

  }

    }

}

static void

_cairo_polygon_add_edge (cairo_polygon_t *polygon,

       const cairo_point_t *p1,

       const cairo_point_t *p2,

       int dir)

{

    /* drop horizontal edges */

    if (p1->y == p2->y)

  return;

    if (p1->y > p2->y) {

  const cairo_point_t *t;

  t = p1, p1 = p2, p2 = t;

  dir = -dir;

    }

    if (polygon->num_limits) {

  if (p2->y <= polygon->limit.p1.y)

      return;

  if (p1->y >= polygon->limit.p2.y)

      return;

  _add_clipped_edge (polygon, p1, p2, p1->y, p2->y, dir);

    } else

  _add_edge (polygon, p1, p2, p1->y, p2->y, dir);

}

cairo_status_t

_cairo_polygon_add_external_edge (void *polygon,

          const cairo_point_t *p1,

          const cairo_point_t *p2)

{

    _cairo_polygon_add_edge (polygon, p1, p2, 1);

    return _cairo_polygon_status (polygon);

}

cairo_status_t

_cairo_polygon_add_line (cairo_polygon_t *polygon,

       const cairo_line_t *line,

       int top, int bottom,

       int dir)

{

    /* drop horizontal edges */

    if (line->p1.y == line->p2.y)

  return CAIRO_STATUS_SUCCESS;

    if (bottom <= top)

  return CAIRO_STATUS_SUCCESS;

    if (polygon->num_limits) {

  if (line->p2.y <= polygon->limit.p1.y)

      return CAIRO_STATUS_SUCCESS;

  if (line->p1.y >= polygon->limit.p2.y)

      return CAIRO_STATUS_SUCCESS;

  _add_clipped_edge (polygon, &line->p1, &line->p2, top, bottom, dir);

    } else

  _add_edge (polygon, &line->p1, &line->p2, top, bottom, dir);

    return polygon->status;

}

cairo_status_t

_cairo_polygon_add_contour (cairo_polygon_t *polygon,

          const cairo_contour_t *contour)

{

    const struct _cairo_contour_chain *chain;

    const cairo_point_t *prev = NULL;

    int i;

    if (contour->chain.num_points <= 1)

  return CAIRO_INT_STATUS_SUCCESS;

    prev = &contour->chain.points[0];

    for (chain = &contour->chain; chain; chain = chain->next) {

  for (i = 0; i < chain->num_points; i++) {

      _cairo_polygon_add_edge (polygon, prev, &chain->points[i],

             contour->direction);

      prev = &chain->points[i];

  }

    }

    return polygon->status;

}

void

_cairo_polygon_translate (cairo_polygon_t *polygon, int dx, int dy)

{

    int n;

    dx = _cairo_fixed_from_int (dx);

    dy = _cairo_fixed_from_int (dy);

    polygon->extents.p1.x += dx;

    polygon->extents.p2.x += dx;

    polygon->extents.p1.y += dy;

    polygon->extents.p2.y += dy;

    for (n = 0; n < polygon->num_edges; n++) {

  cairo_edge_t *e = &polygon->edges[n];

  e->top += dy;

  e->bottom += dy;

  e->line.p1.x += dx;

  e->line.p2.x += dx;

  e->line.p1.y += dy;

  e->line.p2.y += dy;

    }

}