#include "toonz/rasterbrush.h"

#include <vector></vector>

#include <tcurves.h></tcurves.h>

using namespace std;

namespace {

//=====================================================================================

// forward declaration

void lightPixel(const TRasterCM32P &ras, const TPoint &pix, double distance,

                int styleId, bool checkAntialiasedPixel);

class ConeSubVolume {

  const static double m_values[21];

public:

  ConeSubVolume() {}

  // calcola il sottovolume di un cono di raggio e volume unitario in base

  static double compute(double cover) {

    double x = (10 * tcrop(cover, -1.0, 1.0)) + 10;

    assert(0 <= x && x <= 20);

    int i = tfloor(x);

    if (i == 20)

      return m_values[i];

    else

      // Interpolazione lineare.

      return (-(x - (i + 1)) * m_values[i]) - (-(x - i) * m_values[i + 1]);

  }

};

const double ConeSubVolume::m_values[] = {

    1.0,      0.99778,  0.987779,  0.967282,  0.934874,  0.889929,   0.832457,

    0.763067, 0.683002, 0.594266,  0.5,       0.405734,  0.316998,   0.236933,

    0.167543, 0.110071, 0.0651259, 0.0327182, 0.0122208, 0.00221986, 0.0};

//========================================================================================================

// Permette di disegnare un disco pieno

class Disk {

  TPointD m_centre;

  double m_radius;

  bool m_doAntialias;

  // Traccia una linea tra due punti per riempire il disco

  void fill(const TRasterCM32P &ras, const TPoint &p1, const TPoint &p2,

            int styleId) const {

    if (p1.y == p2.y) {

      int xMax, xMin;

      if (p1.x > p2.x) {

        xMax = p1.x;

        xMin = p2.x;

      } else {

        xMin = p1.x;

        xMax = p2.x;

      }

      TPixelCM32 color(styleId, 0, 0);

      for (int i = xMin; i <= xMax; i++) ras->pixels(p1.y)[i] = color;

    } else {

      int yMax, yMin;

      if (p1.y > p2.y) {

        yMax = p1.y;

        yMin = p2.y;

      } else {

        yMin = p1.y;

        yMax = p2.y;

      }

      TPixelCM32 color(styleId, 0, 0);

      for (int i = yMin; i <= yMax; i++) ras->pixels(i)[p1.x] = color;

    }

  }

  // Calcola la distanza di un pixel dal centro "reale"

  inline double distancePointToCentre(const TPoint &point) const {

    double d = sqrt((point.x - m_centre.x) * (point.x - m_centre.x) +

                    (point.y - m_centre.y) * (point.y - m_centre.y));

    return d;

  }

  // Calcola la distanza tra un punto della cironferenza ideale ed il pixel che

  // la approssima

  // Inoltre calcola la distanza dei pixel vicini per gestire l'antialias

  void computeDistances(double distances[3], const TPoint &point,

                        const TPoint ¢re, bool upperPoint) const {

    TPoint p = point - centre;

    if (upperPoint) {

      double d     = distancePointToCentre(point);

      distances[0] = d - m_radius;

      d            = distancePointToCentre(TPoint(point.x, point.y + 1));

      distances[1] = d - m_radius;

      d            = distancePointToCentre(TPoint(point.x, point.y - 1));

      distances[2] = d - m_radius;

      return;

    } else {

      double d     = distancePointToCentre(point);

      distances[0] = d - m_radius;

      d            = distancePointToCentre(TPoint(point.x + 1, point.y));

      distances[1] = d - m_radius;

      d            = distancePointToCentre(TPoint(point.x - 1, point.y));

      distances[2] = d - m_radius;

      return;

    }

  }

  // Illumina otto pixel symmetrici per ottenere la circonferenza del disco.

  // Per ogni pixel illumina anche i vicini per gestire l'antialias.

  void makeAntiAliasedDiskBorder(const TRasterCM32P &ras, const TPoint ¢re,

                                 const TPoint &point,

                                 double distancesAntialias[3], int styleId,

                                 int maxPointToFill) const {

    computeDistances(distancesAntialias, point + centre, centre, true);

    lightPixel(ras, point + centre, distancesAntialias[0], styleId, false);

    lightPixel(ras, TPoint(point.x, point.y + 1) + centre,

               distancesAntialias[1], styleId, false);

    lightPixel(ras, TPoint(point.x, point.y - 1) + centre,

               distancesAntialias[2], styleId, false);

    computeDistances(distancesAntialias, TPoint(point.y, point.x) + centre,

                     centre, false);

    lightPixel(ras, TPoint(point.y, point.x) + centre, distancesAntialias[0],

               styleId, false);

    lightPixel(ras, TPoint(point.y + 1, point.x) + centre,

               distancesAntialias[1], styleId, false);

    lightPixel(ras, TPoint(point.y - 1, point.x) + centre,

               distancesAntialias[2], styleId, false);

    computeDistances(distancesAntialias, TPoint(-point.x, -point.y) + centre,

                     centre, true);

    lightPixel(ras, TPoint(-point.x, -point.y) + centre, distancesAntialias[0],

               styleId, false);

    lightPixel(ras, TPoint(-point.x, -point.y - 1) + centre,

               distancesAntialias[2], styleId, false);

    lightPixel(ras, TPoint(-point.x, -point.y + 1) + centre,

               distancesAntialias[1], styleId, false);

    computeDistances(distancesAntialias,

                     TPoint(-point.y + centre.x, point.x + centre.y), centre,

                     false);

    lightPixel(ras, TPoint(-point.y, point.x) + centre, distancesAntialias[0],

               styleId, false);

    lightPixel(ras, TPoint(-point.y - 1, point.x) + centre,

               distancesAntialias[2], styleId, false);

    lightPixel(ras, TPoint(-point.y + 1, point.x) + centre,

               distancesAntialias[1], styleId, false);

    if ((point.x + centre.x) != (centre.x)) {

      computeDistances(distancesAntialias, TPoint(point.y, -point.x) + centre,

                       centre, false);

      lightPixel(ras, TPoint(point.y, -point.x) + centre, distancesAntialias[0],

                 styleId, false);

      lightPixel(ras, TPoint(point.y + 1, -point.x) + centre,

                 distancesAntialias[1], styleId, false);

      lightPixel(ras, TPoint(point.y - 1, -point.x) + centre,

                 distancesAntialias[2], styleId, false);

      computeDistances(distancesAntialias, TPoint(point.x, -point.y) + centre,

                       centre, true);

      lightPixel(ras, TPoint(point.x, -point.y) + centre, distancesAntialias[0],

                 styleId, false);

      lightPixel(ras, TPoint(point.x, -point.y - 1) + centre,

                 distancesAntialias[2], styleId, false);

      lightPixel(ras, TPoint(point.x, -point.y + 1) + centre,

                 distancesAntialias[1], styleId, false);

      computeDistances(distancesAntialias, TPoint(-point.y, -point.x) + centre,

                       centre, false);

      lightPixel(ras, TPoint(-point.y, -point.x) + centre,

                 distancesAntialias[0], styleId, false);

      lightPixel(ras, TPoint(-point.y - 1, -point.x) + centre,

                 distancesAntialias[2], styleId, false);

      lightPixel(ras, TPoint(-point.y + 1, -point.x) + centre,

                 distancesAntialias[1], styleId, false);

      computeDistances(distancesAntialias, TPoint(-point.x, point.y) + centre,

                       centre, true);

      lightPixel(ras, TPoint(-point.x, point.y) + centre, distancesAntialias[0],

                 styleId, false);

      lightPixel(ras, TPoint(-point.x, point.y + 1) + centre,

                 distancesAntialias[1], styleId, false);

      lightPixel(ras, TPoint(-point.x, point.y - 1) + centre,

                 distancesAntialias[2], styleId, false);

    }

    if (maxPointToFill <= point.y - 2) {

      fill(ras, TPoint(point.x, maxPointToFill) + centre,

           TPoint(point.x, point.y - 2) + centre, styleId);

      fill(ras, TPoint(maxPointToFill, point.x) + centre,

           TPoint(point.y - 2, point.x) + centre, styleId);

      fill(ras, TPoint(point.x, -maxPointToFill) + centre,

           TPoint(point.x, -point.y + 2) + centre, styleId);

      fill(ras, TPoint(-maxPointToFill, -point.x) + centre,

           TPoint(-point.y + 2, -point.x) + centre, styleId);

      if (point.x != 0) {

        fill(ras, TPoint(maxPointToFill, -point.x) + centre,

             TPoint(point.y - 2, -point.x) + centre, styleId);

        fill(ras, TPoint(-point.x, -maxPointToFill) + centre,

             TPoint(-point.x, -point.y + 2) + centre, styleId);

        fill(ras, TPoint(-maxPointToFill, point.x) + centre,

             TPoint(-point.y + 2, point.x) + centre, styleId);

        fill(ras, TPoint(-point.x, maxPointToFill) + centre,

             TPoint(-point.x, point.y - 2) + centre, styleId);

      }

    }

  }

  void makeNoAntiAliasedDiskBorder(const TRasterCM32P &ras,

                                   const TPoint ¢re, const TPoint &point,

                                   int styleId, int maxPointToFill) const {

    if (((int)(m_radius * 2)) % 2 == 0) {

      lightPixel(ras, TPoint(point.x - 1, point.y) + centre, -1, styleId,

                 false);

      lightPixel(ras, TPoint(point.y - 1, -point.x + 1) + centre, -1, styleId,

                 false);

      lightPixel(ras, TPoint(-point.x, -point.y + 1) + centre, -1, styleId,

                 false);

      lightPixel(ras, TPoint(-point.y, point.x) + centre, -1, styleId, false);

    } else {

      lightPixel(ras, point + centre, -1, styleId, false);

      lightPixel(ras, TPoint(point.y, -point.x) + centre, -1, styleId, false);

      lightPixel(ras, TPoint(-point.x, -point.y) + centre, -1, styleId, false);

      lightPixel(ras, TPoint(-point.y, point.x) + centre, -1, styleId, false);

    }

    if ((point.x + centre.x) != (centre.x)) {

      if (((int)(m_radius * 2)) % 2 == 0) {

        lightPixel(ras, TPoint(point.y - 1, point.x) + centre, -1, styleId,

                   false);

        lightPixel(ras, TPoint(point.x - 1, -point.y + 1) + centre, -1, styleId,

                   false);

        lightPixel(ras, TPoint(-point.y, -point.x + 1) + centre, -1, styleId,

                   false);

        lightPixel(ras, TPoint(-point.x, point.y) + centre, -1, styleId, false);

      } else {

        lightPixel(ras, TPoint(point.y, point.x) + centre, -1, styleId, false);

        lightPixel(ras, TPoint(point.x, -point.y) + centre, -1, styleId, false);

        lightPixel(ras, TPoint(-point.y, -point.x) + centre, -1, styleId,

                   false);

        lightPixel(ras, TPoint(-point.x, point.y) + centre, -1, styleId, false);

      }

    }

    if (maxPointToFill <= point.y - 1) {

      if (((int)(m_radius * 2)) % 2 == 0) {

        fill(ras, TPoint(point.x - 1, 0) + centre,

             TPoint(point.x - 1, point.y) + centre, styleId);

        fill(ras, TPoint(0, -point.x + 1) + centre,

             TPoint(point.y - 1, -point.x + 1) + centre, styleId);

        fill(ras, TPoint(-point.x, 0) + centre,

             TPoint(-point.x, -point.y + 1) + centre, styleId);

        fill(ras, TPoint(0, point.x) + centre,

             TPoint(-point.y, point.x) + centre, styleId);

        if (point.x != 0) {

          fill(ras, TPoint(0, point.x) + centre,

               TPoint(point.y - 1, point.x) + centre, styleId);

          fill(ras, TPoint(point.x - 1, 0) + centre,

               TPoint(point.x - 1, -point.y + 1) + centre, styleId);

          fill(ras, TPoint(0, -point.x + 1) + centre,

               TPoint(-point.y, -point.x + 1) + centre, styleId);

          fill(ras, TPoint(-point.x, 0) + centre,

               TPoint(-point.x, point.y) + centre, styleId);

        }

      } else {

        fill(ras, TPoint(point.x, 0) + centre,

             TPoint(point.x, point.y) + centre, styleId);

        fill(ras, TPoint(0, -point.x) + centre,

             TPoint(point.y, -point.x) + centre, styleId);

        fill(ras, TPoint(-point.x, 0) + centre,

             TPoint(-point.x, -point.y) + centre, styleId);

        fill(ras, TPoint(0, point.x) + centre,

             TPoint(-point.y, point.x) + centre, styleId);

        if (point.x != 0) {

          fill(ras, TPoint(0, point.x) + centre,

               TPoint(point.y, point.x) + centre, styleId);

          fill(ras, TPoint(point.x, -0) + centre,

               TPoint(point.x, -point.y) + centre, styleId);

          fill(ras, TPoint(0, -point.x) + centre,

               TPoint(-point.y, -point.x) + centre, styleId);

          fill(ras, TPoint(-point.x, 0) + centre,

               TPoint(-point.x, point.y) + centre, styleId);

        }

      }

    }

  }

public:

  Disk() : m_centre(0.0, 0.0), m_radius(1.0), m_doAntialias(true) {}

  Disk(const TThickPoint &p, bool doAntialias) : m_doAntialias(doAntialias) {

    if (m_doAntialias) {

      m_centre = TPointD(p.x, p.y);

      m_radius = p.thick * 0.5;

    } else {

      m_centre = TPointD(tround(p.x), tround(p.y));

      m_radius = tround(p.thick) * 0.5;

    }

  }

  // Disegna un disco

  void draw(const TRasterCM32P &ras, int styleId) const {

    double distances[3];

    int maxPointToFill = 0;

    TPoint centre      = convert(m_centre);

    TPoint point;

    if (m_doAntialias) {

      TPoint point(0, tround(m_radius));

      int d = 1 - tround(m_radius), dE = 3, dSE = -2 * tround(m_radius) + 5;

      makeAntiAliasedDiskBorder(ras, centre, point, distances, styleId,

                                maxPointToFill);

      while (point.y > point.x) {

        if (d < 0) {

          d += dE;

          dE += 2;

          dSE += 2;

          point.x++;

          maxPointToFill++;

        } else {

          d += dSE;

          dE += 2;

          dSE += 4;

          point.x++;

          point.y--;

          maxPointToFill++;

        }

        makeAntiAliasedDiskBorder(ras, centre, point, distances, styleId,

                                  maxPointToFill);

      }

    } else {

      TPoint point(0, tround(m_radius - 0.5));

      int d = 3 - 2 * (int)m_radius;

      while (point.y > point.x) {

        makeNoAntiAliasedDiskBorder(ras, centre, point, styleId,

                                    maxPointToFill);

        if (d < 0)

          d = d + 4 * point.x + 6;

        else {

          d = d + 4 * (point.x - point.y) + 10;

          point.y--;

        }

        maxPointToFill++;

        point.x++;

      }

      if (point.x == point.y)

        makeNoAntiAliasedDiskBorder(ras, centre, point, styleId,

                                    maxPointToFill);

    }

  }

  TPointD getCentre() const { return m_centre; }

  double getRadius() const { return m_radius; }

  void setCentre(const TPointD ¢re) { m_centre = centre; }

  void setCentre(double x, double y) {

    m_centre.x = x;

    m_centre.y = y;

  }

  void setRadius(double radius) { m_radius = radius; }

};

//===============================================================================================================

// Calcola la distanza di una curva al parametro "t" dalla retta che unisce i

// punti nei parametri "tPrevious" e "tNext"

double findChordalDeviation(const TQuadratic &quadratic, double t,

                            double tPrevious, double tNext) {

  TPointD pPrevious = quadratic.getPoint(tPrevious);

  TPointD pNext     = quadratic.getPoint(tNext);

  TPointD p         = quadratic.getPoint(t);

  TPointD P(p - pPrevious), Q(pNext - pPrevious);

  double PQ = P * Q;

  double QQ = Q * Q;

  return norm(P - (PQ / QQ) * Q);

}

// Accende un pixel calcolandone l'intensita'

/*-- 筆先のアンチエイリアス部分の描画 --*/

void lightPixel(const TRasterCM32P &ras, const TPoint &pix, double distance,

                int styleId, bool checkAntialiasedPixel) {

  TPixelCM32 pixel      = ras->pixels(pix.y)[pix.x];

  double volumeParziale = ConeSubVolume::compute(distance);

  /*- 現在のToneに乗算していく -*/

  int newTone = tround((double)pixel.getTone() * (1.0 - volumeParziale));

  assert(newTone >= 0 && newTone <= 255);

  ras->pixels(pix.y)[pix.x] = TPixelCM32(styleId, pixel.getPaint(), newTone);

}

bool isDiskNecessaryNecessary(const TQuadratic &quadratic, double tCurrent,

                              double lastT, double nextT,

                              bool &idLastDiskDrown) {

  TPoint currentPoint = convert(quadratic.getPoint(tCurrent));

  TPoint lastPoint    = convert(quadratic.getPoint(lastT));

  TPoint nextPoint    = convert(quadratic.getPoint(nextT));

  TPoint lastDiff     = currentPoint - lastPoint;

  TPoint nextDiff     = currentPoint - nextPoint;

  TPoint otherDiff    = lastPoint - nextPoint;

  bool isLastNear     = lastDiff.x == 0 || lastDiff.y == 0;

  bool isNextNear     = nextDiff.x == 0 || nextDiff.y == 0;

  bool isOtheNear     = otherDiff.x == 0 || otherDiff.y == 0;

  if (isLastNear && isNextNear && !isOtheNear && idLastDiskDrown) {

    idLastDiskDrown = false;

    return false;

  }

  return true;

}

// Disegna una porzione di del tratto (un piccolo arco)

void makeLittleArch(const TRasterCM32P &ras, const Disk &disk1,

                    const Disk &disk2, const Disk &disk3, int styleId,

                    bool doAntialias) {

  TPointD center1 = disk1.getCentre();

  TPointD center2 = disk2.getCentre();

  TPointD center3 = disk3.getCentre();

  TQuadratic quadratic(center1, center2, center3);

  disk1.draw(ras, styleId);

  double length = quadratic.getLength();

  if (length < 2) return;

  double t = 0, step = 1 / (length * 1.5), t2 = quadratic.getT(center2);

  bool idLastDiskDrown = true;

  for (t = step; t < 1; t += step) {

    TPointD center = quadratic.getPoint(t);

    double radius =

        disk1.getRadius() + (disk3.getRadius() - disk1.getRadius()) * t;

    Disk disk(TThickPoint(center, radius * 2), doAntialias);

    bool drawDisk = true;

    if (!doAntialias)

      drawDisk = isDiskNecessaryNecessary(

          quadratic, t, t - step, t + step > 1 ? 1 : t + step, idLastDiskDrown);

    if (t != 1 && drawDisk) {

      disk.draw(ras, styleId);

      idLastDiskDrown = true;

    }

  }

  disk3.draw(ras, styleId);

}

// Disegna un piccolo segmento, invece di un archetto.

void makeLittleSegment(const TRasterCM32P &ras, const Disk &disk1,

                       const Disk &disk2, int styleId, bool doAntialias) {

  TPointD center1 = disk1.getCentre();

  TPointD center2 = disk2.getCentre();

  TPointD middle  = (center1 + center2) * 0.5;

  double raius    = (disk1.getRadius() + disk2.getRadius()) * 0.5;

  Disk disk(TThickPoint(middle, raius * 2), doAntialias);

  makeLittleArch(ras, disk1, disk, disk2, styleId, doAntialias);

}

//=============================================================================

}  // namespace

//=============================================================================

// Preso un vettore di punti, disegna una pennelata che li approssima in un

// raster trasparente

void rasterBrush(const TRasterCM32P &rasBuffer,

                 const vector<tthickpoint> &points, int styleId,</tthickpoint>

                 bool doAntialias) {

  int i, n = points.size();

  if (n == 0)

    return;

  else if (n == 1) {

    Disk disk(points[0], doAntialias);

    disk.draw(rasBuffer, styleId);

    return;

  } else if (n == 2) {

    makeLittleSegment(rasBuffer, Disk(points[0], doAntialias),

                      Disk(points[1], doAntialias), styleId, doAntialias);

    return;

  } else if (n == 4) {

    makeLittleArch(rasBuffer, Disk(points[0], doAntialias),

                   Disk(points[1], doAntialias), Disk(points[2], doAntialias),

                   styleId, doAntialias);

    makeLittleSegment(rasBuffer, Disk(points[2], doAntialias),

                      Disk(points[3], doAntialias), styleId, doAntialias);

    return;

  } else {

    for (i = 0; i + 2 < n; i += 2)

      makeLittleArch(rasBuffer, Disk(points[i], doAntialias),

                     Disk(points[i + 1], doAntialias),

                     Disk(points[i + 2], doAntialias), styleId, doAntialias);

  }

}