Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- set_i={}
- vector_i={}
- vector_b={}
- sort={}
- set_i[10]={'n':[],'s':[]}
- set_i[10]['n'].append(0.0244605000)
- set_i[10]['s'].append(0.0198990000)
- vector_i[10]={'n':[],'s':[]}
- vector_i[10]['n'].append(687.9790485000)
- vector_i[10]['s'].append(693.6321805000)
- vector_b[10]={'n':[],'s':[]}
- vector_b[10]['n'].append(86.5708590000)
- vector_b[10]['s'].append(86.0124780000)
- sort[10]={'n':[],'s':[]}
- sort[10]['n'].append(0.0003610000)
- sort[10]['s'].append(0.0003230000)
- set_i[30]={'n':[],'s':[]}
- set_i[30]['n'].append(0.0636460000)
- set_i[30]['s'].append(0.0629870000)
- vector_i[30]={'n':[],'s':[]}
- vector_i[30]['n'].append(691.1723425000)
- vector_i[30]['s'].append(699.8295890000)
- vector_b[30]={'n':[],'s':[]}
- vector_b[30]['n'].append(84.9870685000)
- vector_b[30]['s'].append(83.3700055000)
- sort[30]={'n':[],'s':[]}
- sort[30]['n'].append(0.0012905000)
- sort[30]['s'].append(0.0013700000)
- set_i[90]={'n':[],'s':[]}
- set_i[90]['n'].append(0.1830740000)
- set_i[90]['s'].append(0.1802080000)
- vector_i[90]={'n':[],'s':[]}
- vector_i[90]['n'].append(699.1795235000)
- vector_i[90]['s'].append(718.8977727500)
- vector_b[90]={'n':[],'s':[]}
- vector_b[90]['n'].append(88.8957545000)
- vector_b[90]['s'].append(90.2301260000)
- sort[90]={'n':[],'s':[]}
- sort[90]['n'].append(0.0050220000)
- sort[90]['s'].append(0.0052515000)
- set_i[271]={'n':[],'s':[]}
- set_i[271]['n'].append(0.6951020000)
- set_i[271]['s'].append(0.5922210000)
- vector_i[271]={'n':[],'s':[]}
- vector_i[271]['n'].append(716.6927007500)
- vector_i[271]['s'].append(701.6047875000)
- vector_b[271]={'n':[],'s':[]}
- vector_b[271]['n'].append(86.2946895000)
- vector_b[271]['s'].append(84.1444715000)
- sort[271]={'n':[],'s':[]}
- sort[271]['n'].append(0.0289420000)
- sort[271]['s'].append(0.0169115000)
- set_i[815]={'n':[],'s':[]}
- set_i[815]['n'].append(1.9326585000)
- set_i[815]['s'].append(1.6610015000)
- vector_i[815]={'n':[],'s':[]}
- vector_i[815]['n'].append(709.5025380000)
- vector_i[815]['s'].append(702.2022080000)
- vector_b[815]={'n':[],'s':[]}
- vector_b[815]['n'].append(84.4623145000)
- vector_b[815]['s'].append(84.0353330000)
- sort[815]={'n':[],'s':[]}
- sort[815]['n'].append(0.0606105000)
- sort[815]['s'].append(0.0605120000)
- set_i[2451]={'n':[],'s':[]}
- set_i[2451]['n'].append(5.3122310000)
- set_i[2451]['s'].append(5.7808360000)
- vector_i[2451]={'n':[],'s':[]}
- vector_i[2451]['n'].append(686.6794080000)
- vector_i[2451]['s'].append(683.3603195000)
- vector_b[2451]={'n':[],'s':[]}
- vector_b[2451]['n'].append(85.5694445000)
- vector_b[2451]['s'].append(90.8011540000)
- sort[2451]={'n':[],'s':[]}
- sort[2451]['n'].append(0.2127380000)
- sort[2451]['s'].append(0.2157050000)
- set_i[7368]={'n':[],'s':[]}
- set_i[7368]['n'].append(19.1625330000)
- set_i[7368]['s'].append(18.7387995000)
- vector_i[7368]={'n':[],'s':[]}
- vector_i[7368]['n'].append(742.7923980000)
- vector_i[7368]['s'].append(703.3277805000)
- vector_b[7368]={'n':[],'s':[]}
- vector_b[7368]['n'].append(85.6553090000)
- vector_b[7368]['s'].append(85.7558440000)
- sort[7368]={'n':[],'s':[]}
- sort[7368]['n'].append(0.7025800000)
- sort[7368]['s'].append(0.7094055000)
- set_i[22143]={'n':[],'s':[]}
- set_i[22143]['n'].append(59.7905750000)
- set_i[22143]['s'].append(60.4921435000)
- vector_i[22143]={'n':[],'s':[]}
- vector_i[22143]['n'].append(700.2722935000)
- vector_i[22143]['s'].append(686.1256720000)
- vector_b[22143]={'n':[],'s':[]}
- vector_b[22143]['n'].append(86.0514560000)
- vector_b[22143]['s'].append(87.0559570000)
- sort[22143]={'n':[],'s':[]}
- sort[22143]['n'].append(2.3735980000)
- sort[22143]['s'].append(2.4882290000)
- set_i[66548]={'n':[],'s':[]}
- set_i[66548]['n'].append(197.5712895000)
- set_i[66548]['s'].append(196.5581185000)
- vector_i[66548]={'n':[],'s':[]}
- vector_i[66548]['n'].append(709.2545550000)
- vector_i[66548]['s'].append(683.3994835000)
- vector_b[66548]={'n':[],'s':[]}
- vector_b[66548]['n'].append(86.0151530000)
- vector_b[66548]['s'].append(92.0718050000)
- sort[66548]={'n':[],'s':[]}
- sort[66548]['n'].append(7.8001840000)
- sort[66548]['s'].append(8.0440845000)
- set_i[200000]={'n':[],'s':[]}
- set_i[200000]['n'].append(675.8519965000)
- set_i[200000]['s'].append(673.8654970000)
- vector_i[200000]={'n':[],'s':[]}
- vector_i[200000]['n'].append(693.5815290000)
- vector_i[200000]['s'].append(710.8254870000)
- vector_b[200000]={'n':[],'s':[]}
- vector_b[200000]['n'].append(86.2450735000)
- vector_b[200000]['s'].append(97.1322280000)
- sort[200000]={'n':[],'s':[]}
- sort[200000]['n'].append(25.4888220000)
- sort[200000]['s'].append(27.0257605000)
- from math import *
- import numpy as np
- from matplotlib.pyplot import *
- class gr:
- def __init__(self):
- self.x1=self.x2=self.y1=self.y2=0
- self.sm=[]
- self.nm=[]
- def add(self,*a,**b):
- #print(len(self.sm),end=' - ')
- if len(self.sm)==0:
- self.x1=min(a[0])
- self.x2=max(a[0])
- self.y1=min(a[1])
- self.y2=max(a[1])
- else:
- self.x1=min(self.x1,min(a[0]))
- self.x2=max(self.x2,max(a[0]))
- self.y1=min(self.y1,min(a[1]))
- self.y2=max(self.y2,max(a[1]))
- '''if 'name' in b:
- self.nm+=[b['name']]
- if 'label' in b:
- b['label']=b['label']+':'+b['name']
- else:
- b['label']=b['name']
- del b['name']
- else:
- self.nm+=[None]'''
- self.sm+=[(a,b)]
- def floor(self,x,a):
- return floor(x/a)*a
- def ceil(self,x,a):
- return ceil(x/a)*a
- def show(self,n=0):
- mpl_fig = figure(n)
- subplots_adjust(left=0.06, bottom=0.06, right=0.94, top=0.94, wspace=None, hspace=None)
- ax = mpl_fig.add_subplot(111)
- ax.set_yscale ('log')
- ax.set_xscale ('log')
- def onclick(event):
- print(event)
- #print(event.width,event.height)
- #if (event.width>1000):
- #print("asf")
- #grid(True,which='major', linestyle='-', linewidth=3.5, color='grey',axis='none')
- #grid(False,which='major', linestyle='-', linewidth=2, color='grey',axis='both')
- def ondraw(event):
- print(event)
- ax=gca()
- print(ax.get_xlim(),ax.get_ylim())
- cid = mpl_fig.canvas.mpl_connect('resize_event', onclick)
- cid = mpl_fig.canvas.mpl_connect('draw_event', ondraw)
- def rnd(x):
- c=[0,0.30102999566398125,-0.30102999566398125]
- n=[d+round(x-d) for d in c]
- m=0
- for u in range(1,len(n)):
- if abs(x-n[u])<abs(x-n[m]):
- m=u
- return n[m]
- #print(self.x1,self.x2,self.y1,self.y2)
- st=0.2
- ex=rnd(log(self.x2-self.x1,10)-1.34+st)
- ey=rnd(log(self.y2-self.y1,10)-1.1+st)
- dx=10**ex
- dy=10**ey
- x1=self.floor(self.x1,dx)
- x2=self.ceil(self.x2,dx)
- y1=self.floor(self.y1,dy)
- y2=self.ceil(self.y2,dy)
- #ax.set_aspect(dx/dy)
- #ax.set_xticks(np.arange(x1, x2+dx, dx), minor=False)
- #ax.set_yticks(np.arange(y1, y2+dy, dy), minor=False)
- grid(True,which='major', linestyle='-', linewidth=0.5, color='grey',axis='both')
- #print (dx)
- #print(dy)
- dx*=0.1
- dy*=0.1
- #print(dx)
- #print(dy)
- x1=self.floor(self.x1,dx)
- x2=self.ceil(self.x2,dx)
- y1=self.floor(self.y1,dy)
- y2=self.ceil(self.y2,dy)
- #ax.set_xticks(np.arange(x1, x2+dx, dx), minor=True)
- #ax.set_yticks(np.arange(y1, y2+dy, dy), minor=True)
- grid(True,which='minor', linestyle=':', linewidth=0.4, color=(0.7,0.7,0.7),axis='both')
- for j in range(len(self.sm)):
- c=self.sm[j]
- plot(*c[0],**c[1])
- legend()
- #fig, ax = subplots()
- #show()
- class a_gr:
- ff={}
- def add(self,n,*b,**p):
- if not (n in self.ff):
- self.ff[n]=gr()
- self.ff[n].add(*b,**p)
- def show(self):
- t=1
- for a in self.ff:
- self.ff[a].show(t)
- t+=1
- show()
- A=a_gr()
- xx=np.array(list(set_i))
- for ci,di in zip(['set_i','vector_i','vector_b','sort'],[(0,0.6,0),(1.0,0.6,0),(0.7,0.7,0),(0.5,0,1.0)]):
- c=eval(ci)
- yy1=[]
- yy2=[]
- for x in xx:
- yy1.append(c[x]['n'][0])
- yy2.append(c[x]['s'][0])
- A.add(0,xx,np.array(yy1),'.-',color=di,label=ci)
- A.add(0,xx,np.array(yy2),'.-',color=di,label=ci+'_s')
- A.show()
- '''set_i={}
- vector_i={}
- vector_b={}
- set_i[10]={'n':[],'s':[]}
- set_i[10]['n'].append(0.0000000000)
- set_i[10]['s'].append(0.0000000000)'''
Add Comment
Please, Sign In to add comment