rect.cpp

   1 /* Convert complex z to rectangular form
   2
   3         Input:          push    z
   4
   5         Output:         Result on stack
   6 */
   7
   8 #include "stdafx.h"
   9 #include "defs.h"
  10
  11 void
  12 eval_rect(void)
  13 {
  14         push(cadr(p1));
  15         eval();
  16         rect();
  17 }
  18
  19 void
  20 rect(void)
  21 {
  22         save();
  23         p1 = pop();
  24         if (car(p1) == symbol(ADD)) {
  25                 push_integer(0);
  26                 p1 = cdr(p1);
  27                 while (iscons(p1)) {
  28                         push(car(p1));
  29                         rect();
  30                         add();
  31                         p1 = cdr(p1);
  32                 }
  33         } else {
  34                 push(p1);       // mag(z) * (cos(arg(z)) + i sin(arg(z)))
  35                 mag();
  36                 push(p1);
  37                 arg();
  38                 p1 = pop();
  39                 push(p1);
  40                 cosine();
  41                 push(imaginaryunit);
  42                 push(p1);
  43                 sine();
  44                 multiply();
  45                 add();
  46                 multiply();
  47         }
  48         restore();
  49 }
  50
  51 #if SELFTEST
  52
  53 static char *s[] = {
  54
  55         "rect(a+i*b)",
  56         "a+i*b",
  57
  58         "rect(exp(a+i*b))",
  59         "i*exp(a)*sin(b)+exp(a)*cos(b)",
  60
  61         "rect(1+exp(i*pi/3))",
  62         "3/2+1/2*i*3^(1/2)",
  63
  64         "z=(a+b*i)/(c+d*i)",
  65         "",
  66
  67         "rect(z)-real(z)-i*imag(z)",
  68         "0",
  69
  70         "z=quote(z)",
  71         "",
  72 };
  73
  74 void
  75 test_rect(void)
  76 {
  77         test(__FILE__, s, sizeof s / sizeof (char *));
  78 }
  79
  80 #endif