#AlN_thickness=random.uniform(20,55) AlNH_thickness += 0.2 GivenAgBot = nklib.Ag(Ag_thickness) GivenAlNHBot = nklib.AlNH(AlNH_thickness) GivenAlNHSeed = nklib.AlNH(AlNH_seed_thickness) #GivenAgTop = nklib.Ag(18.9+3.5) #GivenAlNTop = nklib.AlN(60) MyStack = ML([GivenAlNHSeed, GivenAgBot, GivenAlNHBot]) #MyStack = ML([GivenAgBot, GivenAlNBot, GivenAgTop, GivenAlNTop]) MyStack.calculate_TR() MyStack = bright.correct_TR(MyStack) MyStack.A = [1] * len(MyStack.T) - MyStack.T - MyStack.R L,a,b = GenerateColor.color_calc(MyStack.wl, MyStack.T, MyStack.R) index_lowerUV = 0 index_lowervis = 32 index_uppervis = 94 index_upper = 454 photopic_array = np.interp(MyStack.wl[index_lowervis:index_uppervis],P_data[:, 0] , P_data[:, 1]) sol_array = np.interp(MyStack.wl[index_lowerUV:index_upper],S_data[:, 0] , S_data[:, 3]) Tvis = sum(MyStack.T[index_lowervis:index_uppervis]*photopic_array)/(sum(photopic_array)) #TSER = 1 - sum(MyStack.T[index_lowerUV:index_upper]*sol_array)/(sum(sol_array)) - 0.04 TSER = sum(MyStack.R[index_lowerUV:index_upper]*sol_array)/(sum(sol_array))+0.85*sum(MyStack.A[index_lowerUV:index_upper]*sol_array)/(sum(sol_array)) #print 'UV lower bound:', index_lowerUV #print 'Visible lower bound:', index_lowervis
def recalculate_priority(arg): ''' This function returns a tuple containing a list of layers with calculated priority. The structure is given as an argument ''' seed=AlN(arg[1]) first=Ag(arg[2]) second=AlNH(arg[3]) third=Ag(arg[4]) fourth=AlNH(arg[5]) mystruct = ML([seed, first, second, third, fourth], min_wl=MIN_WL, max_wl=MAX_WL) ML.calculate_TR(mystruct) # Doing unchecked attribute mutation here...future person beware. This # only applies because of comment above in correct_TR mystruct = correct_TR(mystruct) index_lowerUV=0 index_lowervis=0 index_uppervis=0 upper_middle=0 index_upper=0 if mystruct.wl[len(mystruct.wl)-1]<3000: index_upper=len(mystruct.wl)-1 for index, i in enumerate(mystruct.wl): if index_lowervis ==0 and i>=myConst.lower_limit: index_lowervis = index elif index_uppervis==0 and i>=myConst.upper_limit: index_uppervis=index elif myConst.upper_limit!= 700 and upper_middle==0: if i>=700: upper_middle=index elif index_upper==0 and i>=3000: index_upper=index break if myConst.upper_limit==700: upper_middle=index_uppervis photopic_array = np.interp(mystruct.wl[index_lowervis:index_uppervis],P_data[:, 0] , P_data[:, 1]) V_normalized_data = np.interp(mystruct.wl[index_lowervis:index_uppervis],V_data[:, 0], V_data[:, 3]) / max(V_data[:, 3]) T_ideal_array = photopic_array/V_normalized_data #T_one_array = [1]*len(mystruct.wl[index_lowervis:index_uppervis]) sol_array = np.interp(mystruct.wl[index_lowerUV:index_upper],S_data[:, 0] , S_data[:, 3]) sol_array_UV = np.interp(mystruct.wl[index_lowerUV:index_lowervis],U_data[:, 0] , U_data[:, 3]) sol_array_vis = np.interp(mystruct.wl[index_lowervis:index_uppervis],V_data[:, 0] , V_data[:, 3]) sol_array_IR = np.interp(mystruct.wl[index_uppervis:index_upper],I_data[:, 0] , I_data[:, 3]) Tvis = sum(mystruct.T[index_lowervis:index_uppervis]*photopic_array)/(sum(photopic_array)) #Tpriority = sum(mystruct.T[index_lowervis:index_uppervis]*T_ideal_array)/(sum(T_ideal_array)) #Tpriority = sum(mystruct.T[index_lowervis:index_uppervis]*sol_array_vis*photopic_array)/(sum(sol_array_vis*photopic_array)) #Tpriority = 1/np.sqrt(sum((mystruct.T[index_lowervis:index_uppervis]*sol_array_vis/max(sol_array_vis)-photopic_array)**2)) #Tpriority = 1/np.sqrt(sum((mystruct.T[index_lowervis:index_uppervis]-sol_array_vis/max(sol_array_vis))**2)) #Tpriority = 1/np.sqrt(sum((mystruct.T[index_lowervis:index_uppervis]-T_one_array)**2)) Tpriority = Tvis TSER_UV = 1 - sum(mystruct.T[index_lowerUV:index_lowervis]*sol_array_UV)/(sum(sol_array_UV)) TSER_vis = 1 - sum(mystruct.T[index_lowervis:index_uppervis]*sol_array_vis)/(sum(sol_array_vis)) TSER_IR = 1 - sum(mystruct.T[index_uppervis:index_upper]*sol_array_IR)/(sum(sol_array_IR)) # 0.0342 = % of AM1.5 energy in UV # 0.426 = % of AM1.5 energy in visible # 0.54 = % of AM1.5 energy in IR (700-2500) #TSER = TSER_UV * 0.0342 + TSER_vis * 0.426 + TSER_IR * 0.54 - .04 TSER = 1 - sum(mystruct.T[index_lowerUV:index_upper]*sol_array)/(sum(sol_array)) - 0.04 #TSER_priority = TSER_UV * 0.0342 / (0.0342 + 0.54) + TSER_IR * 0.54 / (0.0342 + 0.54) TSER_priority = TSER_IR L,a,b = GenerateColor.color_calc(mystruct.wl, mystruct.T, mystruct.R) Cpriority = 1/ (abs(a) + abs(b)) #priority = myConst.p1 * Tpriority + myConst.p2 * (TSER_IR - myConst.wanted_TSER) priority = myConst.p1 * Tpriority + myConst.p2 * TSER_priority + myConst.p3 * Cpriority #priority = myConst.p1 * Tpriority * Tvis + myConst.p2 * TSER_priority #priority = -1 / (Tvis - myConst.p1) - \ # abs(myConst.p2 * (TSER - myConst.wanted_TSER)) return (priority, arg[1], arg[2], arg[3], arg[4], arg[5], Tvis, TSER)
def plot_structure(arg, text_only=False): ''' Basically, this function plots the graph of reflectance and transmittance of a given three layer structure versus wavelength. The argument is given as a tuple or a list of 3 elements. In addition, it provides the values of reflectance, transmittance, R_color and T_color. ''' seed = AlN(arg[1]) first = Ag(arg[2]) second = AlNH(arg[3]) third = Ag(arg[4]) fourth = AlNH(arg[5]) mystruct = ML([seed, first, second, third, fourth], min_wl=MIN_WL, max_wl=MAX_WL) TR(mystruct) show() #### Uncomment for 6-layer structure ## fifth = Ag(arg[4]) ## sixth = DLC80WA(arg[5]) ## ## mystruct = ML([first, second, third, fourth, ## fifth, sixth, AlN(20)], ## min_wl=230, max_wl=2476) print mystruct # Doing unchecked attribute mutation here...future person beware. This # only applies because of comment above in correct_TR mystruct = correct_TR(mystruct) index_lowerUV=0 index_lowervis=0 index_uppervis=0 upper_middle=0 index_upper=0 if mystruct.wl[len(mystruct.wl)-1]<3000: index_upper=len(mystruct.wl)-1 for index, i in enumerate(mystruct.wl): if index_lowervis ==0 and i>=myConst.lower_limit: index_lowervis = index elif index_uppervis==0 and i>=myConst.upper_limit: index_uppervis=index elif myConst.upper_limit!= 700 and upper_middle==0: if i>=700: upper_middle=index elif index_upper==0 and i>=3000: index_upper=index break if myConst.upper_limit==700: upper_middle=index_uppervis photopic_array = np.interp(mystruct.wl[index_lowervis:index_uppervis],P_data[:, 0] , P_data[:, 1]) V_normalized_data = np.interp(mystruct.wl[index_lowervis:index_uppervis],V_data[:, 0], V_data[:, 3]) / max(V_data[:, 3]) T_ideal_array = photopic_array/V_normalized_data #T_one_array = [1]*len(mystruct.wl[index_lowervis:index_uppervis]) sol_array = np.interp(mystruct.wl[index_lowerUV:index_upper],S_data[:, 0] , S_data[:, 3]) sol_array_UV = np.interp(mystruct.wl[index_lowerUV:index_lowervis],U_data[:, 0] , U_data[:, 3]) sol_array_vis = np.interp(mystruct.wl[index_lowervis:index_uppervis],V_data[:, 0] , V_data[:, 3]) sol_array_IR = np.interp(mystruct.wl[index_uppervis:index_upper],I_data[:, 0] , I_data[:, 3]) Tvis = sum(mystruct.T[index_lowervis:index_uppervis]*photopic_array)/(sum(photopic_array)) #Tpriority = sum(mystruct.T[index_lowervis:index_uppervis]*T_ideal_array)/(sum(T_ideal_array)) #Tpriority = sum(mystruct.T[index_lowervis:index_uppervis]*sol_array_vis*photopic_array)/(sum(sol_array_vis*photopic_array)) #Tpriority = 1/np.sqrt(sum((mystruct.T[index_lowervis:index_uppervis]*sol_array_vis/max(sol_array_vis)-photopic_array)**2)) #Tpriority = 1/np.sqrt(sum((mystruct.T[index_lowervis:index_uppervis]-sol_array_vis/max(sol_array_vis))**2)) #Tpriority = 1/np.sqrt(sum((mystruct.T[index_lowervis:index_uppervis]-T_one_array)**2)) Tpriority = Tvis TSER_UV = 1 - sum(mystruct.T[index_lowerUV:index_lowervis]*sol_array_UV)/(sum(sol_array_UV)) TSER_vis = 1 - sum(mystruct.T[index_lowervis:index_uppervis]*sol_array_vis)/(sum(sol_array_vis)) TSER_IR = 1 - sum(mystruct.T[index_uppervis:index_upper]*sol_array_IR)/(sum(sol_array_IR)) # 0.0342 = % of AM1.5 energy in UV # 0.426 = % of AM1.5 energy in visible # 0.54 = % of AM1.5 energy in IR (700-2500) #TSER = TSER_UV * 0.0342 + TSER_vis * 0.426 + TSER_IR * 0.54 - .04 TSER = 1 - sum(mystruct.T[index_lowerUV:index_upper]*sol_array)/(sum(sol_array)) - 0.04 #TSER_priority = TSER_UV * 0.0342 / (0.0342 + 0.54) + TSER_IR * 0.54 / (0.0342 + 0.54) TSER_priority = TSER_IR L,a,b = GenerateColor.color_calc(mystruct.wl, mystruct.T, mystruct.R) Cpriority = 1/ (abs(a) + abs(b)) #priority = myConst.p1 * Tpriority + myConst.p2 * (TSER_IR - myConst.wanted_TSER) priority = myConst.p1 * Tpriority + myConst.p2 * TSER_priority + myConst.p3 * Cpriority #priority = myConst.p1 * Tpriority * Tvis + myConst.p2 * TSER_priority #priority = -1 / (Tvis - myConst.p1) - \ # abs(myConst.p2 * (TSER - myConst.wanted_TSER)) print 'priorities are: ',myConst.p1, myConst.p2, myConst.p3 print 'TSER = %s' % float("{0:.4f}".format(TSER)) print 'TSER in UV = %s' % float("{0:.4f}".format(TSER_UV)) print 'TSER in visible = %s' % float("{0:.4f}".format(TSER_vis)) print 'TSER in IR = %s' % float("{0:.4f}".format(TSER_IR)) # print 'Tpriority = %s' % float("{0:.4f}".format(Tpriority)) print 'Tvis = %s' % float("{0:.4f}".format(Tvis)) print 'T_color:', L,a,b # print 'UV lower bound:', index_lowerUV # print 'Visible lower bound:', index_lowervis # print 'Visible upper bound:', index_uppervis # print 'IR upper bound:', index_upper #print mystruct.wl[:] if not text_only: TR(mystruct, show_solar=True, min_wl=200, max_wl=2500) plot.plt.grid("on") plot.plt.plot([400,400], [0,1], "k--") plot.plt.plot([700,700], [0,1], "k--") plot.plt.text(425, 0.70, "visible", fontsize=16) show() return mystruct
def create_structure(): ''' This function creates a structure composed of randomly generated thicknesses. It returns a tuple with a list of layer thicknesses and priority value. ''' L0_thickness = 15 seedlayer = AlN(L0_thickness) L1_thickness=random.uniform(12, 20.) L1_thickness=float("{0:.1f}".format(L1_thickness)) layer_one=Ag(L1_thickness) L2_thickness=random.uniform(20.0, myConst.max_DLC_thickness) L2_thickness=float("{0:.1f}".format(L2_thickness)) layer_two=AlNH(L2_thickness) L3_thickness=random.uniform(12, 20.) L3_thickness=float("{0:.1f}".format(L3_thickness)) layer_three=Ag(L3_thickness) L4_thickness=random.uniform(20.0, myConst.max_DLC_thickness) ## L4_thickness=random.uniform(20.0, 60) L4_thickness=float("{0:.1f}".format(L4_thickness)) layer_four=AlNH(L4_thickness) mystruct = ML([seedlayer, layer_one, layer_two, layer_three, layer_four], min_wl=MIN_WL, max_wl=MAX_WL) ML.calculate_TR(mystruct) # Doing unchecked attribute mutation here...future person beware. This # only applies because of comment above in correct_TR mystruct = correct_TR(mystruct) index_lowerUV=0 index_lowervis=0 index_uppervis=0 upper_middle=0 index_upper=0 if mystruct.wl[len(mystruct.wl)-1]<3000: index_upper=len(mystruct.wl)-1 for index, i in enumerate(mystruct.wl): if index_lowervis ==0 and i>=myConst.lower_limit: index_lowervis = index elif index_uppervis==0 and i>=myConst.upper_limit: index_uppervis=index elif myConst.upper_limit!= 700 and upper_middle==0: if i>=700: upper_middle=index elif index_upper==0 and i>=3000: index_upper=index break if myConst.upper_limit==700: upper_middle=index_uppervis photopic_array = np.interp(mystruct.wl[index_lowervis:index_uppervis],P_data[:, 0] , P_data[:, 1]) V_normalized_data = np.interp(mystruct.wl[index_lowervis:index_uppervis],V_data[:, 0], V_data[:, 3]) / max(V_data[:, 3]) T_ideal_array = photopic_array/V_normalized_data #T_one_array = [1]*len(mystruct.wl[index_lowervis:index_uppervis]) sol_array = np.interp(mystruct.wl[index_lowerUV:index_upper],S_data[:, 0] , S_data[:, 3]) sol_array_UV = np.interp(mystruct.wl[index_lowerUV:index_lowervis],U_data[:, 0] , U_data[:, 3]) sol_array_vis = np.interp(mystruct.wl[index_lowervis:index_uppervis],V_data[:, 0] , V_data[:, 3]) sol_array_IR = np.interp(mystruct.wl[index_uppervis:index_upper],I_data[:, 0] , I_data[:, 3]) Tvis = sum(mystruct.T[index_lowervis:index_uppervis]*photopic_array)/(sum(photopic_array)) #Tpriority = sum(mystruct.T[index_lowervis:index_uppervis]*T_ideal_array)/(sum(T_ideal_array)) #Tpriority = sum(mystruct.T[index_lowervis:index_uppervis]*sol_array_vis*photopic_array)/(sum(sol_array_vis*photopic_array)) #Tpriority = 1/np.sqrt(sum((mystruct.T[index_lowervis:index_uppervis]*sol_array_vis/max(sol_array_vis)-photopic_array)**2)) #Tpriority = 1/np.sqrt(sum((mystruct.T[index_lowervis:index_uppervis]-sol_array_vis/max(sol_array_vis))**2)) #Tpriority = 1/np.sqrt(sum((mystruct.T[index_lowervis:index_uppervis]-T_one_array)**2)) Tpriority = Tvis TSER_UV = 1 - sum(mystruct.T[index_lowerUV:index_lowervis]*sol_array_UV)/(sum(sol_array_UV)) TSER_vis = 1 - sum(mystruct.T[index_lowervis:index_uppervis]*sol_array_vis)/(sum(sol_array_vis)) TSER_IR = 1 - sum(mystruct.T[index_uppervis:index_upper]*sol_array_IR)/(sum(sol_array_IR)) # 0.0342 = % of AM1.5 energy in UV # 0.426 = % of AM1.5 energy in visible # 0.54 = % of AM1.5 energy in IR (700-2500) #TSER = TSER_UV * 0.0342 + TSER_vis * 0.426 + TSER_IR * 0.54 - .04 TSER = 1 - sum(mystruct.T[index_lowerUV:index_upper]*sol_array)/(sum(sol_array)) - 0.04 #TSER_priority = TSER_UV * 0.0342 / (0.0342 + 0.54) + TSER_IR * 0.54 / (0.0342 + 0.54) TSER_priority = TSER_IR L,a,b = GenerateColor.color_calc(mystruct.wl, mystruct.T, mystruct.R) Cpriority = 1/ (abs(a) + abs(b)) #priority = myConst.p1 * Tpriority + myConst.p2 * (TSER_IR - myConst.wanted_TSER) priority = myConst.p1 * Tpriority + myConst.p2 * TSER_priority + myConst.p3 * Cpriority #priority = myConst.p1 * Tpriority * Tvis + myConst.p2 * TSER_priority #priority = -1 / (Tvis - myConst.p1) - \ # abs(myConst.p2 * (TSER - myConst.wanted_TSER)) return (priority, L0_thickness, L1_thickness, L2_thickness, L3_thickness, L4_thickness, Tvis, TSER)