📄 fdtd_3d_pml_hy.cpp
字号:
{
Hy[i][j][k] = K_Hy_a_2[ii]*Hy[i][j][k] + K_Hy_b_2[ii]/Mat[Ind[i][j][k]][1]*
( (Ez[i+1][j][k] - Ez[i][j][k])*TwoOREp0DIVMu0ORdtDIVdx -
(Ex[i][j][k+1] - Ex[i][j][k])*TwoOREp0DIVMu0ORdtDIVdz );
}
}
ii--;
}
}
void CFDTD_3D_PML::calc_Hy_16()
{
long i, j, k, ja, ka, jj;
double By_r;
//16
for (i = 0; i < nPML_x_1; i++)
{
ja = 0;
jj = nPML_y_2_MIN_1;
for (j = ny_MIN_nPML_y_2; j < ny; j++)
{
ka = 0;
for (k = nPML_z_1; k < nz_MIN_1_MIN_nPML_z_2; k++)
{
By_r = By_16[i][ja][ka];
By_16[i][ja][ka] += ( (Ez[i+1][j][k] - Ez[i][j][k])*dtDIVdx -
(Ex[i][j][k+1] - Ex[i][j][k])*dtDIVdz );
Hy[i][j][k] = K_Hy_a_1[i]*Hy[i][j][k] +
K_Hy_b_1[i]*( K_Hy_c_2[jj]*By_16[i][ja][ka] +
K_Hy_d_2[jj]*By_r )/Mat[Ind[i][j][k]][1];
ka++;
}
ja++;
jj--;
}
}
}
void CFDTD_3D_PML::calc_Hy_17()
{
long i, j, k, ia, ja, ka, jj;
double By_r;
//17
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_nPML_y_2; j < ny; j++)
{
ka = 0;
for (k = nPML_z_1; k < nz_MIN_1_MIN_nPML_z_2; k++)
{
By_r = By_17[ia][ja][ka];
By_17[ia][ja][ka] += ( (Ez[i+1][j][k] - Ez[i][j][k])*dtDIVdx -
(Ex[i][j][k+1] - Ex[i][j][k])*dtDIVdz );
Hy[i][j][k] += inv_TwoOREp0*( K_Hy_c_2[jj]*By_17[ia][ja][ka] +
K_Hy_d_2[jj]*By_r )/Mat[Ind[i][j][k]][1];
ka++;
}
ja++;
jj--;
}
ia++;
}
}
void CFDTD_3D_PML::calc_Hy_18()
{
long i, j, k, ia, ja, ka, ii, jj;
double By_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_nPML_y_2; j < ny; j++)
{
ka = 0;
for (k = nPML_z_1; k < nz_MIN_1_MIN_nPML_z_2; k++)
{
By_r = By_18[ia][ja][ka];
By_18[ia][ja][ka] += ( (Ez[i+1][j][k] - Ez[i][j][k])*dtDIVdx -
(Ex[i][j][k+1] - Ex[i][j][k])*dtDIVdz );
Hy[i][j][k] = K_Hy_a_2[ii]*Hy[i][j][k] +
K_Hy_b_2[ii]*( K_Hy_c_2[jj]*By_18[ia][ja][ka] +
K_Hy_d_2[jj]*By_r )/Mat[Ind[i][j][k]][1];
ka++;
}
ja++;
jj--;
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Hy_19()
{
long i, j, k, ka, kk;
double By_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_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
By_r = By_19[i][j][ka];
By_19[i][j][ka] = K_By_a_2[kk]*By_19[i][j][ka] +
K_By_b_2[kk]*( (Ez[i+1][j][k] - Ez[i][j][k])*inv_dx -
(Ex[i][j][k+1] - Ex[i][j][k])*inv_dz );
Hy[i][j][k] = K_Hy_a_1[i]*Hy[i][j][k] +
K_Hy_b_1[i]*( K_Hy_c_1[j]*By_19[i][j][ka] +
K_Hy_d_1[j]*By_r )/Mat[Ind[i][j][k]][1];
ka++;
kk--;
}
}
}
}
void CFDTD_3D_PML::calc_Hy_20()
{
long i, j, k, ia, ka, kk;
double By_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_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
By_r = By_20[ia][j][ka];
By_20[ia][j][ka] = K_By_a_2[kk]*By_20[ia][j][ka] +
K_By_b_2[kk]*( (Ez[i+1][j][k] - Ez[i][j][k])*inv_dx -
(Ex[i][j][k+1] - Ex[i][j][k])*inv_dz );
Hy[i][j][k] += inv_TwoOREp0*( K_Hy_c_1[j]*By_20[ia][j][ka] +
K_Hy_d_1[j]*By_r )/Mat[Ind[i][j][k]][1];
ka++;
kk--;
}
}
ia++;
}
}
void CFDTD_3D_PML::calc_Hy_21()
{
long i, j, k, ia, ka, ii, kk;
double By_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_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
By_r = By_21[ia][j][ka];
By_21[ia][j][ka] = K_By_a_2[kk]*By_21[ia][j][ka] +
K_By_b_2[kk]*( (Ez[i+1][j][k] - Ez[i][j][k])*inv_dx -
(Ex[i][j][k+1] - Ex[i][j][k])*inv_dz );
Hy[i][j][k] = K_Hy_a_2[ii]*Hy[i][j][k] +
K_Hy_b_2[ii]*( K_Hy_c_1[j]*By_21[ia][j][ka] +
K_Hy_d_1[j]*By_r )/Mat[Ind[i][j][k]][1];
ka++;
kk--;
}
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Hy_22()
{
long i, j, k, ja, ka, kk;
double By_r;
//22
for (i = 0; i < nPML_x_1; i++)
{
ja = 0;
for (j = nPML_y_1; j < ny_MIN_nPML_y_2; j++)
{
ka = 0;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
By_r = By_22[i][ja][ka];
By_22[i][ja][ka] = K_By_a_2[kk]*By_22[i][ja][ka] +
K_By_b_2[kk]*( (Ez[i+1][j][k] - Ez[i][j][k])*inv_dx -
(Ex[i][j][k+1] - Ex[i][j][k])*inv_dz );
Hy[i][j][k] = K_Hy_a_1[i]*Hy[i][j][k] +
K_Hy_b_1[i]*TwoOREp0DIVMu0*( By_22[i][ja][ka] - By_r )/Mat[Ind[i][j][k]][1];
ka++;
kk--;
}
ja++;
}
}
}
void CFDTD_3D_PML::calc_Hy_23()
{
long i, j, k, kk;
//23
for (i = nPML_x_1; i < nx_MIN_1_MIN_nPML_x_2; i++)
{
for (j = nPML_y_1; j < ny_MIN_nPML_y_2; j++)
{
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
Hy[i][j][k] = K_By_a_2[kk]*Hy[i][j][k] + K_By_b_2[kk]/Mat[Ind[i][j][k]][1]*
( (Ez[i+1][j][k] - Ez[i][j][k])*inv_dxORmu_0 -
(Ex[i][j][k+1] - Ex[i][j][k])*inv_dzORmu_0 );
kk--;
}
}
}
}
void CFDTD_3D_PML::calc_Hy_24()
{
long i, j, k, ia, ja, ka, ii, kk;
double By_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_nPML_y_2; j++)
{
ka = 0;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
By_r = By_24[ia][ja][ka];
By_24[ia][ja][ka] = K_By_a_2[kk]*By_24[ia][ja][ka] +
K_By_b_2[kk]*( (Ez[i+1][j][k] - Ez[i][j][k])*inv_dx -
(Ex[i][j][k+1] - Ex[i][j][k])*inv_dz );
Hy[i][j][k] = K_Hy_a_2[ii]*Hy[i][j][k] +
K_Hy_b_2[ii]*TwoOREp0DIVMu0*( By_24[ia][ja][ka] - By_r )/Mat[Ind[i][j][k]][1];
ka++;
kk--;
}
ja++;
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Hy_25()
{
long i, j, k, ja, ka, jj, kk;
double By_r;
//25
for (i = 0; i < nPML_x_1; i++)
{
ja = 0;
jj = nPML_y_2_MIN_1;
for (j = ny_MIN_nPML_y_2; j < ny; j++)
{
ka = 0;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
By_r = By_25[i][ja][ka];
By_25[i][ja][ka] = K_By_a_2[kk]*By_25[i][ja][ka] +
K_By_b_2[kk]*( (Ez[i+1][j][k] - Ez[i][j][k])*inv_dx -
(Ex[i][j][k+1] - Ex[i][j][k])*inv_dz );
Hy[i][j][k] = K_Hy_a_1[i]*Hy[i][j][k] +
K_Hy_b_1[i]*( K_Hy_c_2[jj]*By_25[i][ja][ka] +
K_Hy_d_2[jj]*By_r )/Mat[Ind[i][j][k]][1];
ka++;
kk--;
}
ja++;
jj--;
}
}
}
void CFDTD_3D_PML::calc_Hy_26()
{
long i, j, k, ia, ja, ka, jj, kk;
double By_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_nPML_y_2; j < ny; j++)
{
ka = 0;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
By_r = By_26[ia][ja][ka];
By_26[ia][ja][ka] = K_By_a_2[kk]*By_26[ia][ja][ka] +
K_By_b_2[kk]*( (Ez[i+1][j][k] - Ez[i][j][k])*inv_dx -
(Ex[i][j][k+1] - Ex[i][j][k])*inv_dz );
Hy[i][j][k] += inv_TwoOREp0*( K_Hy_c_2[jj]*By_26[ia][ja][ka] +
K_Hy_d_2[jj]*By_r )/Mat[Ind[i][j][k]][1];
ka++;
kk--;
}
ja++;
jj--;
}
ia++;
}
}
void CFDTD_3D_PML::calc_Hy_27()
{
long i, j, k, ia, ja, ka, ii, jj, kk;
double By_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_nPML_y_2; j < ny; j++)
{
ka = 0;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
By_r = By_27[ia][ja][ka];
By_27[ia][ja][ka] = K_By_a_2[kk]*By_27[ia][ja][ka] +
K_By_b_2[kk]*( (Ez[i+1][j][k] - Ez[i][j][k])*inv_dx -
(Ex[i][j][k+1] - Ex[i][j][k])*inv_dz );
Hy[i][j][k] = K_Hy_a_2[ii]*Hy[i][j][k] +
K_Hy_b_2[ii]*( K_Hy_c_2[jj]*By_27[ia][ja][ka] +
K_Hy_d_2[jj]*By_r )/Mat[Ind[i][j][k]][1];
ka++;
kk--;
}
ja++;
jj--;
}
ia++;
ii--;
}
}⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -