1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677 |
- #include "ehead.h"
- void eexp( x, y )
- unsigned short *x, *y;
- {
- unsigned short num[NE], den[NE], x2[NE];
- long i;
- unsigned short sign, expchk;
- if( ecmp(x, ezero) == 0 )
- {
- emov( eone, y );
- return;
- }
- emov(x, x2);
- expchk = x2[NE-1];
- sign = expchk & 0x8000;
- x2[NE-1] &= 0x7fff;
- expchk &= 0x7fff;
- if( expchk > (EXONE + 15) )
- {
- eclear( y );
- if( sign == 0 )
- einfin( y );
- return;
- }
- eifrac( x2, &i, num );
- if( i != 0 )
- {
- ltoe( &i, den );
- ediv( elog2, den, den );
- eifrac( den, &i, den );
- emul( elog2, den, den );
- eadd( den, num, x2 );
- }
- eldexp( x2, -1L, x2 );
- etanh( x2, x2 );
- eadd( x2, eone, num );
- eneg( x2 );
- eadd( x2, eone, den );
- ediv( den, num, y );
- if( sign )
- {
- ediv( y, eone, y );
- i = -i;
- }
- eldexp( y, i, y );
- }
|