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