octave_packages/fpl-1.2.0/FPL2pdesurf.m

   1 ## Copyright (C) 2004-2008,2009  Carlo de Falco, Massimiliano Culpo
   2 ##
   3 ##  This file is part of
   4 ##
   5 ##                   FPL - Fem PLotting package for octave
   6 ##
   7 ##  FPL is free software; you can redistribute it and/or modify
   8 ##  it under the terms of the GNU General Public License as published by
   9 ##  the Free Software Foundation; either version 2 of the License, or
  10 ##  (at your option) any later version.
  11 ##
  12 ##  FPL is distributed in the hope that it will be useful,
  13 ##  but WITHOUT ANY WARRANTY; without even the implied warranty of
  14 ##  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15 ##  GNU General Public License for more details.
  16 ##
  17 ##  You should have received a copy of the GNU General Public License
  18 ##  along with FPL; If not, see <http://www.gnu.org/licenses/>.
  19 ##
  20 ##
  21 ##  AUTHORS:
  22 ##  Carlo de Falco <cdf _AT_ users.sourceforge.net>
  23 ##
  24 ##  Culpo Massimiliano
  25 ##  Bergische Universitaet Wuppertal
  26 ##  Fachbereich C - Mathematik und Naturwissenschaften
  27 ##  Arbeitsgruppe fuer Angewandte MathematD-42119 Wuppertal  Gaussstr. 20
  28 ##  D-42119 Wuppertal, Germany
  29
  30 ## -*- texinfo -*-
  31 ## @deftypefn {Function File} {} FPL2pdesurf (@var{mesh}, @
  32 ## @var{u} [ @var{property}, @var{value} ...])
  33 ##
  34 ## plots the scalar field @var{u}
  35 ## defined on the triangulation @var{mesh} using opendx.
  36 ##
  37 ## options (default value):
  38 ## @itemize @minus
  39 ## @item data_dep ("positions") defines wether data depends on
  40 ## positions or connections
  41 ## @item plot_field ("scalar") defines wether to plot the scalar field
  42 ## itself or its gradient
  43 ## @end itemize
  44 ##
  45 ## @seealso{MSH2Mgmsh, MSH2Mstructmesh}
  46 ## @end deftypefn
  47
  48 function FPL2pdesurf(varargin)
  49
  50   data_dep   = "positions";
  51   plot_field = "scalar";
  52
  53   if nargin == 1
  54     FPL2showmesh(varargin{1});
  55   elseif nargin == 2
  56     mesh = varargin{1};
  57     u = varargin{2};
  58   elseif ( (nargin > 2) && (rem(nargin,2)==0) )
  59     mesh = varargin{1};
  60     u = varargin{2};
  61     for ii=3:2:nargin
  62       eval([ varargin{ii} " = """ varargin{ii+1} """;" ]);
  63     endfor
  64   else
  65     keyboard ,error(["wrong number of parameters " num2str (nargin)])
  66   endif
  67
  68   dataname = mktemp("/tmp",".dx");
  69   scriptname = mktemp("/tmp",".net");
  70
  71   FPL2dxoutputdata(dataname,mesh.p,mesh.t,u,'u',0,1,1);
  72
  73   switch plot_field
  74     case {"scalar","scal"}
  75       showmesh = file_in_path(path,"FPL2coloredrubbersheet.net");
  76     case {"gradient","grad"}
  77       showmesh = file_in_path(path,"FPL2coloredgradient.net");
  78     otherwise
  79       error ([ "incorrect value " plot_field " for option plot_field "])
  80   endswitch
  81
  82   system (["cp " showmesh " " scriptname]);
  83   system (["sed -i \'s|__FILE__DX__|" dataname "|g\' " scriptname]);
  84
  85   switch data_dep
  86     case {"positions","continuous","interpolate","P1"}
  87       system (["sed -i \'s|__DATA_DEPENDENCY__|positions|g\' " scriptname]);
  88     case {"connections","discontinuous","P0"}
  89       system (["sed -i \'s|__DATA_DEPENDENCY__|positions|g\' " scriptname]);
  90     otherwise
  91       error ([ "incorrect value " data_dep " for option data_dep "])
  92   endswitch
  93
  94   ##command = ["dx -noConfirmedQuit  -noImageRWNetFile -program " scriptname " -execute -image  >& /dev/null & "];
  95   command  = ["dx -noConfirmedQuit  -program " scriptname " -execute -image  >& /dev/null & "];
  96
  97   system(command);
  98
  99 endfunction
 100
 101 function filename = mktemp (direct,ext);
 102
 103   if (~exist(direct,"dir"))
 104     error("trying to save temporary file to non existing directory")
 105   endif
 106
 107   done=false;
 108
 109   while ~done
 110     filename = [direct,"/FPL.",num2str(floor(rand*1e7)),ext];
 111     if ~exist(filename,"file")
 112       done =true;
 113     endif
 114   endwhile
 115
 116 endfunction
 117
 118 %!shared msh
 119 %!demo
 120 %! msh.p = [0 0; 1 0; 1 1; 0 1].';
 121 %! msh.t = [1 2 3 1; 1 3 4 1].';
 122 %! msh.e = [1 2 0 0 1 0 1; 2 3 0 0 2 0 1; 3 4 0 0 3 0 1; 4 1 0 0 4 0 1].';
 123 %! u = [0 1].';
 124 %! FPL2pdesurf (msh, u, "data_dep", "connections");
 125 %! s = input ("do you see a square divided into two triangles with a diagonal red-to-blue gradient color (if you see an empty plot try ctrl-F)? (y/n): " ,"s");
 126 %! assert(s, "y")
 127 %!demo
 128 %! v = [0 0 1 1].';
 129 %! FPL2pdesurf (msh, v);
 130 %! s = input ("do you see a square divided into two triangles with a vertical red-to-blue gradient color (if you see an empty plot try ctrl-F)? (y/n): " ,"s");