📄 fdtd_3d_pml_hx.cpp
字号:
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_15[ia][ja][ka];
Bx_15[ia][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_2[ii]*Bx_15[ia][ja][ka] +
K_Hx_d_2[ii]*Bx_r )/mu_r[Ind[i][j][k]];
ka++;
}
ja++;
}
ia++;
ii--;
}
}
void calc_Hx_16(double ***Hx, double ***Ey, double ***Ez, double ***Bx_16, long ***Ind,
double *mu_r, double *K_Bx_a_2, double *K_Bx_b_2, double *K_Hx_c_1,
double *K_Hx_d_1, double inv_dy, double inv_dz, long nPML_x_1,
long ny_MIN_1_MIN_nPML_y_2, long ny_MIN_1, long nPML_z_1,
long nz_MIN_1_MIN_nPML_z_2, long nPML_y_2_MIN_1, double inv_TwoOREp0)
{
long i, j, k, jj, ja, ka;
double Bx_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_1_MIN_nPML_z_2; k++)
{
Bx_r = Bx_16[i][ja][ka];
Bx_16[i][ja][ka] = K_Bx_a_2[jj]*Bx_16[i][ja][ka] +
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] += inv_TwoOREp0*( K_Hx_c_1[i]*Bx_16[i][ja][ka] +
K_Hx_d_1[i]*Bx_r )/mu_r[Ind[i][j][k]];
ka++;
}
ja++;
jj--;
}
}
}
void calc_Hx_17(double ***Hx, double ***Ey, double ***Ez, long ***Ind, double *mu_r,
double *K_Bx_a_2, double *K_Bx_b_2, double inv_dyORmu_0, double inv_dzORmu_0,
long nPML_x_1, long nx_MIN_nPML_x_2, long ny_MIN_1_MIN_nPML_y_2, long ny_MIN_1,
long nPML_z_1, long nz_MIN_1_MIN_nPML_z_2, long nPML_y_2_MIN_1)
{
long i, j, k, jj;
//17
for (i = nPML_x_1; i < nx_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_1_MIN_nPML_z_2; k++)
{
Hx[i][j][k] = K_Bx_a_2[jj]*Hx[i][j][k] + K_Bx_b_2[jj]/mu_r[Ind[i][j][k]]*
( (Ey[i][j][k+1] - Ey[i][j][k])*inv_dzORmu_0 -
(Ez[i][j+1][k] - Ez[i][j][k])*inv_dyORmu_0);
}
jj--;
}
}
}
void calc_Hx_18(double ***Hx, double ***Ey, double ***Ez, double ***Bx_18, long ***Ind,
double *mu_r, double *K_Bx_a_2, double *K_Bx_b_2, double *K_Hx_c_2,
double *K_Hx_d_2, double inv_dy, double inv_dz, long nx_MIN_nPML_x_2, long nx,
long ny_MIN_1_MIN_nPML_y_2, long ny_MIN_1, long nPML_z_1,
long nz_MIN_1_MIN_nPML_z_2, long nPML_x_2_MIN_1, long nPML_y_2_MIN_1,
double inv_TwoOREp0)
{
long i, j, k, ii, jj, ia, ja, ka;
double Bx_r;
//18
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++)
{
ka = 0;
for (k = nPML_z_1; k < nz_MIN_1_MIN_nPML_z_2; k++)
{
Bx_r = Bx_18[ia][ja][ka];
Bx_18[ia][ja][ka] = K_Bx_a_2[jj]*Bx_18[ia][ja][ka] +
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] += inv_TwoOREp0*( K_Hx_c_2[ii]*Bx_18[ia][ja][ka] +
K_Hx_d_2[ii]*Bx_r )/mu_r[Ind[i][j][k]];
ka++;
}
ja++;
jj--;
}
ia++;
ii--;
}
}
void calc_Hx_19(double ***Hx, double ***Ey, double ***Ez, double ***Bx_19, long ***Ind,
double *mu_r, double *K_Bx_a_1, double *K_Bx_b_1, double *K_Hx_a_2,
double *K_Hx_b_2, double *K_Hx_c_1, double *K_Hx_d_1, double inv_dy,
double inv_dz, long nPML_x_1, long nPML_y_1, long nz_MIN_1_MIN_nPML_z_2,
long nz_MIN_1, long nPML_z_2_MIN_1)
{
long i, j, k, kk, ka;
double Bx_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++)
{
Bx_r = Bx_19[i][j][ka];
Bx_19[i][j][ka] = K_Bx_a_1[j]*Bx_19[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] = K_Hx_a_2[kk]*Hx[i][j][k] +
K_Hx_b_2[kk]*( K_Hx_c_1[i]*Bx_19[i][j][ka] +
K_Hx_d_1[i]*Bx_r )/mu_r[Ind[i][j][k]];
ka++;
kk--;
}
}
}
}
void calc_Hx_20(double ***Hx, double ***Ey, double ***Ez, double ***Bx_20, long ***Ind,
double *mu_r, double *K_Bx_a_1, double *K_Bx_b_1, double *K_Hx_a_2,
double *K_Hx_b_2, double inv_dy, double inv_dz, long nPML_x_1,
long nx_MIN_nPML_x_2, long nPML_y_1, long nz_MIN_1_MIN_nPML_z_2,
long nz_MIN_1, long nPML_z_2_MIN_1, double TwoOREp0DIVMu0)
{
long i, j, k, kk, ia, ka;
double Bx_r;
//20
ia = 0;
for (i = nPML_x_1; i < nx_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++)
{
Bx_r = Bx_20[ia][j][ka];
Bx_20[ia][j][ka] = K_Bx_a_1[j]*Bx_20[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] = K_Hx_a_2[kk]*Hx[i][j][k] + K_Hx_b_2[kk]*TwoOREp0DIVMu0*
( Bx_20[ia][j][ka] - Bx_r )/mu_r[Ind[i][j][k]];
ka++;
kk--;
}
}
ia++;
}
}
void calc_Hx_21(double ***Hx, double ***Ey, double ***Ez, double ***Bx_21, long ***Ind,
double *mu_r, double *K_Bx_a_1, double *K_Bx_b_1, double *K_Hx_a_2,
double *K_Hx_b_2, double *K_Hx_c_2, double *K_Hx_d_2, double inv_dy,
double inv_dz, long nx_MIN_nPML_x_2, long nx, long nPML_y_1,
long nz_MIN_1_MIN_nPML_z_2, long nz_MIN_1, long nPML_x_2_MIN_1,
long nPML_z_2_MIN_1)
{
long i, j, k, ii, kk, ia, ka;
double Bx_r;
//21
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;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
Bx_r = Bx_21[ia][j][ka];
Bx_21[ia][j][ka] = K_Bx_a_1[j]*Bx_21[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] = K_Hx_a_2[kk]*Hx[i][j][k] +
K_Hx_b_2[kk]*( K_Hx_c_2[ii]*Bx_21[ia][j][ka] +
K_Hx_d_2[ii]*Bx_r )/mu_r[Ind[i][j][k]];
ka++;
kk--;
}
}
ia++;
ii--;
}
}
void calc_Hx_22(double ***Hx, double ***Ey, double ***Ez, double ***Bx_22, long ***Ind,
double *mu_r, double *K_Hx_a_2, double *K_Hx_b_2, double *K_Hx_c_1,
double *K_Hx_d_1, double dtDIVdy, double dtDIVdz, long nPML_x_1,
long nPML_y_1, long ny_MIN_1_MIN_nPML_y_2, long nz_MIN_1_MIN_nPML_z_2,
long nz_MIN_1, long nPML_z_2_MIN_1)
{
long i, j, k, kk, ja, ka;
double Bx_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_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
Bx_r = Bx_22[i][ja][ka];
Bx_22[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] = K_Hx_a_2[kk]*Hx[i][j][k] +
K_Hx_b_2[kk]*( K_Hx_c_1[i]*Bx_22[i][ja][ka] +
K_Hx_d_1[i]*Bx_r )/mu_r[Ind[i][j][k]];
ka++;
kk--;
}
ja++;
}
}
}
void calc_Hx_23(double ***Hx, double ***Ey, double ***Ez, long ***Ind, double *mu_r,
double *K_Hx_a_2, double *K_Hx_b_2, double TwoOREp0DIVMu0ORdtDIVdy,
double TwoOREp0DIVMu0ORdtDIVdz, long nPML_x_1, long nx_MIN_nPML_x_2,
long nPML_y_1, long ny_MIN_1_MIN_nPML_y_2, long nz_MIN_1_MIN_nPML_z_2,
long nz_MIN_1, long nPML_z_2_MIN_1)
{
long i, j, k, kk;
//23
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++)
{
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
Hx[i][j][k] = K_Hx_a_2[kk]*Hx[i][j][k] + K_Hx_b_2[kk]/mu_r[Ind[i][j][k]]*
( (Ey[i][j][k+1] - Ey[i][j][k])*TwoOREp0DIVMu0ORdtDIVdz -
(Ez[i][j+1][k] - Ez[i][j][k])*TwoOREp0DIVMu0ORdtDIVdy );
kk--;
}
}
}
}
void calc_Hx_24(double ***Hx, double ***Ey, double ***Ez, double ***Bx_24, long ***Ind,
double *mu_r, double *K_Hx_a_2, double *K_Hx_b_2, double *K_Hx_c_2,
double *K_Hx_d_2, double dtDIVdy, double dtDIVdz, long nx_MIN_nPML_x_2,
long nx, long nPML_y_1, long ny_MIN_1_MIN_nPML_y_2, long nz_MIN_1_MIN_nPML_z_2,
long nz_MIN_1, long nPML_x_2_MIN_1, long nPML_z_2_MIN_1)
{
long i, j, k, ii, kk, ia, ja, ka;
double Bx_r;
//24
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++)
{
ka = 0;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
Bx_r = Bx_24[ia][ja][ka];
Bx_24[ia][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] = K_Hx_a_2[kk]*Hx[i][j][k] +
K_Hx_b_2[kk]*( K_Hx_c_2[ii]*Bx_24[ia][ja][ka] +
K_Hx_d_2[ii]*Bx_r )/mu_r[Ind[i][j][k]];
ka++;
kk--;
}
ja++;
}
ia++;
ii--;
}
}
void calc_Hx_25(double ***Hx, double ***Ey, double ***Ez, double ***Bx_25, long ***Ind,
double *mu_r, double *K_Bx_a_2, double *K_Bx_b_2, double *K_Hx_a_2,
double *K_Hx_b_2, double *K_Hx_c_1, double *K_Hx_d_1, double inv_dy,
double inv_dz, long nPML_x_1, long ny_MIN_1_MIN_nPML_y_2, long ny_MIN_1,
long nz_MIN_1_MIN_nPML_z_2, long nz_MIN_1, long nPML_y_2_MIN_1,
long nPML_z_2_MIN_1)
{
long i, j, k, jj, kk, ja, ka;
double Bx_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_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
Bx_r = Bx_25[i][ja][ka];
Bx_25[i][ja][ka] = K_Bx_a_2[jj]*Bx_25[i][ja][ka] +
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_2[kk]*Hx[i][j][k] +
K_Hx_b_2[kk]*( K_Hx_c_1[i]*Bx_25[i][ja][ka] +
K_Hx_d_1[i]*Bx_r )/mu_r[Ind[i][j][k]];
ka++;
kk--;
}
ja++;
jj--;
}
}
}
void calc_Hx_26(double ***Hx, double ***Ey, double ***Ez, double ***Bx_26, long ***Ind,
double *mu_r, double *K_Bx_a_2, double *K_Bx_b_2, double *K_Hx_a_2,
double *K_Hx_b_2, double inv_dy, double inv_dz, long nPML_x_1,
long nx_MIN_nPML_x_2, long ny_MIN_1_MIN_nPML_y_2, long ny_MIN_1,
long nz_MIN_1_MIN_nPML_z_2, long nz_MIN_1, long nPML_y_2_MIN_1,
long nPML_z_2_MIN_1, double TwoOREp0DIVMu0)
{
long i, j, k, jj, kk, ia, ja, ka;
double Bx_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_1_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++)
{
Bx_r = Bx_26[ia][ja][ka];
Bx_26[ia][ja][ka] = K_Bx_a_2[jj]*Bx_26[ia][ja][ka] +
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_2[kk]*Hx[i][j][k] + K_Hx_b_2[kk]*TwoOREp0DIVMu0*
( Bx_26[ia][ja][ka] - Bx_r )/mu_r[Ind[i][j][k]];
ka++;
kk--;
}
ja++;
jj--;
}
ia++;
}
}
void calc_Hx_27(double ***Hx, double ***Ey, double ***Ez, double ***Bx_27, long ***Ind,
double *mu_r, double *K_Bx_a_2, double *K_Bx_b_2, double *K_Hx_a_2,
double *K_Hx_b_2, double *K_Hx_c_2, double *K_Hx_d_2, double inv_dy,
double inv_dz, long nx_MIN_nPML_x_2, long nx, long ny_MIN_1_MIN_nPML_y_2,
long ny_MIN_1, long nz_MIN_1_MIN_nPML_z_2, long nz_MIN_1, long nPML_x_2_MIN_1,
long nPML_y_2_MIN_1, long nPML_z_2_MIN_1)
{
long i, j, k, ii, jj, kk, ia, ja, ka;
double Bx_r;
//27
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++)
{
ka = 0;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_1_MIN_nPML_z_2; k < nz_MIN_1; k++)
{
Bx_r = Bx_27[ia][ja][ka];
Bx_27[ia][ja][ka] = K_Bx_a_2[jj]*Bx_27[ia][ja][ka] +
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_2[kk]*Hx[i][j][k] +
K_Hx_b_2[kk]*( K_Hx_c_2[ii]*Bx_27[ia][ja][ka] +
K_Hx_d_2[ii]*Bx_r )/mu_r[Ind[i][j][k]];
ka++;
kk--;
}
ja++;
jj--;
}
ia++;
ii--;
}
}⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -