Mercurial > fem-fenics-eugenio
annotate inst/ufl.m @ 268:61830a4f9ab9
Improve formatting
author | Eugenio Gianniti <eugenio.gianniti@mail.polimi.it> |
---|---|
date | Thu, 14 Aug 2014 12:26:55 +0200 |
parents | 072aee55bb75 |
children | f4d6ae912a08 |
rev | line source |
---|---|
220 | 1 ## Copyright (C) 2014 Eugenio Gianniti <eugenio.gianniti@mail.polimi.it> |
2 ## | |
3 ## This program is free software; you can redistribute it and/or modify it under | |
4 ## the terms of the GNU General Public License as published by the Free Software | |
5 ## Foundation; either version 3 of the License, or (at your option) any later | |
6 ## version. | |
7 ## | |
8 ## This program is distributed in the hope that it will be useful, but WITHOUT | |
9 ## ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or | |
10 ## FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more | |
11 ## details. | |
12 ## | |
13 ## You should have received a copy of the GNU General Public License along with | |
14 ## this program; if not, see <http://www.gnu.org/licenses/>. | |
15 | |
268
61830a4f9ab9
Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
255
diff
changeset
|
16 |
220 | 17 ## -*- texinfo -*- |
18 ## @deftypefn {Function File} ufl [command] line | |
19 ## | |
20 ## Writes the given line to a ufl file. Optional commands are: | |
21 ## @itemize @bullet | |
22 ## @item start: opens an ufl file. It is possible to specify a name for it, | |
23 ## otherwise 'default.ufl' will be used. | |
225
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
24 ## @item end: closes the ufl file and compiles it. It is possible to |
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
25 ## specify its @var{content}, otherwise a complete problem is assumed. |
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
26 ## Possible options are: BilinearForm, Functional, FunctionSpace, LinearForm |
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
27 ## and Problem. |
220 | 28 ## @end itemize |
29 ## | |
225
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
30 ## @example |
226
043120f071de
Correct help message for ufl
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
225
diff
changeset
|
31 ## @code{ufl start linform} |
043120f071de
Correct help message for ufl
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
225
diff
changeset
|
32 ## @code{ufl element = FiniteElement ("Lagrange", triangle, 1)} |
043120f071de
Correct help message for ufl
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
225
diff
changeset
|
33 ## @code{ufl v = TestFunction (element)} |
043120f071de
Correct help message for ufl
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
225
diff
changeset
|
34 ## @code{ufl f = Coefficient (element)} |
043120f071de
Correct help message for ufl
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
225
diff
changeset
|
35 ## @code{ufl L = f*v*dx} |
043120f071de
Correct help message for ufl
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
225
diff
changeset
|
36 ## @code{ufl end LinearForm} |
225
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
37 ## @end example |
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
38 ## |
226
043120f071de
Correct help message for ufl
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
225
diff
changeset
|
39 ## This way you can write a ufl file and compile it obtaining a LinearForm. |
043120f071de
Correct help message for ufl
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
225
diff
changeset
|
40 ## |
228
c74309df3489
State in the ufl.m help message that leading brackets cause an error
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
227
diff
changeset
|
41 ## Be aware that a ufl line starting with brackets will result in a syntax |
c74309df3489
State in the ufl.m help message that leading brackets cause an error
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
227
diff
changeset
|
42 ## error. To avoid this, it is necessary to enclose them and their content |
c74309df3489
State in the ufl.m help message that leading brackets cause an error
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
227
diff
changeset
|
43 ## in quotes or apices: |
c74309df3489
State in the ufl.m help message that leading brackets cause an error
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
227
diff
changeset
|
44 ## @example |
c74309df3489
State in the ufl.m help message that leading brackets cause an error
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
227
diff
changeset
|
45 ## @code{ufl "(sigma, u)" = TrialFunctions(W)} |
c74309df3489
State in the ufl.m help message that leading brackets cause an error
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
227
diff
changeset
|
46 ## @end example |
c74309df3489
State in the ufl.m help message that leading brackets cause an error
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
227
diff
changeset
|
47 ## |
225
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
48 ## @seealso{import_ufl_Problem, import_ufl_BilinearForm, import_ufl_Functional, |
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
49 ## import_ufl_FunctionSpace, import_ufl_LinearForm} |
220 | 50 ## @end deftypefn |
51 | |
52 function ufl (varargin) | |
53 | |
54 persistent fid = -1; | |
225
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
55 persistent filename = "default.ufl"; |
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
56 persistent opts = {"BilinearForm", |
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
57 "Functional", |
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
58 "FunctionSpace", |
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
59 "LinearForm", |
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
60 "Problem"}; |
220 | 61 |
62 if (numel (varargin) < 1) | |
255
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
63 if (is_master_node () && ! is_valid_file_id (fid)) |
220 | 64 print_usage (); |
65 endif | |
66 elseif (! all (cellfun ("ischar", varargin))) | |
255
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
67 if (is_master_node ()) |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
68 error ("ufl: all arguments should be strings"); |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
69 endif |
220 | 70 elseif (strcmpi (varargin{1}, "start")) |
255
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
71 if (is_master_node ()) |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
72 if (is_valid_file_id (fid)) |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
73 error ("ufl: a file is already open"); |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
74 else |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
75 if (numel (varargin) > 1) |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
76 filename = varargin{2}; |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
77 if (isempty (regexpi (filename, ".ufl$", "match"))) |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
78 filename = [filename, ".ufl"]; |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
79 else |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
80 filename = [filename(1:end-4), ".ufl"]; |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
81 endif |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
82 endif |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
83 [~, err, ~] = stat (filename); |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
84 if (err == 0) |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
85 error (["ufl: a file named ", filename, " already exists"]); |
224
a108a0bfe201
Ensure that ufl opens a file with .ufl extension
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
223
diff
changeset
|
86 else |
255
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
87 fid = fopen (filename, "w"); |
220 | 88 endif |
89 endif | |
255
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
90 if (! is_valid_file_id (fid)) |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
91 error (["ufl: could not open file ", filename]); |
223
8b159045242f
ufl errors out if filename exists when start command is issued
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
222
diff
changeset
|
92 endif |
220 | 93 endif |
94 elseif (strcmpi (varargin{1}, "end")) | |
255
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
95 if (is_master_node ()) |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
96 if (! is_valid_file_id (fid)) |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
97 error ("ufl: no open file"); |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
98 else |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
99 fclose (fid); |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
100 fid = -1; |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
101 endif |
220 | 102 endif |
225
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
103 compile = "Problem"; |
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
104 if (numel (varargin) > 1) |
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
105 compile = validatestring (varargin{2}, opts, "ufl", "content"); |
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
106 endif |
255
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
107 # This line is executed by all the nodes since it internally enforces |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
108 # that just one is involved |
225
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
109 eval (["import_ufl_", compile, "(\"", filename(1:end-4), "\");"]); |
255
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
110 if (is_master_node ()) |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
111 delete (filename); |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
112 endif |
225
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
113 filename = "default.ufl"; |
222
6671cb83a2dd
Make use of is_valid_file_id for fid checks
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
220
diff
changeset
|
114 elseif (! is_valid_file_id (fid)) |
255
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
115 if (is_master_node ()) |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
116 error ("ufl: no open file"); |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
117 endif |
220 | 118 else |
255
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
119 if (is_master_node ()) |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
120 if (is_valid_file_id (fid)) |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
121 output = strtrim (sprintf ("%s ", varargin{:})); |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
122 fputs (fid, output); |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
123 fprintf (fid, "\n"); |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
124 else |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
125 error ("ufl: no open file"); |
072aee55bb75
Import UFL in a parallel-safe fashion
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
232
diff
changeset
|
126 endif |
222
6671cb83a2dd
Make use of is_valid_file_id for fid checks
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
220
diff
changeset
|
127 endif |
220 | 128 endif |
129 | |
130 endfunction |