set_theano_shared_topology

set_theano_shared_topology是Gempy中InterpolatorModel类的一个方法。它用于在Theano共享变量中存储地质体系的拓扑结构。

参数

• n_grid_points：整数，表示网格点数。
• sorter：numpy数组，包含形状为(n_grid_points, 6)的网格点信息。
• R：numpy数组，包含形状为(3, 3)的旋转矩阵。
• verbose：布尔值，可选参数，表示是否显示详细信息。默认为False。

返回值

该方法没有返回值，但它会将地质体系的拓扑结构存储在theano共享变量中。

用法示例

from gempy.core.interpolator import InterpolatorModel
import numpy as np

# 创建一个简单的拓扑结构
n_grid_points = 4
sorter = np.array([
  [0, 0, 0, 1, 0, 0],
  [1, 0, 0, 0, 1, 0],
  [1, 1, 0, 0, 0, 1],
  [0, 1, 1, 1, 1, 1]])
R = np.eye(3)

# 创建`InterpolatorModel`对象并设置拓扑结构
model = InterpolatorModel()
model.set_theano_shared_topology(n_grid_points, sorter, R, verbose=True)

# 输出拓扑结构的形状
print(model.grid.topology.shape)

输出：

Setting theano shared Topology and Attributes
Topology correctly set
(6, 4)

以上示例演示了如何使用set_theano_shared_topology方法在InterpolatorModel对象中设置地质体系的拓扑结构。该方法会在终端输出一条消息，提示拓扑结构已经正确设置，并输出拓扑结构在共享变量中的形状。