calculation.c

利用不均一分割的FDTD解析程序

				bbz = hy/delta_x-hx/delta_y;
				ezi[i][j][k] = ezi[i][j][k]+rh2*bbz;
			}
		}
	}


	/*****  Absorbing Boundary Condition (Mur's 2nd Approx.)  *****/

	for(i=1;i<=(xs-2);i++){                           /*EXI*/
		uy41 = uy4;
		if((i >= (Fx-x_sub_length))
		    &&(i <= (Fx+x_sub_length))){
			uy41 = uy4_sub;
		}

		for(k=1;k<=(zs-1);k++){
			uy31 = uy3;
			if(k <= (Fz+z_sub_length)){
				uy31 = uy3_sub;
			}

			r1 = exi[i][1][k]+exiya2[i][0][k];
			r2 = exiya1[i][0][k]+exiya1[i][1][k];
			r3 = exiya1[i][0][k+1]-2*exiya1[i][0][k]
			    +exiya1[i][0][k-1]+exiya1[i][1][k+1]
			    -2*exiya1[i][1][k]+exiya1[i][1][k-1];
			r4 = exiya1[i+1][0][k]-2*exiya1[i][0][k]
			    +exiya1[i-1][0][k]+exiya1[i+1][1][k]
			    -2*exiya1[i][1][k]+exiya1[i-1][1][k];

			exi[i][0][k] = -exiya2[i][1][k]+uy1*r1
			    +uy2*r2+uy31*r3+uy41*r4;

			r1 = exi[i][ys-1][k]+exiyb2[i][ys][k];
			r2 = exiyb1[i][ys][k]+exiyb1[i][ys-1][k];
			r3 = exiyb1[i][ys][k+1]-2*exiyb1[i][ys][k]
			    +exiyb1[i][ys][k-1]+exiyb1[i][ys-1][k+1]
			    -2*exiyb1[i][ys-1][k]+exiyb1[i][ys-1][k-1];
			r4 = exiyb1[i+1][ys][k]-2*exiyb1[i][ys][k]
			    +exiyb1[i-1][ys][k]+exiyb1[i+1][ys-1][k]
			    -2*exiyb1[i][ys-1][k]+exiyb1[i-1][ys-1][k];

			exi[i][ys][k] = -exiyb2[i][ys-1][k]+uy1*r1
			    +uy2*r2+uy31*r3+uy41*r4;
		}
	}

	for(k=1;k<=(zs-1);k++){
		r1 = exi[0][1][k]-exiya1[0][0][k];
		exi[0][0][k] = exiya1[0][1][k]+uy1*r1;
		r1 = exi[xs-1][1][k]-exiya1[xs-1][0][k];
		exi[xs-1][0][k] = exiya1[xs-1][1][k]+uy1*r1;
		r1 = exi[0][ys-1][k]-exiyb1[0][ys][k];
		exi[0][ys][k] = exiyb1[0][ys-1][k]+uy1*r1;
		r1 = exi[xs-1][ys-1][k]-exiyb1[xs-1][ys][k];
		exi[xs-1][ys][k] = exiyb1[xs-1][ys-1][k]+uy1*r1;
	}

	for(i=1;i<=(xs-2);i++){
		uz31 = uz3;
		if((i >= (Fx-x_sub_length))
		    && (i <= (Fx+x_sub_length))){
			uz31 = uz3_sub;
		}

		for(j=1;j<=(ys-1);j++){
			uz41 = uz4;
			if((j >= (Fy-y_sub_length))
			    &&(j <= (Fy+y_sub_length))){
				uz41 = uz4_sub;
			}
			         
			r1 = exi[i][j][zs-1]+exizb2[i][j][zs];
			r2 = exizb1[i][j][zs]+exizb1[i][j][zs-1];
			r3 = exizb1[i+1][j][zs]-2*exizb1[i][j][zs]
			    +exizb1[i-1][j][zs]+exizb1[i+1][j][zs-1]
			    -2*exizb1[i][j][zs-1]+exizb1[i-1][j][zs-1];
			r4 = exizb1[i][j+1][zs]-2*exizb1[i][j][zs]
			    +exizb1[i][j-1][zs]+exizb1[i][j+1][zs-1]
			    -2*exizb1[i][j][zs-1]+exizb1[i][j-1][zs-1];

			exi[i][j][zs] = -exizb2[i][j][zs-1]+uz1*r1
			    +uz2*r2+uz31*r3+uz41*r4;
		}
	}

	for(j=1;j<=(ys-1);j++){
		r1 = exi[0][j][1]-exiza1[0][j][0];
		exi[0][j][0] = exiza1[0][j][1]+uz1*r1;
		r1 = exi[xs-1][j][1]-exiza1[xs-1][j][0];
		exi[xs-1][j][0] = exiza1[xs-1][j][1]+uz1*r1;
		r1 = exi[0][j][zs-1]-exizb1[0][j][zs];
		exi[0][j][zs] = exizb1[0][j][zs-1]+uz1*r1;
		r1 = exi[xs-1][j][zs-1]-exizb1[xs-1][j][zs];
		exi[xs-1][j][zs] = exizb1[xs-1][j][zs-1]+uz1*r1;
	}

	for(i=0;i<=(xs-1);i++){
		r1 = exi[i][0][1]-exiza1[i][0][0];
		exi[i][0][0] = exiza1[i][0][1]+uz1*r1;
		r1 = exi[i][ys][1]-exiza1[i][ys][0];
		exi[i][ys][0] = exiza1[i][ys][1]+uz1*r1;
		r1 = exi[i][0][zs-1]-exizb1[i][0][zs];
		exi[i][0][zs] = exizb1[i][0][zs-1]+uz1*r1;
		r1 = exi[i][ys][zs-1]-exizb1[i][ys][zs];
		exi[i][ys][zs] = exizb1[i][ys][zs-1]+uz1*r1;
	}

	for(j=1;j<=(ys-2);j++){                      /*EYI*/
		uz41 = uz4;
		if((j >= (Fy-y_sub_length))
		    &&(j <= (Fy+y_sub_length))){
			uz41 = uz4_sub;
		}

		for(i=1;i<=(xs-1);i++){
			uz31 = uz3;
			if((i >= (Fx-x_sub_length))
			    && (i <= (Fx+x_sub_length))){
				uz31 = uz3_sub;
			}

			r1 = eyi[i][j][zs-1]+eyizb2[i][j][zs];
			r2 = eyizb1[i][j][zs]+eyizb1[i][j][zs-1];
			r3 = eyizb1[i+1][j][zs]-2*eyizb1[i][j][zs]
			    +eyizb1[i-1][j][zs]+eyizb1[i+1][j][zs-1]
			    -2*eyizb1[i][j][zs-1]+eyizb1[i-1][j][zs-1];
			r4 = eyizb1[i][j+1][zs]-2*eyizb1[i][j][zs]
			    +eyizb1[i][j-1][zs]+eyizb1[i][j+1][zs-1]
			    -2*eyizb1[i][j][zs-1]+eyizb1[i][j-1][zs-1];
			eyi[i][j][zs] = -eyizb2[i][j][zs-1]+uz1*r1
			    +uz2*r2+uz31*r3+uz41*r4;
		}
	}

	for(i=1;i<=(xs-1);i++){
		r1 = eyi[i][0][1]-eyiza1[i][0][0];
		eyi[i][0][0] = eyiza1[i][0][1]+uz1*r1;
		r1 = eyi[i][ys-1][1]-eyiza1[i][ys-1][0];
		eyi[i][ys-1][0] = eyiza1[i][ys-1][1]+uz1*r1;
		r1 = eyi[i][0][zs-1]-eyizb1[i][0][zs];
		eyi[i][0][zs] = eyizb1[i][0][zs-1]+uz1*r1;
		r1 = eyi[i][ys-1][zs-1]-eyizb1[i][ys-1][zs];
		eyi[i][ys-1][zs] = eyizb1[i][ys-1][zs-1]+uz1*r1;
	}

	for(j=1;j<=(ys-2);j++){
		ux31 = ux3;
		if((j >= (Fy-y_sub_length))
		    &&(j <= (Fy+y_sub_length))){
			ux31 = ux3_sub;
		}

		for(k=1;k<=(zs-1);k++){
			ux41 = ux4;
			if(k <= (Fz+z_sub_length)){
				ux41 = ux4_sub;
			}

			r1 = eyi[1][j][k]+eyixa2[0][j][k];
			r2 = eyixa1[0][j][k]+eyixa1[1][j][k];
			r3 = eyixa1[0][j+1][k]-2*eyixa1[0][j][k]
			    +eyixa1[0][j-1][k]+eyixa1[1][j+1][k]
			    -2*eyixa1[1][j][k]+eyixa1[1][j-1][k];
			r4 = eyixa1[0][j][k+1]-2*eyixa1[0][j][k]
			    +eyixa1[0][j][k-1]+eyixa1[1][j][k+1]
			    -2*eyixa1[1][j][k]+eyixa1[1][j][k-1];

			eyi[0][j][k] = -eyixa2[1][j][k]+ux1*r1
			    +ux2*r2+ux31*r3+ux41*r4;

			r1 = eyi[xs-1][j][k]+eyixb2[xs][j][k];
			r2 = eyixb1[xs][j][k]+eyixb1[xs-1][j][k];
			r3 = eyixb1[xs][j+1][k]-2*eyixb1[xs][j][k]
			    +eyixb1[xs][j-1][k]+eyixb1[xs-1][j+1][k]
			    -2*eyixb1[xs-1][j][k]+eyixb1[xs-1][j-1][k];
			r4 = eyixb1[xs][j][k+1]-2*eyixb1[xs][j][k]
			    +eyixb1[xs][j][k-1]+eyixb1[xs-1][j][k+1]
			    -2*eyixb1[xs-1][j][k]+eyixb1[xs-1][j][k-1];

			eyi[xs][j][k] = -eyixb2[xs-1][j][k]+ux1*r1
			    +ux2*r2+ux31*r3+ux41*r4;
		}
	}

	for(k=1;k<=(zs-1);k++){
		r1 = eyi[1][0][k]-eyixa1[0][0][k];
		eyi[0][0][k] = eyixa1[1][0][k]+ux1*r1;
		r1 = eyi[1][ys-1][k]-eyixa1[0][ys-1][k];
		eyi[0][ys-1][k] = eyixa1[1][ys-1][k]+ux1*r1;
		r1 = eyi[xs-1][0][k]-eyixb1[xs][0][k];
		eyi[xs][0][k] = eyixb1[xs-1][0][k]+ux1*r1;
		r1 = eyi[xs-1][ys-1][k]-eyixb1[xs][ys-1][k];
		eyi[xs][ys-1][k] = eyixb1[xs-1][ys-1][k]+ux1*r1;
	}

	for(j=0;j<=(ys-1);j++){
		r1 = eyi[0][j][1]-eyiza1[0][j][0];
		eyi[0][j][0] = eyiza1[0][j][1]+uz1*r1;
		r1 = eyi[xs][j][1]-eyiza1[xs][j][0];
		eyi[xs][j][0] = eyiza1[xs][j][1]+uz1*r1;
		r1 = eyi[0][j][zs-1]-eyizb1[0][j][zs];
		eyi[0][j][zs] = eyizb1[0][j][zs-1]+uz1*r1;
		r1 = eyi[xs][j][zs-1]-eyizb1[xs][j][zs];
		eyi[xs][j][zs] = eyizb1[xs][j][zs-1]+uz1*r1;
	}

	for(k=1;k<=(zs-2);k++){                      /*EZI*/
		ux41 = ux4;
		if(k <= (Fz+z_sub_length)){
			ux41 = ux4_sub;
		}

		for(j=1;j<=(ys-1);j++){
			ux31 = ux3;
			if((j >= (Fy-y_sub_length))
			    &&(j <= (Fy+y_sub_length))){
				ux31 = ux3_sub;
			}

			r1 = ezi[1][j][k]+ezixa2[0][j][k];
			r2 = ezixa1[0][j][k]+ezixa1[1][j][k];
			r3 = ezixa1[0][j+1][k]-2*ezixa1[0][j][k]
			    +ezixa1[0][j-1][k]+ezixa1[1][j+1][k]
			    -2*ezixa1[1][j][k]+ezixa1[1][j-1][k];
			r4 = ezixa1[0][j][k+1]-2*ezixa1[0][j][k]
			    +ezixa1[0][j][k-1]+ezixa1[1][j][k+1]
			    -2*ezixa1[1][j][k]+ezixa1[1][j][k-1];

			ezi[0][j][k] = -ezixa2[1][j][k]+ux1*r1
			    +ux2*r2+ux31*r3+ux41*r4;

			r1 = ezi[xs-1][j][k]+ezixb2[xs][j][k];
			r2 = ezixb1[xs][j][k]+ezixb1[xs-1][j][k];
			r3 = ezixb1[xs][j+1][k]-2*ezixb1[xs][j][k]
			    +ezixb1[xs][j-1][k]+ezixb1[xs-1][j+1][k]
			    -2*ezixb1[xs-1][j][k]+ezixb1[xs-1][j-1][k];
			r4 = ezixb1[xs][j][k+1]-2*ezixb1[xs][j][k]
			    +ezixb1[xs][j][k-1]+ezixb1[xs-1][j][k+1]
			    -2*ezixb1[xs-1][j][k]+ezixb1[xs-1][j][k-1];

			ezi[xs][j][k] = -ezixb2[xs-1][j][k]+ux1*r1
			    +ux2*r2+ux31*r3+ux41*r4;
		}
	}

	for(j=1;j<=(ys-1);j++){
		r1 = ezi[1][j][0]-ezixa1[0][j][0];
		ezi[0][j][0] = ezixa1[1][j][0]+ux1*r1;
		r1 = ezi[1][j][zs-1]-ezixa1[0][j][zs-1];
		ezi[0][j][zs-1] = ezixa1[1][j][zs-1]+ux1*r1;
		r1 = ezi[xs-1][j][0]-ezixb1[xs][j][0];
		ezi[xs][j][0] = ezixb1[xs-1][j][0]+ux1*r1;
		r1 = ezi[xs-1][j][zs-1]-ezixb1[xs][j][zs-1];
		ezi[xs][j][zs-1] = ezixb1[xs-1][j][zs-1]+ux1*r1;
	}

	for(k=1;k<=(zs-2);k++){
		uy31 = uy3;
		if(k <= (Fz+z_sub_length)){
			uy31 = uy3_sub;
		}

		for(i=1;i<=(xs-1);i++){
			uy41 = uy4;
			if((i >= (Fx-x_sub_length))
			    && (i <= (Fx+x_sub_length))){
				uy41 = uy4_sub;
			}

			r1 = ezi[i][1][k]+eziya2[i][0][k];
			r2 = eziya1[i][0][k]+eziya1[i][1][k];
			r3 = eziya1[i][0][k+1]-2*eziya1[i][0][k]
			    +eziya1[i][0][k-1]+eziya1[i][1][k+1]
			    -2*eziya1[i][1][k]+eziya1[i][1][k-1];
			r4 = eziya1[i+1][0][k]-2*eziya1[i][0][k]
			    +eziya1[i-1][0][k]+eziya1[i+1][1][k]
			    -2*eziya1[i][1][k]+eziya1[i-1][1][k];

			ezi[i][0][k] = -eziya2[i][1][k]+uy1*r1
			    +uy2*r2+uy31*r3+uy41*r4;

			r1 = ezi[i][ys-1][k]+eziyb2[i][ys][k];
			r2 = eziyb1[i][ys][k]+eziyb1[i][ys-1][k];
			r3 = eziyb1[i][ys][k+1]-2*eziyb1[i][ys][k]
			    +eziyb1[i][ys][k-1]+eziyb1[i][ys-1][k+1]
			    -2*eziyb1[i][ys-1][k]+eziyb1[i][ys-1][k-1];
			r4 = eziyb1[i+1][ys][k]-2*eziyb1[i][ys][k]
			    +eziyb1[i-1][ys][k]+eziyb1[i+1][ys-1][k]
			    -2*eziyb1[i][ys-1][k]+eziyb1[i-1][ys-1][k];

			ezi[i][ys][k] = -eziyb2[i][ys-1][k]+uy1*r1
			    +uy2*r2+uy31*r3+uy41*r4;
		}
	}

	for(i=1;i<=(xs-1);i++){
		r1 = ezi[i][1][0]-eziya1[i][0][0];
		ezi[i][0][0] = eziya1[i][1][0]+uy1*r1;
		r1 = ezi[i][1][zs-1]-eziya1[i][0][zs-1];
		ezi[i][0][zs-1] = eziya1[i][1][zs-1]+uy1*r1;
		r1 = ezi[i][ys-1][0]-eziyb1[i][ys][0];
		ezi[i][ys][0] = eziyb1[i][ys-1][0]+uy1*r1;
		r1 = ezi[i][ys-1][zs-1]-eziyb1[i][ys][zs-1];
		ezi[i][ys][zs-1] = eziyb1[i][ys-1][zs-1]+uy1*r1;
	}

	for(k=0;k<=(zs-1);k++){
		r1 = ezi[0][1][k]-eziya1[0][0][k];
		ezi[0][0][k] = eziya1[0][1][k]+uy1*r1;
		r1 = ezi[xs][1][k]-eziya1[xs][0][k];
		ezi[xs][0][k] = eziya1[xs][1][k]+uy1*r1;
		r1 = ezi[0][ys-1][k]-eziyb1[0][ys][k];
		ezi[0][ys][k] = eziyb1[0][ys-1][k]+uy1*r1;
		r1 = ezi[xs][ys-1][k]-eziyb1[xs][ys][k];
		ezi[xs][ys][k] = eziyb1[xs][ys-1][k]+uy1*r1;
	}
}