📄 fdtd_3d_pml_ey.cpp
字号:
double *K_Gy_b_2, double *K_Ey_c_1, double *K_Ey_d_1, double inv_dx,
double inv_dz, long nPML_x_1, long nx_MIN_nPML_x_2, long nPML_y_1,
long nz_MIN_nPML_z_2, long nz_MIN_1, long nPML_z_2_MIN_1,
double inv_TwoOREp0)
{
long i, j, k, kk, ia, ka;
double Fy_r, Gy_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_nPML_z_2; k < nz_MIN_1 ; k++)
{
Fy_r = Fy_20[ia][j][ka];
Fy_20[ia][j][ka] = K_a[Ind[i][j][k]]*Fy_20[ia][j][ka] +
K_b[Ind[i][j][k]]*( (Hx[i][j][k] - Hx[i][j][k-1])*inv_dz -
(Hz[i][j][k] - Hz[i-1][j][k])*inv_dx );
Gy_r = Gy_20[ia][j][ka];
Gy_20[ia][j][ka] = K_Gy_a_2[kk]*Gy_20[ia][j][ka] +
K_Gy_b_2[kk]*(Fy_20[ia][j][ka] - Fy_r);
Ey[i][j][k] += inv_TwoOREp0*( K_Ey_c_1[j]*Gy_20[ia][j][ka] +
K_Ey_d_1[j]*Gy_r);
ka++;
kk--;
}
}
ia++;
}
}
void calc_Ey_21(double ***Ey, double ***Hx, double ***Hz, double ***Fy_21,
double ***Gy_21, long ***Ind, double *K_a, double *K_b, double *K_Gy_a_2,
double *K_Gy_b_2, double *K_Ey_a_2, double *K_Ey_b_2, double *K_Ey_c_1,
double *K_Ey_d_1, double inv_dx, double inv_dz, long nx_MIN_nPML_x_2,
long nx_MIN_1, long nPML_y_1, long nz_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 Fy_r, Gy_r;
//21
ia = 0;
ii = nPML_x_2_MIN_1;
for (i = nx_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_nPML_z_2; k < nz_MIN_1 ; k++)
{
Fy_r = Fy_21[ia][j][ka];
Fy_21[ia][j][ka] = K_a[Ind[i][j][k]]*Fy_21[ia][j][ka] +
K_b[Ind[i][j][k]]*( (Hx[i][j][k] - Hx[i][j][k-1])*inv_dz -
(Hz[i][j][k] - Hz[i-1][j][k])*inv_dx );
Gy_r = Gy_21[ia][j][ka];
Gy_21[ia][j][ka] = K_Gy_a_2[kk]*Gy_21[ia][j][ka] +
K_Gy_b_2[kk]*(Fy_21[ia][j][ka] - Fy_r);
Ey[i][j][k] = K_Ey_a_2[ii]*Ey[i][j][k] +
K_Ey_b_2[ii]*( K_Ey_c_1[j]*Gy_21[ia][j][ka] + K_Ey_d_1[j]*Gy_r);
ka++;
kk--;
}
}
ia++;
ii--;
}
}
void calc_Ey_22(double ***Ey, double ***Hx, double ***Hz, double ***Fy_22,
double ***Gy_22, long ***Ind, double *K_a, double *K_b, double *K_Gy_a_2,
double *K_Gy_b_2, double *K_Ey_a_1, double *K_Ey_b_1, double inv_dx,
double inv_dz, long nPML_x_1, long nPML_y_1, long ny_MIN_1_MIN_nPML_y_2,
long nz_MIN_nPML_z_2, long nz_MIN_1, long nPML_z_2_MIN_1, double TwoOREp0)
{
long i, j, k, kk, ja, ka;
double Fy_r, Gy_r;
//22
for (i = 1; 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_nPML_z_2; k < nz_MIN_1 ; k++)
{
Fy_r = Fy_22[i][ja][ka];
Fy_22[i][ja][ka] = K_a[Ind[i][j][k]]*Fy_22[i][ja][ka] +
K_b[Ind[i][j][k]]*( (Hx[i][j][k] - Hx[i][j][k-1])*inv_dz -
(Hz[i][j][k] - Hz[i-1][j][k])*inv_dx );
Gy_r = Gy_22[i][ja][ka];
Gy_22[i][ja][ka] = K_Gy_a_2[kk]*Gy_22[i][ja][ka] +
K_Gy_b_2[kk]*(Fy_22[i][ja][ka] - Fy_r);
Ey[i][j][k] = K_Ey_a_1[i]*Ey[i][j][k] +
K_Ey_b_1[i]*TwoOREp0*( Gy_22[i][ja][ka] - Gy_r);
ka++;
kk--;
}
ja++;
}
}
}
void calc_Ey_23(double ***Ey, double ***Hx, double ***Hz, double ***Fy_23, long ***Ind,
double *K_a, double *K_b, double *K_Gy_a_2, double *K_Gy_b_2,
double inv_dx, double inv_dz, 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_nPML_z_2,
long nz_MIN_1, long nPML_z_2_MIN_1)
{
long i, j, k, kk, ia, ja, ka;
double Fy_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_1_MIN_nPML_y_2; j++)
{
ka = 0;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_nPML_z_2; k < nz_MIN_1 ; k++)
{
Fy_r = Fy_23[ia][ja][ka];
Fy_23[ia][ja][ka] = K_a[Ind[i][j][k]]*Fy_23[ia][ja][ka] +
K_b[Ind[i][j][k]]*( (Hx[i][j][k] - Hx[i][j][k-1])*inv_dz -
(Hz[i][j][k] - Hz[i-1][j][k])*inv_dx );
Ey[i][j][k] = K_Gy_a_2[kk]*Ey[i][j][k] +
K_Gy_b_2[kk]*(Fy_23[ia][ja][ka] - Fy_r);
ka++;
kk--;
}
ja++;
}
ia++;
}
}
void calc_Ey_24(double ***Ey, double ***Hx, double ***Hz, double ***Fy_24,
double ***Gy_24, long ***Ind, double *K_a, double *K_b, double *K_Gy_a_2,
double *K_Gy_b_2, double *K_Ey_a_2, double *K_Ey_b_2, double inv_dx,
double inv_dz, long nx_MIN_nPML_x_2, long nx_MIN_1, long nPML_y_1,
long ny_MIN_1_MIN_nPML_y_2, long nz_MIN_nPML_z_2, long nz_MIN_1,
long nPML_x_2_MIN_1, long nPML_z_2_MIN_1, double TwoOREp0)
{
long i, j, k, ii, kk, ia, ja, ka;
double Fy_r, Gy_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_1_MIN_nPML_y_2; j++)
{
ka = 0;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_nPML_z_2; k < nz_MIN_1 ; k++)
{
Fy_r = Fy_24[ia][ja][ka];
Fy_24[ia][ja][ka] = K_a[Ind[i][j][k]]*Fy_24[ia][ja][ka] +
K_b[Ind[i][j][k]]*( (Hx[i][j][k] - Hx[i][j][k-1])*inv_dz -
(Hz[i][j][k] - Hz[i-1][j][k])*inv_dx );
Gy_r = Gy_24[ia][ja][ka];
Gy_24[ia][ja][ka] = K_Gy_a_2[kk]*Gy_24[ia][ja][ka] +
K_Gy_b_2[kk]*(Fy_24[ia][ja][ka] - Fy_r);
Ey[i][j][k] = K_Ey_a_2[ii]*Ey[i][j][k] +
K_Ey_b_2[ii]*TwoOREp0*( Gy_24[ia][ja][ka] - Gy_r);
ka++;
kk--;
}
ja++;
}
ia++;
ii--;
}
}
void calc_Ey_25(double ***Ey, double ***Hx, double ***Hz, double ***Fy_25,
double ***Gy_25, long ***Ind, double *K_a, double *K_b, double *K_Gy_a_2,
double *K_Gy_b_2, double *K_Ey_a_1, double *K_Ey_b_1, double *K_Ey_c_2,
double *K_Ey_d_2, double inv_dx, double inv_dz, long nPML_x_1,
long ny_MIN_1_MIN_nPML_y_2, long ny_MIN_1, long nz_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 Fy_r, Gy_r;
//25
for (i = 1; 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_nPML_z_2; k < nz_MIN_1 ; k++)
{
Fy_r = Fy_25[i][ja][ka];
Fy_25[i][ja][ka] = K_a[Ind[i][j][k]]*Fy_25[i][ja][ka] +
K_b[Ind[i][j][k]]*( (Hx[i][j][k] - Hx[i][j][k-1])*inv_dz -
(Hz[i][j][k] - Hz[i-1][j][k])*inv_dx );
Gy_r = Gy_25[i][ja][ka];
Gy_25[i][ja][ka] = K_Gy_a_2[kk]*Gy_25[i][ja][ka] +
K_Gy_b_2[kk]*(Fy_25[i][ja][ka] - Fy_r);
Ey[i][j][k] = K_Ey_a_1[i]*Ey[i][j][k] +
K_Ey_b_1[i]*( K_Ey_c_2[jj]*Gy_25[i][ja][ka] + K_Ey_d_2[jj]*Gy_r);
ka++;
kk--;
}
ja++;
jj--;
}
}
}
void calc_Ey_26(double ***Ey, double ***Hx, double ***Hz, double ***Fy_26,
double ***Gy_26, long ***Ind, double *K_a, double *K_b, double *K_Gy_a_2,
double *K_Gy_b_2, double *K_Ey_c_2, double *K_Ey_d_2, double inv_dx,
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_nPML_z_2,
long nz_MIN_1, long nPML_y_2_MIN_1, long nPML_z_2_MIN_1,
double inv_TwoOREp0)
{
long i, j, k, jj, kk, ia, ja, ka;
double Fy_r, Gy_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_nPML_z_2; k < nz_MIN_1 ; k++)
{
Fy_r = Fy_26[ia][ja][ka];
Fy_26[ia][ja][ka] = K_a[Ind[i][j][k]]*Fy_26[ia][ja][ka] +
K_b[Ind[i][j][k]]*( (Hx[i][j][k] - Hx[i][j][k-1])*inv_dz -
(Hz[i][j][k] - Hz[i-1][j][k])*inv_dx );
Gy_r = Gy_26[ia][ja][ka];
Gy_26[ia][ja][ka] = K_Gy_a_2[kk]*Gy_26[ia][ja][ka] +
K_Gy_b_2[kk]*(Fy_26[ia][ja][ka] - Fy_r);
Ey[i][j][k] += inv_TwoOREp0*( K_Ey_c_2[jj]*Gy_26[ia][ja][ka] +
K_Ey_d_2[jj]*Gy_r);
ka++;
kk--;
}
ja++;
jj--;
}
ia++;
}
}
void calc_Ey_27(double ***Ey, double ***Hx, double ***Hz, double ***Fy_27,
double ***Gy_27, long ***Ind, double *K_a, double *K_b, double *K_Gy_a_2,
double *K_Gy_b_2, double *K_Ey_a_2, double *K_Ey_b_2, double *K_Ey_c_2,
double *K_Ey_d_2, double inv_dx, double inv_dz, long nx_MIN_nPML_x_2,
long nx_MIN_1, long ny_MIN_1_MIN_nPML_y_2, long ny_MIN_1,
long nz_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 Fy_r, Gy_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_1_MIN_nPML_y_2; j < ny_MIN_1; j++)
{
ka = 0;
kk = nPML_z_2_MIN_1;
for (k = nz_MIN_nPML_z_2; k < nz_MIN_1 ; k++)
{
Fy_r = Fy_27[ia][ja][ka];
Fy_27[ia][ja][ka] = K_a[Ind[i][j][k]]*Fy_27[ia][ja][ka] +
K_b[Ind[i][j][k]]*( (Hx[i][j][k] - Hx[i][j][k-1])*inv_dz -
(Hz[i][j][k] - Hz[i-1][j][k])*inv_dx );
Gy_r = Gy_27[ia][ja][ka];
Gy_27[ia][ja][ka] = K_Gy_a_2[kk]*Gy_27[ia][ja][ka] +
K_Gy_b_2[kk]*(Fy_27[ia][ja][ka] - Fy_r);
Ey[i][j][k] = K_Ey_a_2[ii]*Ey[i][j][k] +
K_Ey_b_2[ii]*( K_Ey_c_2[jj]*Gy_27[ia][ja][ka] + K_Ey_d_2[jj]*Gy_r);
ka++;
kk--;
}
ja++;
jj--;
}
ia++;
ii--;
}
}⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -