|
|
Power 6/32 Unix version 1.21
/* THIS IS BASED ON TAHOE FLOATING POINT REPR. ! */
#include <FP.h>
double d_sign(a,b)
double *a, *b;
{
double x;
/*
x = (*a >= 0 ? *a : - *a);
return( *b >= 0 ? x : -x);
*/
x = *a;
if ((*a < 0) || (*b < 0))
*(unsigned long *)&x ^= SIGN_BIT;
return(x);
}
This archive runs on limited infrastructure. Preserving old code on modern bandwidth. Automated agents are requested to crawl responsibly.