mycomplex.cpp

《陈必红算法》一书的附带的源代码

// MyComplex.cpp: implementation of the CMyComplex class.
//
//////////////////////////////////////////////////////////////////////

#include "stdafx.h"
#include "MyComplex.h"

#ifdef _DEBUG
#undef THIS_FILE
static char THIS_FILE[]=__FILE__;
#define new DEBUG_NEW
#endif

//////////////////////////////////////////////////////////////////////
// Construction/Destruction
//////////////////////////////////////////////////////////////////////

CMyComplex::CMyComplex()
{
m_real=0.0;
m_image=0.0;
}
CMyComplex::CMyComplex(CMyComplex &value)
{
m_real=value.m_real;
m_image=value.m_image;
}
CMyComplex::CMyComplex(double myreal,double myimage)
{
m_real=myreal;
m_image=myimage;
}
CMyComplex::CMyComplex(double myreal)
{
m_real=myreal;
m_image=0.0;
}
CMyComplex::~CMyComplex()
{
m_real=0.0;
m_image=0.0;
}
CMyComplex& CMyComplex::operator=(double a)
{
m_real=a;
m_image=0.0;
return (*this);
}
CMyComplex  CMyComplex::operator*(double a)
{
complex<double>temp(m_real,m_image);
CMyComplex result;
temp=temp*a;
result.m_real=temp.real();
result.m_image=temp.imag();
return result;
}
CMyComplex  CMyComplex::operator*(CMyComplex & a)
{
complex<double>temp(m_real,m_image);
CMyComplex result;
complex<double>temp2(a.m_real,a.m_image);
temp=temp*temp2;
result.m_real=temp.real();
result.m_image=temp.imag();
return result;
}
CMyComplex& CMyComplex::operator*=(double a)
{
complex<double>temp(m_real,m_image);
temp*=a;
m_real=temp.real();
m_image=temp.imag();
return (*this);
}
CMyComplex& CMyComplex::operator*=(CMyComplex &a)
{
complex<double>temp(m_real,m_image);
complex<double>temp2(a.m_real,a.m_image);
temp=temp*temp2;
m_real=temp.real();
m_image=temp.imag();
return (*this);
}
CMyComplex CMyComplex::operator+(double a)
{
CMyComplex result;
complex<double>temp(m_real,m_image);
temp=temp+a;
result.m_real=temp.real();
result.m_image=temp.imag();
return result;
}
CMyComplex CMyComplex::operator+(CMyComplex &a)
{
complex<double>temp(m_real,m_image);
CMyComplex result;
complex<double>temp2(a.m_real,a.m_image);
temp=temp+temp2;
result.m_real=temp.real();
result.m_image=temp.imag();
return result;
}
CMyComplex CMyComplex::operator+=(CMyComplex &a)
{
complex<double>temp(m_real,m_image);
complex<double>temp2(a.m_real,a.m_image);
temp=temp+temp2;
m_real=temp.real();
m_image=temp.imag();
return (*this);
}
CMyComplex &CMyComplex::operator+=(double a)
{
complex<double>temp(m_real,m_image);
complex<double>temp2;
temp2=a;
temp=temp+temp2;
m_real=temp.real();
m_image=temp.imag();
return (*this);
}
CMyComplex CMyComplex::operator-(double a)
{
CMyComplex result;
complex<double>temp(m_real,m_image);
temp=temp-a;
result.m_real=temp.real();
result.m_image=temp.imag();
return result;
}
CMyComplex &CMyComplex::operator-=(double a)
{
complex<double>temp(m_real,m_image);
complex<double>temp2;
temp2=a;
temp=temp-temp2;
m_real=temp.real();
m_image=temp.imag();
return (*this);
}
CMyComplex CMyComplex::operator-(CMyComplex &a)
{
complex<double>temp(m_real,m_image);
CMyComplex result;
complex<double>temp2(a.m_real,a.m_image);
temp=temp-temp2;
result.m_real=temp.real();
result.m_image=temp.imag();
return result;
}
CMyComplex &CMyComplex::operator-=(CMyComplex & a)
{
complex<double>temp(m_real,m_image);
complex<double>temp2(a.m_real,a.m_image);
temp=temp-temp2;
m_real=temp.real();
m_image=temp.imag();
return (*this);
}
CMyComplex CMyComplex::operator/(double a)
{
CMyComplex result;
complex<double>temp(m_real,m_image);
temp=temp/a;
result.m_real=temp.real();
result.m_image=temp.imag();
return result;
}
CMyComplex CMyComplex ::operator/(CMyComplex& a)
{
complex<double>temp(m_real,m_image);
CMyComplex result;
complex<double>temp2(a.m_real,a.m_image);
temp=temp/temp2;
result.m_real=temp.real();
result.m_image=temp.imag();
return result;
}
CMyComplex& CMyComplex::operator/=(double a)
{
complex<double>temp(m_real,m_image);
complex<double>temp2;
temp2=a;
temp=temp/temp2;
m_real=temp.real();
m_image=temp.imag();
return (*this);
}
CMyComplex& CMyComplex::operator/=(CMyComplex &a)
{
complex<double>temp(m_real,m_image);
complex<double>temp2(a.m_real,a.m_image);
temp=temp/temp2;
m_real=temp.real();
m_image=temp.imag();
return (*this);
}
CMyComplex CMyComplex::operator-()
{
CMyComplex result;
result.m_real=m_real*-1;
result.m_image=m_image*-1;
return result;
}

CMyComplex  operator/(double in,CMyComplex value)// 全局函数，负责倒数的计算1/CMyComplex
{
CMyComplex result;
complex<double> temp(value.m_real,value.m_image);
temp= in/temp;
result.m_real=temp.real();
result.m_image=temp.imag();
return result;
}

CMyComplex CMyComplex::pow(double b)
{
 CMyComplex result;
 complex<double> temp(m_real,m_image);
 complex<double> temp2;
 temp2=std::pow(temp,b);
 double real=temp2.real();
 double image=temp2.imag();
 if(fabs(real)<=1e-10) real=0;
 if (fabs(image)<=1e-10) image=0;
 result.m_real=real;
 result.m_image=image;
return result;
}
CMyComplex CMyComplex::pow(CMyComplex b)
{
complex<double> a(m_real,m_image);
complex<double> c(b.m_real,b.m_image);
complex<double> rt;
rt=std::pow(a,c);
CMyComplex result;
double real=rt.real();
double image=rt.imag();
 if(fabs(real)<=1e-10) real=0;
 if (fabs(image)<=1e-10) image=0;
 result.m_real=real;
 result.m_image=image;
 return result;
}
CMyComplex CMyComplex::conj()
{
CMyComplex result;
result.m_real=m_real;
result.m_image=-m_image;
return result;
}
std::ostream & operator<<(std::ostream& o, CMyComplex& b)
{
    o<<'('<<b.m_real<<','<<b.m_image<<')';
	return o;
}