61 lines
773 B
C++
61 lines
773 B
C++
|
// clock(z) = mag(z) * (-1) ^ (arg(z) / pi)
|
||
|
|
|
||
|
|
// For example clock(exp(i pi/3)) gives the result (-1)^(1/3)
|
||
|
|
|
||
|
|
#include "stdafx.h"
|
||
|
|
#include "defs.h"
|
||
|
|
|
||
|
|
void
|
||
|
|
eval_clock(void)
|
||
|
|
{
|
||
|
|
#if 1
|
||
|
|
push(cadr(p1));
|
||
|
|
eval();
|
||
|
|
p1 = pop();
|
||
|
|
push(p1);
|
||
|
|
mag();
|
||
|
|
push_integer(-1);
|
||
|
|
push(p1);
|
||
|
|
arg();
|
||
|
|
push(symbol(PI));
|
||
|
|
divide();
|
||
|
|
power();
|
||
|
|
multiply();
|
||
|
|
#else
|
||
|
|
push(cadr(p1));
|
||
|
|
eval();
|
||
|
|
p1 = pop();
|
||
|
|
push(p1);
|
||
|
|
mag();
|
||
|
|
push(symbol(E));
|
||
|
|
push(p1);
|
||
|
|
arg();
|
||
|
|
push(imaginaryunit);
|
||
|
|
multiply();
|
||
|
|
power();
|
||
|
|
multiply();
|
||
|
|
#endif
|
||
|
|
}
|
||
|
|
|
||
|
|
#if SELFTEST
|
||
|
|
|
||
|
|
static char *s[] = {
|
||
|
|
|
||
|
|
"clock(exp(i pi/3))",
|
||
|
|
"(-1)^(1/3)",
|
||
|
|
|
||
|
|
"clock(exp(-i pi/3))",
|
||
|
|
"-(-1)^(2/3)",
|
||
|
|
|
||
|
|
"rect(clock(3+4*i))", // needs sin(arctan(x)) and cos(arctan(x))
|
||
|
|
"3+4*i",
|
||
|
|
};
|
||
|
|
|
||
|
|
void
|
||
|
|
test_clock(void)
|
||
|
|
{
|
||
|
|
test(__FILE__, s, sizeof s / sizeof (char *));
|
||
|
|
}
|
||
|
|
|
||
|
|
#endif
|