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

Updated to account for new PML boundary ID naming.

浏览代码
这个提交包含在：
Craig Warren
2017-01-30 18:16:49 +00:00
父节点 8e419b26e6
当前提交 22fe57e3da
共有 1 个文件被更改，包括 13 次插入 和 13 次删除
下载 Patch 文件 下载 Diff 文件 展开所有文件 折叠所有文件

26
gprMax/input_cmds_singleuse.py
查看文件

@@ -133,18 +133,18 @@ def process_singlecmds(singlecmds, G):
if G.nx == 1: if G.nx == 1:
G.dt = 1 / (c * np.sqrt((1 / G.dy) * (1 / G.dy) + (1 / G.dz) * (1 / G.dz))) G.dt = 1 / (c * np.sqrt((1 / G.dy) * (1 / G.dy) + (1 / G.dz) * (1 / G.dz)))
G.dimension = '2D' G.dimension = '2D'
G.pmlthickness['xminus'] = 0 G.pmlthickness['x0'] = 0
G.pmlthickness['xplus'] = 0 G.pmlthickness['xmax'] = 0
elif G.ny == 1: elif G.ny == 1:
G.dt = 1 / (c * np.sqrt((1 / G.dx) * (1 / G.dx) + (1 / G.dz) * (1 / G.dz))) G.dt = 1 / (c * np.sqrt((1 / G.dx) * (1 / G.dx) + (1 / G.dz) * (1 / G.dz)))
G.dimension = '2D' G.dimension = '2D'
G.pmlthickness['yminus'] = 0 G.pmlthickness['y0'] = 0
G.pmlthickness['yplus'] = 0 G.pmlthickness['ymax'] = 0
elif G.nz == 1: elif G.nz == 1:
G.dt = 1 / (c * np.sqrt((1 / G.dx) * (1 / G.dx) + (1 / G.dy) * (1 / G.dy))) G.dt = 1 / (c * np.sqrt((1 / G.dx) * (1 / G.dx) + (1 / G.dy) * (1 / G.dy)))
G.dimension = '2D' G.dimension = '2D'
G.pmlthickness['zminus'] = 0 G.pmlthickness['z0'] = 0
G.pmlthickness['zplus'] = 0 G.pmlthickness['zmax'] = 0
else: else:
G.dt = 1 / (c * np.sqrt((1 / G.dx) * (1 / G.dx) + (1 / G.dy) * (1 / G.dy) + (1 / G.dz) * (1 / G.dz))) G.dt = 1 / (c * np.sqrt((1 / G.dx) * (1 / G.dx) + (1 / G.dy) * (1 / G.dy) + (1 / G.dz) * (1 / G.dz)))
G.dimension = '3D' G.dimension = '3D'
@@ -200,13 +200,13 @@ def process_singlecmds(singlecmds, G):
for key in G.pmlthickness.keys(): for key in G.pmlthickness.keys():
G.pmlthickness[key] = int(tmp[0]) G.pmlthickness[key] = int(tmp[0])
else: else:
G.pmlthickness['xminus'] = int(tmp[0]) G.pmlthickness['x0'] = int(tmp[0])
G.pmlthickness['yminus'] = int(tmp[1]) G.pmlthickness['y0'] = int(tmp[1])
G.pmlthickness['zminus'] = int(tmp[2]) G.pmlthickness['z0'] = int(tmp[2])
G.pmlthickness['xplus'] = int(tmp[3]) G.pmlthickness['xmax'] = int(tmp[3])
G.pmlthickness['yplus'] = int(tmp[4]) G.pmlthickness['ymax'] = int(tmp[4])
G.pmlthickness['zplus'] = int(tmp[5]) G.pmlthickness['zmax'] = int(tmp[5])
if 2 * G.pmlthickness['xminus'] >= G.nx or 2 * G.pmlthickness['yminus'] >= G.ny or 2 * G.pmlthickness['zminus'] >= G.nz or 2 * G.pmlthickness['xplus'] >= G.nx or 2 * G.pmlthickness['yplus'] >= G.ny or 2 * G.pmlthickness['zplus'] >= G.nz: if 2 * G.pmlthickness['x0'] >= G.nx or 2 * G.pmlthickness['y0'] >= G.ny or 2 * G.pmlthickness['z0'] >= G.nz or 2 * G.pmlthickness['xmax'] >= G.nx or 2 * G.pmlthickness['ymax'] >= G.ny or 2 * G.pmlthickness['zmax'] >= G.nz:
raise CmdInputError(cmd + ' has too many cells for the domain size') raise CmdInputError(cmd + ' has too many cells for the domain size')
# src_steps # src_steps