function
   petsc-3.7.7 2017-09-25
    Set = PetscReadBinaryMatlab(filename)
  matlabHeader = ['%$$ '; '#$$ ']; 
  matlabHeaderLen = size(matlabHeader,2);
  if (isempty(strfind(filename,'.info')))
     filename = [filename,'.info'];
  end
  fid=fopen(filename,'r');
  if (fid == -1)
     error(sprintf('PetscReadBinaryMatlab: cannot load file %s',filename))
  end
  str=fgets(fid);
  while (ischar(str))
     
     if strncmp(str,matlabHeader(1,:),matlabHeaderLen) || strncmp(str,matlabHeader(2,:),matlabHeaderLen)
	 str = str(1+matlabHeaderLen:end);
     	 
	 if strfind(str,'fopen(Set.filename')
	    str = 'fd = PetscOpenFile(Set.filename);';
	 elseif strfind(str,'if (fd < 0) error'); 
	    str = '%';
	 elseif strfind(str,'fclose(fd)');
	    str = 'close(fd);';
	 end
	 eval(str);
     end
     str=fgets(fid);
  end
  fclose(fid);
  return