diff options
author | wkj <devnull@localhost> | 2004-04-21 02:16:43 +0000 |
---|---|---|
committer | wkj <devnull@localhost> | 2004-04-21 02:16:43 +0000 |
commit | cd5bae7871bc0f0bc68b4d2a84703929a7a3c9d1 (patch) | |
tree | 336db54785d2b77113a6e570574be715c7eb7d1d /src/cmd/astro/cosadd.c | |
parent | 95f57b01e21feb457e79eaf52d593422c318024f (diff) | |
download | plan9port-cd5bae7871bc0f0bc68b4d2a84703929a7a3c9d1.tar.gz plan9port-cd5bae7871bc0f0bc68b4d2a84703929a7a3c9d1.tar.bz2 plan9port-cd5bae7871bc0f0bc68b4d2a84703929a7a3c9d1.zip |
Astro with some minor changes to placate Unix.
Diffstat (limited to 'src/cmd/astro/cosadd.c')
-rw-r--r-- | src/cmd/astro/cosadd.c | 64 |
1 files changed, 64 insertions, 0 deletions
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; +} |