octave_packages/secs2d-0.0.8/Utilities/Updemesh.m

   1 function Updemesh(varargin);
   2
   3 % Udrawedge(mesh);
   4
   5
   6 % This file is part of
   7 %
   8 %            SECS2D - A 2-D Drift--Diffusion Semiconductor Device Simulator
   9 %         -------------------------------------------------------------------
  10 %            Copyright (C) 2004-2006  Carlo de Falco
  11 %
  12 %
  13 %
  14 %  SECS2D is free software; you can redistribute it and/or modify
  15 %  it under the terms of the GNU General Public License as published by
  16 %  the Free Software Foundation; either version 2 of the License, or
  17 %  (at your option) any later version.
  18 %
  19 %  SECS2D is distributed in the hope that it will be useful,
  20 %  but WITHOUT ANY WARRANTY; without even the implied warranty of
  21 %  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  22 %  GNU General Public License for more details.
  23 %
  24 %  You should have received a copy of the GNU General Public License
  25 %  along with SECS2D; If not, see <http://www.gnu.org/licenses/>.
  26
  27
  28
  29
  30 if nargin == 1
  31
  32         Udrawedge(varargin{1});
  33
  34 elseif nargin == 2
  35
  36         mesh = varargin{1};
  37         u = varargin{2};
  38
  39         dataname = mktemp("/tmp",".dx");
  40         scriptname = mktemp("/tmp",".net");
  41
  42         UDXoutput2Ddata(dataname,mesh.p,mesh.t,u,'u',0,1,1);
  43
  44
  45         showmesh = file_in_path(path,"Urubbersheet.net");
  46
  47         system (["cp " showmesh " " scriptname]);
  48         system (["sed -i \'s|__FILE__DX__|" dataname "|g\' " scriptname]);
  49
  50         command = ["dx -program " scriptname " -execute -image >& /dev/null &"];
  51         system(command);
  52
  53 else
  54
  55         fprintf(1,"wrong number of parameters\n\n")
  56
  57 end
  58
  59 endfunction
  60
  61 function filename = mktemp (direct,ext);
  62
  63 if (~exist(direct,"dir"))
  64   error("trying to save temporary file to non existing directory")
  65 end
  66
  67 done=false;
  68
  69 while ~done
  70   filename = [direct,"/SECS2D.",num2str(floor(rand*1e7)),ext];
  71   if ~exist(filename,"file")
  72     done =true;
  73   end
  74 end
  75 endfunction