//PolyValueTwoDim.cpp		二維實(復)多項式求值示例

#include <iostream>			//輸入輸出流頭文件
#include <polynomials.h>	//多項式及連分式頭文件
using namespace std;		//名字空間

//求多項式p(x,y)在多點處的函數值
void main(void)
{	
	cout << "PolyValueTwoDim()" << endl << endl;

	//多項式系數與自變量都為實數
	cout << "Coefficient and variable are real. Result is real: " << endl;
	double dCoff[4][5] = 
	{
		{ 1.0,  2.0,  3.0,  4.0,  5.0},
		{ 6.0,  7.0,  8.0,  9.0, 10.0},
		{11.0, 12.0, 13.0, 14.0, 15.0},
		{16.0, 17.0, 18.0, 19.0, 20.0}
	};
	
	matrix<double> mdCoff(&dCoff[0][0],4,5);

	double dX = 0.6081, dY = -1.3054;
	
	cout.setf(ios::fixed|ios::left);	//輸出數據為定點法，靠左對齊
	cout.precision(4);					//精度4位
	
	for(size_t st = 0; st < 6; st++)	//輸出多項式值
	{
		double dValue = PolyValueTwoDim(mdCoff, dX, dY);
		cout << "p(" << dX << ", " << dY << ") = " << dValue << endl;
		dX = dX + 0.1101;
		dY = dY - 0.1101;
	}
	cout << endl;

	//多項式系數與自變量都為復數
	cout << "Coefficient and variable are complex. Result is complex: " << endl;
	complex<double> cdCoff[5][3] = 
	{
		{complex<double>(1.2, 3.4),complex<double>(-3.2, -4.8), complex<double>(2.7, -6.4)},
		{complex<double>(3.1, -0.9),complex<double>(-3.5, 1.9),complex<double>(2.1, -5.4)},
		{complex<double>(-2.5, 1.8),complex<double>(-2.4, 2.4),complex<double>(-1.5, 0.4)},
		{complex<double>(-2.5, 2.1),complex<double>(-2.5, 8.4),complex<double>(2.3, -1.4)},
		{complex<double>(1.1, -0.2),complex<double>(-0.4, 3.4),complex<double>(0.5, -2.4)},
	};
	
	matrix<complex<double> > mcdCoff(&cdCoff[0][0],5,3);

	complex<double> cdX = complex<double>( 0.11, -1.99);
	complex<double> cdY = complex<double>(-0.11,  1.99);
	
	for(st = 0; st < 6; st++)	//輸出多項式值
	{
		complex<double> cp = PolyValueTwoDim(mcdCoff, cdX, cdY);
		cout << "cp(cdX, cdY) = cp(" << cdX << ", " << cdY << ") = " << cp << endl;
		cdX = cdX + complex<double>(0.11, 0.11);
		cdY = cdY - complex<double>(0.11, 0.11);
	}	
	cout << endl;

	//多項式系數為實數，自變量為復數
	cout << "Coefficient is real. Variable is complex. Result is complex: " << endl;
	double dcCoff[5][3] =
	{
		{ 1.0,  2.0,  3.0},
		{ 4.0,  5.0,  6.0},
		{ 7.0,  8.0,  9.0},
		{10.0, 11.0, 12.0},
		{13.0, 14.0, 15.0} 
	};

	matrix<double> mdcCoff(&dcCoff[0][0], 5,3);

	cdX = complex<double>( 0.11, -1.99);
	cdY = complex<double>(-0.11,  1.99);
		
	for(st = 0; st < 6; st++)	//輸出多項式值
	{
		complex<double> dp = PolyValueTwoDim(mdcCoff, cdX, cdY);
		cout << "dp(cdX, cdY) = dp(" << cdX << ", " << cdY << ") = " << dp << endl;
		cdX = cdX + complex<double>(0.11, 0.11);
		cdY = cdY - complex<double>(0.11, 0.11);
	}
}