def drawmins(G,n,lb,ub,SUPPORT=300,mode = [ESutils.SUPPORT_SLICELCB],SLICELCB_PARA=1.): #draw support points W = sp.vstack([ESutils.draw_support(G, lb,ub,SUPPORT/len(mode),m, para = SLICELCB_PARA) for m in mode]) if False: print 'how did I get here' plt.figure() plt.plot(W[:,0],W[:,1],'g.') plt.show() #draw in samples on the support print "drawing mins from support" R = ESutils.draw_min(G,W,n) if False: plt.figure() plt.plot(R[:,0],R[:,1],'r.') plt.show() #plt.show() return R
lb = sp.array([-2.,-2.,-2.]) ub = sp.array([2.,2.,2.]) #MLEH = GPdc.searchMLEhyp(X,Y,S,D,lb,ub,GPdc.SQUEXP,mx=10000) G = GPdc.GPcore(X,Y,S,D,GPdc.kernel(GPdc.SQUEXP,2,sp.array([1.5,0.15,0.15]))) #np=180 #sup = sp.linspace(-1,1,np) #Dp = [[sp.NaN]]*np #Xp = sp.vstack([sp.array([i]) for i in sup]) #[m,v] = G.infer_diag(Xp,Dp) #a0.plot(sup,m) #sq = sp.sqrt(v) #a0.fill_between(sup, sp.array(m-2.*sq).flatten(), sp.array(m+2.*sq).flatten(), facecolor='lightblue',edgecolor='lightblue') Z = ESutils.draw_support(G, sp.array([-1.,-1.]),sp.array([1.,1.]),500,ESutils.SUPPORT_SLICEEI) np=30 R = ESutils.draw_min(G,Z,500) Z2 = ESutils.draw_support(G, sp.array([-1.,-1.]),sp.array([1.,1.]),500,ESutils.SUPPORT_LAPAPR,para=np) R2 = ESutils.draw_min(G,Z2,500) plt.figure() for i in xrange(Z.shape[0]): plt.plot(Z[i,0],Z[i,1],'r.') for i in xrange(Z2.shape[0]): plt.plot(Z2[i,0],Z2[i,1],'b.') for i in xrange(R.shape[0]): plt.plot(R[i,0],R[i,1],'gx') plt.axis([-1,1,-1,1]) ng = 30 A = sp.empty([ng,ng]) for i in xrange(ng): for j in xrange(ng):
def drawmins_inplane(G,n,lb,ub,axis,value, SUPPORT=300, mode=ESutils.SUPPORT_SLICELCB, SLICELCB_PARA=1.): W = sp.vstack([ESutils.draw_support_inplane(G, lb,ub,SUPPORT/len(mode),m, axis,value,para = SLICELCB_PARA) for m in mode]) #draw in samples on the support R = ESutils.draw_min(G,W,n) return R
#MLEH = GPdc.searchMLEhyp(X,Y,S,D,lb,ub,GPdc.SQUEXP,mx=10000) G = GPdc.GPcore(X, Y, S, D, GPdc.kernel(GPdc.SQUEXP, 2, sp.array([1.5, 0.15, 0.15]))) #np=180 #sup = sp.linspace(-1,1,np) #Dp = [[sp.NaN]]*np #Xp = sp.vstack([sp.array([i]) for i in sup]) #[m,v] = G.infer_diag(Xp,Dp) #a0.plot(sup,m) #sq = sp.sqrt(v) #a0.fill_between(sup, sp.array(m-2.*sq).flatten(), sp.array(m+2.*sq).flatten(), facecolor='lightblue',edgecolor='lightblue') plt.figure() Z = ESutils.draw_support_inplane(G, sp.array([-1., -1.]), sp.array([1., 1.]), 400, ESutils.SUPPORT_LAPAPR, 1, 0.5) R = ESutils.draw_min(G, Z, 100) for i in xrange(Z.shape[0]): plt.plot(Z[i, 0], Z[i, 1], 'r.') for i in xrange(R.shape[0]): plt.plot(R[i, 0], R[i, 1], 'gx') Z = ESutils.draw_support_inplane(G, sp.array([-1., -1.]), sp.array([1., 1.]), 400, ESutils.SUPPORT_LAPAPR, 0, 0.5) R = ESutils.draw_min(G, Z, 100) print Z.shape print R.shape for i in xrange(Z.shape[0]): plt.plot(Z[i, 0], Z[i, 1], 'r.') for i in xrange(R.shape[0]): plt.plot(R[i, 0], R[i, 1], 'gx') plt.axis([-1, 1, -1, 1])
lb = sp.array([-2.,-2.,-2.]) ub = sp.array([2.,2.,2.]) #MLEH = GPdc.searchMLEhyp(X,Y,S,D,lb,ub,GPdc.SQUEXP,mx=10000) G = GPdc.GPcore(X,Y,S,D,GPdc.kernel(GPdc.SQUEXP,2,sp.array([1.5,0.15,0.15]))) #np=180 #sup = sp.linspace(-1,1,np) #Dp = [[sp.NaN]]*np #Xp = sp.vstack([sp.array([i]) for i in sup]) #[m,v] = G.infer_diag(Xp,Dp) #a0.plot(sup,m) #sq = sp.sqrt(v) #a0.fill_between(sup, sp.array(m-2.*sq).flatten(), sp.array(m+2.*sq).flatten(), facecolor='lightblue',edgecolor='lightblue') plt.figure() Z = ESutils.draw_support_inplane(G, sp.array([-1.,-1.]),sp.array([1.,1.]),400,ESutils.SUPPORT_LAPAPR,1,0.5) R = ESutils.draw_min(G,Z,100) for i in xrange(Z.shape[0]): plt.plot(Z[i,0],Z[i,1],'r.') for i in xrange(R.shape[0]): plt.plot(R[i,0],R[i,1],'gx') Z = ESutils.draw_support_inplane(G, sp.array([-1.,-1.]),sp.array([1.,1.]),400,ESutils.SUPPORT_LAPAPR,0,0.5) R = ESutils.draw_min(G,Z,100) print Z.shape print R.shape for i in xrange(Z.shape[0]): plt.plot(Z[i,0],Z[i,1],'r.') for i in xrange(R.shape[0]): plt.plot(R[i,0],R[i,1],'gx') plt.axis([-1,1,-1,1]) ng = 30