fork download
copy 
  1. #include "VectorCpx.h"
  2. #include <cstdlib>
  3. #include <cstdio>
  4. #include <climits>
  5. #include <limits>
  6. #include <iostream>
  7. #include <fstream>
  8.  
  9. using namespace std;
  10.  
  11. /*
  12. 第一段是Complex class，除了Complex與VectorCpx之間乘法的operator函式以外，都是用助教提供的
  13. */
  14. /***********************************constructor***********************************/
  15. Complex::Complex(double _real, double _imag): mReal(_real), mImag(_imag)
  16. {
  17.  
  18. }
  19. /***********************************constructor***********************************/
  20.  
  21. /*************************************功能函式************************************/
  22. const double& Complex::real() const
  23. {
  24.     return mReal;
  25. }
  26.  
  27. const double& Complex::imag() const
  28. {
  29. 	return mImag;
  30. }
  31.  
  32. void Complex::real(double value)
  33. {
  34. 	mReal = value;
  35. }
  36.  
  37. void Complex::imag(double value)
  38. { 
  39. 	mImag = value;
  40. }
  41. /*************************************功能函式************************************/
  42.  
  43. /************************************重載運算子***********************************/
  44. Complex	Complex::operator + (const Complex& cpx)const
  45. {
  46. 	Complex cpOut;
  47. 	cpOut.real(mReal + cpx.real());
  48. 	cpOut.imag(mImag + cpx.imag());
  49. 	return cpOut;
  50. }
  51.  
  52. Complex	Complex::operator - (const Complex& cpx)const
  53. {
  54. 	Complex cpOut;
  55. 	cpOut.real(mReal - cpx.real());
  56. 	cpOut.imag(mImag - cpx.imag());
  57. 	return cpOut;
  58. }
  59.  
  60. Complex	Complex::operator * (const Complex& cpx)const
  61. {
  62. 	Complex cpOut;
  63. 	cpOut.real( mReal * cpx.real() - mImag * cpx.imag() );
  64. 	cpOut.imag( mReal * cpx.imag() + mImag * cpx.real() );
  65. 	return cpOut;
  66. }
  67.  
  68. VectorCpx Complex::operator * (const VectorCpx& vc)const	//Complex*vector的自定義乘法，把每個vector裡的資料都再乘以Complex值。這段要自己寫
  69. {
  70. 	VectorCpx outVec;
  71.  
  72. 	for(unsigned int i = 0; i < vc.size(); i++)
  73. 	{
  74. 		Complex cpOut;
  75. 		cpOut.real( mReal * vc[i].real() - mImag * vc[i].imag() );
  76. 		cpOut.imag( mReal * vc[i].imag() + mImag * vc[i].real() );
  77. 		outVec.push_back(cpOut);
  78. 	}
  79.  
  80. 	return outVec;
  81. }
  82.  
  83. Complex	Complex::operator / (const Complex& cpx)const
  84. {
  85. 	Complex cpOut;
  86. 	double rDiv = cpx.real() * cpx.real() + cpx.imag() * cpx.imag();
  87. 	if( rDiv < DBL_MIN) return cpOut;
  88.     cpOut.real( (mReal * cpx.real() + mImag * cpx.imag()) / rDiv );
  89.     cpOut.imag( (mImag * cpx.real() - mReal * cpx.imag()) / rDiv );
  90. 	return cpOut;
  91. }
  92. /************************************重載運算子***********************************/
  93.  
  94. /************************************自定義輸出***********************************/
  95. ostream& operator<<(ostream& ostr, const Complex& cp)
  96. {
  97.     ostr << "R:" << cp.real() << ", I:" << cp.imag();
  98.     return ostr;
  99. }
  100. /************************************自定義輸出***********************************/
  101.  
  102.  
  103. //////////////////////////////////////////////////////VectorCpx//////////////////////////////////////////////////////////////////////
  104.  
  105.  
  106. /***********************************constructor***********************************/
  107. VectorCpx::VectorCpx():pCpx(NULL),mSize(0),mCapacity(0)	//default constructor，不給預設值的話會造成程式崩潰
  108. {
  109. 	//耶不用寫
  110. }
  111.  
  112. VectorCpx::VectorCpx(unsigned size, const Complex& cpx)	//constructor
  113. {
  114. 	int capacity = 1;
  115. 	while(size > capacity)	//一直乘到capacity比size大，因為capacity一定是2的次方因此乘以2
  116. 		capacity *= 2;
  117.  
  118. 	pCpx = new Complex[capacity];	//最後new出一個該大小的新陣列，把給定的Complex值填進去
  119. 	for(int i = 0; i < size; i++)
  120. 		pCpx[i] = cpx;
  121.  
  122. 	mSize = size;
  123. 	mCapacity = capacity;
  124. }
  125.  
  126. VectorCpx::VectorCpx(unsigned capacity)	//constructor，因為實在很不爽一定要給資料數量和資料內容才能做vector，所以自己寫了一個只要給capatity就可以做出一個空vector
  127. {
  128. 	pCpx = new Complex[capacity];
  129. 	mSize = 0;
  130. 	mCapacity = capacity;
  131. }
  132.  
  133. VectorCpx::VectorCpx(const VectorCpx& vc) //copy constructor
  134. {
  135. 	mSize = vc.mSize;
  136. 	mCapacity = vc.mCapacity;
  137.  
  138. 	pCpx = new Complex[mCapacity];
  139.  
  140. 	for(unsigned int i = 0; i < mSize; i++)
  141. 		pCpx[i] = vc.pCpx[i];
  142. }
  143.  
  144. VectorCpx::~VectorCpx()	//destructor
  145. {
  146. 	delete[] pCpx;
  147. }
  148. /***********************************constructor***********************************/
  149.  
  150. /*************************************功能函式************************************/
  151. VectorCpx& VectorCpx::push_back(const Complex& cpx)	//這個函式是要把一個值存到vector裡
  152. {
  153. 	if(mSize == mCapacity)	//如果原本的vector已經滿了，那就做一個新的vector，大小會是原本vector兩倍
  154. 	{
  155. 		VectorCpx newVec(mCapacity + 1, cpx);	//這裡的mCapacity+1，事實上等於newVec的mSize。為了省事全部填入cpx
  156. 		for(unsigned int i = 0; i < mSize; i++)	//然後除了位於位置mCapacity+1上的cpx，前面全部換回原本pCpx vector裡的所有值
  157. 		{
  158. 			const Complex reg = pCpx[i];
  159. 			newVec.pCpx[i] = reg;
  160. 		}
  161. 		delete[] pCpx;	//別忘記砍掉它，它已經沒有利用價值了
  162.  
  163. 		return newVec;
  164. 	}
  165. 	else
  166. 	{
  167. 		pCpx[mSize] = cpx;	//沒滿的話就直接扔進去就好
  168. 		return *this;
  169. 	}
  170. }
  171.  
  172. VectorCpx& VectorCpx::resize(unsigned place)
  173. {
  174. 	VectorCpx newVec(place, Complex(0.0,0.0));
  175. 	delete[] pCpx;
  176. 	return newVec;
  177. }
  178.  
  179. VectorCpx& VectorCpx::resize(unsigned place, const Complex& cpx)
  180. {
  181. 	/*
  182. 	這個函式會傳入一個整數place，可能還會有一個Complex cpx。
  183. 	如果傳入的整數比vector的大小還大，那就把vector擴張到適合的大小，並且從原本的資料之後一直到place的位置都填入cpx，如果當初沒有傳入cpx則填0
  184. 	如果傳入的整數比vector裡面資料個數還少，那vector大小不變，但原本的資料從尾開始砍到剩下符合的數量
  185. 	*/
  186. 	if(place >= mSize)
  187. 	{
  188. 		VectorCpx newVecBig(place, cpx);
  189.  
  190. 		for(unsigned int i = 0; i < mSize; i++)
  191. 			newVecBig.pCpx[i] = pCpx[i];
  192.  
  193. 		delete[] pCpx;
  194. 		return newVecBig;
  195. 	}
  196. 	else
  197. 	{
  198. 		VectorCpx newVecSmall(mCapacity);
  199.  
  200. 		for(unsigned int i = 0; i < place; i++)
  201. 			newVecSmall.pCpx[i] = pCpx[i];
  202.  
  203. 		delete[] pCpx;
  204. 		newVecSmall.mSize = place;
  205. 		return newVecSmall;
  206. 	}
  207. }
  208.  
  209. unsigned VectorCpx::size() const	//回傳size大小
  210. {
  211. 	return mSize;
  212. }
  213.  
  214. unsigned VectorCpx::capacity() const	//回傳capacity大小
  215. {
  216. 	return mCapacity;
  217. }
  218.  
  219. /*Complex VectorCpx::output(int place) const	//回傳vector裡特定位置的值
  220. {
  221. 	return pCpx[place];
  222. }*/
  223. /*************************************功能函式************************************/
  224.  
  225. /************************************重載運算子***********************************/
  226. Complex VectorCpx::operator*(const VectorCpx& vc) const	//vector*vector的自定義乘法，把每個香對應位置的Complex值乘起來之後全部加起來回傳，有多出來的不用管他
  227. {
  228. 	Complex result;
  229.  
  230. 	if(mSize >= vc.mSize)
  231. 	{
  232. 		for(unsigned int i = 0; i < vc.mSize; i++)
  233. 			result = result + pCpx[i] * vc.pCpx[i];
  234. 	}
  235. 	else
  236. 	{
  237. 		for(unsigned int i = 0; i < mSize; i++)
  238. 			result = result + pCpx[i] * vc.pCpx[i];
  239. 	}
  240.  
  241. 	return result;
  242. }
  243.  
  244. VectorCpx VectorCpx::operator*(const Complex& cpx) const	//vector*Complex的自定義乘法，把每個vector裡的資料都再乘以Complex值
  245. {
  246. 	for(unsigned int i = 0; i < mSize; i++)
  247. 		pCpx[i] = pCpx[i] * cpx;
  248.  
  249. 	return *this;
  250. }
  251.  
  252. Complex& VectorCpx::operator[](const int& place) const	//vector的自定義中括號，vc[i]即可得到位於vector vc中i位置的值
  253. {
  254. 	if(place < 0)	//如果輸入位置小於0，則自動變成0
  255. 		return pCpx[0];
  256. 	else if(place >= mSize)	//如果輸入位置超過資料數目，則扔出最後一筆
  257. 		return pCpx[mSize - 1];
  258. 	else
  259. 		return pCpx[place];
  260. }
  261.  
  262. VectorCpx VectorCpx::operator=(const VectorCpx vc) const	//vector的自定義=符號
  263. {
  264. 	VectorCpx newVec(vc.mSize, Complex(0.0,0.0));
  265.  
  266. 	for(unsigned int i = 0; i < mSize; i++)
  267. 		newVec.pCpx[i] = vc.pCpx[i];
  268.  
  269. 	delete[] pCpx;
  270. 	return newVec;
  271. }
  272. /************************************重載運算子***********************************/
  273.  
  274. /************************************自定義輸出***********************************/
  275. ostream& operator<<(ostream& ostr, const VectorCpx vc)
  276. {
  277. 	unsigned int size = vc.size();
  278.  
  279. 	if(size==0)	//由於size和vc.size()都是unsigned型態，如果直接在下面的for迴圈0-1，會直接跳到最大值去，造成程式無法結束，因此要另外設一個size==0的if判斷式
  280. 	{
  281. 		ostr << endl;
  282. 		return ostr;
  283. 	}
  284.  
  285. 	for(unsigned int i = 0; i < size - 1; i++)
  286. 		ostr << vc[i] << "; ";
  287.  
  288. 	ostr << vc[size - 1] << endl;
  289.  
  290. 	return ostr;
  291. }
  292. /************************************自定義輸出***********************************/
Compilation error #stdin compilation error #stdout 0s 0KB
comments (?)
stdin
copy 
Standard input is empty
compilation info 
prog.cpp:1:23: fatal error: VectorCpx.h: No such file or directory
compilation terminated.
stdout
copy 
Standard output is empty
https://ideone.com/u8voXB
language:
C++ (gcc 8.3)
created:
12 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!