#%% if 1: # plot mesh fig = plt.figure(8); fig.clf() tets = bba.get_convergences_max_tets() varsz = bba.get_variable_vs(swpvar) Y = {} for key in tets.keys(): Y[varsz[key]] = tets[key] y = pd.Series(Y ) #.values(), index = varsz.values()) y.plot(marker = '*', ms = 20) ax7t = ax7.twinx() ax7t.plot(y, marker = '*', ms = 10, c = 'g') #%% if 1: fig = plt.figure(21); fig.clf() tts = bba.get_convergences_Tets_vs_pass() for key, x in tts.iteritems(): #np.log10(x).plot(label = varsz[key]) x.plot(label = varsz[key]) plt.legend(loc = 0) #%% if 0: variation = '0'; pJ_method = 'J_surf_mag'; #pJ_mj_series = bbp.calc_Pjs_from_I_for_mode(variation, bbp.U_H,bbp.U_E, bbp.LJs, junc_rects, junc_lens, method = pJ_method) # to be implemented res = bbp.calc_avg_current_J_surf_mag(variation,junc_rects[0], junc_lens[0]) if 0: # for debug variation = '0'; junc_rect = 'juncV'; print_color(' Setup: ' + bbp.setup.name) lv = bbp.get_lv(variation)