git.maemo.org Git - opencv/blob - 3rdparty/lapack/pow_ri.c

   1 #include "clapack.h"
   2
   3 double pow_ri(real *ap, integer *bp)
   4 {
   5     double pow, x;
   6     integer n;
   7     unsigned long u;
   8
   9     pow = 1;
  10     x = *ap;
  11     n = *bp;
  12
  13     if(n != 0)
  14     {
  15         if(n < 0)
  16         {
  17             n = -n;
  18             x = 1/x;
  19         }
  20         for(u = n; ; )
  21         {
  22             if(u & 01)
  23                 pow *= x;
  24             if(u >>= 1)
  25                 x *= x;
  26             else
  27                 break;
  28         }
  29     }
  30     return(pow);
  31 }