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