#ifndef FIELD_H
#define FIELD_H
#include "Def.h"

//起始中止點
int Beginx=0;
int Endx=0;
int Beginy=0;
int Endy=0;

//初始化場
void InitialField()
{
	for (int i=0; i<DOMAINX; i++)
	{
		for (int j=0; j<DOMAINY; j++)
		{				
				Ez[i][j]=0.0;
				Hx[i][j]=0.0;
				Hy[i][j]=0.0;		
		}
	}
}

//Hx,Hy推導場Ez
void  EzField(const double Coefficient1,const double Coefficient2)
{
	Beginx=NPML+1;
	Endx=DOMAINX-NPML-1;
	Beginy=NPML+1;
	Endy=DOMAINY-NPML-1;

	for (int i=Beginx; i<Endx; i++)
	{
		for (int j=Beginy; j<Endy; j++)
		{						
			 Ez[i][j]=Ez[i][j]+(
							Coefficient1*((Hy[i][j]-Hy[i-1][j])-\
							              (Hx[i][j]-Hx[i][j-1]))+
							Coefficient2*((Hy[i+1][j]-Hy[i-2][j])-\
							              (Hx[i][j+1]-Hx[i][j-2])));						
		}
	}
}

//Ez推導場Hx
void  HxField(const double Coefficient1,const double Coefficient2)
{
	Beginx=NPML+1;
	Endx=DOMAINX-NPML-1;
	Beginy=NPML;
	Endy=DOMAINY-NPML-1;

	for (int i=Beginx; i<Endx; i++)
	{
		for (int j=Beginy; j<Endy; j++)
		{				
			Hx[i][j]=Hx[i][j]+(
							Coefficient1*(Ez[i][j]-Ez[i][j+1])+\
							Coefficient2*(Ez[i][j-1]-Ez[i][j+2]));									 
		}
	}
}

//Ez推導場Hy
void  HyField(const double Coefficient1,const double Coefficient2)
{
    Beginx=NPML;
	Endx=DOMAINX-NPML-1;
	Beginy=NPML+1;
	Endy=DOMAINY-NPML-1;

	for (int i=Beginx; i<Endx; i++)
	{
		for (int j=Beginy; j<Endy; j++)
		{
			Hy[i][j]=Hy[i][j]+(
						Coefficient1*(Ez[i+1][j]-Ez[i][j])+\
						Coefficient2*(Ez[i+2][j]-Ez[i-1][j]));
						
		}
	}
}
	
#endif
/*
#ifndef FIELD_H
#define FIELD_H
#include "Definition.h"

//起始中止點
int Beginx=0;
int Endx=0;
int Beginy=0;
int Endy=0;

//初始化場
void InitialField()
{
	for (int i=0; i<DOMAINX; i++)
	{
		for (int j=0; j<DOMAINY; j++)
		{				
				Ez[i][j]=0.0;
				Hx[i][j]=0.0;
				Hy[i][j]=0.0;		
		}
	}
}

//Hx,Hy推導場Ez
void  EzField(const double Coefficient1,const double Coefficient2)
{
	Beginx=NPML+1;
	Endx=DOMAINX-NPML-1;
	Beginy=NPML+1;
	Endy=DOMAINY-NPML-1;

	for (int i=Beginx; i<Endx; i++)
	{
		for (int j=Beginy; j<Endy; j++)
		{						
			 Ez[i][j]=Ez[i][j]+(
							Coefficient1*((Hy[i][j]-Hy[i-1][j])-\
							              (Hx[i][j]-Hx[i][j-1]))+
							Coefficient2*((Hy[i+1][j]-Hy[i-2][j])-\
							              (Hx[i][j+1]-Hx[i][j-2])));						
		}
	}
}

//Ez推導場Hx
void  HxField(const double Coefficient1,const double Coefficient2)
{
	Beginx=NPML+1;
	Endx=DOMAINX-NPML-1;
	Beginy=NPML;
	Endy=DOMAINY-NPML-1;

	for (int i=Beginx; i<Endx; i++)
	{
		for (int j=Beginy; j<Endy; j++)
		{				
			Hx[i][j]=Hx[i][j]+(
							Coefficient1*(Ez[i][j]-Ez[i][j+1])+\
							Coefficient2*(Ez[i][j-1]-Ez[i][j+2]));									 
		}
	}
}

//Ez推導場Hy
void  HyField(const double Coefficient1,const double Coefficient2)
{
    Beginx=NPML;
	Endx=DOMAINX-NPML-1;
	Beginy=NPML+1;
	Endy=DOMAINY-NPML-1;

	for (int i=Beginx; i<Endx; i++)
	{
		for (int j=Beginy; j<Endy; j++)
		{
			Hy[i][j]=Hy[i][j]+(
						Coefficient1*(Ez[i+1][j]-Ez[i][j])+\
						Coefficient2*(Ez[i+2][j]-Ez[i-1][j]));
						
		}
	}
}
	
#endif
	  */