/*

 * Copyright © 2004 Carl Worth

 * Copyright © 2006 Red Hat, Inc.

 * Copyright © 2008 Chris Wilson

 *

 * 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 Carl Worth

 *

 * Contributor(s):

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

 *  Chris Wilson <chris@chris-wilson.co.uk>

 */

/* Provide definitions for standalone compilation */

#include "cairoint.h"

#include "cairo-boxes-private.h"

#include "cairo-combsort-inline.h"

#include "cairo-error-private.h"

#include "cairo-traps-private.h"

typedef struct _cairo_bo_edge cairo_bo_edge_t;

typedef struct _cairo_bo_trap cairo_bo_trap_t;

/* A deferred trapezoid of an edge */

struct _cairo_bo_trap {

    cairo_bo_edge_t *right;

    int32_t top;

};

struct _cairo_bo_edge {

    cairo_edge_t edge;

    cairo_bo_edge_t *prev;

    cairo_bo_edge_t *next;

    cairo_bo_trap_t deferred_trap;

};

typedef enum {

    CAIRO_BO_EVENT_TYPE_START,

    CAIRO_BO_EVENT_TYPE_STOP

} cairo_bo_event_type_t;

typedef struct _cairo_bo_event {

    cairo_bo_event_type_t type;

    cairo_point_t point;

    cairo_bo_edge_t *edge;

} cairo_bo_event_t;

typedef struct _cairo_bo_sweep_line {

    cairo_bo_event_t **events;

    cairo_bo_edge_t *head;

    cairo_bo_edge_t *stopped;

    int32_t current_y;

    cairo_bo_edge_t *current_edge;

} cairo_bo_sweep_line_t;

static inline int

_cairo_point_compare (const cairo_point_t *a,

          const cairo_point_t *b)

{

    int cmp;

    cmp = a->y - b->y;

    if (likely (cmp))

  return cmp;

    return a->x - b->x;

}

static inline int

_cairo_bo_edge_compare (const cairo_bo_edge_t *a,

      const cairo_bo_edge_t *b)

{

    int cmp;

    cmp = a->edge.line.p1.x - b->edge.line.p1.x;

    if (likely (cmp))

  return cmp;

    return b->edge.bottom - a->edge.bottom;

}

static inline int

cairo_bo_event_compare (const cairo_bo_event_t *a,

      const cairo_bo_event_t *b)

{

    int cmp;

    cmp = _cairo_point_compare (&a->point, &b->point);

    if (likely (cmp))

  return cmp;

    cmp = a->type - b->type;

    if (cmp)

  return cmp;

    return a - b;

}

static inline cairo_bo_event_t *

_cairo_bo_event_dequeue (cairo_bo_sweep_line_t *sweep_line)

{

    return *sweep_line->events++;

}

CAIRO_COMBSORT_DECLARE (_cairo_bo_event_queue_sort,

      cairo_bo_event_t *,

      cairo_bo_event_compare)

static void

_cairo_bo_sweep_line_init (cairo_bo_sweep_line_t *sweep_line,

         cairo_bo_event_t **events,

         int        num_events)

{

    _cairo_bo_event_queue_sort (events, num_events);

    events[num_events] = NULL;

    sweep_line->events = events;

    sweep_line->head = NULL;

    sweep_line->current_y = INT32_MIN;

    sweep_line->current_edge = NULL;

}

static void

_cairo_bo_sweep_line_insert (cairo_bo_sweep_line_t  *sweep_line,

           cairo_bo_edge_t    *edge)

{

    if (sweep_line->current_edge != NULL) {

  cairo_bo_edge_t *prev, *next;

  int cmp;

  cmp = _cairo_bo_edge_compare (sweep_line->current_edge, edge);

  if (cmp < 0) {

      prev = sweep_line->current_edge;

      next = prev->next;

      while (next != NULL && _cairo_bo_edge_compare (next, edge) < 0)

    prev = next, next = prev->next;

      prev->next = edge;

      edge->prev = prev;

      edge->next = next;

      if (next != NULL)

    next->prev = edge;

  } else if (cmp > 0) {

      next = sweep_line->current_edge;

      prev = next->prev;

      while (prev != NULL && _cairo_bo_edge_compare (prev, edge) > 0)

    next = prev, prev = next->prev;

      next->prev = edge;

      edge->next = next;

      edge->prev = prev;

      if (prev != NULL)

    prev->next = edge;

      else

    sweep_line->head = edge;

  } else {

      prev = sweep_line->current_edge;

      edge->prev = prev;

      edge->next = prev->next;

      if (prev->next != NULL)

    prev->next->prev = edge;

      prev->next = edge;

  }

    } else {

  sweep_line->head = edge;

    }

    sweep_line->current_edge = edge;

}

static void

_cairo_bo_sweep_line_delete (cairo_bo_sweep_line_t  *sweep_line,

           cairo_bo_edge_t  *edge)

{

    if (edge->prev != NULL)

  edge->prev->next = edge->next;

    else

  sweep_line->head = edge->next;

    if (edge->next != NULL)

  edge->next->prev = edge->prev;

    if (sweep_line->current_edge == edge)

  sweep_line->current_edge = edge->prev ? edge->prev : edge->next;

}

static inline cairo_bool_t

edges_collinear (const cairo_bo_edge_t *a, const cairo_bo_edge_t *b)

{

    return a->edge.line.p1.x == b->edge.line.p1.x;

}

static cairo_status_t

_cairo_bo_edge_end_trap (cairo_bo_edge_t  *left,

       int32_t     bot,

       cairo_bool_t    do_traps,

       void     *container)

{

    cairo_bo_trap_t *trap = &left->deferred_trap;

    cairo_status_t status = CAIRO_STATUS_SUCCESS;

    /* Only emit (trivial) non-degenerate trapezoids with positive height. */

    if (likely (trap->top < bot)) {

  if (do_traps) {

      _cairo_traps_add_trap (container,

           trap->top, bot,

           &left->edge.line, &trap->right->edge.line);

      status =  _cairo_traps_status ((cairo_traps_t *) container);

  } else {

      cairo_box_t box;

      box.p1.x = left->edge.line.p1.x;

      box.p1.y = trap->top;

      box.p2.x = trap->right->edge.line.p1.x;

      box.p2.y = bot;

      status = _cairo_boxes_add (container, CAIRO_ANTIALIAS_DEFAULT, &box);

  }

    }

    trap->right = NULL;

    return status;

}

/* Start a new trapezoid at the given top y coordinate, whose edges

 * are `edge' and `edge->next'. If `edge' already has a trapezoid,

 * then either add it to the traps in `traps', if the trapezoid's

 * right edge differs from `edge->next', or do nothing if the new

 * trapezoid would be a continuation of the existing one. */

static inline cairo_status_t

_cairo_bo_edge_start_or_continue_trap (cairo_bo_edge_t  *left,

               cairo_bo_edge_t  *right,

               int               top,

               cairo_bool_t  do_traps,

               void   *container)

{

    cairo_status_t status;

    if (left->deferred_trap.right == right)

  return CAIRO_STATUS_SUCCESS;

    if (left->deferred_trap.right != NULL) {

  if (right != NULL && edges_collinear (left->deferred_trap.right, right))

  {

      /* continuation on right, so just swap edges */

      left->deferred_trap.right = right;

      return CAIRO_STATUS_SUCCESS;

  }

  status = _cairo_bo_edge_end_trap (left, top, do_traps, container);

  if (unlikely (status))

      return status;

    }

    if (right != NULL && ! edges_collinear (left, right)) {

  left->deferred_trap.top = top;

  left->deferred_trap.right = right;

    }

    return CAIRO_STATUS_SUCCESS;

}

static inline cairo_status_t

_active_edges_to_traps (cairo_bo_edge_t   *left,

      int32_t      top,

      cairo_fill_rule_t  fill_rule,

      cairo_bool_t     do_traps,

      void      *container)

{

    cairo_bo_edge_t *right;

    cairo_status_t status;

    if (fill_rule == CAIRO_FILL_RULE_WINDING) {

  while (left != NULL) {

      int in_out;

      /* Greedily search for the closing edge, so that we generate the

       * maximal span width with the minimal number of trapezoids.

       */

      in_out = left->edge.dir;

      /* Check if there is a co-linear edge with an existing trap */

      right = left->next;

      if (left->deferred_trap.right == NULL) {

    while (right != NULL && right->deferred_trap.right == NULL)

        right = right->next;

    if (right != NULL && edges_collinear (left, right)) {

        /* continuation on left */

        left->deferred_trap = right->deferred_trap;

        right->deferred_trap.right = NULL;

    }

      }

      /* End all subsumed traps */

      right = left->next;

      while (right != NULL) {

    if (right->deferred_trap.right != NULL) {

        status = _cairo_bo_edge_end_trap (right, top, do_traps, container);

        if (unlikely (status))

      return status;

    }

    in_out += right->edge.dir;

    if (in_out == 0) {

        /* skip co-linear edges */

        if (right->next == NULL ||

      ! edges_collinear (right, right->next))

        {

      break;

        }

    }

    right = right->next;

      }

      status = _cairo_bo_edge_start_or_continue_trap (left, right, top,

                  do_traps, container);

      if (unlikely (status))

    return status;

      left = right;

      if (left != NULL)

    left = left->next;

  }

    } else {

  while (left != NULL) {

      int in_out = 0;

      right = left->next;

      while (right != NULL) {

    if (right->deferred_trap.right != NULL) {

        status = _cairo_bo_edge_end_trap (right, top, do_traps, container);

        if (unlikely (status))

      return status;

    }

    if ((in_out++ & 1) == 0) {

        cairo_bo_edge_t *next;

        cairo_bool_t skip = FALSE;

        /* skip co-linear edges */

        next = right->next;

        if (next != NULL)

      skip = edges_collinear (right, next);

        if (! skip)

      break;

    }

    right = right->next;

      }

      status = _cairo_bo_edge_start_or_continue_trap (left, right, top,

                  do_traps, container);

      if (unlikely (status))

    return status;

      left = right;

      if (left != NULL)

    left = left->next;

  }

    }

    return CAIRO_STATUS_SUCCESS;

}

static cairo_status_t

_cairo_bentley_ottmann_tessellate_rectilinear (cairo_bo_event_t   **start_events,

                 int       num_events,

                 cairo_fill_rule_t   fill_rule,

                 cairo_bool_t    do_traps,

                 void     *container)

{

    cairo_bo_sweep_line_t sweep_line;

    cairo_bo_event_t *event;

    cairo_status_t status;

    _cairo_bo_sweep_line_init (&sweep_line, start_events, num_events);

    while ((event = _cairo_bo_event_dequeue (&sweep_line))) {

  if (event->point.y != sweep_line.current_y) {

      status = _active_edges_to_traps (sweep_line.head,

               sweep_line.current_y,

               fill_rule, do_traps, container);

      if (unlikely (status))

    return status;

      sweep_line.current_y = event->point.y;

  }

  switch (event->type) {

  case CAIRO_BO_EVENT_TYPE_START:

      _cairo_bo_sweep_line_insert (&sweep_line, event->edge);

      break;

  case CAIRO_BO_EVENT_TYPE_STOP:

      _cairo_bo_sweep_line_delete (&sweep_line, event->edge);

      if (event->edge->deferred_trap.right != NULL) {

    status = _cairo_bo_edge_end_trap (event->edge,

              sweep_line.current_y,

              do_traps, container);

    if (unlikely (status))

        return status;

      }

      break;

  }

    }

    return CAIRO_STATUS_SUCCESS;

}

cairo_status_t

_cairo_bentley_ottmann_tessellate_rectilinear_polygon_to_boxes (const cairo_polygon_t *polygon,

                cairo_fill_rule_t   fill_rule,

                cairo_boxes_t *boxes)

{

    cairo_status_t status;

    cairo_bo_event_t stack_events[CAIRO_STACK_ARRAY_LENGTH (cairo_bo_event_t)];

    cairo_bo_event_t *events;

    cairo_bo_event_t *stack_event_ptrs[ARRAY_LENGTH (stack_events) + 1];

    cairo_bo_event_t **event_ptrs;

    cairo_bo_edge_t stack_edges[ARRAY_LENGTH (stack_events)];

    cairo_bo_edge_t *edges;

    int num_events;

    int i, j;

    if (unlikely (polygon->num_edges == 0))

  return CAIRO_STATUS_SUCCESS;

    num_events = 2 * polygon->num_edges;

    events = stack_events;

    event_ptrs = stack_event_ptrs;

    edges = stack_edges;

    if (num_events > ARRAY_LENGTH (stack_events)) {

  events = _cairo_malloc_ab_plus_c (num_events,

            sizeof (cairo_bo_event_t) +

            sizeof (cairo_bo_edge_t) +

            sizeof (cairo_bo_event_t *),

            sizeof (cairo_bo_event_t *));

  if (unlikely (events == NULL))

      return _cairo_error (CAIRO_STATUS_NO_MEMORY);

  event_ptrs = (cairo_bo_event_t **) (events + num_events);

  edges = (cairo_bo_edge_t *) (event_ptrs + num_events + 1);

    }

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

  edges[i].edge = polygon->edges[i];

  edges[i].deferred_trap.right = NULL;

  edges[i].prev = NULL;

  edges[i].next = NULL;

  event_ptrs[j] = &events[j];

  events[j].type = CAIRO_BO_EVENT_TYPE_START;

  events[j].point.y = polygon->edges[i].top;

  events[j].point.x = polygon->edges[i].line.p1.x;

  events[j].edge = &edges[i];

  j++;

  event_ptrs[j] = &events[j];

  events[j].type = CAIRO_BO_EVENT_TYPE_STOP;

  events[j].point.y = polygon->edges[i].bottom;

  events[j].point.x = polygon->edges[i].line.p1.x;

  events[j].edge = &edges[i];

  j++;

    }

    status = _cairo_bentley_ottmann_tessellate_rectilinear (event_ptrs, j,

                  fill_rule,

                  FALSE, boxes);

    if (events != stack_events)

  free (events);

    return status;

}

cairo_status_t

_cairo_bentley_ottmann_tessellate_rectilinear_traps (cairo_traps_t *traps,

                 cairo_fill_rule_t fill_rule)

{

    cairo_bo_event_t stack_events[CAIRO_STACK_ARRAY_LENGTH (cairo_bo_event_t)];

    cairo_bo_event_t *events;

    cairo_bo_event_t *stack_event_ptrs[ARRAY_LENGTH (stack_events) + 1];

    cairo_bo_event_t **event_ptrs;

    cairo_bo_edge_t stack_edges[ARRAY_LENGTH (stack_events)];

    cairo_bo_edge_t *edges;

    cairo_status_t status;

    int i, j, k;

    if (unlikely (traps->num_traps == 0))

  return CAIRO_STATUS_SUCCESS;

    assert (traps->is_rectilinear);

    i = 4 * traps->num_traps;

    events = stack_events;

    event_ptrs = stack_event_ptrs;

    edges = stack_edges;

    if (i > ARRAY_LENGTH (stack_events)) {

  events = _cairo_malloc_ab_plus_c (i,

            sizeof (cairo_bo_event_t) +

            sizeof (cairo_bo_edge_t) +

            sizeof (cairo_bo_event_t *),

            sizeof (cairo_bo_event_t *));

  if (unlikely (events == NULL))

      return _cairo_error (CAIRO_STATUS_NO_MEMORY);

  event_ptrs = (cairo_bo_event_t **) (events + i);

  edges = (cairo_bo_edge_t *) (event_ptrs + i + 1);

    }

    for (i = j = k = 0; i < traps->num_traps; i++) {

  edges[k].edge.top = traps->traps[i].top;

  edges[k].edge.bottom = traps->traps[i].bottom;

  edges[k].edge.line = traps->traps[i].left;

  edges[k].edge.dir = 1;

  edges[k].deferred_trap.right = NULL;

  edges[k].prev = NULL;

  edges[k].next = NULL;

  event_ptrs[j] = &events[j];

  events[j].type = CAIRO_BO_EVENT_TYPE_START;

  events[j].point.y = traps->traps[i].top;

  events[j].point.x = traps->traps[i].left.p1.x;

  events[j].edge = &edges[k];

  j++;

  event_ptrs[j] = &events[j];

  events[j].type = CAIRO_BO_EVENT_TYPE_STOP;

  events[j].point.y = traps->traps[i].bottom;

  events[j].point.x = traps->traps[i].left.p1.x;

  events[j].edge = &edges[k];

  j++;

  k++;

  edges[k].edge.top = traps->traps[i].top;

  edges[k].edge.bottom = traps->traps[i].bottom;

  edges[k].edge.line = traps->traps[i].right;

  edges[k].edge.dir = -1;

  edges[k].deferred_trap.right = NULL;

  edges[k].prev = NULL;

  edges[k].next = NULL;

  event_ptrs[j] = &events[j];

  events[j].type = CAIRO_BO_EVENT_TYPE_START;

  events[j].point.y = traps->traps[i].top;

  events[j].point.x = traps->traps[i].right.p1.x;

  events[j].edge = &edges[k];

  j++;

  event_ptrs[j] = &events[j];

  events[j].type = CAIRO_BO_EVENT_TYPE_STOP;

  events[j].point.y = traps->traps[i].bottom;

  events[j].point.x = traps->traps[i].right.p1.x;

  events[j].edge = &edges[k];

  j++;

  k++;

    }

    _cairo_traps_clear (traps);

    status = _cairo_bentley_ottmann_tessellate_rectilinear (event_ptrs, j,

                  fill_rule,

                  TRUE, traps);

    traps->is_rectilinear = TRUE;

    if (events != stack_events)

  free (events);

    return status;

}