1 #include "f2c.h"
2 
3 #ifdef KR_headers
r_sign(a,b)4 double r_sign(a,b) real *a, *b;
5 #else
6 double r_sign(real *a, real *b)
7 #endif
8 {
9 double x;
10 x = (*a >= 0 ? *a : - *a);
11 return( *b >= 0 ? x : -x);
12 }
13