static char SCCSID[]="@(#)softarc.c	1.2";
/* <: t-5 d :> */
#include <stdio.h>
#include "../include/dev.h"
#define FAC   	(0.01745)
#define RADIAN	(57.2957)
static int cnt,ccw,circ,*ptr;
static double x1,y1,x2,y2,x3,y3,xc,yc,rad,x1x2,y1y2,x1x3,y1y3,xpt,ypt;
double atan2(),sin(),cos(),sqrt();

softarc(cdptr)
struct command *cdptr;
{	int x0,y0,points,ang,ang1,ang2,ainc,wt,st,col;
/*
 * softarc decomposes arcs of a circle into chord points using the
 * sin/cos function.
 */
	if(cdptr->cnt !=8) {displines(cdptr);return(FAIL);}
	ptr=cdptr->aptr;
	x1 = *ptr++;  y1 = *ptr++;
	x2 = *ptr++;  y2 = *ptr++;
	x3 = *ptr++;  y3 = *ptr++;
	x1x2 = x2-x1; y1y2 = y2-y1; x1x3 = x3-x1; y1y3 = y3-y1;
	if(!((x1x2 *y1y3)-(y1y2 * x1x3)) &&   /*  straight line  */
	  !(x1==x3 && y1==y3)){
		displines(cdptr);
		return(FAIL);
	}
	arcrad(); /*  sets xc,yc,cco,and rad */

	wt=cdptr->weight; st=cdptr->style; col=cdptr->color;
	ang1 = (atan2(y1-yc,x1-xc) * RADIAN) + 0.5;
	ang2 = (atan2(y3-yc,x3-xc) * RADIAN) + 0.5;
	if(ang1 < 0)ang1 += 360;  if(ang2 < 0)ang2 += 360;
	if(circ)ang2 = ang1 + 360;
	if(ccw)ang = ang2 - ang1;
	else	ang = ang1 - ang2;
	if(ang < 0)ang += 360;
	points=ang / (ainc = 6);  /*  chord pts every 6 degrees */
	if(!ccw)ainc *= -1;
	x0=(int)x1;  y0=(int)y1;
	for(ang = ang1 + ainc; (--points > 0); ang += ainc) {
		xpt = gslclipu((long)(xc + cos(ang * FAC) * rad));
		ypt = gslclipu((long)(yc + sin(ang * FAC) * rad));
		clipseg(x0,y0,(int)xpt,(int)ypt,wt,st,col); x0=(int)xpt; y0=(int)ypt;
	}
	clipseg(x0,y0,(int)x3,(int)y3,wt,st,col);
	return(SUCCESS);
}

static
arcrad()
{
	double den,fxc,fyc,v1,v2,dxx,dyy;
/*
 *  arcrad computes the center point
 *  radius and arc orientation
 */
	circ = 0;  ccw = 1;
	v1 = (x1x2 * x1x2) + (y1y2 * y1y2);
	v2 = (x1x3 * x1x3) + (y1y3 * y1y3);
	if(v2 == 0.0) {	  /*  process full circle  */
		circ = 1;
		xc = (x1 + x2)/2.0;
		yc = (y1 + y2)/2.0;
		dxx = x1 -xc;   dyy = y1 - yc;
		rad = sqrt((dxx * dxx) + (dyy * dyy));
	}
	else  {
		den = 2.0 * ((x1x2 * y1y3) - (y1y2 * x1x3));
		fxc = ((y1y3 * v1) - (y1y2 * v2))/den;
		fyc = ((x1x2 * v2) - (x1x3 * v1))/den;
		rad = sqrt((fxc * fxc) + (fyc * fyc));
		xc = fxc + x1;
		yc = fyc + y1;
	}
	if(circ == 0) {
		if(((x2 * y1y3)-(x3 *y1y2)-(x1 * (y3-y2))) <0) ccw=0;
	}
	return;
}

