#ifndef _COMPLEX_H
#define _COMPLEX_H

#include <iostream>
using namespace std;

class Complex {
  private:
	  float re_;
	  float im_;
  public:
	  Complex(float r=0, float i=0) { re_=r; im_=i;}
	  int operator==( const Complex& c );
	  bool operator!=( const Complex& c );
	  Complex operator+( const Complex& c );
	  Complex operator-( const Complex& c );
	  Complex operator/( const Complex& c );
	  Complex operator*( const Complex& c );
	  float real( void ) { return( re_ ); }
	  float imag( void ) { return( im_ ); }
};

//std::ostream& operator<< (std::ostream& s, Complex& z);
ostream& operator<<(ostream& s, Complex& z);


inline int Complex::operator==( const Complex& c )
{
  return( re_==c.re_ && im_==c.im_ );
}

inline bool Complex::operator!=( const Complex& c )
{
  return !operator==(c);
}

inline Complex Complex::operator+( const Complex& c )
{
  return( Complex (re_+c.re_, im_+c.im_) );  // Cria Objeto temporario
}

inline Complex Complex::operator-( const Complex& c )
{
  return( Complex (re_-c.re_, im_-c.im_) );  // Cria Objeto temporario
}

inline Complex Complex::operator*( const Complex& c )
{
  return( Complex (re_*c.re_-im_*c.im_, re_*c.im_ + im_*c.re_) );  // Cria Objeto temporario
}

inline Complex Complex::operator/( const Complex& c )
{
  if( c.re_ > c.im_ )
  {
	 float d1 = c.im_/c.re_,      // (a+bi)/(c-di) = ((a+bd')+i(b-ad'))/(c+dd')
								  //                 onde d' = d/c
	 div = c.re_ + c.im_*d1;

	 return( Complex ((re_+im_*d1)/div, (im_ - re_*d1)/div) );
  }
  else
  {
	 float c1 = c.re_/c.im_,  // (a+bi)/(c-di) = ((a+bc')+i(b-ac'))/(cc'+d)
								  //                 onde c' = c/d
 	 div = c.re_*c1 + c.im_;

	 return( Complex ((re_*c1+im_)/div, (im_*c1 - re_)/div) );
  }
}

#endif