BSD 4_3_Tahoe release

[unix-history] / usr / src / usr.lib / libF77 / pow_ci.c

diff --git a/usr/src/usr.lib/libF77/pow_ci.c b/usr/src/usr.lib/libF77/pow_ci.c
index 60d7488..44874ed 100644 (file)
--- a/usr/src/usr.lib/libF77/pow_ci.c
+++ b/usr/src/usr.lib/libF77/pow_ci.c
@@ -3,22 +3,64 @@
  * All rights reserved.  The Berkeley software License Agreement
  * specifies the terms and conditions for redistribution.
  *
  * All rights reserved.  The Berkeley software License Agreement
  * specifies the terms and conditions for redistribution.
  *

- *     @(#)pow_ci.c    5.1     6/7/85
+ *     @(#)pow_ci.c    5.3     1/24/88

  */
 
 #include "complex"
 
  */
 
 #include "complex"
 

+#ifdef tahoe
+
+#define        C_MULEQ(A,B)    \
+       t = (A).real * (B).real - (A).imag * (B).imag,\
+       (A).imag = (A).real * (B).imag + (A).imag * (B).real,\
+       (A).real = t    /* A *= B */
+
+void
+pow_ci(p, a, b)        /* p = a**b  */
+       complex *p, *a;
+       long *b;
+{
+       register long n = *b;
+       register float t;
+       complex x;
+
+       x = *a;
+       p->real = (float)1, p->imag = (float)0;
+       if (!n)
+               return;
+       if (n < 0) {
+               c_div(&x, p, a);
+               n = -n;
+       }
+       while (!(n&1)) {
+               C_MULEQ(x, x);
+               n >>= 1;
+       }
+       for (*p = x; --n > 0; C_MULEQ(*p, x))
+               while (!(n&1)) {
+                       C_MULEQ(x, x);
+                       n >>= 1;
+               }
+}
+
+#else /* !tahoe */
+
+extern void pow_zi();
+
+void

 pow_ci(p, a, b)        /* p = a**b  */
 pow_ci(p, a, b)        /* p = a**b  */

-complex *p, *a;
-long int *b;
+       complex *p, *a;
+       long *b;

 {
 {

-dcomplex p1, a1;
+       dcomplex p1, a1;

 
 

-a1.dreal = a->real;
-a1.dimag = a->imag;
+       a1.dreal = a->real;
+       a1.dimag = a->imag;

 
 

-pow_zi(&p1, &a1, b);
+       pow_zi(&p1, &a1, b);

 
 

-p->real = p1.dreal;
-p->imag = p1.dimag;
+       p->real = p1.dreal;
+       p->imag = p1.dimag;

 }
 }

+
+#endif /* tahoe */