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