From e5888a1ffdae813d7575f5fb02275c6bb07e5199 Mon Sep 17 00:00:00 2001 From: Taru Karttunen Date: Wed, 30 Mar 2011 15:46:40 +0300 Subject: Import sources from 2011-03-30 iso image --- sys/src/cmd/astro/cosadd.c | 64 ++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 64 insertions(+) create mode 100755 sys/src/cmd/astro/cosadd.c (limited to 'sys/src/cmd/astro/cosadd.c') diff --git a/sys/src/cmd/astro/cosadd.c b/sys/src/cmd/astro/cosadd.c new file mode 100755 index 000000000..a632149a0 --- /dev/null +++ b/sys/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