% eml_sps_test.m% compare aspire and matlab versions of E-ML-SPS% Copyright Mar 2000, Jeff Fessler, The University of Michigan%% generate data%if ~isvar('yi'),	disp 'setup eml_sps_test'	f.dir	= test_dir;	f.wtf	= [f.dir 't,g.wtf'];	f.wtr	= strrep(f.wtf, 'wtf', 'wtr');	f.yi	= [f.dir 'yi.fld'];	f.ci	= [f.dir 'ci.fld'];	f.ri	= [f.dir 'ri.fld'];	f.mask	= [f.dir 'mask.fld'];	em_test_setuppromptend%% matlab iterations%if ~isvar('xmat'),	disp 'matlab E-ML-SPS'	f.niter = 9;	xinit = ones(size(xtrue));	xmat = eml_sps(xinit(mask(:)), G', yi, ci, ri, f.niter, 'oc');	xmat = embed(xmat, mask);	clf, im(xmat, 'matlab')promptendif ~has_aspire, return, end%% aspire iterations%if ~isvar('xasp'),	disp 'aspire E-ML-SPS'	f.init	= [f.dir 't,init.fld'];	f.out	= [f.dir 't,out.fld'];	fld_write(f.init, xinit, '-nocheck')	if (exist(f.out) == 2), delete(f.out), end	f.scaleinit = 0;	f.saver = 'stack,1';	f.alg = 'ospsc,oc,1,1,1,0';	f.penal = sprintf('%g,quad,0,-', -100);	f.method = sprintf('@%d@%s@%s', f.niter-1, f.alg, f.penal);	f.fitype = ['2z@' f.wtr '@-'];	f.com = sprintf(['i -chat 0 empl3 %s %s  %s %s 1 %s 1 %s -' ...			' %s %s 0 1 1e30 0 -'], ...		f.out, f.init, f.yi, f.ci, f.ri, f.fitype, ...		f.method, f.saver);	os_run(f.com)	xasp = double(fld_read(f.out));endif 1	clf, im(221, xmat, 'xhat matlab'), cbar	im(222, xasp, 'xhat aspire'), cbar	im(223, (xasp-xmat)/max(xmat(:)), 'aspire-matlab'), cbar	t = vcorrcoef(xasp, xmat);	printf('corr. %g,%g', t, t-1)	t1 = eql_obj(xmat, G, yi(:), ci(:), ri(:), [], mask);	t2 = eql_obj(xasp, G, yi(:), ci(:), ri(:), [], mask);	subplot(224)	plot(0:f.niter-1, t1-t1(1), '-o', 0:f.niter-1, t2-t1(1), '-x')	legend('mat', 'asp', 4)promptend