#include <tools assistant.h=""></tools>

#include <tgl.h></tgl.h>

#include <tproperty.h></tproperty.h>

#include <limits></limits>

#include <cassert></cassert>

#ifdef MACOSX

const double line_width_scale = 1.5;

#else

const double line_width_scale = 1.0;

#endif

//************************************************************************

//    TGuideline implementation

//************************************************************************

void

TGuideline::drawSegment(

  const TPointD &p0,

  const TPointD &p1,

  double pixelSize,

  bool active,

  bool enabled ) const

{

  double colorBlack[4] = { 0.0, 0.0, 0.0, 0.25 };

  double colorWhite[4] = { 1.0, 1.0, 1.0, 0.25 };

  if (!this->enabled || !enabled)

    colorBlack[3] = (colorWhite[3] = 0.125);

  else if (active)

    colorBlack[3] = (colorWhite[3] = 0.5);

  glPushAttrib(GL_ALL_ATTRIB_BITS);

  tglEnableBlending();

  tglEnableLineSmooth(true, 1.0 * line_width_scale);

  TPointD d = p1 - p0;

  double k = norm2(d);

  if (k > TConsts::epsilon*TConsts::epsilon) {

    k = 0.5*pixelSize*line_width_scale/sqrt(k);

    d = TPointD(-k*d.y, k*d.x);

    glColor4dv(colorWhite);

    tglDrawSegment(p0 - d, p1 - d);

    glColor4dv(colorBlack);

    tglDrawSegment(p0 + d, p1 + d);

  }

  glPopAttrib();

}

//---------------------------------------------------------------------------------------------------

double

TGuideline::calcTrackWeight(const TTrack &track, const TAffine &toScreen, bool &outLongEnough) const {

  outLongEnough = false;

  if (!enabled || track.size() < 2)

    return std::numeric_limits<double>::infinity();</double>

  const double snapLenght = 20.0;

  const double snapScale = 1.0;

  const double maxLength = 2.0*snapLenght*snapScale;

  double sumWeight = 0.0;

  double sumLength = 0.0;

  double sumDeviation = 0.0;

  TPointD prev = toScreen*track[0].position;

  for(int i = 1; i < track.size(); ++i) {

    const TTrackPoint &tp = track[i];

    TTrackPoint mid = TTrack::interpolationLinear(track[i-1], track[i], 0.5);

    TPointD p = toScreen*tp.position;

    double length = tdistance(p, prev);

    sumLength += length;

    double midStepLength = sumLength - 0.5*length;

    if (midStepLength > TConsts::epsilon) {

      double weight = length*logNormalDistribuitionUnscaled(midStepLength, snapLenght, snapScale);

      sumWeight += weight;

      TPointD ntp = nearestPoint(mid.position);

      double deviation = tdistance(toScreen*ntp, p);

      sumDeviation += weight*deviation;

    }

    prev = p;

    if (sumLength >= maxLength)

      { outLongEnough = i < track.fixedSize(); break; }

  }

  return sumWeight > TConsts::epsilon

       ? sumDeviation/sumWeight

       : std::numeric_limits<double>::infinity();</double>

}

//---------------------------------------------------------------------------------------------------

TGuidelineP

TGuideline::findBest(const TGuidelineList &guidelines, const TTrack &track, const TAffine &toScreen, bool &outLongEnough) {

  outLongEnough = true;

  double bestWeight = 0.0;

  TGuidelineP best;

  for(TGuidelineList::const_iterator i = guidelines.begin(); i != guidelines.end(); ++i) {

    double weight = (*i)->calcTrackWeight(track, toScreen, outLongEnough);

    if (!best || weight < bestWeight)

      { bestWeight = weight; best = *i; }

  }

  return best;

}

//************************************************************************

//    TAssistantPoint implementation

//************************************************************************

TAssistantPoint::TAssistantPoint(const TStringId &name, const TPointD &defPosition):

  name(name),

  defPosition(defPosition),

  type(Circle),

  position(defPosition),

  radius(10.0),

  visible(true),

  selected() { }

//************************************************************************

//    TAssistantType implementation

//************************************************************************

TMetaObjectHandler*

TAssistantType::createHandler(TMetaObject &obj) const

  { return createAssistant(obj); }

//************************************************************************

//    TAssistant implementation

//************************************************************************

TAssistant::TAssistant(TMetaObject &object):

  TMetaObjectHandler(object),

  m_idEnabled("enabled"),

  m_idPoints("points"),

  m_idX("x"),

  m_idY("y"),

  m_idMagnetism("magnetism"),

  m_basePoint()

{

  addProperty( new TBoolProperty(m_idEnabled.str(), getEnabled()) );

  addProperty( new TDoubleProperty(m_idMagnetism.str(), 0.0, 1.0, getMagnetism()) );

}

//---------------------------------------------------------------------------------------------------

TAssistantPoint&

TAssistant::addPoint(

  const TStringId &name,

  TAssistantPoint::Type type,

  const TPointD &defPosition,

  bool visible,

  double radius )

{

  assert(!m_points.count(name));

  TAssistantPoint &p = m_points.insert(

    TAssistantPointMap::value_type(name, TAssistantPoint(name, defPosition)) ).first->second;

  m_pointsOrder.push_back(&p);

  p.type     = type;

  p.radius   = radius;

  p.visible  = visible;

  if (!m_basePoint) m_basePoint = &p;

  return p;

}

//---------------------------------------------------------------------------------------------------

TAssistantPoint&

TAssistant::addPoint(

  const TStringId &name,

  TAssistantPoint::Type type,

  const TPointD &defPosition,

  bool visible )

    { return addPoint(name, type, defPosition, visible, 10.0); }

//---------------------------------------------------------------------------------------------------

const TAssistantPoint&

TAssistant::getBasePoint() const

  { assert(m_basePoint); return *m_basePoint; }

//---------------------------------------------------------------------------------------------------

void

TAssistant::addProperty(TProperty *p)

  { m_properties.add(p); }

//---------------------------------------------------------------------------------------------------

void

TAssistant::setTranslation(const TStringId &name, const QString &localName) const

  { m_properties.getProperty(name)->setQStringName( localName ); }

//---------------------------------------------------------------------------------------------------

void

TAssistant::updateTranslation() const {

  setTranslation(m_idEnabled, tr("Enabled"));

  setTranslation(m_idMagnetism, tr("Magnetism"));

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::onSetDefaults() {

  setEnabled(true);

  setMagnetism(1.0);

  for(TAssistantPointMap::iterator i = m_points.begin(); i != m_points.end(); ++i)

    i->second.position = i->second.defPosition;

  fixPoints();

  fixData();

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::fixPoints()

  { onFixPoints(); }

//---------------------------------------------------------------------------------------------------

bool

TAssistant::move(const TPointD &position) {

  TPointD d = position - getBasePoint().position;

  if (d != TPointD()) {

    for(TAssistantPointMap::iterator i = m_points.begin(); i != m_points.end(); ++i)

      i->second.position += d;

    fixPoints();

    return true;

  }

  return false;

}

//---------------------------------------------------------------------------------------------------

bool

TAssistant::movePoint(const TStringId &name, const TPointD &position) {

  TAssistantPointMap::iterator i = m_points.find(name);

  if (i != m_points.end() && i->second.position != position) {

    onMovePoint(i->second, position);

    return true;

  }

  return false;

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::setPointSelection(const TStringId &name, bool selected)  const {

  if (const TAssistantPoint *p = findPoint(name))

    p->selected = selected;

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::setAllPointsSelection(bool selected) const {

  for(TAssistantPointMap::const_iterator i = points().begin(); i != points().end(); ++i)

    i->second.selected = selected;

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::onDataChanged(const TVariant &value) {

  const TVariant& pointsData = data()[m_idPoints];

  TVariantPathEntry entry;

  if (&value == &data() || &value == &pointsData)

    onAllDataChanged();

  else

  if (pointsData.getChildPathEntry(value, entry) && entry.isField()) {

    const TVariant& pointData = pointsData[entry];

    TPointD position = TPointD(

      pointData[m_idX].getDouble(),

      pointData[m_idY].getDouble() );

    movePoint(entry.field(), position);

  } else

  if (data().getChildPathEntry(value, entry) && entry.isField()) {

    onDataFieldChanged(entry.field(), data()[entry.field()]);

  }

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::onDataFieldChanged(const TStringId &name, const TVariant &value) {

  updateProperty(name, value);

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::onAllDataChanged() {

  const TVariant& pointsData = data()[m_idPoints];

  for(TAssistantPointMap::iterator i = m_points.begin(); i != m_points.end(); ++i) {

    const TVariant& pointData = pointsData[i->first];

    i->second.position = TPointD(

      pointData[m_idX].getDouble(),

      pointData[m_idY].getDouble() );

  }

  fixPoints();

  updateProperties();

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::onFixPoints()

  { }

//---------------------------------------------------------------------------------------------------

void

TAssistant::onMovePoint(TAssistantPoint &point, const TPointD &position)

  { point.position = position; }

//---------------------------------------------------------------------------------------------------

void

TAssistant::onFixData() {

  TVariant& pointsData = data()[m_idPoints];

  for(TAssistantPointMap::const_iterator i = points().begin(); i != points().end(); ++i) {

    TVariant& pointData = pointsData[i->first];

    pointData[m_idX].setDouble( i->second.position.x );

    pointData[m_idY].setDouble( i->second.position.y );

  }

  setMagnetism( std::max(0.0, std::min(1.0, getMagnetism())) );

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::updateProperties() {

  const TVariantMap &map = data().getMap();

  for(TVariantMap::const_iterator i = map.begin(); i != map.end(); ++i)

    if (i->first != m_idPoints)

      updateProperty(i->first, i->second);

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::updateProperty(const TStringId &name, const TVariant &value) {

  TProperty *property = m_properties.getProperty(name);

  if (!property)

    return;

  if (TBoolProperty *boolProperty = dynamic_cast<tboolproperty*>(property)) {</tboolproperty*>

    boolProperty->setValue( value.getBool() );

  } else

  if (TDoubleProperty *doubleProperty = dynamic_cast<tdoubleproperty*>(property)) {</tdoubleproperty*>

    doubleProperty->setValue( value.getDouble() );

  } else

  if (TStringProperty *stringProperty = dynamic_cast<tstringproperty*>(property)) {</tstringproperty*>

    stringProperty->setValue( to_wstring(value.getString()) );

  } else

  if (TEnumProperty *enumProperty = dynamic_cast<tenumproperty*>(property)) {</tenumproperty*>

    enumProperty->setValue( to_wstring(value.getString()) );

  }

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::onPropertyChanged(const TStringId &name) {

  TProperty *property = m_properties.getProperty(name);

  if (!property)

    return;

  if (name == m_idPoints)

    return;

  if (TBoolProperty *boolProperty = dynamic_cast<tboolproperty*>(property)) {</tboolproperty*>

    data()[name].setBool( boolProperty->getValue() );

  } else

  if (TDoubleProperty *doubleProperty = dynamic_cast<tdoubleproperty*>(property)) {</tdoubleproperty*>

    data()[name].setDouble( doubleProperty->getValue() );

  } else

  if (TStringProperty *stringProperty = dynamic_cast<tstringproperty*>(property)) {</tstringproperty*>

    data()[name].setString( to_string(stringProperty->getValue()) );

  } else

  if (TEnumProperty *enumProperty = dynamic_cast<tenumproperty*>(property)) {</tenumproperty*>

    data()[name].setString( to_string(enumProperty->getValue()) );

  }

}

//---------------------------------------------------------------------------------------------------

double

TAssistant::getDrawingAlpha(bool enabled) const

  { return enabled && this->getEnabled() ? 0.5 : 0.25; }

//---------------------------------------------------------------------------------------------------

double

TAssistant::getDrawingGridAlpha() const

  { return 0.2; }

//---------------------------------------------------------------------------------------------------

void

TAssistant::drawSegment(const TPointD &p0, const TPointD &p1, double pixelSize, double alpha) const {

  double colorBlack[4] = { 0.0, 0.0, 0.0, alpha };

  double colorWhite[4] = { 1.0, 1.0, 1.0, alpha };

  glPushAttrib(GL_ALL_ATTRIB_BITS);

  tglEnableBlending();

  tglEnableLineSmooth(true, 1.0 * line_width_scale);

  TPointD d = p1 - p0;

  double k = norm2(d);

  if (k > TConsts::epsilon*TConsts::epsilon) {

    k = 0.5*pixelSize*line_width_scale/sqrt(k);

    d = TPointD(-k*d.y, k*d.x);

    glColor4dv(colorWhite);

    tglDrawSegment(p0 - d, p1 - d);

    glColor4dv(colorBlack);

    tglDrawSegment(p0 + d, p1 + d);

  }

  glPopAttrib();

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::drawDot(const TPointD &p, double alpha) const {

  double colorBlack[4] = { 0.0, 0.0, 0.0, alpha };

  double colorWhite[4] = { 1.0, 1.0, 1.0, alpha };

  glPushAttrib(GL_ALL_ATTRIB_BITS);

  tglEnableBlending();

  glColor4dv(colorWhite);

  tglEnablePointSmooth(6.0);

  glBegin(GL_POINTS);

  glVertex2d(p.x, p.y);

  glEnd();

  glColor4dv(colorBlack);

  tglEnablePointSmooth(3.0);

  glBegin(GL_POINTS);

  glVertex2d(p.x, p.y);

  glEnd();

  glPopAttrib();

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::drawPoint(const TAssistantPoint &point, double pixelSize) const {

  if (!point.visible) return;

  double radius = point.radius;

  double crossSize = 1.2*radius;

  double alpha = 0.5;

  double colorBlack[4] = { 0.0, 0.0, 0.0, alpha };

  double colorGray[4]  = { 0.5, 0.5, 0.5, alpha };

  double colorWhite[4] = { 1.0, 1.0, 1.0, alpha };

  double width = 1.5;

  if (point.selected) {

    colorBlack[2] = 1.0;

    colorGray[2] = 1.0;

    width = 2.0;

  }

  glPushAttrib(GL_ALL_ATTRIB_BITS);

  // fill

  tglEnableBlending();

  if (point.type == TAssistantPoint::CircleFill) {

    glColor4dv(colorGray);

    tglDrawDisk(point.position, radius*pixelSize);

  }

  TPointD crossDx(pixelSize*crossSize, 0.0);

  TPointD crossDy(0.0, pixelSize*crossSize);

  TPointD gridDx(pixelSize*radius, 0.0);

  TPointD gridDy(0.0, pixelSize*radius);

  // back line

  tglEnableLineSmooth(true, 2.0*width*line_width_scale);

  glColor4dv(colorWhite);

  if (point.type == TAssistantPoint::CircleCross) {

    tglDrawSegment(point.position - crossDx, point.position + crossDx);

    tglDrawSegment(point.position - crossDy, point.position + crossDy);

  }

  tglDrawCircle(point.position, radius*pixelSize);

  // front line

  glLineWidth(width * line_width_scale);

  glColor4dv(colorBlack);

  if (point.type == TAssistantPoint::CircleCross) {

    tglDrawSegment(point.position - crossDx, point.position + crossDx);

    tglDrawSegment(point.position - crossDy, point.position + crossDy);

  }

  tglDrawCircle(point.position, radius*pixelSize);

  // dots

  switch(point.type) {

  case TAssistantPoint::CircleDoubleDots:

    drawDot(point.position - gridDx*0.5, alpha);

    drawDot(point.position + gridDx*0.5, alpha);

    drawDot(point.position - gridDy*0.5, alpha);

    drawDot(point.position + gridDy*0.5, alpha);

    //no break

  case TAssistantPoint::CircleDots:

    drawDot(point.position - gridDx, alpha);

    drawDot(point.position + gridDx, alpha);

    drawDot(point.position - gridDy, alpha);

    drawDot(point.position + gridDy, alpha);

    //no break

  case TAssistantPoint::Circle:

    drawDot(point.position, alpha);

    break;

  default:

    break;

  }

  glPopAttrib();

}

//---------------------------------------------------------------------------------------------------

void

TAssistant::getGuidelines(const TPointD &position, const TAffine &toTool, TGuidelineList &outGuidelines) const

  { }

//---------------------------------------------------------------------------------------------------

void

TAssistant::draw(TToolViewer *viewer, bool enabled) const

  { }

//---------------------------------------------------------------------------------------------------

void

TAssistant::drawEdit(TToolViewer *viewer) const {

  // paint all points

  draw(viewer);

  double pixelSize = sqrt(tglGetPixelSize2());

  for(TAssistantPointMap::const_iterator i = points().begin(); i != points().end(); ++i)

    drawPoint(i->second, pixelSize);

}

//---------------------------------------------------------------------------------------------------

bool

TAssistant::calcPerspectiveStep(

  double minStep,

  double minX,

  double maxX,

  double x0,

  double x1,

  double x2,

  double &outK,

  double &outMin,

  double &outMax )

{

  outK = outMin = outMax = 0.0;

  double dx1 = x1 - x0;

  double dx2 = x2 - x0;

  if (fabs(dx1) <= TConsts::epsilon) return false;

  if (fabs(dx2) <= TConsts::epsilon) return false;

  if ((dx1 < 0.0) != (dx2 < 0.0)) dx2 = -dx2;

  if (fabs(dx2 - dx1) <= minStep) return false;

  if (fabs(dx2) < fabs(dx1)) std::swap(dx1, dx2);

  if (x0 <= minX + TConsts::epsilon && dx1 < 0.0) return false;

  if (x0 >= maxX - TConsts::epsilon && dx1 > 0.0) return false;

  outK = dx2/dx1;

  double minI = log(minStep/fabs(dx1*(1.0 - 1.0/outK)))/log(outK);

  outMin = dx1*pow(outK, floor(minI - TConsts::epsilon));

  if (fabs(outMin) < TConsts::epsilon) return false;

  outMax = (dx1 > 0.0 ? maxX : minX) - x0;

  return true;

}