📄 fdtd_3d_pml_hx.cpp
字号:
#include "FDTD_3D_PML_Hx.h"
///////////////////////////////////////////////////////////////////////////////////////
//Calculate Hx field
///////////////////////////////////////////////////////////////////////////////////////
void CFDTD_3D_PML::calc_Hx_1()
{
long i, j, k;
double Bx_r;
//1
for (i = 0; i < nPML_x_1; i++)
{
for (j = 0; j < nPML_y_1; j++)
{
for (k = 0; k < nPML_z_1; k++)
{
Bx_r = Bx_1[i][j][k];
Bx_1[i][j][k] = K_Bx_a_1[j]*Bx_1[i][j][k] +
K_Bx_b_1[j]*( (Ey[i][j][k+1] - Ey[i][j][k])*inv_dz -
(Ez[i][j+1][k] - Ez[i][j][k])*inv_dy);
Hx[i][j][k] = K_Hx_a_1[k]*Hx[i][j][k] +
K_Hx_b_1[k]*( K_Hx_c_1[i]*Bx_1[i][j][k] +
K_Hx_d_1[i]*Bx_r )/Mat[Ind[i][j][k]][1];
}
}
}
}
void CFDTD_3D_PML::calc_Hx_2()
{
long i, j, k, ia;
double Bx_r;
//2
ia = 0;
for (i = nPML_x_1; i < nx_MIN_nPML_x_2; i++)
{
for (j = 0; j < nPML_y_1; j++)
{
for (k = 0; k < nPML_z_1; k++)
{
Bx_r = Bx_2[ia][j][k];
Bx_2[ia][j][k] = K_Bx_a_1[j]*Bx_2[ia][j][k] +
K_Bx_b_1[j]*( (Ey[i][j][k+1] - Ey[i][j][k])*inv_dz -
(Ez[i][j+1][k] - Ez[i][j][k])*inv_dy);
Hx[i][j][k] = K_Hx_a_1[k]*Hx[i][j][k] + K_Hx_b_1[k]*TwoOREp0DIVMu0*
( Bx_2[ia][j][k] - Bx_r )/Mat[Ind[i][j][k]][1];
}
}
ia++;
}
}
void CFDTD_3D_PML::calc_Hx_3()
{
long i, j, k, ii, ia;
double Bx_r;
//3
ia = 0;
ii = nPML_x_2_MIN_1;
for (i = nx_MIN_nPML_x_2; i < nx; i++)
{
for (j = 0; j < nPML_y_1; j++)
{
for (k = 0; k < nPML_z_1; k++)
{
Bx_r = Bx_3[ia][j][k];
Bx_3[ia][j][k] = K_Bx_a_1[j]*Bx_3[ia][j][k] +
K_Bx_b_1[j]*( (Ey[i][j][k+1] - Ey[i][j][k])*inv_dz -
(Ez[i][j+1][k] - Ez[i][j][k])*inv_dy);
Hx[i][j][k] = K_Hx_a_1[k]*Hx[i][j][k] +
K_Hx_b_1[k]*( K_Hx_c_2[ii]*Bx_3[ia][j][k] +
K_Hx_d_2[ii]*Bx_r )/Mat[Ind[i][j][k]][1];
}
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Hx_4()
{
long i, j, k, ja;
double Bx_r;
//4
for (i = 0; i < nPML_x_1; i++)
{
ja = 0;
for (j = nPML_y_1; j < ny_MIN_1_MIN_nPML_y_2; j++)
{
for (k = 0; k < nPML_z_1; k++)
{
Bx_r = Bx_4[i][ja][k];
Bx_4[i][ja][k] += ( (Ey[i][j][k+1] - Ey[i][j][k])*dtDIVdz -
(Ez[i][j+1][k] - Ez[i][j][k])*dtDIVdy );
Hx[i][j][k] = K_Hx_a_1[k]*Hx[i][j][k] +
K_Hx_b_1[k]*( K_Hx_c_1[i]*Bx_4[i][ja][k] +
K_Hx_d_1[i]*Bx_r )/Mat[Ind[i][j][k]][1];
}
ja++;
}
}
}
void CFDTD_3D_PML::calc_Hx_5()
{
long i, j, k;
//5
for (i = nPML_x_1; i < nx_MIN_nPML_x_2; i++)
{
for (j = nPML_y_1; j < ny_MIN_1_MIN_nPML_y_2; j++)
{
for (k = 0; k < nPML_z_1; k++)
{
Hx[i][j][k] = K_Hx_a_1[k]*Hx[i][j][k] + K_Hx_b_1[k]/Mat[Ind[i][j][k]][1]*
( (Ey[i][j][k+1] - Ey[i][j][k])*TwoOREp0DIVMu0ORdtDIVdz -
(Ez[i][j+1][k] - Ez[i][j][k])*TwoOREp0DIVMu0ORdtDIVdy );
}
}
}
}
void CFDTD_3D_PML::calc_Hx_6()
{
long i, j, k, ii, ia, ja;
double Bx_r;
//6
ia = 0;
ii = nPML_x_2_MIN_1;
for (i = nx_MIN_nPML_x_2; i < nx; i++)
{
ja = 0;
for (j = nPML_y_1; j < ny_MIN_1_MIN_nPML_y_2; j++)
{
for (k = 0; k < nPML_z_1; k++)
{
Bx_r = Bx_6[ia][ja][k];
Bx_6[ia][ja][k] += ( (Ey[i][j][k+1] - Ey[i][j][k])*dtDIVdz -
(Ez[i][j+1][k] - Ez[i][j][k])*dtDIVdy );
Hx[i][j][k] = K_Hx_a_1[k]*Hx[i][j][k] +
K_Hx_b_1[k]*( K_Hx_c_2[ii]*Bx_6[ia][ja][k] +
K_Hx_d_2[ii]*Bx_r )/Mat[Ind[i][j][k]][1];
}
ja++;
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Hx_7()
{
long i, j, k, jj, ja;
double Bx_r;
//7
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++)
{
for (k = 0; k < nPML_z_1; k++)
{
Bx_r = Bx_7[i][ja][k];
Bx_7[i][ja][k] = K_Bx_a_2[jj]*Bx_7[i][ja][k] +
K_Bx_b_2[jj]*( (Ey[i][j][k+1] - Ey[i][j][k])*inv_dz -
(Ez[i][j+1][k] - Ez[i][j][k])*inv_dy);
Hx[i][j][k] = K_Hx_a_1[k]*Hx[i][j][k] +
K_Hx_b_1[k]*( K_Hx_c_1[i]*Bx_7[i][ja][k] +
K_Hx_d_1[i]*Bx_r )/Mat[Ind[i][j][k]][1];
}
ja++;
jj--;
}
}
}
void CFDTD_3D_PML::calc_Hx_8()
{
long i, j, k, jj, ia, ja;
double Bx_r;
//8
ia = 0;
for (i = nPML_x_1; i < nx_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++)
{
for (k = 0; k < nPML_z_1; k++)
{
Bx_r = Bx_8[ia][ja][k];
Bx_8[ia][ja][k] = K_Bx_a_2[jj]*Bx_8[ia][ja][k] +
K_Bx_b_2[jj]*( (Ey[i][j][k+1] - Ey[i][j][k])*inv_dz -
(Ez[i][j+1][k] - Ez[i][j][k])*inv_dy);
Hx[i][j][k] = K_Hx_a_1[k]*Hx[i][j][k] + K_Hx_b_1[k]*TwoOREp0DIVMu0*
( Bx_8[ia][ja][k] - Bx_r )/Mat[Ind[i][j][k]][1];
}
ja++;
jj--;
}
ia++;
}
}
void CFDTD_3D_PML::calc_Hx_9()
{
long i, j, k, ii, jj, ia, ja;
double Bx_r;
//9
ia = 0;
ii = nPML_x_2_MIN_1;
for (i = nx_MIN_nPML_x_2; i < nx; i++)
{
ja = 0;
jj = nPML_y_2_MIN_1;
for (j = ny_MIN_1_MIN_nPML_y_2; j < ny_MIN_1; j++)
{
for (k = 0; k < nPML_z_1; k++)
{
Bx_r = Bx_9[ia][ja][k];
Bx_9[ia][ja][k] = K_Bx_a_2[jj]*Bx_9[ia][ja][k] +
K_Bx_b_2[jj]*( (Ey[i][j][k+1] - Ey[i][j][k])*inv_dz -
(Ez[i][j+1][k] - Ez[i][j][k])*inv_dy);
Hx[i][j][k] = K_Hx_a_1[k]*Hx[i][j][k] +
K_Hx_b_1[k]*( K_Hx_c_2[ii]*Bx_9[ia][ja][k] +
K_Hx_d_2[ii]*Bx_r )/Mat[Ind[i][j][k]][1];
}
ja++;
jj--;
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Hx_10()
{
long i, j, k, ka;
double Bx_r;
//10
for (i = 0; i < nPML_x_1; i++)
{
for (j = 0; j < nPML_y_1; j++)
{
ka = 0;
for (k = nPML_z_1; k < nz_MIN_1_MIN_nPML_z_2; k++)
{
Bx_r = Bx_10[i][j][ka];
Bx_10[i][j][ka] = K_Bx_a_1[j]*Bx_10[i][j][ka] +
K_Bx_b_1[j]*( (Ey[i][j][k+1] - Ey[i][j][k])*inv_dz -
(Ez[i][j+1][k] - Ez[i][j][k])*inv_dy);
Hx[i][j][k] = Hx[i][j][k] +
inv_TwoOREp0*( K_Hx_c_1[i]*Bx_10[i][j][ka] +
K_Hx_d_1[i]*Bx_r )/Mat[Ind[i][j][k]][1];
ka++;
}
}
}
}
void CFDTD_3D_PML::calc_Hx_11()
{
long i, j, k;
//11
for (i = nPML_x_1; i < nx_MIN_nPML_x_2; i++)
{
for (j = 0; j < nPML_y_1; j++)
{
for (k = nPML_z_1; k < nz_MIN_1_MIN_nPML_z_2; k++)
{
Hx[i][j][k] = K_Bx_a_1[j]*Hx[i][j][k] + K_Bx_b_1[j]/Mat[Ind[i][j][k]][1]*
( (Ey[i][j][k+1] - Ey[i][j][k])*inv_dzORmu_0 -
(Ez[i][j+1][k] - Ez[i][j][k])*inv_dyORmu_0);
}
}
}
}
void CFDTD_3D_PML::calc_Hx_12()
{
long i, j, k, ii, ia, ka;
double Bx_r;
//12
ia = 0;
ii = nPML_x_2_MIN_1;
for (i = nx_MIN_nPML_x_2; i < nx; i++)
{
for (j = 0; j < nPML_y_1; j++)
{
ka = 0;
for (k = nPML_z_1; k < nz_MIN_1_MIN_nPML_z_2; k++)
{
Bx_r = Bx_12[ia][j][ka];
Bx_12[ia][j][ka] = K_Bx_a_1[j]*Bx_12[ia][j][ka] +
K_Bx_b_1[j]*( (Ey[i][j][k+1] - Ey[i][j][k])*inv_dz -
(Ez[i][j+1][k] - Ez[i][j][k])*inv_dy);
Hx[i][j][k] += inv_TwoOREp0*( K_Hx_c_2[ii]*Bx_12[ia][j][ka] +
K_Hx_d_2[ii]*Bx_r )/Mat[Ind[i][j][k]][1];
ka++;
}
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Hx_13()
{
long i, j, k, ja, ka;
double Bx_r;
//13
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;
for (k = nPML_z_1; k < nz_MIN_1_MIN_nPML_z_2; k++)
{
Bx_r = Bx_13[i][ja][ka];
Bx_13[i][ja][ka] += ( (Ey[i][j][k+1] - Ey[i][j][k])*dtDIVdz -
(Ez[i][j+1][k] - Ez[i][j][k])*dtDIVdy );
Hx[i][j][k] += inv_TwoOREp0*( K_Hx_c_1[i]*Bx_13[i][ja][ka] +
K_Hx_d_1[i]*Bx_r )/Mat[Ind[i][j][k]][1];
ka++;
}
ja++;
}
}
}
void CFDTD_3D_PML::calc_Hx_14(long nx_a, long nx_b, long ny_a, long ny_b, long nz_a, long nz_b)
{
long i, j, k;
//14 - Outside of the PML region
for (i = nx_a; i < nx_b; i++)
{
for (j = ny_a; j < ny_b; j++)
{
for (k = nz_a; k < nz_b; k++)
{
Hx[i][j][k] += ( (Ey[i][j][k+1] - Ey[i][j][k])*dtDIVMu0DIVdz -
(Ez[i][j+1][k] - Ez[i][j][k])*dtDIVMu0DIVdy )/
Mat[Ind[i][j][k]][1];
}
}
}
}
void CFDTD_3D_PML::calc_Hx_15()
{
long i, j, k, ii, ia, ja, ka;
double Bx_r;
//15
ia = 0;
ii = nPML_x_2_MIN_1;
for (i = nx_MIN_nPML_x_2; i < nx; i++)
{
ja = 0;
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -