eigenmath/rect.cpp

81 lines
946 B
C++

/* Convert complex z to rectangular form
Input: push z
Output: Result on stack
*/
#include "stdafx.h"
#include "defs.h"
void
eval_rect(void)
{
push(cadr(p1));
eval();
rect();
}
void
rect(void)
{
save();
p1 = pop();
if (car(p1) == symbol(ADD)) {
push_integer(0);
p1 = cdr(p1);
while (iscons(p1)) {
push(car(p1));
rect();
add();
p1 = cdr(p1);
}
} else {
push(p1); // mag(z) * (cos(arg(z)) + i sin(arg(z)))
mag();
push(p1);
arg();
p1 = pop();
push(p1);
cosine();
push(imaginaryunit);
push(p1);
sine();
multiply();
add();
multiply();
}
restore();
}
#if SELFTEST
static char *s[] = {
"rect(a+i*b)",
"a+i*b",
"rect(exp(a+i*b))",
"i*exp(a)*sin(b)+exp(a)*cos(b)",
"rect(1+exp(i*pi/3))",
"3/2+1/2*i*3^(1/2)",
"z=(a+b*i)/(c+d*i)",
"",
"rect(z)-real(z)-i*imag(z)",
"0",
"z=quote(z)",
"",
};
void
test_rect(void)
{
test(__FILE__, s, sizeof s / sizeof (char *));
}
#endif