gmevallog.m

matlab有限元分析工具,比经较全面的一个手册,请大家下载呀

function gmevallog(GELinfile, GELoutfile)
% gmevallog(infile, outfile);
% This routine evaluates a script file (first argument) and
% places the result in the output file (second argument). 
% This command will usually correctly handle block statements
% like for and if, but there may be some squirrel cases that escape.


GELifid = fopen(GELinfile,'rt');
GELallcmds = setstr(fread(GELifid))';
fclose(GELifid);

GELofid = fopen(GELoutfile,'wt');
GELversion = version;
GELvnum = str2num(GELversion(1:3));

if GELvnum < 5.3
  error('gmevallog not implemented for Matlab version < 5.3');
end

while 1
  [GELcmd, GELrest] = strtok(GELallcmds, setstr(10));
  GELorigcmd = GELcmd;
  fprintf(GELofid, '%s\n', gm_evalclock);
  if length(GELcmd) > 0
    if GELcmd(end) == 13
      GELcmd = GELcmd(1:end - 1);
    end
  end
  GELexec = 0;
  while ~GELexec
    GELtryit = 1;
    if length(GELcmd) >= 3 
      if all(GELcmd(end - 2 : end) == '...')
        if length(GELrest) < 2
          disp('errcall 1')
          GELerr1(GELcmd, 'Continuation mark ''...'' at end of file')
        end
        [GELf,GELrest2] = strtok(GELrest, setstr(10));
        if length(GELf) > 0
          if GELf(end) == 13
            GELf = GELf(1:end - 1);
          end
        end
        GELcmd = [GELcmd(1:end - 3), ' ', GELf];
        GELorigcmd = [GELorigcmd, setstr(10), GELf];
        GELrest = GELrest2;
        GELtryit = 0;
      end
    end
    if GELtryit
      try
        if length(GELcmd) > 0
          GELrval = evalc(GELcmd);
        else
          GELrval = '';
        end
        GELexec = 1;
      catch
        GELlasterr = lasterr;
        if length(findstr(GELlasterr, 'End of file reached inside')) > 0
          if length(GELrest) < 2
            disp('errcall 2')
            GELerr1(GELcmd, GELlasterr)
          end
          [GELf,GELrest2] = strtok(GELrest, setstr(10));
          if length(GELf) > 0
            if GELf(end) == 13
              GELf = GELf(1:end - 1);
            end
          end
          GELcmd = [GELcmd, ', ', GELf];
          GELorigcmd = [GELorigcmd, setstr(10), GELf];
          GELrest = GELrest2;
        elseif length(findstr(GELlasterr, 'Unterminated bracketed')) > 0
          if length(GELrest) < 2
            disp('errcall 3')
            GELerr1(GELcmd, GELlasterr)
          end
          [GELf, GELrest2] = strtok(GELrest, setstr(10));
          if length(GELf) > 0
            if GELf(end) == 13
              GELf = GELf(1:end - 1);
            end
          end
          GELcmd = [GELcmd, '; ', GELf];
          GELorigcmd = [GELorigcmd, setstr(10), GELf];
          GELrest = GELrest2;
        else
          disp('errcall 4')
          GELcmd
          GELlasterr
          GELerr1(GELcmd, GELlasterr)
        end
      end
    end
  end
  fprintf(GELofid, '-- >> %s\n', GELorigcmd);
  fprintf(GELofid, '%s\n', GELrval);
  if length(GELrest) > 0
    GELallcmds = GELrest;
  else
    break
  end
end
fclose(GELofid);




function GELerr1(cmd, errmsg)
disp(sprintf('during execution of command: %s', cmd))
disp(sprintf('encountered error: %s',errmsg))
error('error encountered during gmevallog')




% ------------------------------------------------------------------
% Copyright (c) 1999 by Cornell University.  All rights reserved.
% See the accompanying file 'Copyright' for authorship information,
% the terms of the license governing this software, and disclaimers
% concerning this software.
% ------------------------------------------------------------------
% This file is part of the QMG software.  
% Version 2.0 of QMG, release date September 3, 1999
% ------------------------------------------------------------------