本文整理汇总了Python中hyperion.model.Model.set_voronoi_grid方法的典型用法代码示例。如果您正苦于以下问题:Python Model.set_voronoi_grid方法的具体用法?Python Model.set_voronoi_grid怎么用?Python Model.set_voronoi_grid使用的例子?那么, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类hyperion.model.Model
的用法示例。
在下文中一共展示了Model.set_voronoi_grid方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: random_in_sphere
# 需要导入模块: from hyperion.model import Model [as 别名]
# 或者: from hyperion.model.Model import set_voronoi_grid [as 别名]
# Inner edge
x3, y3, z3 = random_in_sphere(RSUB * 0.95, RSUB * 1.05, n3)
x3 = np.abs(x3)
y3 = np.abs(y3)
z3 = np.abs(z3)
# Concatenate
x = np.hstack([x0, x1, x2, x3])
y = np.hstack([y0, y1, y2, y3])
z = np.hstack([z0, z1, z2, z3])
# np.savetxt('points.txt', list(zip(x, y, z)))
m.set_voronoi_grid(x, y, z,
xmin=-10 * au, xmax=60 * au,
ymin=-10 * au, ymax=60 * au,
zmin=-10 * au, zmax=60 * au)
density = np.zeros(len(x))
# Ambient medium
in_box = (x > 0.) & (x < 60 * au) & (y > 0.) & (y < 60 * au) & (z > 0.) & (z < 60 * au)
density[in_box] = rho_0
# Set up sphere 1
in_sphere_1 = (x - 10 * au) ** 2 + (y - 15 * au) ** 2 + (z - 20 * au) ** 2 < r_1 ** 2
density[in_sphere_1] = rho_1
# Set up sphere 2
in_sphere_2 = (x - 26.666667 * au) ** 2 + (y - 31.666667 * au) ** 2 + (z - 28.333333 * au) ** 2 < r_2 ** 2
density[in_sphere_2] = rho_2