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