home *** CD-ROM | disk | FTP | other *** search
/ Liren Large Software Subsidy 7 / 07.iso / c / c129 / 1.ddi / COMPLEX.C < prev    next >
Encoding:
C/C++ Source or Header  |  1989-02-23  |  2.9 KB  |  118 lines
  1. #include <math.h>
  2. #include "miscio.h"
  3.  
  4. struct complext{
  5.  float x,y;
  6.  
  7. };
  8.  
  9. float CMag(struct complext *CNum1)
  10. {
  11.    float _fret;
  12.  
  13.       _fret = sqrt( (*CNum1).x * (*CNum1).x +
  14.               (*CNum1).y * (*CNum1).y );
  15.    return(_fret);
  16. }
  17.  
  18.  
  19. void ComplexMath(struct complext c1,char op,
  20.                  struct complext c2,struct complext *result)
  21. {
  22.    struct complext temp1;
  23.    struct complext temp2;
  24.    float       tempr;
  25.  
  26.    switch ( op ) {
  27.       case '+':
  28.          {
  29.  
  30.             (*result).x = c1.x + c2.x;
  31.             (*result).y = c1.y + c2.y;
  32.  
  33.         }
  34.          break;
  35.       case '-':
  36.          {
  37.             (*result).x = c1.x - c2.x;
  38.             (*result).y = c1.y - c2.y;
  39.          }
  40.          break;
  41.       case '*':
  42.          {
  43.             (*result).x = c1.x * c2.x - c1.y * c2.y;
  44.             (*result).y = c1.x * c2.y + c1.y * c2.x;
  45.          }
  46.          break;
  47.       case '/':
  48.          {
  49.         (*result).x = (c1.x * c2.x + c1.y * c2.y) / ((c2.x * c2.x) + (c2.y * c2.y));
  50.             (*result).y = (c1.y * c2.x - c1.x * c2.y) / ((c2.x * c2.x) + (c2.y * c2.y));
  51.          }
  52.          break;
  53.       case 'i':
  54.          {
  55.             tempr = CMag(&c1);
  56.             (*result).y = -(c1.y / tempr);
  57.             (*result).x = c1.x / tempr;
  58.          }
  59.          break;
  60.    }
  61. }
  62.  
  63. void CExp(struct complext c1,struct complext *c2)
  64. {
  65.  
  66.    (*c2).x = exp(c1.x) * cos(c1.y);
  67.    (*c2).y = exp(c1.x) * sin(c1.y);
  68. }
  69.  
  70.  
  71. float CAngle(struct complext CNum1)
  72. {
  73.    float _fret;
  74.    float temp;
  75.    # define NearZero 1.0e-20
  76.  
  77.       {
  78.          if ( fabs(CNum1.x) < NearZero ) {
  79.             if ( fabs(CNum1.y) < NearZero ) {
  80.                _fret = 0.0;
  81.             }
  82.             if ( (CNum1.y < -(NearZero)) ) {
  83.                _fret = 3.0 * M_PI / 2.0;
  84.             }
  85.             if ( CNum1.y > NearZero ) {
  86.                _fret = M_PI / 2.0;
  87.             }
  88.             return(_fret);
  89.          }
  90.          else {
  91.             if ( fabs(CNum1.y) < NearZero ) {
  92.                if ( (CNum1.x < -(NearZero)) ) {
  93.                   _fret = M_PI;
  94.                }
  95.                if ( CNum1.x > NearZero ) {
  96.                   _fret = 0.0;
  97.                }
  98.                return(_fret);
  99.             }
  100.             else {
  101.                if ( (CNum1.x > NearZero) && (CNum1.y > NearZero) ) {
  102.                   _fret = atan(CNum1.y / CNum1.x);
  103.                }
  104.                if ( (CNum1.x < -NearZero) && (CNum1.y > NearZero) ) {
  105.                   _fret = M_PI - atan(-CNum1.y / CNum1.x);
  106.                }
  107.                if ( (CNum1.x < -NearZero) && (CNum1.y < -NearZero) ) {
  108.                   _fret = M_PI + atan(CNum1.y / CNum1.x);
  109.                }
  110.                if ( (CNum1.x > NearZero) && (CNum1.y < -NearZero) ) {
  111.                   _fret = 2.0 * M_PI - atan(-CNum1.y / CNum1.x);
  112.                }
  113.             }
  114.          }
  115.       }
  116.    return(_fret);
  117. }
  118.