1 #include "v3p_f2c.h"
2 #ifdef __cplusplus
3 extern "C" {
4 #endif
5 
6 #ifdef KR_headers
r_sign(a,b)7 double r_sign(a,b) real *a, *b;
8 #else
9 double r_sign(real *a, real *b)
10 #endif
11 {
12 double x;
13 x = (*a >= 0 ? *a : - *a);
14 return( *b >= 0 ? x : -x);
15 }
16 #ifdef __cplusplus
17 }
18 #endif
19