#include "f2c.h" #undef abs #include shortint h_dnnt (doublereal * x) { return (shortint) (*x >= 0. ? floor (*x + .5) : -floor (.5 - *x)); }