00001 #include "f2c.h" 00002 00003 #ifdef KR_headers 00004 double d_sign(a,b) doublereal *a, *b; 00005 #else 00006 double d_sign(doublereal *a, doublereal *b) 00007 #endif 00008 { 00009 double x; 00010 x = (*a >= 0 ? *a : - *a); 00011 return( *b >= 0 ? x : -x); 00012 }