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
