publicImage/gprMax
1
0
派生 0
镜像自地址 https://gitee.com/sunhf/gprMax.git 已同步 2025-08-08 07:24:19 +08:00
代码 工单 软件包 项目 版本发布 百科 动态

Update knl_fields_updates.py

浏览代码
这个提交包含在：
rania-pa
2023-07-25 12:40:32 +01:00
提交者 GitHub
父节点 988128e3d2
当前提交 79f319abbf
共有 1 个文件被更改，包括 4 次插入 和 4 次删除
下载 Patch 文件 下载 Diff 文件 展开所有文件 折叠所有文件

4
gprMax/cuda_opencl/knl_fields_updates.py
查看文件

@@ -222,6 +222,7 @@ update_electric_dispersive_A = {
// dispersive, ID and field // dispersive, ID and field
// component arrays. // component arrays.
$CUDA_IDX
// Convert the linear index to subscripts for 3D field arrays // Convert the linear index to subscripts for 3D field arrays
int x = i / ($NY_FIELDS * $NZ_FIELDS); int x = i / ($NY_FIELDS * $NZ_FIELDS);
@@ -238,7 +239,6 @@ update_electric_dispersive_A = {
int y_T = ((i % ($NX_T * $NY_T * $NZ_T)) % ($NY_T * $NZ_T)) / $NZ_T; int y_T = ((i % ($NX_T * $NY_T * $NZ_T)) % ($NY_T * $NZ_T)) / $NZ_T;
int z_T = ((i % ($NX_T * $NY_T * $NZ_T)) % ($NY_T * $NZ_T)) % $NZ_T; int z_T = ((i % ($NX_T * $NY_T * $NZ_T)) % ($NY_T * $NZ_T)) % $NZ_T;
$CUDA_IDX
// Ex component // Ex component
if ((NY != 1 || NZ != 1) && x >= 0 && x < NX && y > 0 && y < NY && z > 0 && z < NZ) { if ((NY != 1 || NZ != 1) && x >= 0 && x < NX && y > 0 && y < NY && z > 0 && z < NZ) {
@@ -334,6 +334,7 @@ update_electric_dispersive_B = {
// dispersive, ID and field // dispersive, ID and field
// component arrays. // component arrays.
$CUDA_IDX
// Convert the linear index to subscripts for 3D field arrays // Convert the linear index to subscripts for 3D field arrays
int x = i / ($NY_FIELDS * $NZ_FIELDS); int x = i / ($NY_FIELDS * $NZ_FIELDS);
@@ -350,7 +351,6 @@ update_electric_dispersive_B = {
int y_T = ((i % ($NX_T * $NY_T * $NZ_T)) % ($NY_T * $NZ_T)) / $NZ_T; int y_T = ((i % ($NX_T * $NY_T * $NZ_T)) % ($NY_T * $NZ_T)) / $NZ_T;
int z_T = ((i % ($NX_T * $NY_T * $NZ_T)) % ($NY_T * $NZ_T)) % $NZ_T; int z_T = ((i % ($NX_T * $NY_T * $NZ_T)) % ($NY_T * $NZ_T)) % $NZ_T;
$CUDA_IDX
// Ex component // Ex component
if ((NY != 1 || NZ != 1) && x >= 0 && x < NX && y > 0 && y < NY && z > 0 && z < NZ) { if ((NY != 1 || NZ != 1) && x >= 0 && x < NX && y > 0 && y < NY && z > 0 && z < NZ) {