/*	qasin.c		*/
/* arc sine */

#include "qhead.h"

extern QELT qhalf[], qone[], qpi[];

int qasin( x, y )
QELT x[], y[];
{
QELT a[NQ], z[NQ], zz[NQ], temp[NQ];
int sign, flg;

sign = x[0];
qmov( x, a );
a[0] = 0;
if( qcmp( a, qone ) > 0 )
	{
	mtherr( "qasin", DOMAIN );
	qclear(y);
	return 0;
	}
if( qcmp( a, qhalf) > 0 )
	{
	qsub( a, qhalf, zz );	/*zz = 0.5 -a;*/
	qadd( qhalf, zz, zz );	/*zz = (zz + 0.5)/2.0;*/
	if( zz[1] > 0 )
		zz[1] -= 1;
	qsqrt( zz, z );		/* z = sqrt( zz );*/
	flg = 1;
	}
else
	{
	qmov( a, z );		/*z = a;*/
	qmul( z, z, zz );	/* zz = z * z;*/
	flg = 0;
	}

qsub( zz, qone, temp );	/* 1 - x**2 */
qsqrt( temp, temp );
qdiv( temp, z, temp );
qatn( temp, y );
if( flg != 0 )
	{
	y[1] += 1;
	qmov( qpi, a );
	a[1] -= 1;
	qsub( y, a, y );	/* z = PIO2 - z;*/
	}
y[0] = sign;
return 0;
}



int qacos( x, y )
QELT x[], y[];
{
QELT r[NQ], temp[NQ];

qmov( x, temp );
temp[0] = 0;
if( qcmp( temp, qone ) > 0 )
	{
	mtherr( "qacos", DOMAIN );
	qclear(y);
	return 0;
	}
qasin( x, r );
qmov( qpi, temp );
temp[1] -= 1;
qsub( r, temp, y );
return 0;
}