import post_processing
import pandas as pd
import gdf
pd.options.display.float_format = '{:.3f}'.format
pd.set_option('display.width', 1000)
try:
#from main import *
pass
except Exception:
pass
#total = data_extraction.collect_from_dir('results/profiled_gamma_constant/total')
#non_turbine = data_extraction.collect_from_dir('results/profiled_gamma_constant/non_turbine')
print('started pi_trend')
pi_trend = data_extraction.collect_from_dir('results/profiled_gamma_constant/pi_stag_trend')
print('finished pi_trend')
print('started perfect')
perfect = pi_trend.ix[(pi_trend.root_c_u >= 0) & (pi_trend.D_out_1 <= 0.95)]
print('finished perfect')
#df_list = [compressor.get_stage_info_data_frame() for compressor in perfect.compressor]
#mach_list = [df.M_w_1 for df in df_list]
#mach_df = pd.DataFrame.from_records(mach_list)
compressor = perfect.ix[perfect.D_out_1 == perfect.D_out_1.min()].compressor.values[0]
#plt.scatter(total.D_out_1, total.root_c_u)
#plt.grid()
#rotor_profiler = compressor.first_stage.rotor_profiler
import data_extraction
from post_processing import *
from start_functions import calculate_mean_radius, do_profiling
import blade_strength
import geometry_results
import post_processing
import materials
import supporting_blade_profilers
pd.options.display.float_format = "{:.3f}".format
pd.set_option("display.width", 1000)
result = data_extraction.collect_from_dir("results/fit_dir/optimal")
compressor = result.compressor.values[0]
lattice_density = 1.2
blade_elongation = 3
blade_windage = 1
profiler = compressor.first_stage.rotor_profiler
D_out_inlet = profiler.blading_geometry.D_out_inlet
d_rel_inlet = profiler.blading_geometry.d_rel_inlet
temp = supporting_blade_profilers.A40InletStatorProfiler()
temp.stage_profiler = profiler
temp.mean_lattice_density = lattice_density
temp.blade_elongation = blade_elongation
temp.blade_windage = blade_windage
temp.D_out_inlet = D_out_inlet
temp.d_rel_inlet = d_rel_inlet
