Differential D87 Diff 389 src/EvtGenBase/EvtComplex.cpp

Changeset View

Standalone View

	Show All 25 Lines
	#include <math.h>			#include <math.h>
	using std::ostream;			using std::ostream;

	ostream& operator<<( ostream& s, const EvtComplex& c )			ostream& operator<<( ostream& s, const EvtComplex& c )
	{			{
	s << "(" << c._rpart << "," << c._ipart << ")";			s << "(" << c._rpart << "," << c._ipart << ")";
	return s;			return s;
	}			}

	EvtComplex& EvtComplex::operator*=( EvtComplex c )
	{
	double r = _rpart * c._rpart - _ipart * c._ipart;
	double i = _rpart * c._ipart + _ipart * c._rpart;

	_rpart = r;
	_ipart = i;

	return *this;
	}

	EvtComplex& EvtComplex::operator/=( EvtComplex c )
	{
	double inv = 1.0 / ( c._rpart * c._rpart + c._ipart * c._ipart );

	double r = inv * ( _rpart * c._rpart + _ipart * c._ipart );
	double i = inv * ( _ipart * c._rpart - _rpart * c._ipart );

	_rpart = r;
	_ipart = i;

	return *this;
	}