nbuffer = 10;
ndraws=3*nbuffer;

n=3;
yhatw = zeros(n^2,ndraws);
seldraw = [3 7];

[fid,message]=fopen('outyhat.bin')
%[xd,count]=fread(fid,[n^2,length(seldraw)],'double');   % (1) working
[xd,count]=fread(fid,inf,'single');    % (2) working.  with 'single' or 'double'
fid
message
count
%xdd=reshape(xd,n^2,length(seldraw))   % (1) working
xdd=reshape(xd,n^2,ndraws)  % (2) working