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

Updated to use IDlookup method.

浏览代码
这个提交包含在：
Craig Warren
2016-05-06 12:30:39 +01:00
父节点 a39cd5853e
当前提交 c91c9ba583
共有 1 个文件被更改，包括 8 次插入 和 8 次删除
下载 Patch 文件 下载 Diff 文件 展开所有文件 折叠所有文件

16
gprMax/yee_cell_build.pyx
查看文件

@@ -108,10 +108,10 @@ cpdef void build_electric_components(np.uint32_t[:, :, ::1] solid, np.int8_t[:,
""" """
cdef Py_ssize_t i, j, k cdef Py_ssize_t i, j, k
cdef int numID1, numID2, numID3, numID4 cdef int numID1, numID2, numID3, numID4, componentID
# Ex component # Ex component
componentID = 0 componentID = G.IDlookup['Ex']
for i in range(0, G.nx): for i in range(0, G.nx):
for j in range(1, G.ny): for j in range(1, G.ny):
for k in range(1, G.nz): for k in range(1, G.nz):
@@ -133,7 +133,7 @@ cpdef void build_electric_components(np.uint32_t[:, :, ::1] solid, np.int8_t[:,
create_electric_average(i, j, k, numID1, numID2, numID3, numID4, componentID, G) create_electric_average(i, j, k, numID1, numID2, numID3, numID4, componentID, G)
# Ey component # Ey component
componentID = 1 componentID = G.IDlookup['Ey']
for i in range(1, G.nx): for i in range(1, G.nx):
for j in range(0, G.ny): for j in range(0, G.ny):
for k in range(1, G.nz): for k in range(1, G.nz):
@@ -155,7 +155,7 @@ cpdef void build_electric_components(np.uint32_t[:, :, ::1] solid, np.int8_t[:,
create_electric_average(i, j, k, numID1, numID2, numID3, numID4, componentID, G) create_electric_average(i, j, k, numID1, numID2, numID3, numID4, componentID, G)
# Ez component # Ez component
componentID = 2 componentID = G.IDlookup['Ez']
for i in range(1, G.nx): for i in range(1, G.nx):
for j in range(1, G.ny): for j in range(1, G.ny):
for k in range(0, G.nz): for k in range(0, G.nz):
@@ -185,10 +185,10 @@ cpdef void build_magnetic_components(np.uint32_t[:, :, ::1] solid, np.int8_t[:,
""" """
cdef Py_ssize_t i, j, k cdef Py_ssize_t i, j, k
cdef int numID1, numID2 cdef int numID1, numID2, componentID
# Hx component # Hx component
componentID = 3 componentID = G.IDlookup['Hx']
for i in range(1, G.nx): for i in range(1, G.nx):
for j in range(0, G.ny): for j in range(0, G.ny):
for k in range(0, G.nz): for k in range(0, G.nz):
@@ -208,7 +208,7 @@ cpdef void build_magnetic_components(np.uint32_t[:, :, ::1] solid, np.int8_t[:,
create_magnetic_average(i, j, k, numID1, numID2, componentID, G) create_magnetic_average(i, j, k, numID1, numID2, componentID, G)
# Hy component # Hy component
componentID = 4 componentID = G.IDlookup['Hy']
for i in range(0, G.nx): for i in range(0, G.nx):
for j in range(1, G.ny): for j in range(1, G.ny):
for k in range(0, G.nz): for k in range(0, G.nz):
@@ -228,7 +228,7 @@ cpdef void build_magnetic_components(np.uint32_t[:, :, ::1] solid, np.int8_t[:,
create_magnetic_average(i, j, k, numID1, numID2, componentID, G) create_magnetic_average(i, j, k, numID1, numID2, componentID, G)
# Hz component # Hz component
componentID = 5 componentID = G.IDlookup['Hz']
for i in range(0, G.nx): for i in range(0, G.nx):
for j in range(0, G.ny): for j in range(0, G.ny):
for k in range(1, G.nz): for k in range(1, G.nz):