polar.cpp

   1 /* Convert complex z to polar form
   2
   3         Input:          push    z
   4
   5         Output:         Result on stack
   6
   7         polar(z) = mag(z) * exp(i * arg(z))
   8 */
   9
  10 #include "stdafx.h"
  11 #include "defs.h"
  12
  13 void
  14 eval_polar(void)
  15 {
  16         push(cadr(p1));
  17         eval();
  18         polar();
  19 }
  20
  21 void
  22 polar(void)
  23 {
  24         save();
  25         p1 = pop();
  26         push(p1);
  27         mag();
  28         push(imaginaryunit);
  29         push(p1);
  30         arg();
  31         multiply();
  32         exponential();
  33         multiply();
  34         restore();
  35 }
  36
  37 #if SELFTEST
  38
  39 static char *s[] = {
  40
  41         "polar(1+i)",
  42         "2^(1/2)*exp(1/4*i*pi)",
  43
  44         "polar(-1+i)",
  45         "2^(1/2)*exp(3/4*i*pi)",
  46
  47         "polar(-1-i)",
  48         "2^(1/2)*exp(-3/4*i*pi)",
  49
  50         "polar(1-i)",
  51         "2^(1/2)*exp(-1/4*i*pi)",
  52
  53         "rect(polar(3+4*i))",
  54         "3+4*i",
  55
  56         "rect(polar(-3+4*i))",
  57         "-3+4*i",
  58
  59         "rect(polar(3-4*i))",
  60         "3-4*i",
  61
  62         "rect(polar(-3-4*i))",
  63         "-3-4*i",
  64 };
  65
  66 void
  67 test_polar(void)
  68 {
  69         test(__FILE__, s, sizeof s / sizeof (char *));
  70 }
  71
  72 #endif