f2clib.c

Go to the documentation of this file.

/*

$Log$
Revision 1.14  2003/02/12 13:59:15  matteo
mer feb 12 14:56:57 CET 2003

Revision 1.1.1.1  2003/02/12 13:59:15  matteo
mer feb 12 14:56:57 CET 2003

Revision 1.2  2000/01/05 08:20:39  markster
Some OSS fixes and a few lpc changes to make it actually work

 * Revision 1.1  1996/08/19  22:32:10  jaf
 * Initial revision
 *

*/

/*
 * f2clib.c
 *
 * SCCS ID:  @(#)f2clib.c 1.2 96/05/19
 */

#include "f2c.h"

#ifdef KR_headers
integer pow_ii(ap, bp) integer *ap, *bp;
#else
integer pow_ii(integer *ap, integer *bp)
#endif
{
   integer pow, x, n;
   unsigned long u;

   x = *ap;
   n = *bp;

   if (n <= 0) {
      if (n == 0 || x == 1)
         return 1;
      if (x != -1)
         return x == 0 ? 0 : 1/x;
      n = -n;
      }
   u = n;
   for(pow = 1; ; )
      {
      if(u & 01)
         pow *= x;
      if(u >>= 1)
         x *= x;
      else
         break;
      }
   return(pow);
   }



#ifdef KR_headers
double r_sign(a,b) real *a, *b;
#else
double r_sign(real *a, real *b)
#endif
{
double x;
x = (*a >= 0 ? *a : - *a);
return( *b >= 0 ? x : -x);
}



#ifdef KR_headers
double floor();
integer i_nint(x) real *x;
#else
#undef abs
#include "math.h"
integer i_nint(real *x)
#endif
{
return( (integer)((*x)>=0 ?
   floor(*x + .5) : -floor(.5 - *x)) );
}