murabc.h

基于电磁波方程

void TSMurEz_x(int i,int j1,int j2)
{
	const float d1=-0.33333333;
	const float d2=0.16666667;
	float rotEz;
	Ez[i][j1]=Ez[i][j2]-d1*Ez[i][j1];
	rotEz=Hy[i][j1]-Hy[i-1][j1]+Hy[i][j2]-Hy[i-1][j2];
	Ez[i][j1]=Ez[i][j1]+d2*rotEz;
	rotEz=Hy[i][j2]-Hy[i-1][j2]+Hx[i][j2-1]-Hx[i][j2];
	Ez[i][j1]=Ez[i][j1]+d1*(Ez[i][j2]+0.5*rotEz);
}
void TSMurEz_y(int j,int i1,int i2)
{
	const float d1=-0.33333333;
	const float d2=0.16666667;
	float rotEz;
	Ez[i1][j]=Ez[i2][j]-d1*Ez[i1][j];
	rotEz=Hx[i1][j]-Hx[i1][j-1]+Hx[i2][j]-Hx[i2][j-1];
	Ez[i1][j]=Ez[i1][j]-d2*rotEz;
	rotEz=Hy[i2][j]-Hy[i2-1][j]+Hx[i2][j-1]-Hx[i2][j];
	Ez[i1][j]=Ez[i1][j]+d1*(Ez[i2][j]+0.5*rotEz);
}
void TSMurEz_Point(int L,int i,int j,int f1,int f2)
{
	const float d3=0.2928932;
    const float d4=0.7071068;
	Ez[i][j]=EAC[L][0];
	EAC[L][0]=EAC[L][1];
	EAC[L][1]=d3*Ez[i][j]+d4*Ez[i+f1][j+f2];
}