// CallCircle.cpp: implementation of the CCallCircle class.

//

//////////////////////////////////////////////////////////////////////

#include <math.h></math.h>

#include <stdlib.h></stdlib.h>

#include <search.h></search.h>

//#include "myetc.h"

#include "CallCircle.h"

//////////////////////////////////////////////////////////////////////

// Construction/Destruction

//////////////////////////////////////////////////////////////////////

int callcircle_xydwCompare(const void *a, const void *b) {

  SXYDW *aa = (SXYDW *)a;

  SXYDW *bb = (SXYDW *)b;

  if (aa->w < bb->w) return -1;

  if (aa->w > bb->w) return 1;

  return 0;

}

CCallCircle::CCallCircle(const double r) : m_r(r), m_nb(0) {

  int rr = (int)r + 1;

  rr *= 2;

  int dd  = rr * 2 + 1;

  int dd2 = dd * dd;

  if (dd2 == 0) {

    null();

    return;

  }

  m_c.reset(new SXYDW[dd2]);

  if (!m_c) throw SMemAllocError("in callCircle");

  for (int y = -rr; y <= rr; y++)

    for (int x = -rr; x <= rr; x++) {

      double d = sqrt((double)(x * x + y * y));

      if (d <= r && m_nb < dd2) {

        m_c[m_nb].x = x;

        m_c[m_nb].y = y;

        m_c[m_nb].w = d;

        m_nb++;

      }

    }

  qsort(m_c.get(), m_nb, sizeof(SXYDW), callcircle_xydwCompare);

}

void CCallCircle::null() {

  m_nb = 0;

  m_r  = 0.0;

  m_c.reset();

}

CCallCircle::~CCallCircle() { null(); }

void CCallCircle::print() {

  /*smsg_info(" --- CCallCircle ---");

  smsg_info(" m_nb=%d", m_nb);

  smsg_info(" m_r=%f", m_r);

  for( int i=0; i

          smsg_info(" %d. = (%d,%d,%f)", i, m_c[i].x,m_c[i].y,m_c[i].w);

  smsg_info(" --- ----------- ---\n");

*/

}

void CCallCircle::draw(UCHAR *drawB, const int lX, const int lY, const int xx,

                       const int yy, const double r) {

  double aa = 2.0 * r / 3.0;  // Size of antialising

  for (int i = 0; i < m_nb && m_c[i].w <= r; i++) {

    double w = m_c[i].w;

    int x    = xx + m_c[i].x;

    int y    = yy + m_c[i].y;

    if (x >= 0 && y >= 0 && x < lX && y < lY) {

      UCHAR *pDB = drawB + y * lX + x;

      if (w <= aa) {

        *pDB = (UCHAR)255;

      } else {

        double q  = 255.0 * (r - w) / (r - aa);

        q         = D_CUT_0_255(q);

        UCHAR ucq = UC_ROUND(q);

        (*pDB)    = (*pDB) < ucq ? ucq : (*pDB);

      }

    }

  }

}