Mercurial > forge

comparison extra/fpl/deprecated/FPL2pdesurf.m @ 12671:20e8aca47b2c octave-forge

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
prepare for release
author cdf
date Mon, 17 Aug 2015 10:19:39 +0000
parents extra/fpl/inst/deprecated/FPL2pdesurf.m@2748190086ad
children
comparison
equal deleted inserted replaced
12670:d68da2f2417b 12671:20e8aca47b2c
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 3 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 %!test
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 %!test
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");