trig.u

操作系统SunOS 4.1.3版本的源码

|*
|*	Single Precision Sin, Cos, Tan, Sincos
|*
routine xs.sin w01000000000 las ssub
	;		las csla;		reg;		;	;	ptr5!c_spio4
	;		;			;		;	rcssp ptr5;	lpreg!0
	call strig;	lbsf scmpm csla;	reg;		;	;
routine xs.cos w01000000010 las ssub
	;		las csla;		reg;		;	;	ptr5!c_spio4
	;		;			;		;	rcssp ptr5;	lpreg!1
	call strig;	lbsf scmpm csla;	reg;		;	;
#ifdef STAN
routine xs.tan w01000000100 las ssub
	;		las csla;		reg;		;	;	ptr5!c_spio4
	;		;			;		;	rcssp ptr5;	lpreg!2
	call strig;	lbsf scmpm csla;	reg;		;	;
#endif
routine xs.sincos w01001100000 las ssub
	;		las csla;		reg;		;	;	ptr5!c_spio4
	;		;			;		;	rcssp ptr5;	lpreg!3
	call strig;	lbsf scmpm csla;	reg;		;	;


|*
|*	Double Precision Sin, Cos, Tan, Sincos
|*
routine xd.sin w01000000001 las ssub
	;		lam csla;		reg;		recclk;	rcslsw ptr2;	lpreg!0
	;		lal csla;		reg;		;	;	ptr5!c_dpio4
	call dtrig;	;			;		;	rcsmsw ptr5;
routine xd.cos w01000000011 las ssub
	;		lam csla;		reg;		recclk;	rcslsw ptr2;	lpreg!1
	;		lal csla;		reg;		;	;	ptr5!c_dpio4
	call dtrig;	;			;		;	rcsmsw ptr5;
#ifdef DTAN
routine xd.tan w01000000101 las ssub
	;		lam csla;		reg;		recclk;	rcslsw ptr2;	lpreg!2
	;		lal csla;		reg;		;	;	ptr5!c_dpio4
	call dtrig;	;			;		;	rcsmsw ptr5;
#endif
routine xd.sincos w01001100001 las ssub
	;		lam csla;		reg;		recclk;	rcslsw ptr2;	lpreg!3
	;		lal csla;		reg;		;	;	ptr5!c_dpio4
	call dtrig;	;			;		;	rcsmsw ptr5;


|*
|*	Single Precision Sin, Cos, Sincos
|*
|*  enter with:
|*	comparison of x and spio4 just starting in ALU
|*	lpreg is 0 for sin, 1 for cos, 2 for tan, 3 for sincos
|*
strig:
	;		;			;		recclk;	;
	;		;			;		;	;	ptr5!c_sthreshold
	;		;			;		;	rcssp ptr5;	ptr5!c_dSCRATCH_0
	;		lbsf scmpm csla;	reg;		;	ptr5;
	;		csuxa usp;		;		recoe;	ramwe rcssp ptr5;
	;		;			;		recoe;	ramwe ptr5;	ptr5!c_spi3o4
	jgt strig.1;	woea;			;		recclk;	rcssp ptr1;
	;		;			;		recclk;	rcssp ptr5;	ptr5!c_dSCRATCH_4
|*		save x in SCRATCH_4 for restoration
	;		csuxa usp;		;		;	ptr5;
	;		csuxa usp;		;		recoe;	ramwe rcssp ptr5;
	jle strig.2;	woea;			;		recoe;	ramwe ptr5;
	jwerr soberr;	woea;			;		;	;
	call strig.ex0;	;			;		;	;
strig.a:
	jloop0 ssinapprox1;;			;		;	;	lpreg-
	jloop0 scosapprox1;;			;		;	;	lpreg-
	jloop0 stanapprox1;;			;		;	;	lpreg-
	call scosapprox;;			;		;	;
	call swrt4;	;			;		;	;
	call ssinapprox;;			dtor;		;	ramwe ptr4;
	call swrt1.4;	;			;		;	;

swrt1.4:
|* this is for sincos in case we get an error and have to restore reg4
	;		;			;		;	;
	;		;			;		;	;
	;		csuxa usp;		;		;	;
	;		;			;		;	ptr1;
	;		woea;			dtor;		;	ramwe rcssp ptr1;
	call sreturninexact.4;;			dtor;		;	ramwe ptr1;

swrt4:
|		for sincos we need to save ptr4 for possible later restoration
	;		;			;		recclk;	rcssp ptr4;	ptr5!c_dSCRATCH_11
	;		;			;		recclk;	ptr5;
	;		;			;		recoe;	ramwe rcssp ptr5;
	;		csuxa usp;		;		recoe;	ramwe ptr5;
	;		;			;		;	ptr4;
	rtn;		woea;			dtor;		;	ramwe rcssp ptr4;


strig.2:
|*  abs(x) <= threshhold
	jloop0 ssin0;	;			;		;	;	lpreg-
	jloop0 scos0;	;			;		;	;	lpreg-
	;		;			;		;	;	lpreg-
stan0:
ssin0:
	;		lf scmp0 csla;		;		recclk;	rcssp ptr2;
	;		;			;		recclk;	ptr1;
	;		;			;		recoe;	ramwe rcssp ptr1;
	jloop0 ssincos0;;			;		recoe;	ramwe ptr1;
	call wait1;	csuxa usp;		;		;	;
	jeq sreturnzerostat;woea;		;		;	;
	call sreturninexact.4;;			;		;	;

ssincos0:
|		for sincos we need to save ptr4 for possible later restoration
	;		lf scmp0 csla;		;		recclk;	rcssp ptr4;	ptr5!c_dSCRATCH_11
	;		;			;		recclk;	ptr5;
	;		;			;		recoe;	ramwe rcssp ptr5;
	;		;			;		recoe;	ramwe ptr5;	ptr5!sone
	;		csuxa usp;		;		recclk;	rcssp ptr5;
	;		;			;		recclk;	ptr4;
	jeq sreturnzerostat;woea;		;		recoe;	ramwe rcssp ptr4;
	call sreturninexact.4;;			;		recoe;	ramwe ptr4;

scos0:
	;		lf scmp0 csla;		;		;	;	ptr5!sone
	;		;			;		recclk;	rcssp ptr5;
	;		;			;		recclk;	ptr1;
	;		;			;		recoe;	ramwe rcssp ptr1;
	call wait1;	csuxa usp;		;		recoe;	ramwe ptr1;
	jeq sreturnexact;woea;			;		;	;
	call sreturninexact;;			;		;	;

strig.1:
|*  pi/4 < abs(x)
	call wait1;	lbsf scmpm csla;	reg;		;	;
	;		;			;		;	;	ptr5!c_spi5o4
	;		;			;		;	rcssp ptr5;	ptr5!c_dSCRATCH_4
|*		save x in SCRATCH_4 for restoration
	;		lbs csla;		reg;		;	ptr5;
	;		lf scmp0 csla csuxa usp;;		recoe;	ramwe rcssp ptr5;
	;		;			;		recoe;	ramwe ptr5;
	jgt strig.3;	woea;			;		;	;
	;		;			;		;	;
|*  pi/4 < abs(x) <= 3*pi/4
	;		csuxa usp;		;		;	;
	;		;			;		;	;	ptr5!c_dpio2a
	jge strig.11;	woea;			;		;	;
	;		;			;		;	rcsmsw ptr5;

|*  -pi/4 > x > -3*pi/4
|* x = (dpio2b + (dx + dpio2a))
|* ex = (pio2b - (dble(x) - (dx + dpio2a)))
	;		lbm csla;		reg;		;	rcslsw ptr5;
	call wait1;	lblf dadd csla;		reg;		;	;
	;		;			;		;	;	ptr5!c_dpio2b
	;		;			;		;	rcsmsw ptr5;
	;		lam csla;		reg;		;	rcslsw ptr5;
	;		lal csla;		reg;		;	;
	;		lbm csla csuxa umsw;	;		;	;
	call wait3;	lblf dadd csla csuxa ulsw;;		;	;
	;		lam csla csuxa umsw;	;		;	;
	call wait3;	lalf dcvts csla csuxa ulsw;;		;	;
|		*start dble(x) - (dx+PIO2A)
	;		las csla csuxa usp;	;		;	;	ptr5!c_dSCRATCH_0
	;		lf dsub csla;		;		;	ptr5;
|		*save x in c_dSCRATCH_0
	;		woea;			dtor;		;	ramwe rcssp ptr5;
	call wait1;	;			dtor;		;	ramwe ptr5;
	;		lbm csla csuxa umsw;	;		;	;	ptr5!c_dpio2b
	;		lbl csla csuxa ulsw;	;		;	rcsmsw ptr5;
	;		lam csla;		reg;		;	rcslsw ptr5;
|		*start PIO2B - (dble(x) - (dx+PIO2A))
	call wait3;	lalf dsub csla;		reg;		;	;
	;		lam csla csuxa umsw;	;		;	;
	call wait3;	lalf dcvts csla csuxa ulsw;;		;	;
	;		csuxa usp;		;		;	;	ptr5!c_dSCRATCH_1
	;		;			;		;	ptr5;
	;		woea;			dtor;		;	ramwe rcssp ptr5;
strig.e:
|		*Note that here (and elsewhere) we can enter and generate a 
|		*write enable without a chip select
	jloop0 scosnapprox1;;			dtor;		;	ramwe ptr5;	lpreg-
	jloop0 ssinapprox1;;			;		;	;	lpreg-
	jloop0 scotnapprox1;;			;		;	;	lpreg-

|*	sincos
	call ssinapprox;;			;		;	;
	call swrt4;	;			;		;	;
	call scosapprox;;			dtor;		;	ramwe ptr4;
	call snegx;	;			;		;	;
	call swrt1.4;	;			;		;	;

strig.11:
|*  pi/4 < x < 3*pi/4
|*  x = TOREAL(-PIO2B - (PIO2A - dx))
|*  ex = TOREAL(-PIO2B - (dble(x) + (PIO2A - dx)))
	;		lam csla;		reg;		;	rcslsw ptr5;	ptr5!c_dSCRATCH_0
	;		lal csla;		reg;		;	rcssp ptr5;
	;		lbs csla;		reg;		;	;
	call wait3;	lf dsub csla;		;		;	;
	;		lbm csla csuxa umsw;	;		;	;	ptr5!c_dmpio2b
	;		lbl csla csuxa ulsw;	;		;	rcsmsw ptr5;
	;		lam csla;		reg;		;	rcslsw ptr5;
	call wait3;	lalf dsub csla;		reg;		;	;
	;		lam csla csuxa umsw;	;		;	;
	call wait3;	lalf dcvts csla csuxa ulsw;;		;	;
	;		las csla csuxa usp;	;		;	;	ptr5!c_dSCRATCH_0
	;		lf dadd csla;		;		;	ptr5;
	;		woea;			dtor;		;	ramwe rcssp ptr5;
	call wait1;	;			dtor;		;	ramwe ptr5;
	;		lbm csla csuxa umsw;	;		;	;	ptr5!c_dmpio2b
	;		lbl csla csuxa ulsw;	;		;	rcsmsw ptr5;
	;		lam csla;		reg;		;	rcslsw ptr5;
	call wait3;	lalf dsub csla;		reg;		;	;
	;		lam csla csuxa umsw;	;		;	;
	call wait3;	lalf dcvts csla csuxa ulsw;;		;	;
	;		csuxa usp;		;		;	;	ptr5!c_dSCRATCH_1
	;		;			;		;	ptr5;
	;		woea;			dtor;		;	ramwe rcssp ptr5;
strig.c:
	jloop0 scosapprox1;;			dtor;		;	ramwe ptr5;	lpreg-
	jloop0 ssinnapprox1;;			;		;	;	lpreg-
	jloop0 scotnapprox1;;			;		;	;	lpreg-

|*	sincos
	call ssinapprox;;			;		;	;
	call snegx;	;			;		;	;
	call swrt4;	;			;		;	;
	call scosapprox;;			dtor;		;	ramwe ptr4;
	call swrt1.4;	;			;		;	;

strig.3:
|*  abs(x) > 3*pi/4
	call wait1;	lf scmpm csla;		;		;	;
	;		;			;		;	;	ptr5!c_stwo19pi
	;		;			;		;	rcssp ptr5;
	;		lbsf scmp0 csla;	reg;		;	;
	call wait1;	csuxa usp;		;		;	;
	jgt strig.4;	woea;			;		;	;

|*  3*pi/4 < abs(x) <= 5*pi/4
|*  x = TOREAL(PIB + (dx + PIA))
|*  ex = TOREAL(PIB - (dble(x) - (dx + PIA)))
	;		csuxa usp;		;		;	;
	;		;			;		;	;	ptr5!c_dpia
	jge strig.31;	woea;			;		;	;
	;		;			;		;	rcsmsw ptr5;

|*  -3*pi/4 > x > -5*pi/4
	;		lbm csla;		reg;		;	rcslsw ptr5;
	call wait1;	lblf dadd csla;		reg;		;	;
	;		;			;		;	;	ptr5!c_dpib
	;		;			;		;	rcsmsw ptr5;
	;		lam csla;		reg;		;	rcslsw ptr5;
	;		lal csla;		reg;		;	;
	;		lbm csla csuxa umsw;	;		;	;
	call wait3;	lblf dadd csla csuxa ulsw;;		;	;
	;		lam csla csuxa umsw;	;		;	;
	call wait3;	lalf dcvts csla csuxa ulsw;;		;	;
|		*start dble(x) - (dx+PIA)
	;		las csla csuxa usp;	;		;	;	ptr5!c_dSCRATCH_0
	;		lf dsub csla;		;		;	ptr5;
	;		woea;			dtor;		;	ramwe rcssp ptr5;
	call wait1;	;			dtor;		;	ramwe ptr5;
	;		lbm csla csuxa umsw;	;		;	;	ptr5!c_dpib
	;		lbl csla csuxa ulsw;	;		;	rcsmsw ptr5;
	;		lam csla;		reg;		;	rcslsw ptr5;
|		*start PIB - (dble(x) - (dx+PIA))
	;		lalf dsub csla;		reg;		;	;
strig.32:
	call wait2;	;			;		;	;
	;		lam csla csuxa umsw;	;		;	;
	call wait3;	lalf dcvts csla csuxa ulsw;;		;	;
	;		csuxa usp;		;		;	;	ptr5!c_dSCRATCH_1
	;		;			;		;	ptr5;
	;		woea;			dtor;		;	ramwe rcssp ptr5;
strig.b:
strig.d:
	jloop0 ssinnapprox1;;			dtor;		;	ramwe ptr5;	lpreg-
	jloop0 scosnapprox1;;			;		;	;	lpreg-
	jloop0 stanapprox1;;			;		;	;	lpreg-

|*	sincos
	call scosapprox;;			;		;	;
	call snegx;	;			;		;	;
	call swrt4;	;			;		;	;
	call ssinapprox;;			dtor;		;	ramwe ptr4;
	call snegx;	;			;		;	;
	call swrt1.4;	;			;		;	;

strig.31:
|*  3*pi/4 < x < 5*pi/4
|*  x = TOREAL(-PIB + (dx - PIA))
|*  ex = TOREAL(-PIB - (dble(x) - (dx - PIA)))
	;		lbm csla;		reg;		;	rcslsw ptr5;
	call wait3;	lblf dsub csla;		reg;		;	;
	;		lbm csla csuxa umsw;	;		;	;	ptr5!c_dmpib
	;		lbl csla csuxa ulsw;	;		;	rcsmsw ptr5;
	;		lam csla;		reg;		;	rcslsw ptr5;
	call wait3;	lalf dadd csla;		reg;		;	;
	;		lam csla csuxa umsw;	;		;	;
	call wait3;	lalf dcvts csla csuxa ulsw;;		;	;
	;		las csla csuxa usp;	;		;	;	ptr5!c_dSCRATCH_0
	;		lf dsub csla;		;		;	ptr5;
	;		woea;			dtor;		;	ramwe rcssp ptr5;
	call wait1;	;			dtor;		;	ramwe ptr5;
	;		lbm csla csuxa umsw;	;		;	;	ptr5!c_dmpib
	;		lbl csla csuxa ulsw;	;		;	rcsmsw ptr5;
	jmp strig.32;	lam csla;		reg;		;	rcslsw ptr5;
	;		lalf dsub csla;		reg;		;	;

strig.4:
|*  abs(x) > 5*pi/4
	call wait3;	lf scmpm csla;		;		;	;
	call wait1;	csuxa usp;		;		;	;
	jgt soberr;	woea;			;		;	;
|*  5*pi/4 < abs(x) <= pi * 2**19
	call ckmode;	;			;		;	;
	;		;			;		;	;	ptr5!c_stwoopi
	;		lmode3.0 0x0 csla;	;		;	rcssp ptr5;	ptr5!c_dSCRATCH_0
	;		las cslm;		reg;		;	rcssp ptr5;
	call wait3;	lbsf smul cslm;		reg;		;	;
	call wait1;	csuxm usp;		;		;	;
	call wait3;	lasf scvtint csla woem;	;		;	;
	;		lalf sfloat csla csuxa ulsw;;		;	;
	;		;			;		;	;	ptr5!c_sfourth
	;		;			;		;	rcssp ptr5;
	;		las cslm;		reg;		;	;
	;		lbs csla csuxa usp;	;		;	;	ptr5!c_dSCRATCH_2
	;		;			;		;	ptr5;
	;		lbsf smul cslm woea;	dtor;		;	ramwe rcssp ptr5;
	call wait2;	;			dtor;		;	ramwe ptr5;
	call wait1;	csuxm usp;		;		;	;
	call wait3;	lasf scvtint csla woem;	;		;	;
	;		lalf sfloat csla csuxa ulsw;;		;	;
	;		;			;		;	;	ptr5!sfour
	;		;			;		;	rcssp ptr5;
	;		las cslm;		reg;		;	;
	call wait1;	csuxa usp;		;		;	;
	;		lbsf smul cslm woea;	;		;	;
	;		;			;		;	;	ptr5!c_dSCRATCH_2
	;		;			;		;	rcssp ptr5;
	;		las csla;		reg;		;	;
	call wait1;	csuxm usp;		;		;	;
	call wait3;	lbsf ssub csla woem;	;		;	;
	;		csuxa usp;		;		;	;	ptr5!c_dSCRATCH_3
	;		;			;		;	ptr5;
|*  save qrem
	;		woea;			dtor;		;	ramwe rcssp ptr5;
	;		;			dtor;		;	ramwe ptr5;	ptr5!c_dSCRATCH_2
	;		lmode3.0 0x2 csla;	;		;	rcssp ptr5;	ptr5!c_dpio2a
	;		las cslm;		reg;		;	rcsmsw ptr5;
	;		lbm cslm;		reg;		;	rcslsw ptr5;
	call wait3;	lblf dmul cslm;		reg;		;	;
	;		;			;		;	;	ptr5!c_dSCRATCH_0
	;		;			;		;	rcssp ptr5;
	;		las csla csuxm umsw;	reg;		;	;
	;		csuxm ulsw;		;		;	;
	;		lbm csla woem;		;		;	;
	;		lblf dsub csla woem;	;		;	;	ptr5!c_dpio2b
	;		;			;		;	rcsmsw ptr5;
	;		lbm cslm;		reg;		;	rcslsw ptr5;
	;		lblf dmul cslm;		reg;		;	;
	;		lam csla csuxa umsw;	;		;	;	ptr5!c_dSCRATCH_5
	;		lal csla csuxa ulsw;	;		;	ptr5;
|		*store dx - q * pio2a
	;		woea;			dtor;		;	ramwe rcsmsw ptr5;
	;		woea;			dtor;		;	ramwe rcslsw ptr5;
	;		;			dtor;		;	ramwe ptr5;
	;		csuxm umsw;		;		;	;	ptr5!c_dSCRATCH_6
	;		csuxm ulsw;		;		;	ptr5;
|		*store q * pio2b
	;		lbm csla woem;		dtor;		;	ramwe rcsmsw ptr5;
	;		lblf dsub csla woem;	dtor;		;	ramwe rcslsw ptr5;
	call wait2;	;			dtor;		;	ramwe ptr5;
	;		lam csla csuxa umsw;	;		;	;
	call wait2;	lalf dcvts csla csuxa ulsw;;		;	;
	;		;			;		;	;	ptr5!c_dSCRATCH_5
	;		lbs csla csuxa usp;	;		;	rcsmsw ptr5;
	;		lam csla;		reg;		;	rcslsw ptr5;
	;		lalf dsub csla;		reg;		;	;	ptr5!c_dSCRATCH_0
	;		csuxa usp;		;		;	;
	;		;			;		;	ptr5;
|		*store x
	;		woea;			dtor;		;	ramwe rcssp ptr5;
	;		;			dtor;		;	ramwe ptr5;
	;		lam csla csuxa umsw;	;		;	;	ptr5!c_dSCRATCH_6
	;		lal csla csuxa ulsw;	;		;	rcsmsw ptr5;
	;		lbm csla;		reg;		;	rcslsw ptr5;
	call wait3;	lblf dsub csla;		reg;		;	;
	;		lam csla csuxa umsw;	;		;	;
	call wait3;	lalf dcvts csla csuxa ulsw;;		;	;
	;		csuxa usp;		;		;	;	ptr5!c_dSCRATCH_1
	;		;			;		;	ptr5;
|		*store ex
	;		woea;			dtor;		;	ramwe rcssp ptr5;
	;		;			dtor;		;	ramwe ptr5;	ptr5!c_dSCRATCH_3
	;		;			;		;	rcssp ptr5;	ptr5!sone
	;		lasf scmp0 csla;	reg;		;	rcssp ptr5;
	call wait2;	lbs csla;		reg;		;	;
	;		lf sadd csla csuxa usp;	;		;	;
	;		csuxa usp;		;		;	;
	jlt strig.41;	woea;			;		;	;
	jeq strig.a;	woea;			;		;	;
	;		;			;		;	;
|*  qrem > 0
	call wait3;	lf ssub csla;		;		;	;
	;		csuxa usp;		;		;	;
	;		;			;		;	;
	jeq strig.c;	woea;			;		;	;
	call strig.b;	;			;		;	;

strig.41:
	call wait1;	csuxa usp;		;		;	;
	jeq strig.e;	woea;			;		;	;
	call strig.d;	;			;		;	;

snegx:
	;		;			;		;	;
	;		;			;		;	;
	rtn;		lasf sneg csla csuxa usp;;		;	;

|*	set "ex" to 0 for call to approx routines
strig.ex0:
	;		;			;		;	;	ptr5!szero
	;		;			;		recclk;	rcssp ptr5;	ptr5!c_dSCRATCH_1
	;		;			;		recclk;	ptr5;
	;		;			;		recoe;	ramwe rcssp ptr5;
	rtn;		;			;		recoe;	ramwe ptr5;

|*
|*	Single Precision Cos Approximation
|*
|*  enter with:
|*	x in c_dSCRATCH_0
|*	ex in c_dSCRATCH_1
|*  exit with:
|*	ALU starting
|*
scosapprox:
	;		;			;		;	;	ptr5!c_dSCRATCH_0
	;		;			;		;	rcssp ptr5;
	;		las cslm;		reg;		;	rcssp ptr5;
	;		lbsf smul cslm;		reg;		;	;
	;		;			;		;	;
	;		;			;		;	;
	;		;			;		;	;	ptr5!c_scos2
	;		las cslm csuxm usp;	;		;	rcssp ptr5;
	;		lbsf smul cslm;		reg;		;	;	lpreg!2
scos.app.lp:
	;		;			;		;	;	ptr5-
	;		;			;		;	;
	;		;			;		;	rcssp ptr5;
	;		las csla csuxm usp;	reg;		;	;	lpreg-
	;		;			;		;	;
	;		lbsf sadd csla woem;	;		;	;
	;		;			;		;	;
	;		;			;		;	;
	;		;			;		;	;
	;		csuxa usp;		;		;	;
	jloop scos.app.lp;;			;		;	;
	;		lbsf smul cslm woea;	;		;	;
	;		;			;		;	;
	;		;			;		;	;
	;		;			;		;	;
	;		lbsf smul cslm csuxm usp;;		;	;
	;		;			;		;	;
	;		;			;		;	;	ptr5!shalf
	;		;			;		;	rcssp ptr5;
|		*start z * .5
	;		lbsf smul cslm;		reg;		;	;
	;		csuxm usp;		;		;	;
	;		;			;		;	;	ptr5!c_dSCRATCH_0
|		* put ca into a reg of ALU
	;		las csla woem;		;		;	;
	;		;			;		;	rcssp ptr5;	ptr5!c_dSCRATCH_1
	;		las cslm;		reg;		;	rcssp ptr5;
	;		lbsf smul cslm;		reg;		;	;
	;		csuxm usp;		;		;	;	ptr5!c_dSCRATCH_2
	;		;			;		;	ptr5;
|		*hold z*HALF in c_dSCRATCH_2 for now
	;		woem;			dtor;		;	ramwe rcssp ptr5;
	;		csuxm usp;		dtor;		;	ramwe ptr5;
	;		;			;		;	;
|		*start ca - ex*x
	;		lbsf ssub csla woem;	;		;	;
	;		;			;		;	;
	;		;			;		;	;	ptr5!c_dSCRATCH_2
	;		;			;		;	rcssp ptr5;	ptr5!c_scosthresh
	;		las csla;		reg;		;	rcssp ptr5;
	;		lbsf scmp csla;		reg;		;	;
	;		csuxa usp;		;		;	;	ptr5!c_dSCRATCH_3
	;		;			;		;	ptr5;
|		*hold ca in c_dSCRATCH_3 for now
	;		woea;			dtor;		;	ramwe rcssp ptr5;
	;		csuxa usp;		dtor;		;	ramwe ptr5;
	;		;			;		;	;	ptr5!shalf
	jlt scos.1;	woea;			;		;	;
	;		;			;		;	rcssp ptr5;	ptr5!c_dSCRATCH_3

|*  z >= scosthresh
	;		lbsf ssub csla;		reg;		;	;
	;		;			;		;	;
	;		;			;		;	;
	;		;			;		;	;
	;		las csla csuxa usp;	;		;	rcssp ptr5;
	;		lbsf ssub csla;		reg;		;	;
	;		;			;		;	;
	;		;			;		;	;
	;		;			;		;	;	ptr5!shalf
	;		lbs csla csuxa usp;	;		;	rcssp ptr5;
	rtn;		lasf ssub csla;		reg;		;	;


scos.1:
|*  z < scosthresh
	;		;			;		;	rcssp ptr5;
	;		lbsf ssub csla;		reg;		;	;
	;		;			;		;	;
	;		;			;		;	;
	;		;			;		;	;	ptr5!sone
	;		lbs csla csuxa usp;	;		;	rcssp ptr5;
	rtn;		lasf ssub csla;		reg;		;	;


|*
|*	Single Precision Sin Approximation