plt.imshow(true_mod, **ima)
plt.subplot(1,2,2)
plt.imshow(tim.data, **ima)
plt.title('S/N %f, r_e %f' % (sn, re))
ps.savefig()
## run one_blob code? Or shortcut?
src = PointSource(RaDecPos(0., 0.),
NanoMaggies(**{band: flux}))
nblob,iblob,Isrcs = 0, 1, np.array([0])
brickwcs = wcs
bx0, by0, blobw, blobh = 0, 0, sz, sz
blobmask = np.ones((sz,sz), bool)
timargs = [(tim.data, tim.getInvError(), tim.wcs, tim.wcs.wcs,
tim.getPhotoCal(), tim.getSky(), tim.psf,
'tim', 0, sz, 0, sz, band, sig1,
tim.modelMinval, None)]
srcs = [src]
bands = band
plots,psx = False, None
simul_opt, use_ceres, hastycho = False, False, False
X = (nblob, iblob, Isrcs, brickwcs, bx0, by0, blobw, blobh,
blobmask, timargs, srcs, bands, plots, psx, simul_opt,
use_ceres, hastycho)
R = one_blob(X)
#print('Got:', R)
print('Sources:', R.sources)
print('Dchisqs:', R.dchisqs)
plt.imshow(true_mod, **ima)
plt.subplot(1,2,2)
plt.imshow(tim.data, **ima)
plt.title('S/N %f, r_e %f' % (sn, re))
ps.savefig()
## run one_blob code? Or shortcut?
src = PointSource(RaDecPos(0., 0.),
NanoMaggies(**{band: flux}))
nblob,iblob,Isrcs = 0, 1, np.array([0])
brickwcs = wcs
bx0, by0, blobw, blobh = 0, 0, sz, sz
blobmask = np.ones((sz,sz), bool)
timargs = [(tim.data, tim.getInvError(), tim.wcs, tim.wcs.wcs,
tim.getPhotoCal(), tim.getSky(), tim.psf,
'tim', 0, sz, 0, sz, band, sig1,
tim.modelMinval, None)]
srcs = [src]
bands = band
plots,psx = False, None
simul_opt, use_ceres, hastycho = False, False, False
X = (nblob, iblob, Isrcs, brickwcs, bx0, by0, blobw, blobh,
blobmask, timargs, srcs, bands, plots, psx, simul_opt,
use_ceres, hastycho)
R = one_blob(X)
#print('Got:', R)
print('Sources:', R.sources)
print('Dchisqs:', R.dchisqs)
