d_sign X
1: double d_sign(a,b) 2: double *a, *b; 3: { 4: double x; 5: x = (*a >= 0 ? *a : - *a); 6: return( *b >= 0 ? x : -x); 7: }