📄 fdtd_3d_pml_ex.cpp
字号:
#include "FDTD_3D_PML_Ex.h"
///////////////////////////////////////////////////////////////////////////////////////
//Calculate Ex field - decomposed in 27 blocks
///////////////////////////////////////////////////////////////////////////////////////
void CFDTD_3D_PML::calc_Ex_1()
{
double Fx_r, Gx_r;
long i, j, k;
//1
for (i = 0; i < nPML_x_1; i++)
{
for (j = 1; j < nPML_y_1; j++)
{
for (k = 1; k < nPML_z_1; k++)
{
Fx_r = Fx_1[i][j][k];
Fx_1[i][j][k] = K_a[Ind[i][j][k]]*Fx_1[i][j][k] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
Gx_r = Gx_1[i][j][k];
Gx_1[i][j][k] = K_Gx_a_1[j]*Gx_1[i][j][k] +
K_Gx_b_1[j]*(Fx_1[i][j][k] - Fx_r);
Ex[i][j][k] = K_Ex_a_1[k]*Ex[i][j][k] +
K_Ex_b_1[k]*( K_Ex_c_1[i]*Gx_1[i][j][k] + K_Ex_d_1[i]*Gx_r);
}
}
}
}
void CFDTD_3D_PML::calc_Ex_2()
{
double Fx_r, Gx_r;
long i, j, k, ia;
//2
ia = 0;
for (i = nPML_x_1; i < nx_MIN_1_MIN_nPML_x_2; i++)
{
for (j = 1; j < nPML_y_1; j++)
{
for (k = 1; k < nPML_z_1; k++)
{
Fx_r = Fx_2[ia][j][k];
Fx_2[ia][j][k] = K_a[Ind[i][j][k]]*Fx_2[ia][j][k] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
Gx_r = Gx_2[ia][j][k];
Gx_2[ia][j][k] = K_Gx_a_1[j]*Gx_2[ia][j][k] +
K_Gx_b_1[j]*(Fx_2[ia][j][k] - Fx_r);
Ex[i][j][k] = K_Ex_a_1[k]*Ex[i][j][k] +
K_Ex_b_1[k]*TwoOREp0*( Gx_2[ia][j][k] - Gx_r);
}
}
ia++;
}
}
void CFDTD_3D_PML::calc_Ex_3()
{
double Fx_r, Gx_r;
long i, j, k, ii, ia;
//3
ia = 0;
ii = nPML_x_2_MIN_1;
for (i = nx_MIN_1_MIN_nPML_x_2; i < nx_MIN_1; i++)
{
for (j = 1; j < nPML_y_1; j++)
{
for (k = 1; k < nPML_z_1; k++)
{
Fx_r = Fx_3[ia][j][k];
Fx_3[ia][j][k] = K_a[Ind[i][j][k]]*Fx_3[ia][j][k] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
Gx_r = Gx_3[ia][j][k];
Gx_3[ia][j][k] = K_Gx_a_1[j]*Gx_3[ia][j][k] +
K_Gx_b_1[j]*(Fx_3[ia][j][k] - Fx_r);
Ex[i][j][k] = K_Ex_a_1[k]*Ex[i][j][k] +
K_Ex_b_1[k]*( K_Ex_c_2[ii]*Gx_3[ia][j][k] + K_Ex_d_2[ii]*Gx_r);
}
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Ex_4()
{
double Fx_r;
long i, j, k, ja;
//4
for (i = 0; i < nPML_x_1; i++)
{
ja = 0;
for (j = nPML_y_1; j < ny_MIN_nPML_y_2; j++)
{
for (k = 1; k < nPML_z_1; k++)
{
Fx_r = Fx_4[i][ja][k];
Fx_4[i][ja][k] = K_a[Ind[i][j][k]]*Fx_4[i][ja][k] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
Ex[i][j][k] = K_Ex_a_1[k]*Ex[i][j][k] +
K_Ex_b_1[k]*( K_Ex_c_1[i]*Fx_4[i][ja][k] + K_Ex_d_1[i]*Fx_r);
}
ja++;
}
}
}
void CFDTD_3D_PML::calc_Ex_5()
{
double Fx_r;
long i, j, k, ia, ja;
//5
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_nPML_y_2; j++)
{
for (k = 1; k < nPML_z_1; k++)
{
Fx_r = Fx_5[ia][ja][k];
Fx_5[ia][ja][k] = K_a[Ind[i][j][k]]*Fx_5[ia][ja][k] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
Ex[i][j][k] = K_Ex_a_1[k]*Ex[i][j][k] +
K_Ex_b_1[k]*TwoOREp0*( Fx_5[ia][ja][k] - Fx_r);
}
ja++;
}
ia++;
}
}
void CFDTD_3D_PML::calc_Ex_6()
{
double Fx_r;
long i, j, k, ii, ia, ja;
//6
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++)
{
for (k = 1; k < nPML_z_1; k++)
{
Fx_r = Fx_6[ia][ja][k];
Fx_6[ia][ja][k] = K_a[Ind[i][j][k]]*Fx_6[ia][ja][k] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
Ex[i][j][k] = K_Ex_a_1[k]*Ex[i][j][k] +
K_Ex_b_1[k]*(K_Ex_c_2[ii]*Fx_6[ia][ja][k] + K_Ex_d_2[ii]*Fx_r);
}
ja++;
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Ex_7()
{
double Fx_r, Gx_r;
long i, j, k, jj, ja;
//7
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_MIN_1; j++)
{
for (k = 1; k < nPML_z_1; k++)
{
Fx_r = Fx_7[i][ja][k];
Fx_7[i][ja][k] = K_a[Ind[i][j][k]]*Fx_7[i][ja][k] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
Gx_r = Gx_7[i][ja][k];
Gx_7[i][ja][k] = K_Gx_a_2[jj]*Gx_7[i][ja][k] +
K_Gx_b_2[jj]*(Fx_7[i][ja][k] - Fx_r);
Ex[i][j][k] = K_Ex_a_1[k]*Ex[i][j][k] +
K_Ex_b_1[k]*( K_Ex_c_1[i]*Gx_7[i][ja][k] + K_Ex_d_1[i]*Gx_r);
}
ja++;
jj--;
}
}
}
void CFDTD_3D_PML::calc_Ex_8()
{
double Fx_r, Gx_r;
long i, j, k, jj, ia, ja;
//8
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_MIN_1; j++)
{
for (k = 1; k < nPML_z_1; k++)
{
Fx_r = Fx_8[ia][ja][k];
Fx_8[ia][ja][k] = K_a[Ind[i][j][k]]*Fx_8[ia][ja][k] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
Gx_r = Gx_8[ia][ja][k];
Gx_8[ia][ja][k] = K_Gx_a_2[jj]*Gx_8[ia][ja][k] +
K_Gx_b_2[jj]*(Fx_8[ia][ja][k] - Fx_r);
Ex[i][j][k] = K_Ex_a_1[k]*Ex[i][j][k] +
K_Ex_b_1[k]*TwoOREp0*( Gx_8[ia][ja][k] - Gx_r);
}
ja++;
jj--;
}
ia++;
}
}
void CFDTD_3D_PML::calc_Ex_9()
{
double Fx_r, Gx_r;
long i, j, k, ii, jj, ia, ja;
//9
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_MIN_1; j++)
{
for (k = 1; k < nPML_z_1; k++)
{
Fx_r = Fx_9[ia][ja][k];
Fx_9[ia][ja][k] = K_a[Ind[i][j][k]]*Fx_9[ia][ja][k] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
Gx_r = Gx_9[ia][ja][k];
Gx_9[ia][ja][k] = K_Gx_a_2[jj]*Gx_9[ia][ja][k] +
K_Gx_b_2[jj]*(Fx_9[ia][ja][k] - Fx_r);
Ex[i][j][k] = K_Ex_a_1[k]*Ex[i][j][k] +
K_Ex_b_1[k]*(K_Ex_c_2[ii]*Gx_9[ia][ja][k] + K_Ex_d_2[ii]*Gx_r);
}
ja++;
jj--;
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Ex_10()
{
double Fx_r, Gx_r;
long i, j, k, ka;
//10
for (i = 0; i < nPML_x_1; i++)
{
for (j = 1; j < nPML_y_1; j++)
{
ka = 0;
for (k = nPML_z_1; k < nz_MIN_nPML_z_2; k++)
{
Fx_r = Fx_10[i][j][ka];
Fx_10[i][j][ka] = K_a[Ind[i][j][k]]*Fx_10[i][j][ka] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
Gx_r = Gx_10[i][j][ka];
Gx_10[i][j][ka] = K_Gx_a_1[j]*Gx_10[i][j][ka] +
K_Gx_b_1[j]*(Fx_10[i][j][ka] - Fx_r);
Ex[i][j][k] += inv_TwoOREp0*( K_Ex_c_1[i]*Gx_10[i][j][ka] +
K_Ex_d_1[i]*Gx_r);
ka++;
}
}
}
}
void CFDTD_3D_PML::calc_Ex_11()
{
double Fx_r;
long i, j, k, ia, ka;
//11
ia = 0;
for (i = nPML_x_1; i < nx_MIN_1_MIN_nPML_x_2; i++)
{
for (j = 1; j < nPML_y_1; j++)
{
ka = 0;
for (k = nPML_z_1; k < nz_MIN_nPML_z_2; k++)
{
Fx_r = Fx_11[ia][j][ka];
Fx_11[ia][j][ka] = K_a[Ind[i][j][k]]*Fx_11[ia][j][ka] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
Ex[i][j][k] = K_Gx_a_1[j]*Ex[i][j][k] +
K_Gx_b_1[j]*(Fx_11[ia][j][ka] - Fx_r);
ka++;
}
}
ia++;
}
}
void CFDTD_3D_PML::calc_Ex_12()
{
double Fx_r, Gx_r;
long i, j, k, ii, ia, ka;
//12
ia = 0;
ii = nPML_x_2_MIN_1;
for (i = nx_MIN_1_MIN_nPML_x_2; i < nx_MIN_1; i++)
{
for (j = 1; j < nPML_y_1; j++)
{
ka = 0;
for (k = nPML_z_1; k < nz_MIN_nPML_z_2; k++)
{
Fx_r = Fx_12[ia][j][ka];
Fx_12[ia][j][ka] = K_a[Ind[i][j][k]]*Fx_12[ia][j][ka] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
Gx_r = Gx_12[ia][j][ka];
Gx_12[ia][j][ka] = K_Gx_a_1[j]*Gx_12[ia][j][ka] +
K_Gx_b_1[j]*(Fx_12[ia][j][ka] - Fx_r);
Ex[i][j][k] += inv_TwoOREp0*( K_Ex_c_2[ii]*Gx_12[ia][j][ka] +
K_Ex_d_2[ii]*Gx_r);
ka++;
}
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Ex_13()
{
double Fx_r;
long i, j, k, ja, ka;
//13
for (i = 0; i < nPML_x_1; i++)
{
ja = 0;
for (j = nPML_y_1; j < ny_MIN_nPML_y_2; j++)
{
ka = 0;
for (k = nPML_z_1; k < nz_MIN_nPML_z_2; k++)
{
Fx_r = Fx_13[i][ja][ka];
Fx_13[i][ja][ka] = K_a[Ind[i][j][k]]*Fx_13[i][ja][ka] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
Ex[i][j][k] += inv_TwoOREp0*( K_Ex_c_1[i]*Fx_13[i][ja][ka] +
K_Ex_d_1[i]*Fx_r);
ka++;
}
ja++;
}
}
}
//Outside of the PML zone
void CFDTD_3D_PML::calc_Ex_14(long nx_a, long nx_b, long ny_a, long ny_b, long nz_a, long nz_b)
{
long i, j, k;
//14
for (i = nx_a; i < nx_b; i++)
{
for (j = ny_a; j < ny_b; j++)
{
for (k = nz_a; k < nz_b; k++)
{
Ex[i][j][k] = K_a[Ind[i][j][k]]*Ex[i][j][k] +
K_b[Ind[i][j][k]]*( (Hz[i][j][k] - Hz[i][j-1][k])*inv_dy -
(Hy[i][j][k] - Hy[i][j][k-1])*inv_dz );
}
}
}
}
void CFDTD_3D_PML::calc_Ex_15()
{
double Fx_r;
long i, j, k, ii, ia, ja, ka;
//15
ia = 0;
ii = nPML_x_2_MIN_1;
for (i = nx_MIN_1_MIN_nPML_x_2; i < nx_MIN_1; i++)
{
ja = 0;
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -