fdtd_3d_pml_hz.cpp

采用UMPL边界的FDTD编程

				Hz[i][j][k] = K_Bz_a_2[ii]*Hz[i][j][k] + K_Bz_b_2[ii]/Mat[Ind[i][j][k]][1]*
							      ( (Ex[i][j+1][k] - Ex[i][j][k])*inv_dyORmu_0 - 
				                    (Ey[i+1][j][k] - Ey[i][j][k])*inv_dxORmu_0 );
			}
		}
		ii--;
	}
}

void CFDTD_3D_PML::calc_Hz_16()
{
	long  i, j, k, ja, ka, jj;
	double Bz_r;

	//16
	for (i = 0; i < nPML_x_1; i++)
	{	
		ja = 0;
		jj = nPML_y_2_MIN_1;
		for (j = ny_MIN_1_MIN_nPML_y_2; j < ny_MIN_1; j++)
		{	
			ka = 0;
			for (k = nPML_z_1; k < nz_MIN_nPML_z_2; k++)
			{
				Bz_r = Bz_16[i][ja][ka];

				Bz_16[i][ja][ka] = K_Bz_a_1[i]*Bz_16[i][ja][ka] +  
					               K_Bz_b_1[i]*( (Ex[i][j+1][k] - Ex[i][j][k])*inv_dy - 
					                             (Ey[i+1][j][k] - Ey[i][j][k])*inv_dx );

				Hz[i][j][k] = K_Hz_a_2[jj]*Hz[i][j][k] + 
					          K_Hz_b_2[jj]*TwoOREp0DIVMu0*( Bz_16[i][ja][ka] - Bz_r )/Mat[Ind[i][j][k]][1];
				ka++;
			}
			ja++;
			jj--;
		}
	}
}

void CFDTD_3D_PML::calc_Hz_17()
{
	long  i, j, k, jj;
	
	//17
	for (i = nPML_x_1; i < nx_MIN_1_MIN_nPML_x_2; i++)
	{	
		jj = nPML_y_2_MIN_1;
		for (j = ny_MIN_1_MIN_nPML_y_2; j < ny_MIN_1; j++)
		{	
			for (k = nPML_z_1; k < nz_MIN_nPML_z_2; k++)
			{
				Hz[i][j][k] = K_Hz_a_2[jj]*Hz[i][j][k] +  K_Hz_b_2[jj]/Mat[Ind[i][j][k]][1]*
					          ( (Ex[i][j+1][k] - Ex[i][j][k])*TwoOREp0DIVMu0ORdtDIVdy - 
					            (Ey[i+1][j][k] - Ey[i][j][k])*TwoOREp0DIVMu0ORdtDIVdx );
			}
			jj--;
		}
	}
}

void CFDTD_3D_PML::calc_Hz_18()
{
	long  i, j, k, ia, ja, ka, ii, jj;
	double Bz_r;

	//18
	ia = 0;
	ii = nPML_x_2_MIN_1;
	for (i = nx_MIN_1_MIN_nPML_x_2; i < nx_MIN_1; i++)
	{	
		ja = 0;
		jj = nPML_y_2_MIN_1;
		for (j = ny_MIN_1_MIN_nPML_y_2; j < ny_MIN_1; j++)
		{	
			ka = 0;
			for (k = nPML_z_1; k < nz_MIN_nPML_z_2; k++)
			{
				Bz_r = Bz_18[ia][ja][ka];

				Bz_18[ia][ja][ka] = K_Bz_a_2[ii]*Bz_18[ia][ja][ka] +  
					                K_Bz_b_2[ii]*( (Ex[i][j+1][k] - Ex[i][j][k])*inv_dy - 
					                               (Ey[i+1][j][k] - Ey[i][j][k])*inv_dx );

				Hz[i][j][k] = K_Hz_a_2[jj]*Hz[i][j][k] + K_Hz_b_2[jj]*TwoOREp0DIVMu0*
					          ( Bz_18[ia][ja][ka] - Bz_r )/Mat[Ind[i][j][k]][1];
				ka++;
			}
			ja++;
			jj--;
		}
		ia++;
		ii--;
	}
}

void CFDTD_3D_PML::calc_Hz_19()
{
	long  i, j, k, ka, kk;
	double Bz_r;

	//19
	for (i = 0; i < nPML_x_1; i++)
	{	
		for (j = 0; j < nPML_y_1; j++)
		{	
			ka = 0;
			kk = nPML_z_2_MIN_1;
			for (k = nz_MIN_nPML_z_2; k < nz; k++)
			{
				Bz_r = Bz_19[i][j][ka];

				Bz_19[i][j][ka] = K_Bz_a_1[i]*Bz_19[i][j][ka] +  
					              K_Bz_b_1[i]*( (Ex[i][j+1][k] - Ex[i][j][k])*inv_dy - 
					                            (Ey[i+1][j][k] - Ey[i][j][k])*inv_dx );

				Hz[i][j][k] = K_Hz_a_1[j]*Hz[i][j][k] + 
					          K_Hz_b_1[j]*( K_Hz_c_2[kk]*Bz_19[i][j][ka] + 
							                K_Hz_d_2[kk]*Bz_r )/Mat[Ind[i][j][k]][1];
				ka++;
				kk--;
			}
		}
	}
}

void CFDTD_3D_PML::calc_Hz_20()
{
	long  i, j, k, ia, ka, kk;
	double Bz_r;

	//20
	ia = 0;
	for (i = nPML_x_1; i < nx_MIN_1_MIN_nPML_x_2; i++)
	{	
		for (j = 0; j < nPML_y_1; j++)
		{	
			ka = 0;
			kk = nPML_z_2_MIN_1;
			for (k = nz_MIN_nPML_z_2; k < nz; k++)
			{
				Bz_r = Bz_20[ia][j][ka];

				Bz_20[ia][j][ka] +=  (Ex[i][j+1][k] - Ex[i][j][k])*dtDIVdy - 
					                 (Ey[i+1][j][k] - Ey[i][j][k])*dtDIVdx;

				Hz[i][j][k] = K_Hz_a_1[j]*Hz[i][j][k] + 
					          K_Hz_b_1[j]*( K_Hz_c_2[kk]*Bz_20[ia][j][ka] + 
							                K_Hz_d_2[kk]*Bz_r )/Mat[Ind[i][j][k]][1];
				ka++;
				kk--;
			}
		}
		ia++;
	}
}

void CFDTD_3D_PML::calc_Hz_21()
{
	long  i, j, k, ia, ka, ii, kk;
	double Bz_r;

	//21
	ia = 0;
	ii = nPML_x_2_MIN_1;
	for (i = nx_MIN_1_MIN_nPML_x_2; i < nx_MIN_1; i++)
	{	
		for (j = 0; j < nPML_y_1; j++)
		{	
			ka = 0;
			kk = nPML_z_2_MIN_1;
			for (k = nz_MIN_nPML_z_2; k < nz; k++)
			{
				Bz_r = Bz_21[ia][j][ka];

				Bz_21[ia][j][ka] = K_Bz_a_2[ii]*Bz_21[ia][j][ka] +  
					               K_Bz_b_2[ii]*( (Ex[i][j+1][k] - Ex[i][j][k])*inv_dy - 
					                              (Ey[i+1][j][k] - Ey[i][j][k])*inv_dx );

				Hz[i][j][k] = K_Hz_a_1[j]*Hz[i][j][k] + 
					          K_Hz_b_1[j]*( K_Hz_c_2[kk]*Bz_21[ia][j][ka] + 
							                K_Hz_d_2[kk]*Bz_r )/Mat[Ind[i][j][k]][1];
				ka++;
				kk--;
			}
		}
		ia++;
		ii--;
	}
}

void CFDTD_3D_PML::calc_Hz_22()
{
	long  i, j, k, ja, ka, kk;
	double Bz_r;

	//22
	for (i = 0; i < nPML_x_1; i++)
	{	
		ja = 0;
		for (j = nPML_y_1; j < ny_MIN_1_MIN_nPML_y_2; j++)
		{	
			ka = 0;
			kk = nPML_z_2_MIN_1;
			for (k = nz_MIN_nPML_z_2; k < nz; k++)
			{
				Bz_r = Bz_22[i][ja][ka];

				Bz_22[i][ja][ka] = K_Bz_a_1[i]*Bz_22[i][ja][ka] +  
					               K_Bz_b_1[i]*( (Ex[i][j+1][k] - Ex[i][j][k])*inv_dy - 
					                             (Ey[i+1][j][k] - Ey[i][j][k])*inv_dx );

				Hz[i][j][k] += inv_TwoOREp0*( K_Hz_c_2[kk]*Bz_22[i][ja][ka] + 
							                  K_Hz_d_2[kk]*Bz_r )/Mat[Ind[i][j][k]][1];
				ka++;
				kk--;
			}
			ja++;
		}
	}
}

void CFDTD_3D_PML::calc_Hz_23()
{
	long  i, j, k, ia, ja, ka, kk;
	double Bz_r;

	//23
	ia = 0;
	for (i = nPML_x_1; i < nx_MIN_1_MIN_nPML_x_2; i++)
	{	
		ja = 0;
		for (j = nPML_y_1; j < ny_MIN_1_MIN_nPML_y_2; j++)
		{	
			ka = 0;
			kk = nPML_z_2_MIN_1;
			for (k = nz_MIN_nPML_z_2; k < nz; k++)
			{
				Bz_r = Bz_23[ia][ja][ka];

				Bz_23[ia][ja][ka] += (Ex[i][j+1][k] - Ex[i][j][k])*dtDIVdy - 
					                 (Ey[i+1][j][k] - Ey[i][j][k])*dtDIVdx;

				Hz[i][j][k] += inv_TwoOREp0*( K_Hz_c_2[kk]*Bz_23[ia][ja][ka] + 
							                  K_Hz_d_2[kk]*Bz_r )/Mat[Ind[i][j][k]][1];
				ka++;
				kk--;
			}
			ja++;
		}
		ia++;
	}
}

void CFDTD_3D_PML::calc_Hz_24()
{
	long  i, j, k, ia, ja, ka, ii, kk;
	double Bz_r;

	//24
	ia = 0;
	ii = nPML_x_2_MIN_1;
	for (i = nx_MIN_1_MIN_nPML_x_2; i < nx_MIN_1; i++)
	{	
		ja = 0;
		for (j = nPML_y_1; j < ny_MIN_1_MIN_nPML_y_2; j++)
		{	
			ka = 0;
			kk = nPML_z_2_MIN_1;
			for (k = nz_MIN_nPML_z_2; k < nz; k++)
			{
				Bz_r = Bz_24[ia][ja][ka];

				Bz_24[ia][ja][ka] = K_Bz_a_2[ii]*Bz_24[ia][ja][ka] +  
					                K_Bz_b_2[ii]*( (Ex[i][j+1][k] - Ex[i][j][k])*inv_dy - 
					                               (Ey[i+1][j][k] - Ey[i][j][k])*inv_dx );

				Hz[i][j][k] += inv_TwoOREp0*( K_Hz_c_2[kk]*Bz_24[ia][ja][ka] + 
							                  K_Hz_d_2[kk]*Bz_r )/Mat[Ind[i][j][k]][1];
				ka++;
				kk--;
			}
			ja++;
		}
		ia++;
		ii--;
	}
}

void CFDTD_3D_PML::calc_Hz_25()
{
	long  i, j, k, ja, ka, jj, kk;
	double Bz_r;

	//25
	for (i = 0; i < nPML_x_1; i++)
	{	
		ja = 0;
		jj = nPML_y_2_MIN_1;
		for (j = ny_MIN_1_MIN_nPML_y_2; j < ny_MIN_1; j++)
		{	
			ka = 0;
			kk = nPML_z_2_MIN_1;
			for (k = nz_MIN_nPML_z_2; k < nz; k++)
			{
				Bz_r = Bz_25[i][ja][ka];

				Bz_25[i][ja][ka] = K_Bz_a_1[i]*Bz_25[i][ja][ka] +  
					               K_Bz_b_1[i]*( (Ex[i][j+1][k] - Ex[i][j][k])*inv_dy - 
					                             (Ey[i+1][j][k] - Ey[i][j][k])*inv_dx );

				Hz[i][j][k] = K_Hz_a_2[jj]*Hz[i][j][k] + 
					          K_Hz_b_2[jj]*( K_Hz_c_2[kk]*Bz_25[i][ja][ka] + 
							                 K_Hz_d_2[kk]*Bz_r )/Mat[Ind[i][j][k]][1];
				ka++;
				kk--;
			}
			ja++;
			jj--;
		}
	}
}

void CFDTD_3D_PML::calc_Hz_26()
{
	long  i, j, k, ia, ja, ka, jj, kk;
	double Bz_r;

	//26
	ia = 0;
	for (i = nPML_x_1; i < nx_MIN_1_MIN_nPML_x_2; i++)
	{	
		ja = 0;
		jj = nPML_y_2_MIN_1;
		for (j = ny_MIN_1_MIN_nPML_y_2; j < ny_MIN_1; j++)
		{	
			ka = 0;
			kk = nPML_z_2_MIN_1;
			for (k = nz_MIN_nPML_z_2; k < nz; k++)
			{
				Bz_r = Bz_26[ia][ja][ka];

				Bz_26[ia][ja][ka] += (Ex[i][j+1][k] - Ex[i][j][k])*dtDIVdy - 
					                 (Ey[i+1][j][k] - Ey[i][j][k])*dtDIVdx;

				Hz[i][j][k] = K_Hz_a_2[jj]*Hz[i][j][k] + 
					          K_Hz_b_2[jj]*( K_Hz_c_2[kk]*Bz_26[ia][ja][ka] + 
							                 K_Hz_d_2[kk]*Bz_r )/Mat[Ind[i][j][k]][1];
				ka++;
				kk--;
			}
			ja++;
			jj--;
		}
		ia++;
	}
}

void CFDTD_3D_PML::calc_Hz_27()
{
	long  i, j, k, ia, ja, ka, ii, jj, kk;
	double Bz_r;

	//27
	ia = 0;
	ii = nPML_x_2_MIN_1;
	for (i = nx_MIN_1_MIN_nPML_x_2; i < nx_MIN_1; i++)
	{	
		ja = 0;
		jj = nPML_y_2_MIN_1;
		for (j = ny_MIN_1_MIN_nPML_y_2; j < ny_MIN_1; j++)
		{	
			ka = 0;
			kk = nPML_z_2_MIN_1;
			for (k = nz_MIN_nPML_z_2; k < nz; k++)
			{
				Bz_r = Bz_27[ia][ja][ka];

				Bz_27[ia][ja][ka] = K_Bz_a_2[ii]*Bz_27[ia][ja][ka] +  
					                K_Bz_b_2[ii]*( (Ex[i][j+1][k] - Ex[i][j][k])*inv_dy - 
					                               (Ey[i+1][j][k] - Ey[i][j][k])*inv_dx );

				Hz[i][j][k] = K_Hz_a_2[jj]*Hz[i][j][k] + 
					          K_Hz_b_2[jj]*( K_Hz_c_2[kk]*Bz_27[ia][ja][ka] + 
							                 K_Hz_d_2[kk]*Bz_r )/Mat[Ind[i][j][k]][1];
				ka++;
				kk--;
			}
			ja++;
			jj--;
		}
		ia++;
		ii--;
	}
}