system.create_model()
binary_vertices = system.get_vertices()
binary_norms = system.get_norms()
# dot surface visualization
# -------------------------
if False:
Plt.plot_3d(vertices=[binary_vertices["system"]],
faces=None,
normals_view=False,
points_view=True,
faces_view=False,
point_color="r",
point_size=10.0,
face_color="c",
edge_color="k",
azim=45,
elev=45,
save=False,
s_format="pdf",
x_range=[-0.2, 1.2],
y_range=[-0.5, 0.5],
z_range=[-0.5, 0.5],
dpi=1000)
# application of 3d Delaunay convex hull triangulation to primary component of star system
# ----------------------------------------------------------------------------------------
if False:
primary_faces = d3d(vertices=binary_vertices["primary"])[1]
# vizualisation
# for plotting without vertices, set variable "points_view" in function "plot_3d" to False
# visualization
faces = [model["primary"]["t_object"]]
color = [[colors[-1]] * len(model["primary"]["t_object"])]
for i in model["primary"]["spots"]:
faces.append(model["primary"]["spots"][i])
# noinspection PyTypeChecker
color.append([colors[i]] * len(model["primary"]["spots"][i]))
Plt.plot_3d(vertices=None,
faces=faces,
normals_view=False,
points_view=False,
faces_view=True,
point_color=None,
point_size=5.0,
face_color=color,
edge_color="k",
elev=45,
azim=45,
save=False)
# two spots on the secondary component
# ------------------------------------
if False:
system.create_spots(meta=[{
"lon": 1.5,
"lat": 1.0,
"diameter": 0.4,
"steps_azimuthal": 4,
"steps_radial": 4,
