Alphashapefix2

from ntrfc.cascade_case.domain import DomainParameters
import importlib.resources
import numpy as np
import pyvista as pv
import os

profilepoints_file = importlib.resources.files("ntrfc") / "data/turbine_cascade_2/profilepoints.txt"



points = np.loadtxt(profilepoints_file)
pointcloud = pv.PolyData(points)

ans = []
for i in range(0, 100):

    turbine_domain = DomainParameters()
    turbine_domain.generate_params_by_pointcloud(pointcloud)
    ans.append(turbine_domain.alpha)
    warn = "" if turbine_domain.alpha < 0.015 else "WARNING"
    print(f"{i}" + warn)

# pv.set_jupyter_backend("static")
#
# if os.getenv('DISPLAY') is None:
#     pv.start_xvfb()  # Start X virtual framebuffer (Xvfb)
#
# plotpath = turbine_domain.plot_domainparas("bla.jpg")
#
import matplotlib.pyplot as plt
import numpy as np



plt.hist(ans, bins=10)  # density=False would make counts
plt.ylabel('Count')
plt.xlabel('Alpha')
plt.xlim(0.006,0.016)
plt.grid()
plt.show()