initial.c

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

/*
 *
 *   Initialize  of Electric and Magnetic Field Components
 *
 */

void initial(void)
{
	int i,j,k;

	for(i=0;i<=(xs-1);i++){
		for(j=0;j<=ys;j++){
			for(k=0;k<=zs;k++){
				exi[i][j][k] = 0.0;
			}
		}
	}

	for(i=0;i<=xs;i++){
		for(j=0;j<=(ys-1);j++){
			for(k=0;k<=zs;k++){
				eyi[i][j][k] = 0.0;
			}
		}
	}

	for(i=0;i<=xs;i++){
		for(j=0;j<=ys;j++){
			for(k=0;k<=(zs-1);k++){
				ezi[i][j][k] = 0.0;
			}
		}
	}

	for(i=0;i<=xs;i++){
		for(j=0;j<=(ys-1);j++){
			for(k=0;k<=(zs-1);k++){
				hxi[i][j][k] = 0.0;
			}
		}
	}

	for(i=0;i<=(xs-1);i++){
		for(j=0;j<=ys;j++){
			for(k=0;k<=(zs-1);k++){
				hyi[i][j][k] = 0.0;
			}
		}
	}

	for(i=0;i<=(xs-1);i++){
		for(j=0;j<=(ys-1);j++){
			for(k=0;k<=zs;k++){
				hzi[i][j][k] = 0.0;
			}
		}
	}

	for(i=0;i<=xs;i++){
		for(j=0;j<=1;j++){
			for(k=0;k<=(zs-1);k++){
				hx_y_absorb[i][j][k] = 0.0;
			}
		}
	}

	for(i=0;i<=xs;i++){
		for(j=0;j<=(ys-1);j++){
			for(k=0;k<=1;k++){
				hx_z_absorb[i][j][k] = 0.0;
			}
		}
	}

	for(i=0;i<=1;i++){
		for(j=0;j<=ys;j++){
			for(k=0;k<=(zs-1);k++){
				hy_x_absorb[i][j][k] = 0.0;
			}
		}
	}

	for(i=0;i<=(xs-1);i++){
		for(j=0;j<=ys;j++){
			for(k=0;k<=1;k++){
				hy_z_absorb[i][j][k] = 0.0;
			}
		}
	}

	for(i=0;i<=1;i++){
		for(j=0;j<=(ys-1);j++){
			for(k=0;k<=zs;k++){
				hz_x_absorb[i][j][k] = 0.0;
			}
		}
	}


	for(i=0;i<=(xs-1);i++){
		for(j=0;j<=1;j++){
			for(k=0;k<=zs;k++){
				hz_y_absorb[i][j][k] = 0.0;
			}
		}
	}
}