1 2 3 4 5 6 7 8
#include "f2c.h" #undef abs #include <math.h> shortint h_dnnt(doublereal *x) { return (shortint)(*x >= 0. ? floor(*x + .5) : -floor(.5 - *x)); }