nbuffer = 10;
ndraws=3*nbuffer;

n=3;
yhatw = zeros(n^2,nbuffer);

wdraws=0;
for draws=1:ndraws
	tmp = draws*ones(n);
	yhatw(:,draws-wdraws) = tmp(:);

	if ~mod(draws,nbuffer)
      wdraws=draws

   	fwriteid = fopen('outyhat.bin','a');
      count = fwrite(fwriteid,yhatw,'single');   % 'single' or 'double'
   	status = fclose('all');
	end
end