C++面向对象编程<四>：类comlpex的完整代码

这里给出完整代码

#ifndef __MYCOMPLEX__#define __MYCOMPLEX__class complex; complex&  __doapl (complex* ths, const complex& r);complex&  __doami (complex* ths, const complex& r);complex&  __doaml (complex* ths, const complex& r);class complex{public:  complex (double r = 0, double i = 0): re (r), im (i) { }  complex& operator += (const complex&);  complex& operator -= (const complex&);  complex& operator *= (const complex&);  complex& operator /= (const complex&);  double real () const { return re; }  double imag () const { return im; }private:  double re, im;  friend complex& __doapl (complex *, const complex&);  friend complex& __doami (complex *, const complex&);  friend complex& __doaml (complex *, const complex&);};inline complex&__doapl (complex* ths, const complex& r){  ths->re += r.re;  ths->im += r.im;  return *ths;}inline complex&complex::operator += (const complex& r){  return __doapl (this, r);}inline complex&__doami (complex* ths, const complex& r){  ths->re -= r.re;  ths->im -= r.im;  return *ths;}inline complex&complex::operator -= (const complex& r){  return __doami (this, r);}inline complex&__doaml (complex* ths, const complex& r){  double f = ths->re * r.re - ths->im * r.im;  ths->im = ths->re * r.im + ths->im * r.re;  ths->re = f;  return *ths;}inline complex&complex::operator *= (const complex& r){  return __doaml (this, r);}inline doubleimag (const complex& x){  return x.imag ();}inline doublereal (const complex& x){  return x.real ();}inline complexoperator + (const complex& x, const complex& y){  return complex (real (x) + real (y), imag (x) + imag (y));}inline complexoperator + (const complex& x, double y){  return complex (real (x) + y, imag (x));}inline complexoperator + (double x, const complex& y){  return complex (x + real (y), imag (y));}inline complexoperator - (const complex& x, const complex& y){  return complex (real (x) - real (y), imag (x) - imag (y));}inline complexoperator - (const complex& x, double y){  return complex (real (x) - y, imag (x));}inline complexoperator - (double x, const complex& y){  return complex (x - real (y), - imag (y));}inline complexoperator * (const complex& x, const complex& y){  return complex (real (x) * real (y) - imag (x) * imag (y),               real (x) * imag (y) + imag (x) * real (y));}inline complexoperator * (const complex& x, double y){  return complex (real (x) * y, imag (x) * y);}inline complexoperator * (double x, const complex& y){  return complex (x * real (y), x * imag (y));}complexoperator / (const complex& x, double y){  return complex (real (x) / y, imag (x) / y);}inline complexoperator + (const complex& x){  return x;}inline complexoperator - (const complex& x){  return complex (-real (x), -imag (x));}inline booloperator == (const complex& x, const complex& y){  return real (x) == real (y) && imag (x) == imag (y);}inline booloperator == (const complex& x, double y){  return real (x) == y && imag (x) == 0;}inline booloperator == (double x, const complex& y){  return x == real (y) && imag (y) == 0;}inline booloperator != (const complex& x, const complex& y){  return real (x) != real (y) || imag (x) != imag (y);}inline booloperator != (const complex& x, double y){  return real (x) != y || imag (x) != 0;}inline booloperator != (double x, const complex& y){  return x != real (y) || imag (y) != 0;}#include <cmath>inline complexpolar (double r, double t){  return complex (r * cos (t), r * sin (t));}inline complexconj (const complex& x) {  return complex (real (x), -imag (x));}inline doublenorm (const complex& x){  return real (x) * real (x) + imag (x) * imag (x);}#endif   //__MYCOMPLEX__

下面是测试文件

#include <iostream>#include "complex.h"using namespace std;ostream&operator << (ostream& os, const complex& x){  return os << '(' << real (x) << ',' << imag (x) << ')';}int main(){  complex c1(2, 1);  complex c2(4, 0);  cout << c1 << endl;  cout << c2 << endl;  cout << c1+c2 << endl;  cout << c1-c2 << endl;  cout << c1*c2 << endl;  cout << c1 / 2 << endl;  cout << conj(c1) << endl;  cout << norm(c1) << endl;  cout << polar(10,4) << endl;  cout << (c1 += c2) << endl;  cout << (c1 == c2) << endl;  cout << (c1 != c2) << endl;  cout << +c2 << endl;  cout << -c2 << endl;  cout << (c2 - 2) << endl;  cout << (5 + c2) << endl;  return 0;}