📄 fdtd_3d_pml_ez.cpp
字号:
for (k = nPML_z_1; k < nz_MIN_1_MIN_nPML_z_2; k++)
{
Fz_r = Fz_15[ia][ja][ka];
Fz_15[ia][ja][ka] = K_a[Ind[i][j][k]]*Fz_15[ia][ja][ka] +
K_b[Ind[i][j][k]]*( (Hy[i][j][k] - Hy[i-1][j][k])*inv_dx -
(Hx[i][j][k] - Hx[i][j-1][k])*inv_dy );
Ez[i][j][k] = K_Gz_a_2[ii]*Ez[i][j][k] +
K_Gz_b_2[ii]*(Fz_15[ia][ja][ka] - Fz_r);
ka++;
}
ja++;
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Ez_16()
{
long i, j, k, jj, ja, ka;
double Fz_r, Gz_r;
//16
for (i = 1; 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++)
{
ka = 0;
for (k = nPML_z_1; k < nz_MIN_1_MIN_nPML_z_2; k++)
{
Fz_r = Fz_16[i][ja][ka];
Fz_16[i][ja][ka] = K_a[Ind[i][j][k]]*Fz_16[i][ja][ka] +
K_b[Ind[i][j][k]]*( (Hy[i][j][k] - Hy[i-1][j][k])*inv_dx -
(Hx[i][j][k] - Hx[i][j-1][k])*inv_dy );
Gz_r = Gz_16[i][ja][ka];
Gz_16[i][ja][ka] = K_Gz_a_1[i]*Gz_16[i][ja][ka] +
K_Gz_b_1[i]*(Fz_16[i][ja][ka] - Fz_r);
Ez[i][j][k] = K_Ez_a_2[jj]*Ez[i][j][k] +
K_Ez_b_2[jj]*TwoOREp0*( Gz_16[i][ja][ka] - Gz_r);
ka++;
}
ja++;
jj--;
}
}
}
void CFDTD_3D_PML::calc_Ez_17()
{
long i, j, k, jj, ia, ja, ka;
double Fz_r;
//17
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_nPML_y_2; j < ny_MIN_1; j++)
{
ka = 0;
for (k = nPML_z_1; k < nz_MIN_1_MIN_nPML_z_2; k++)
{
Fz_r = Fz_17[ia][ja][ka];
Fz_17[ia][ja][ka] = K_a[Ind[i][j][k]]*Fz_17[ia][ja][ka] +
K_b[Ind[i][j][k]]*( (Hy[i][j][k] - Hy[i-1][j][k])*inv_dx -
(Hx[i][j][k] - Hx[i][j-1][k])*inv_dy );
Ez[i][j][k] = K_Ez_a_2[jj]*Ez[i][j][k] +
K_Ez_b_2[jj]*TwoOREp0*( Fz_17[ia][ja][ka] - Fz_r);
ka++;
}
ja++;
jj--;
}
ia++;
}
}
void CFDTD_3D_PML::calc_Ez_18()
{
long i, j, k, ii, jj, ia, ja, ka;
double Fz_r, Gz_r;
//18
ia = 0;
ii = nPML_x_2_MIN_1;
for (i = nx_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++)
{
ka = 0;
for (k = nPML_z_1; k < nz_MIN_1_MIN_nPML_z_2; k++)
{
Fz_r = Fz_18[ia][ja][ka];
Fz_18[ia][ja][ka] = K_a[Ind[i][j][k]]*Fz_18[ia][ja][ka] +
K_b[Ind[i][j][k]]*( (Hy[i][j][k] - Hy[i-1][j][k])*inv_dx -
(Hx[i][j][k] - Hx[i][j-1][k])*inv_dy );
Gz_r = Gz_18[ia][ja][ka];
Gz_18[ia][ja][ka] = K_Gz_a_2[ii]*Gz_18[ia][ja][ka] +
K_Gz_b_2[ii]*(Fz_18[ia][ja][ka] - Fz_r);
Ez[i][j][k] = K_Ez_a_2[jj]*Ez[i][j][k] +
K_Ez_b_2[jj]*TwoOREp0*( Gz_18[ia][ja][ka] - Gz_r);
ka++;
}
ja++;
jj--;
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Ez_19()
{
long i, j, k, kk, ka;
double Fz_r, Gz_r;
//19
for (i = 1; i < nPML_x_1; i++)
{
for (j = 1; 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++)
{
Fz_r = Fz_19[i][j][ka];
Fz_19[i][j][ka] = K_a[Ind[i][j][k]]*Fz_19[i][j][ka] +
K_b[Ind[i][j][k]]*( (Hy[i][j][k] - Hy[i-1][j][k])*inv_dx -
(Hx[i][j][k] - Hx[i][j-1][k])*inv_dy );
Gz_r = Gz_19[i][j][ka];
Gz_19[i][j][ka] = K_Gz_a_1[i]*Gz_19[i][j][ka] +
K_Gz_b_1[i]*(Fz_19[i][j][ka] - Fz_r);
Ez[i][j][k] = K_Ez_a_1[j]*Ez[i][j][k] +
K_Ez_b_1[j]*( K_Ez_c_2[kk]*Gz_19[i][j][ka] + K_Ez_d_2[kk]*Gz_r);
ka++;
kk--;
}
}
}
}
void CFDTD_3D_PML::calc_Ez_20()
{
long i, j, k, kk, ia, ka;
double Fz_r;
//20
ia = 0;
for (i = nPML_x_1; i < nx_MIN_nPML_x_2; i++)
{
for (j = 1; 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++)
{
Fz_r = Fz_20[ia][j][ka];
Fz_20[ia][j][ka] = K_a[Ind[i][j][k]]*Fz_20[ia][j][ka] +
K_b[Ind[i][j][k]]*( (Hy[i][j][k] - Hy[i-1][j][k])*inv_dx -
(Hx[i][j][k] - Hx[i][j-1][k])*inv_dy );
Ez[i][j][k] = K_Ez_a_1[j]*Ez[i][j][k] +
K_Ez_b_1[j]*( K_Ez_c_2[kk]*Fz_20[ia][j][ka] + K_Ez_d_2[kk]*Fz_r);
ka++;
kk--;
}
}
ia++;
}
}
void CFDTD_3D_PML::calc_Ez_21()
{
long i, j, k, ii, kk, ia, ka;
double Fz_r, Gz_r;
//21
ia = 0;
ii = nPML_x_2_MIN_1;
for (i = nx_MIN_nPML_x_2; i < nx_MIN_1; i++)
{
for (j = 1; 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++)
{
Fz_r = Fz_21[ia][j][ka];
Fz_21[ia][j][ka] = K_a[Ind[i][j][k]]*Fz_21[ia][j][ka] +
K_b[Ind[i][j][k]]*( (Hy[i][j][k] - Hy[i-1][j][k])*inv_dx -
(Hx[i][j][k] - Hx[i][j-1][k])*inv_dy );
Gz_r = Gz_21[ia][j][ka];
Gz_21[ia][j][ka] = K_Gz_a_2[ii]*Gz_21[ia][j][ka] +
K_Gz_b_2[ii]*(Fz_21[ia][j][ka] - Fz_r);
Ez[i][j][k] = K_Ez_a_1[j]*Ez[i][j][k] +
K_Ez_b_1[j]*( K_Ez_c_2[kk]*Gz_21[ia][j][ka] + K_Ez_d_2[kk]*Gz_r);
ka++;
kk--;
}
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Ez_22()
{
long i, j, k, kk, ja, ka;
double Fz_r, Gz_r;
//22
for (i = 1; 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++)
{
Fz_r = Fz_22[i][ja][ka];
Fz_22[i][ja][ka] = K_a[Ind[i][j][k]]*Fz_22[i][ja][ka] +
K_b[Ind[i][j][k]]*( (Hy[i][j][k] - Hy[i-1][j][k])*inv_dx -
(Hx[i][j][k] - Hx[i][j-1][k])*inv_dy );
Gz_r = Gz_22[i][ja][ka];
Gz_22[i][ja][ka] = K_Gz_a_1[i]*Gz_22[i][ja][ka] +
K_Gz_b_1[i]*(Fz_22[i][ja][ka] - Fz_r);
Ez[i][j][k] += inv_TwoOREp0*( K_Ez_c_2[kk]*Gz_22[i][ja][ka] + K_Ez_d_2[kk]*Gz_r);
ka++;
kk--;
}
ja++;
}
}
}
void CFDTD_3D_PML::calc_Ez_23()
{
long i, j, k, kk, ia, ja, ka;
double Fz_r;
//23
ia = 0;
for (i = nPML_x_1; i < nx_MIN_nPML_x_2; 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++)
{
Fz_r = Fz_23[ia][ja][ka];
Fz_23[ia][ja][ka] = K_a[Ind[i][j][k]]*Fz_23[ia][ja][ka] +
K_b[Ind[i][j][k]]*( (Hy[i][j][k] - Hy[i-1][j][k])*inv_dx -
(Hx[i][j][k] - Hx[i][j-1][k])*inv_dy );
Ez[i][j][k] += inv_TwoOREp0*( K_Ez_c_2[kk]*Fz_23[ia][ja][ka] + K_Ez_d_2[kk]*Fz_r);
ka++;
kk--;
}
ja++;
}
ia++;
}
}
void CFDTD_3D_PML::calc_Ez_24()
{
long i, j, k, ii, kk, ia, ja, ka;
double Fz_r, Gz_r;
//24
ia = 0;
ii = nPML_x_2_MIN_1;
for (i = nx_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++)
{
Fz_r = Fz_24[ia][ja][ka];
Fz_24[ia][ja][ka] = K_a[Ind[i][j][k]]*Fz_24[ia][ja][ka] +
K_b[Ind[i][j][k]]*( (Hy[i][j][k] - Hy[i-1][j][k])*inv_dx -
(Hx[i][j][k] - Hx[i][j-1][k])*inv_dy );
Gz_r = Gz_24[ia][ja][ka];
Gz_24[ia][ja][ka] = K_Gz_a_2[ii]*Gz_24[ia][ja][ka] +
K_Gz_b_2[ii]*(Fz_24[ia][ja][ka] - Fz_r);
Ez[i][j][k] += inv_TwoOREp0*( K_Ez_c_2[kk]*Gz_24[ia][ja][ka] + K_Ez_d_2[kk]*Gz_r);
ka++;
kk--;
}
ja++;
}
ia++;
ii--;
}
}
void CFDTD_3D_PML::calc_Ez_25()
{
long i, j, k, jj, kk, ja, ka;
double Fz_r, Gz_r;
//25
for (i = 1; 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++)
{
ka = 0;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
Fz_r = Fz_25[i][ja][ka];
Fz_25[i][ja][ka] = K_a[Ind[i][j][k]]*Fz_25[i][ja][ka] +
K_b[Ind[i][j][k]]*( (Hy[i][j][k] - Hy[i-1][j][k])*inv_dx -
(Hx[i][j][k] - Hx[i][j-1][k])*inv_dy );
Gz_r = Gz_25[i][ja][ka];
Gz_25[i][ja][ka] = K_Gz_a_1[i]*Gz_25[i][ja][ka] +
K_Gz_b_1[i]*(Fz_25[i][ja][ka] - Fz_r);
Ez[i][j][k] = K_Ez_a_2[jj]*Ez[i][j][k] +
K_Ez_b_2[jj]*( K_Ez_c_2[kk]*Gz_25[i][ja][ka] + K_Ez_d_2[kk]*Gz_r);
ka++;
kk--;
}
ja++;
jj--;
}
}
}
void CFDTD_3D_PML::calc_Ez_26()
{
long i, j, k, jj, kk, ia, ja, ka;
double Fz_r;
//26
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_nPML_y_2; j < ny_MIN_1; j++)
{
ka = 0;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
Fz_r = Fz_26[ia][ja][ka];
Fz_26[ia][ja][ka] = K_a[Ind[i][j][k]]*Fz_26[ia][ja][ka] +
K_b[Ind[i][j][k]]*( (Hy[i][j][k] - Hy[i-1][j][k])*inv_dx -
(Hx[i][j][k] - Hx[i][j-1][k])*inv_dy );
Ez[i][j][k] = K_Ez_a_2[jj]*Ez[i][j][k] +
K_Ez_b_2[jj]*( K_Ez_c_2[kk]*Fz_26[ia][ja][ka] + K_Ez_d_2[kk]*Fz_r);
ka++;
kk--;
}
ja++;
jj--;
}
ia++;
}
}
void CFDTD_3D_PML::calc_Ez_27()
{
long i, j, k, ii, jj, kk, ia, ja, ka;
double Fz_r, Gz_r;
//27
ia = 0;
ii = nPML_x_2_MIN_1;
for (i = nx_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++)
{
ka = 0;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
Fz_r = Fz_27[ia][ja][ka];
Fz_27[ia][ja][ka] = K_a[Ind[i][j][k]]*Fz_27[ia][ja][ka] +
K_b[Ind[i][j][k]]*( (Hy[i][j][k] - Hy[i-1][j][k])*inv_dx -
(Hx[i][j][k] - Hx[i][j-1][k])*inv_dy );
Gz_r = Gz_27[ia][ja][ka];
Gz_27[ia][ja][ka] = K_Gz_a_2[ii]*Gz_27[ia][ja][ka] +
K_Gz_b_2[ii]*(Fz_27[ia][ja][ka] - Fz_r);
Ez[i][j][k] = K_Ez_a_2[jj]*Ez[i][j][k] +
K_Ez_b_2[jj]*( K_Ez_c_2[kk]*Gz_27[ia][ja][ka] + K_Ez_d_2[kk]*Gz_r);
ka++;
kk--;
}
ja++;
jj--;
}
ia++;
ii--;
}
}⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -