Esempio n. 1
0
        ]
    })

#%%==============================================================================
#     Plot results for sweep
#==============================================================================
if 1:
    swpvar = 'pin_shift'
    RES = []
    SWP = []
    for key, s in sol.iteritems():
        varz = hfss_variables[key]
        SWP += [ureg.Quantity(varz['_' + swpvar]).magnitude]
        RES += [
            eBBQ_Pmj_to_H_params(s,
                                 meta_datas[key],
                                 cos_trunc=cos_trunc,
                                 fock_trunc=fock_trunc)
        ]
    import matplotlib.gridspec as gridspec
    #%%
    fig = plt.figure(num=1, figsize=(19, 5))
    gs1 = gridspec.GridSpec(1, 4, width_ratios=[2, 2, 2, 1])
    gs1.update(left=0.05, right=0.98)  # wspace=0.05
    ax1 = plt.subplot(gs1[0])
    ax2 = plt.subplot(gs1[1])
    ax3 = plt.subplot(gs1[2])
    ax3b = plt.subplot(gs1[3])

    ax = ax1
    ID = 1
    args = {'lw': 0, 'marker': 'o', 'ms': 5}
Esempio n. 2
0
    swp_var   = 'inductance_shift'
    bba = BbqAnalysis(file_name)
    hfss_variables = bba.hfss_variables
    sols           = bba.sols
    meta_datas     = bba.meta_data


#%% 
if analyze_BBQ:  
    RES = []; SWP = [];
    for key, s in sols.iteritems():
        print '\r Analyzing ', key,
        try:
            varz  = hfss_variables[key]
            SWP  += [ ureg.Quantity(varz['_'+swp_var]).magnitude ]  
            RES  += [ eBBQ_Pmj_to_H_params(s, meta_datas[key], cos_trunc = cos_trunc, fock_trunc = fock_trunc) ]
        except Exception as e:
            print_color(" !ERROR %s" % e)
#%%==============================================================================
# Plot Chis
#==============================================================================
if plot_Fs:
    #Fs  = bba.get_Fs(swp_var=swp_var)
    Fs = []
    for i in range(len(RES)):
        Fs.append(RES[i][9])
        
    plt.plot(sorted(SWP), np.array(Fs).transpose()[0]*10**-9, label='D mode freq', marker = 'o', lw =1)
    plt.plot(sorted(SWP), np.array(Fs).transpose()[1]*10**-9, label='B mode freq', marker = 'o', lw =1)
    
    plt.title('D and B mode frequencies (nonlinear)');
Esempio n. 3
0
#    print 'Q={:8.1e} {:7.1e} {:6.0f}'.format(*(Qs))
    #print pd.Series({ key:varz[key] for key in ['_join_w','_join_h','_padV_width', '_padV_height','_padH_width', '_padH_height','_scaleV','_scaleH', '_LJ1','_LJ2'] })

#%%==============================================================================
#     Plot results for sweep
#==============================================================================
if 1:
    swpvar='LJ'    
    use_1st_order = True  # use 1st O PT  to identify correct eigenvectors in ND
    RES = []; SWP = [];
    for key, s in sol.iteritems():     
        print '\r Analyzing ', key,
        try:
            varz  = hfss_variables[key]
            SWP  += [ ureg.Quantity(varz['_'+swpvar]).magnitude ]  
            RES  += [ eBBQ_Pmj_to_H_params(s, meta_datas[key], cos_trunc = cos_trunc, fock_trunc = fock_trunc, use_1st_order = False) ]
        except Exception as e:
            print_color(" !ERROR %s" % e)
    import matplotlib.gridspec as gridspec;
    #%%
    fig = plt.figure(num = 1, figsize=(19,5)) 
    gs1 = gridspec.GridSpec(1, 4, width_ratios=[2,2,2,1]); gs1.update(left=0.05, right=0.98)  # wspace=0.05
    ax1 = plt.subplot(gs1[0]); ax2 = plt.subplot(gs1[1]); ax3 = plt.subplot(gs1[2]); ax3b = plt.subplot(gs1[3])
    
    ax = ax1; ID = 1; 
    args = {'lw':0,'marker':'o','ms':5}
    ax.plot(SWP, [r[ID][0,1]for r in RES], label = '$\\chi_{DB}$', **args)
    ax.plot(SWP, [r[ID][0,2]for r in RES], label = '$\\chi_{DC}$', **args)
    ax.plot(SWP, [r[ID][1,2]for r in RES], label = '$\\chi_{BC}$', **args)
    ax.set_ylim([0.01,10**2]); ax.set_xlabel(swpvar); ax.set_title('cross-Kerr'); ax.set_ylabel('$\\chi$ (MHz)'); ax.legend(loc = 0)
    ax.set_yscale('log');   ax.set_ylim(0.1,100)    
Esempio n. 4
0
#    print 'Q={:8.1e} {:7.1e} {:6.0f}'.format(*(Qs))
    #print pd.Series({ key:varz[key] for key in ['_join_w','_join_h','_padV_width', '_padV_height','_padH_width', '_padH_height','_scaleV','_scaleH', '_LJ1','_LJ2'] })

#%%==============================================================================
#     Plot results for sweep
#==============================================================================
if 1:
    swpvar='LJ'    
    RES = []; SWP = [];
    for key, s in sol.iteritems():     
        print '\r Analyzing ', key,
        try:
            varz  = hfss_variables[key]
            SWP  += [ ureg.Quantity(varz['_'+swpvar]).magnitude ]  
            RES  += [ eBBQ_Pmj_to_H_params(s, meta_datas[key], 
                                           cos_trunc = cos_trunc, fock_trunc = fock_trunc,
                                           _renorm_pj = True) ]
        except Exception as e:
            print_color(" !ERROR %s" % e)
    import matplotlib.gridspec as gridspec;
    #%%
    fig = plt.figure(num = 1, figsize=(19,5)) 
    gs1 = gridspec.GridSpec(1, 4, width_ratios=[2,2,2,1]); gs1.update(left=0.05, right=0.98)  # wspace=0.05
    ax1 = plt.subplot(gs1[0]); ax2 = plt.subplot(gs1[1]); ax3 = plt.subplot(gs1[2]); ax3b = plt.subplot(gs1[3])
    
    ax = ax1; ID = 1; 
    ax.set_title('cross-Kerr')
    args = {'lw':0,'marker':'o','ms':5}
    ax.plot(SWP, [r[ID][0,1]for r in RES], label = '$\\chi_{DB}$', **args)
    ax.set_xlabel(swpvar)
    ax.axhline(5)