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

Corrected bug where rough surfaces would not get created.

浏览代码
这个提交包含在：
craig-warren
2015-10-13 21:28:37 +03:00
提交者 Craig Warren
父节点 193cc6d649
当前提交 172278a6f0
共有 1 个文件被更改，包括 3 次插入 和 1 次删除
下载 Patch 文件 下载 Diff 文件 展开所有文件 折叠所有文件

4
gprMax/input_cmds_geometry.py
查看文件

@@ -1408,13 +1408,15 @@ def process_geometrycmds(geometry, G):
for k in range(volume.zs, volume.zf + 1): for k in range(volume.zs, volume.zf + 1):
if volume.mask[i - volume.xs, j - volume.ys, k - volume.zs] == 1: if volume.mask[i - volume.xs, j - volume.ys, k - volume.zs] == 1:
numID = numIDx = numIDy = numIDz = volume.fractalvolume[i - volume.xs, j - volume.ys, k - volume.zs] numID = numIDx = numIDy = numIDz = volume.fractalvolume[i - volume.xs, j - volume.ys, k - volume.zs]
build_voxel(i, j, k, numID, numIDx, numIDy, numIDz, False, G.solid, G.rigidE, G.rigidH, G.ID)
elif volume.mask[i - volume.xs, j - volume.ys, k - volume.zs] == 2: elif volume.mask[i - volume.xs, j - volume.ys, k - volume.zs] == 2:
waternumID = next(x.numID for x in G.materials if x.ID == 'water') waternumID = next(x.numID for x in G.materials if x.ID == 'water')
numID = numIDx = numIDy = numIDz = waternumID numID = numIDx = numIDy = numIDz = waternumID
build_voxel(i, j, k, numID, numIDx, numIDy, numIDz, False, G.solid, G.rigidE, G.rigidH, G.ID)
elif volume.mask[i - volume.xs, j - volume.ys, k - volume.zs] == 3: elif volume.mask[i - volume.xs, j - volume.ys, k - volume.zs] == 3:
grassnumID = next(x.numID for x in G.materials if x.ID == 'grass') grassnumID = next(x.numID for x in G.materials if x.ID == 'grass')
numID = numIDx = numIDy = numIDz = grassnumID numID = numIDx = numIDy = numIDz = grassnumID
build_voxel(i, j, k, numID, numIDx, numIDy, numIDz, False, G.solid, G.rigidE, G.rigidH, G.ID) build_voxel(i, j, k, numID, numIDx, numIDy, numIDz, False, G.solid, G.rigidE, G.rigidH, G.ID)
else: else:
if volume.nbins == 1: if volume.nbins == 1: