28 integer pow_ii(ap, bp) integer *ap, *bp;
30 integer pow_ii(integer *ap, integer *bp)
43 return x == 0 ? 0 : 1/x;
62 double r_sign(a,b) real *a, *b;
64 double r_sign(real *a, real *b)
68 x = (*a >= 0 ? *a : - *a);
69 return( *b >= 0 ? x : -x);
76 integer i_nint(x) real *x;
80 integer i_nint(real *x)
83 return( (integer)((*x)>=0 ?
84 floor(*x + .5) : -floor(.5 - *x)) );