insilmaril@139: #include <math.h>
insilmaril@139: 
insilmaril@0: #include <qregexp.h>
insilmaril@0: #include <qpoint.h>
insilmaril@164: #include <stdlib.h>
insilmaril@0: 
insilmaril@0: #include "misc.h"
insilmaril@0: 
insilmaril@129: QString qpointToString (const QPoint &p)
insilmaril@129: {
insilmaril@129: 	return "(" + QString("%1").arg(p.x()) +","+ QString ("%1").arg (p.y()) +")";
insilmaril@129: }
insilmaril@0: 
insilmaril@408: QString qpointfToString (const QPointF &p)
insilmaril@408: {
insilmaril@408: 	return "(" + QString("%1").arg(p.x()) +","+ QString ("%1").arg (p.y()) +")";
insilmaril@408: }
insilmaril@408: 
insilmaril@0: ostream &operator<< (ostream &stream, QPoint const &p)
insilmaril@0: { 
insilmaril@366: 	stream << "("<<p.x()<<","<<p.y()<<")";
insilmaril@366: 	return stream;
insilmaril@0: }
insilmaril@0: 
insilmaril@408: ostream &operator<< (ostream &stream, QPointF const &p)
insilmaril@408: { 
insilmaril@408: 	stream << "("<<p.x()<<","<<p.y()<<")";
insilmaril@408: 	return stream;
insilmaril@408: }
insilmaril@408: 
insilmaril@408: qreal getAngle(const QPointF &p)
insilmaril@0: {	
insilmaril@0: 	// Calculate angle of vector to y-axis
insilmaril@0: 	if (p.y()==0)
insilmaril@0: 	{
insilmaril@0: 		if (p.x()>=0)
insilmaril@0: 			return M_PI_2;
insilmaril@0: 		else
insilmaril@0: 			return 3* M_PI_2;
insilmaril@0: 	} else
insilmaril@0: 	{
insilmaril@0: 		if (p.y()>0) 
insilmaril@408: 			return (qreal)(M_PI  - atan ( (qreal)(p.x()) / (qreal)(p.y()) ) );
insilmaril@0: 		else	
insilmaril@0: 			if (p.x()<0)
insilmaril@408: 				return (qreal)( 2*M_PI - atan ( (qreal)(p.x()) / (qreal)(p.y()) ) );
insilmaril@0: 			else	
insilmaril@408: 				return (qreal)( - atan ( (qreal)(p.x()) / (qreal)(p.y()) ) );
insilmaril@0: 	}	
insilmaril@0: }
insilmaril@0: 
insilmaril@408: QPointF normalise(const QPointF &p)
insilmaril@0: {	
insilmaril@0: 	// Calculate normalised position (fixed length) 
insilmaril@0: 
insilmaril@408: 	qreal px=p.x();
insilmaril@408: 	qreal py=p.y();
insilmaril@408: 	qreal x;
insilmaril@408: 	qreal y;
insilmaril@408: 	qreal r=150;
insilmaril@0: 
insilmaril@0: 	if (px==0)
insilmaril@0: 	{
insilmaril@0: 		x=0;
insilmaril@0: 		if (py>=0)
insilmaril@0: 			y=r;
insilmaril@0: 		else
insilmaril@0: 			y=-r;
insilmaril@0: 	} else
insilmaril@0: 	{
insilmaril@408: 		qreal sign;
insilmaril@408: 		qreal a;
insilmaril@0: 		if (px>0) 
insilmaril@0: 			sign=1; 
insilmaril@0: 		else 
insilmaril@0: 			sign=-1;
insilmaril@0: 		
insilmaril@0: 		a=atan (py / px);
insilmaril@0: 		x=cos (a) * r *sign;
insilmaril@0: 		y=sin (a) * r *sign;
insilmaril@0: 	}	
insilmaril@0: 	return QPoint ((int) (x),(int) (y));
insilmaril@0: }
insilmaril@0: 
insilmaril@166: 
insilmaril@408: qreal max(qreal a, qreal b)
insilmaril@0: {
insilmaril@0: 	if (a>b)
insilmaril@0: 		return a;
insilmaril@0: 	return b;
insilmaril@0: }
insilmaril@0: