#ifndef DEFINITION_H
#define DEFINITION_H

/****************范圍設置******************/
//總迭代范圍
const int DOMAINX=60;  
const int DOMAINY=60;

//點源位置(中心)
const int IC=DOMAINX/2;
const int JC=DOMAINY/2;

/*************電磁參數**************/
#define PI 3.1415926535897932384626433832795

const double V=3.0*pow(10.0,8.0);         //波速
const double MIU0=4.0*PI*pow(10.0,-7.0);  //真空磁導率
const double EPSILON0=1.0/(4.0*PI*9.0*pow(10.0,9.0));  //真空電導率
const double Z0=sqrt(MIU0/EPSILON0);      //真空波阻抗

const double CFL=0.5;          //*迭代常數(c*deltat/deltax)
const double Coefficient1=CFL;  //總系數1

/******************源*******************/
const int	 T=1500;						    //迭代時間
const double FREQUENCY=3.0*pow(10.0,9.0);       //頻率
const double WAVELENGTH=V/FREQUENCY;            //波長
const double DELTAX=0.005;  //0.05;				//空間步
const double DELTAT=DELTAX*CFL/V;	//時間步
const double TIME_period=1/(FREQUENCY*DELTAT);

//脈沖源設置
const double T0=25;  //2.0*pow(10.0,-9.0); //*中心
const double SPREAD=8.0;				   //*寬度

//PML厚度
const int NPML=10;

//PML測試
double maxA=0.;
double maxB=0.;
double EzA[1000];        //測試向量A
double EzB[1000];        //測試向量B

//場外部變量(x,y,z表示極化);
static double Ez[DOMAINX][DOMAINY];
static double Hx[DOMAINX][DOMAINY];
static double Hy[DOMAINX][DOMAINY];
#endif