From: David Wasley <dlw@ucbvax.Berkeley.EDU>
Date: Sat, 22 Jan 1983 03:16:34 +0000 (-0800)
Subject: date and time created 83/01/21 11:16:34 by dlw
X-Git-Tag: BSD-4_1c_2-Snapshot-Development~888
X-Git-Url: https://git.subgeniuskitty.com/unix-history/.git/commitdiff_plain/69fcff2c7fe80d68ef38181765bb5024dcead357

date and time created 83/01/21 11:16:34 by dlw

SCCS-vsn: usr.bin/f77/libF77/pow_di.c 1.1
---

diff --git a/usr/src/usr.bin/f77/libF77/pow_di.c b/usr/src/usr.bin/f77/libF77/pow_di.c
new file mode 100644
index 0000000000..b51e44c2fd
--- /dev/null
+++ b/usr/src/usr.bin/f77/libF77/pow_di.c
@@ -0,0 +1,38 @@
+/*
+ *	"@(#)pow_di.c	1.1"
+ */
+
+double pow_di(ap, bp)
+double *ap;
+long int *bp;
+{
+double pow, x;
+long int n;
+
+pow = 1;
+x = *ap;
+n = *bp;
+
+if(n != 0)
+	{
+	if(n < 0)
+		{
+		if(x == 0)
+			{
+			return(pow);
+			}
+		n = -n;
+		x = 1/x;
+		}
+	for( ; ; )
+		{
+		if(n & 01)
+			pow *= x;
+		if(n >>= 1)
+			x *= x;
+		else
+			break;
+		}
+	}
+return(pow);
+}