From cd5bae7871bc0f0bc68b4d2a84703929a7a3c9d1 Mon Sep 17 00:00:00 2001
From: wkj <devnull@localhost>
Date: Wed, 21 Apr 2004 02:16:43 +0000
Subject: Astro with some minor changes to placate Unix.

---
 src/cmd/astro/cosadd.c | 64 ++++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 64 insertions(+)
 create mode 100644 src/cmd/astro/cosadd.c

(limited to 'src/cmd/astro/cosadd.c')

diff --git a/src/cmd/astro/cosadd.c b/src/cmd/astro/cosadd.c
new file mode 100644
index 00000000..a632149a
--- /dev/null
+++ b/src/cmd/astro/cosadd.c
@@ -0,0 +1,64 @@
+#include "astro.h"
+
+
+void
+icosadd(double *fp, char *cp)
+{
+
+	cafp = fp;
+	cacp = cp;
+}
+
+double
+cosadd(int n, double coef, ...)
+{
+	double *coefp;
+	char *cp;
+	int i;
+	double sum, a1, a2;
+
+	sum = 0;
+	cp = cacp;
+
+loop:
+	a1 = *cafp++;
+	if(a1 == 0) {
+		cacp = cp;
+		return sum;
+	}
+	a2 = *cafp++;
+	i = n;
+	coefp = &coef;
+	do
+		a2 += *cp++ * *coefp++;
+	while(--i);
+	sum += a1 * cos(a2);
+	goto loop;
+}
+
+double
+sinadd(int n, double coef, ...)
+{
+	double *coefp;
+	char *cp;
+	int i;
+	double sum, a1, a2;
+
+	sum = 0;
+	cp = cacp;
+
+loop:
+	a1 = *cafp++;
+	if(a1 == 0) {
+		cacp = cp;
+		return sum;
+	}
+	a2 = *cafp++;
+	i = n;
+	coefp = &coef;
+	do
+		a2 += *cp++ * *coefp++;
+	while(--i);
+	sum += a1 * sin(a2);
+	goto loop;
+}
-- 
cgit v1.2.3