subplot(1,2,1)
for i in range(0,N_samps):
f=GPSampler.f.trace()[i](x)
plot(x,GPSampler.f.trace()[i](x))
mid_traces.append(f[len(f)/2])
plot(fmesh,GPSampler.d.value,"k.",markersize=16)
axis([x.min(),x.max(),-5.,10.])
title("Some samples of f")
After Change
subplot(1,2,1)
for i in range(0,N_samps):
f=GPSampler.submod.f.trace()[i](x)
plot(x,f)
mid_traces.append(f[len(f)/2])
plot(fmesh,GPSampler.d.value,"k.",markersize=16)
axis([x.min(),x.max(),-5.,10.])
title("Some samples of f")