print "debugging of system class"
start_delta_time = time.time()
# Structural age
service_time = np.arange(START_AGE+TIME_INTERVAL,END_AGE+TIME_INTERVAL,TIME_INTERVAL)
comp_type_list = ['flexure', 'shear', 'deck']
icorr_mean_list = [1, 1, 1]
# construct virgin component
str_yr_list = [0, 0, 0]
cost_list = []
virgin_component_list = []
for comp_type,str_yr,icorr_mean in zip(comp_type_list, str_yr_list, icorr_mean_list):
virgin_component = Component(comp_type, maintain_tag=False, str_yr=str_yr)
resistance_mean,resistance_cov,cost = simpleCorrosionLHS(comp_type, service_time, icorr_mean, str_yr)
virgin_component.setServiceTime(service_time)
virgin_component.setResistanceMean(resistance_mean)
virgin_component.setResistanceCov(resistance_cov)
virgin_component.setCESampling(NUM_COMPONENT, NUM_ADAPTATION, NUM_PRE_SMP, NUM_MAIN_SMP)
virgin_component_list.append(virgin_component)
cost_list.append(cost)
# strengthened component
str_yr_list = [0, 50, 0]
cost_list = []
component_list = []
for comp_type,str_yr,icorr_mean in zip(comp_type_list, str_yr_list, icorr_mean_list):
component = Component(comp_type, maintain_tag=False, str_yr=str_yr)
resistance_mean,resistance_cov,cost = simpleCorrosionLHS(comp_type, service_time, icorr_mean, str_yr)
component.setServiceTime(service_time)
component.setResistanceMean(resistance_mean)
