complex.h

计算电磁学的３维矢量有限元方法通用的原程序

/***********************complex_function*******************************/
struct complex { double x,y; };
typedef struct complex complex;
complex COMplex_Cmplx();
complex COMplex_Conjg();
complex COMplex_Add();
complex COMplex_Add2();
complex COMplex_Sub();
complex COMplex_Mul();
complex Real_Mul();
complex COMplex_Div();
complex Real_Div();
complex COMplex_Null();
double  Real();
double  Aimag();
double  COMplex_Abs();
complex COMplex_Expon();
 
complex COMplex_Cmplx(x,y)
double x,y;
     {complex z; z.x = x; z.y = y; return z;}
complex COMplex_Conjg(z)
complex z;
     {return COMplex_Cmplx(z.x, -z.y);}
complex COMplex_Add(a,b)
complex a,b;
     {return COMplex_Cmplx(a.x + b.x,a.y + b.y);}
complex COMplex_Add2(a,b,c)
complex a,b,c;
     {return COMplex_Cmplx(a.x + b.x + c.x,a.y + b.y + c.y);}
complex COMplex_Sub(a,b)
complex a,b;
     {return COMplex_Cmplx(a.x - b.x,a.y - b.y);}
complex COMplex_Mul(a,b)
complex a,b;
     {return COMplex_Cmplx(a.x*b.x - a.y*b.y,a.x*b.y + a.y*b.x);}
complex Real_Mul(a,z)
double a;complex z;
     {return COMplex_Cmplx(a*z.x,a*z.y);}
complex COMplex_Div(a,b)
complex a,b;
     {double D = b.x*b.x + b.y*b.y;
      return COMplex_Cmplx((a.x*b.x + a.y*b.y)/D,(a.y*b.x - a.x*b.y)/D);}
complex Real_Div(z,a)
double a;complex z;
     {return COMplex_Cmplx(z.x/a,z.y/a);}
double Real(z)
complex z;
     {return z.x;}
double Aimag(z)
complex z;
     {return z.y;}
double COMplex_Abs(z)
complex z;
     {return sqrt(z.x*z.x + z.y*z.y);}
complex COMplex_Expon(a,z)
double a; complex z;
     {double R = exp((a*z.x));
      return COMplex_Cmplx(R*cos(z.y),(a*R)*sin(z.y));}
complex COMplex_Null()
{ complex z; z.x=0.0; z.y=0.0; return z; }
 
/***************************/