using System;
using Science.Mathematics;
using L=Science.Mathematics.PartialDifferentialEquation;

namespace ScienceTest.MathematicsTest.PartialDifferentialEquationTest
{
	/// <summary>
	/// PoissonEquationTwoDimensionTest
	/// </summary>
	public class PoissonEquation2DTest
	{
		public PoissonEquation2DTest()
		{
		}
		private string result;
		public string Result
		{
			get{return result;}
		}
		public void Compute()
		{
            Function.ToLastType<double, double, double>[] cf
                = new Function.ToLastType<double, double, double>[4];
            cf[0] = new Function.ToLastType<double, double, double>(Line1);
            cf[1] = new Function.ToLastType<double, double, double>(Line2);
            cf[2] = new Function.ToLastType<double, double, double>(Line3);
            cf[3] = new Function.ToLastType<double, double, double>(Line4);
			L.Domain2D xd = new L.Domain2D();
            xd.ConditionFunctionsLessThanZero = cf;
			xd.LowerBoundOfX = 0.0;
			xd.LowerBoundOfY = 0.0;
			xd.UpperBoundOfX = 1.0;
			xd.UpperBoundOfY = 1.0;

			L.PoissonEquation2D eq 
				= new L.PoissonEquation2D();

			L.PoissonEquation2D.Source cs 
				= new L.PoissonEquation2D.Source(S);
		
			eq.SourceTerm = cs;

			eq.MeshSize = 1.0/32.0;
			eq.JacobiRadius = Math.Cos(Math.PI/33.0);

			eq.Domain = xd;
			eq.SetBoundaryValue(-2.0,-2.0,2.0,2.0,10.0);
			eq.Solve();

			for(int i = 0; i < eq.Solution.GetLength(0); i++)
			{
				for(int j = 0; j < eq.Solution.GetLength(1)-1; j++)
					result += Convert.ToString(eq.Solution[i,j])+"\t";
				result += Convert.ToString(eq.Solution[i,eq.Solution.GetLength(1)-1]);
				result +="\r\n";
			}
		}

		private double S(double x, double y)
		{
			if(0.5-0.01 < x && x < 0.5+0.01 && 0.5-0.01 < y && y < 0.5+0.01) return -300.0;
			else return 0.0;
		}
		
		private double Line1(double x, double y)
		{
			return -x;
		}
		private double Line2(double x, double y)
		{
			return -y;
		}
		private double Line3(double x, double y)
		{
			return x-1.0;
		}
		private double Line4(double x, double y)
		{
			return y-1.0;
		}
	}
}
/*
10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10
10	10.0003123016431	10.000624602069	10.0009367398064	10.0012482150093	10.0015580458299	10.0018645870864	10.0021654175514	10.0024572098068	10.0027357111942	10.0029957479204	10.0032314171727	10.0034363359479	10.0036041154863	10.0037288636193	10.0038058091245	10.0038318178675	10.0038058091245	10.0037288636193	10.0036041154863	10.0034363359479	10.0032314171727	10.0029957479204	10.0027357111942	10.0024572098068	10.0021654175514	10.0018645870864	10.0015580458299	10.0012482150093	10.0009367398064	10.000624602069	10.0003123016431	10
10	10.000624602069	10.0012493713722	10.001874138535	10.0024980819392	10.0031193780043	10.0037348891373	10.0043398676432	10.0049277179425	10.0054898785389	10.0060158716827	10.0064935762584	10.0069098200113	10.0072512613942	10.0075055396338	10.0076625532836	10.0077156655747	10.0076625532836	10.0075055396338	10.0072512613942	10.0069098200113	10.0064935762584	10.0060158716827	10.0054898785389	10.0049277179425	10.0043398676432	10.0037348891373	10.0031193780043	10.0024980819392	10.001874138535	10.0012493713722	10.000624602069	10
10	10.0009367398064	10.001874138535	10.0028123666718	10.0037505949813	10.0046864958578	10.0056157271057	10.0065314449035	10.0074239146581	10.0082802181707	10.0090842792055	10.0098171946762	10.0104581057833	10.0109855705026	10.0113794795514	10.0116232035869	10.0117057320674	10.0116232035869	10.0113794795514	10.0109855705026	10.0104581057833	10.0098171946762	10.0090842792055	10.0082802181707	10.0074239146581	10.0065314449035	10.0056157271057	10.0046864958578	10.0037505949813	10.0028123666718	10.001874138535	10.0009367398064	10
10	10.0012482150093	10.0024980819392	10.0037505949813	10.0050054415903	10.0062602872168	10.0075100720984	10.0087462780561	10.0099562713821	10.0111228037585	10.0122238294943	10.0132328194469	10.0141198300624	10.0148534459727	10.0154035980384	10.0157450450983	10.015860854666	10.0157450450983	10.0154035980384	10.0148534459727	10.0141198300624	10.0132328194469	10.0122238294943	10.0111228037585	10.0099562713821	10.0087462780561	10.0075100720984	10.0062602872168	10.0050054415903	10.0037505949813	10.0024980819392	10.0012482150093	10
10	10.0015580458299	10.0031193780043	10.0046864958578	10.0062602872168	10.0078391442397	10.0094179995473	10.0109873231488	10.0125320923697	10.0140309044282	10.0154554142689	10.0167704241711	10.0179349613321	10.0189047847193	10.0196364284036	10.0200925307424	10.0202475932543	10.0200925307424	10.0196364284036	10.0189047847193	10.0179349613321	10.0167704241711	10.0154554142689	10.0140309044282	10.0125320923697	10.0109873231488	10.0094179995473	10.0078391442397	10.0062602872168	10.0046864958578	10.0031193780043	10.0015580458299	10
10	10.0018645870864	10.0037348891373	10.0056157271057	10.0075100720984	10.0094179995473	10.0113354608221	10.0132529261073	10.0151538702062	10.0170133026662	10.0187965068318	10.0204584987847	10.0219448032112	10.0231943186244	10.0241447958285	10.0247410559858	10.0249444606711	10.0247410559858	10.0241447958285	10.0231943186244	10.0219448032112	10.0204584987847	10.0187965068318	10.0170133026662	10.0151538702062	10.0132529261073	10.0113354608221	10.0094179995473	10.0075100720984	10.0056157271057	10.0037348891373	10.0018645870864	10
10	10.0021654175514	10.0043398676432	10.0065314449035	10.0087462780561	10.0109873231488	10.0132529261073	10.0155350434051	10.0178171679578	10.0200719327866	10.0222588088875	10.0243222541193	10.0261914478657	10.0277828788494	10.0290073980565	10.0297824287512	10.0300481426808	10.0297824287512	10.0290073980565	10.0277828788494	10.0261914478657	10.0243222541193	10.0222588088875	10.0200719327866	10.0178171679578	10.0155350434051	10.0132529261073	10.0109873231488	10.0087462780561	10.0065314449035	10.0043398676432	10.0021654175514	10
10	10.0024572098068	10.0049277179425	10.0074239146582	10.0099562713821	10.0125320923697	10.0151538702062	10.0178171679578	10.0205078092534	10.0231984552618	10.0258445371115	10.0283802620985	10.0307158444064	10.032738371687	10.0343194741038	10.035333125156	10.0356832459872	10.035333125156	10.0343194741038	10.032738371687	10.0307158444064	10.0283802620985	10.0258445371115	10.0231984552618	10.0205078092534	10.0178171679578	10.0151538702062	10.0125320923697	10.0099562713821	10.0074239146581	10.0049277179425	10.0024572098068	10
10	10.0027357111942	10.0054898785389	10.0082802181707	10.0111228037585	10.0140309044282	10.0170133026662	10.0200719327866	10.0231984552618	10.0263695379343	10.0295406204461	10.0326383990177	10.0355533023271	10.0381352747282	10.0401990109079	10.0415473405705	10.0420185827891	10.0415473405705	10.0401990109079	10.0381352747282	10.0355533023271	10.0326383990177	10.0295406204461	10.0263695379343	10.0231984552618	10.0200719327866	10.0170133026662	10.0140309044282	10.0111228037585	10.0082802181707	10.0054898785389	10.0027357111942	10
10	10.0029957479204	10.0060158716827	10.0090842792055	10.0122238294943	10.0154554142689	10.0187965068318	10.0222588088875	10.0258445371115	10.0295406204461	10.0333100162635	10.0370794063834	10.0407236925418	10.0440504285807	10.0467939533298	10.0486386441378	10.0492964102848	10.0486386441378	10.0467939533298	10.0440504285807	10.0407236925418	10.0370794063834	10.0333100162635	10.0295406204461	10.0258445371115	10.0222588088875	10.0187965068318	10.0154554142689	10.0122238294943	10.0090842792055	10.0060158716827	10.0029957479204	10
10	10.0032314171727	10.0064935762584	10.0098171946762	10.0132328194469	10.0167704241711	10.0204584987847	10.0243222541193	10.0283802620985	10.0326383990177	10.0370794063834	10.0416455232916	10.0462116290941	10.0505487960098	10.0542877254819	10.0569168785695	10.0578897554964	10.0569168785695	10.0542877254819	10.0505487960098	10.0462116290941	10.0416455232916	10.0370794063834	10.0326383990177	10.0283802620985	10.0243222541193	10.0204584987847	10.0167704241711	10.0132328194469	10.0098171946762	10.0064935762584	10.0032314171727	10
10	10.0034363359479	10.0069098200113	10.0104581057833	10.0141198300624	10.0179349613321	10.0219448032112	10.0261914478657	10.0307158444064	10.0355533023271	10.0407236925418	10.0462116290941	10.0519285111078	10.0576453963394	10.0628912769685	10.0668513668327	10.0684288643261	10.0668513668327	10.0628912769685	10.0576453963394	10.0519285111078	10.0462116290941	10.0407236925418	10.0355533023271	10.0307158444064	10.0261914478657	10.0219448032112	10.0179349613321	10.0141198300624	10.0104581057833	10.0069098200113	10.0034363359479	10
10	10.0036041154863	10.0072512613943	10.0109855705026	10.0148534459727	10.0189047847193	10.0231943186244	10.0277828788494	10.032738371687	10.0381352747282	10.0440504285807	10.0505487960098	10.0576453963394	10.065213002698	10.0727806215775	10.079168454279	10.0821229615625	10.079168454279	10.0727806215775	10.065213002698	10.0576453963394	10.0505487960098	10.0440504285807	10.0381352747282	10.032738371687	10.0277828788494	10.0231943186244	10.0189047847193	10.0148534459727	10.0109855705026	10.0072512613942	10.0036041154863	10
10	10.0037288636193	10.0075055396338	10.0113794795514	10.0154035980384	10.0196364284036	10.0241447958285	10.0290073980565	10.0343194741038	10.0401990109079	10.0467939533298	10.0542877254819	10.0628912769685	10.0727806215775	10.0838497466477	10.0949188746994	10.1017260696734	10.0949188746994	10.0838497466477	10.0727806215775	10.0628912769685	10.0542877254819	10.0467939533298	10.0401990109079	10.0343194741038	10.0290073980565	10.0241447958285	10.0196364284036	10.0154035980384	10.0113794795514	10.0075055396338	10.0037288636193	10
10	10.0038058091245	10.0076625532836	10.0116232035869	10.0157450450983	10.0200925307424	10.0247410559858	10.0297824287512	10.035333125156	10.0415473405705	10.0486386441378	10.0569168785695	10.0668513668327	10.079168454279	10.0949188746994	10.1149312209088	10.1349435690865	10.1149312209088	10.0949188746994	10.079168454279	10.0668513668327	10.0569168785695	10.0486386441378	10.0415473405705	10.035333125156	10.0297824287512	10.0247410559858	10.0200925307424	10.0157450450983	10.0116232035869	10.0076625532836	10.0038058091245	10
10	10.0038318178675	10.0077156655747	10.0117057320674	10.015860854666	10.0202475932543	10.0249444606711	10.0300481426808	10.0356832459872	10.0420185827891	10.0492964102848	10.0578897554964	10.0684288643261	10.0821229615625	10.1017260696734	10.1349435690865	10.208185754971	10.1349435690865	10.1017260696734	10.0821229615625	10.0684288643261	10.0578897554964	10.0492964102848	10.0420185827891	10.0356832459872	10.0300481426808	10.0249444606711	10.0202475932543	10.015860854666	10.0117057320674	10.0077156655747	10.0038318178675	10
10	10.0038058091245	10.0076625532836	10.0116232035869	10.0157450450983	10.0200925307424	10.0247410559858	10.0297824287512	10.035333125156	10.0415473405705	10.0486386441378	10.0569168785695	10.0668513668327	10.079168454279	10.0949188746994	10.1149312209088	10.1349435690865	10.1149312209088	10.0949188746994	10.079168454279	10.0668513668327	10.0569168785695	10.0486386441378	10.0415473405705	10.035333125156	10.0297824287512	10.0247410559858	10.0200925307424	10.0157450450983	10.0116232035869	10.0076625532836	10.0038058091245	10
10	10.0037288636193	10.0075055396338	10.0113794795514	10.0154035980384	10.0196364284036	10.0241447958285	10.0290073980565	10.0343194741038	10.0401990109079	10.0467939533298	10.0542877254819	10.0628912769685	10.0727806215775	10.0838497466477	10.0949188746994	10.1017260696734	10.0949188746994	10.0838497466477	10.0727806215775	10.0628912769685	10.0542877254819	10.0467939533298	10.0401990109079	10.0343194741038	10.0290073980565	10.0241447958285	10.0196364284036	10.0154035980384	10.0113794795514	10.0075055396338	10.0037288636193	10
10	10.0036041154863	10.0072512613943	10.0109855705026	10.0148534459727	10.0189047847193	10.0231943186244	10.0277828788494	10.032738371687	10.0381352747282	10.0440504285807	10.0505487960098	10.0576453963394	10.065213002698	10.0727806215775	10.079168454279	10.0821229615625	10.079168454279	10.0727806215775	10.065213002698	10.0576453963394	10.0505487960098	10.0440504285807	10.0381352747282	10.032738371687	10.0277828788494	10.0231943186244	10.0189047847193	10.0148534459727	10.0109855705026	10.0072512613942	10.0036041154863	10
10	10.0034363359479	10.0069098200113	10.0104581057833	10.0141198300624	10.0179349613321	10.0219448032112	10.0261914478657	10.0307158444064	10.0355533023271	10.0407236925418	10.0462116290941	10.0519285111078	10.0576453963394	10.0628912769685	10.0668513668327	10.0684288643261	10.0668513668327	10.0628912769685	10.0576453963394	10.0519285111078	10.0462116290941	10.0407236925418	10.0355533023271	10.0307158444064	10.0261914478657	10.0219448032112	10.0179349613321	10.0141198300624	10.0104581057833	10.0069098200113	10.0034363359479	10
10	10.0032314171727	10.0064935762584	10.0098171946762	10.0132328194469	10.0167704241711	10.0204584987847	10.0243222541193	10.0283802620985	10.0326383990177	10.0370794063834	10.0416455232916	10.0462116290941	10.0505487960098	10.0542877254819	10.0569168785695	10.0578897554964	10.0569168785695	10.0542877254819	10.0505487960098	10.0462116290941	10.0416455232916	10.0370794063834	10.0326383990177	10.0283802620985	10.0243222541193	10.0204584987847	10.0167704241711	10.0132328194469	10.0098171946762	10.0064935762584	10.0032314171727	10
10	10.0029957479204	10.0060158716827	10.0090842792055	10.0122238294943	10.0154554142689	10.0187965068318	10.0222588088875	10.0258445371115	10.0295406204461	10.0333100162635	10.0370794063834	10.0407236925418	10.0440504285807	10.0467939533298	10.0486386441378	10.0492964102848	10.0486386441378	10.0467939533298	10.0440504285807	10.0407236925418	10.0370794063834	10.0333100162635	10.0295406204461	10.0258445371115	10.0222588088875	10.0187965068318	10.0154554142689	10.0122238294943	10.0090842792055	10.0060158716827	10.0029957479204	10
10	10.0027357111942	10.0054898785389	10.0082802181707	10.0111228037585	10.0140309044282	10.0170133026662	10.0200719327866	10.0231984552618	10.0263695379343	10.0295406204461	10.0326383990177	10.0355533023271	10.0381352747282	10.0401990109079	10.0415473405705	10.0420185827891	10.0415473405705	10.0401990109079	10.0381352747282	10.0355533023271	10.0326383990177	10.0295406204461	10.0263695379343	10.0231984552618	10.0200719327866	10.0170133026662	10.0140309044282	10.0111228037585	10.0082802181707	10.0054898785389	10.0027357111942	10
10	10.0024572098068	10.0049277179425	10.0074239146582	10.0099562713821	10.0125320923697	10.0151538702062	10.0178171679578	10.0205078092534	10.0231984552618	10.0258445371115	10.0283802620985	10.0307158444064	10.032738371687	10.0343194741038	10.035333125156	10.0356832459872	10.035333125156	10.0343194741038	10.032738371687	10.0307158444064	10.0283802620985	10.0258445371115	10.0231984552618	10.0205078092534	10.0178171679578	10.0151538702062	10.0125320923697	10.0099562713821	10.0074239146581	10.0049277179425	10.0024572098068	10
10	10.0021654175514	10.0043398676432	10.0065314449035	10.0087462780561	10.0109873231488	10.0132529261073	10.0155350434051	10.0178171679578	10.0200719327866	10.0222588088875	10.0243222541193	10.0261914478657	10.0277828788494	10.0290073980565	10.0297824287512	10.0300481426808	10.0297824287512	10.0290073980565	10.0277828788494	10.0261914478657	10.0243222541193	10.0222588088875	10.0200719327866	10.0178171679578	10.0155350434051	10.0132529261073	10.0109873231488	10.0087462780561	10.0065314449035	10.0043398676432	10.0021654175514	10
10	10.0018645870864	10.0037348891373	10.0056157271057	10.0075100720984	10.0094179995473	10.0113354608221	10.0132529261073	10.0151538702062	10.0170133026662	10.0187965068318	10.0204584987847	10.0219448032112	10.0231943186244	10.0241447958285	10.0247410559858	10.0249444606711	10.0247410559858	10.0241447958285	10.0231943186244	10.0219448032112	10.0204584987847	10.0187965068318	10.0170133026662	10.0151538702062	10.0132529261073	10.0113354608221	10.0094179995473	10.0075100720984	10.0056157271057	10.0037348891373	10.0018645870864	10
10	10.0015580458299	10.0031193780043	10.0046864958578	10.0062602872168	10.0078391442397	10.0094179995473	10.0109873231488	10.0125320923697	10.0140309044282	10.0154554142689	10.0167704241711	10.0179349613321	10.0189047847193	10.0196364284036	10.0200925307424	10.0202475932543	10.0200925307424	10.0196364284036	10.0189047847193	10.0179349613321	10.0167704241711	10.0154554142689	10.0140309044282	10.0125320923697	10.0109873231488	10.0094179995473	10.0078391442397	10.0062602872168	10.0046864958578	10.0031193780043	10.0015580458299	10
10	10.0012482150093	10.0024980819392	10.0037505949813	10.0050054415903	10.0062602872168	10.0075100720984	10.0087462780561	10.0099562713821	10.0111228037585	10.0122238294943	10.0132328194469	10.0141198300624	10.0148534459727	10.0154035980384	10.0157450450983	10.015860854666	10.0157450450983	10.0154035980384	10.0148534459727	10.0141198300624	10.0132328194469	10.0122238294943	10.0111228037585	10.0099562713821	10.0087462780561	10.0075100720984	10.0062602872168	10.0050054415903	10.0037505949813	10.0024980819392	10.0012482150093	10
10	10.0009367398064	10.001874138535	10.0028123666718	10.0037505949813	10.0046864958578	10.0056157271057	10.0065314449035	10.0074239146581	10.0082802181707	10.0090842792055	10.0098171946762	10.0104581057833	10.0109855705026	10.0113794795514	10.0116232035869	10.0117057320674	10.0116232035869	10.0113794795514	10.0109855705026	10.0104581057833	10.0098171946762	10.0090842792055	10.0082802181707	10.0074239146581	10.0065314449035	10.0056157271057	10.0046864958578	10.0037505949813	10.0028123666718	10.001874138535	10.0009367398064	10
10	10.000624602069	10.0012493713722	10.001874138535	10.0024980819392	10.0031193780043	10.0037348891373	10.0043398676432	10.0049277179425	10.0054898785389	10.0060158716827	10.0064935762584	10.0069098200113	10.0072512613942	10.0075055396338	10.0076625532836	10.0077156655747	10.0076625532836	10.0075055396338	10.0072512613942	10.0069098200113	10.0064935762584	10.0060158716827	10.0054898785389	10.0049277179425	10.0043398676432	10.0037348891373	10.0031193780043	10.0024980819392	10.001874138535	10.0012493713722	10.000624602069	10
10	10.0003123016431	10.000624602069	10.0009367398064	10.0012482150093	10.0015580458299	10.0018645870864	10.0021654175514	10.0024572098068	10.0027357111942	10.0029957479204	10.0032314171727	10.0034363359479	10.0036041154863	10.0037288636193	10.0038058091245	10.0038318178675	10.0038058091245	10.0037288636193	10.0036041154863	10.0034363359479	10.0032314171727	10.0029957479204	10.0027357111942	10.0024572098068	10.0021654175514	10.0018645870864	10.0015580458299	10.0012482150093	10.0009367398064	10.000624602069	10.0003123016431	10
10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10	10
*/