Mercurial > fem-fenics-eugenio
annotate inst/ufl.m @ 225:72af6354bf02
Compile right away after ufl end
| author | Eugenio Gianniti <eugenio.gianniti@mail.polimi.it> |
|---|---|
| date | Thu, 22 May 2014 12:34:47 +0200 |
| parents | a108a0bfe201 |
| children | 043120f071de |
| 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 | |
| 16 ## -*- texinfo -*- | |
| 17 ## @deftypefn {Function File} ufl [command] line | |
| 18 ## | |
| 19 ## Writes the given line to a ufl file. Optional commands are: | |
| 20 ## @itemize @bullet | |
| 21 ## @item start: opens an ufl file. It is possible to specify a name for it, | |
| 22 ## 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
|
23 ## @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
|
24 ## 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
|
25 ## 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
|
26 ## and Problem. |
| 220 | 27 ## @end itemize |
| 28 ## | |
|
225
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
29 ## @example |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
30 ## @code{ufl start fs.ufl} |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
31 ## @code{ufl element = FunctionSpace ("Lagrange", triangle, 1)} |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
32 ## @code{ufl end FunctionSpace} |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
33 ## This way you can write a ufl file and compile it obtaining a FunctionSpace. |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
34 ## @end example |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
35 ## |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
36 ## @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
|
37 ## import_ufl_FunctionSpace, import_ufl_LinearForm} |
| 220 | 38 ## @end deftypefn |
| 39 | |
| 40 function ufl (varargin) | |
| 41 | |
| 42 persistent fid = -1; | |
|
225
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
43 persistent filename = "default.ufl"; |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
44 persistent opts = {"BilinearForm", |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
45 "Functional", |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
46 "FunctionSpace", |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
47 "LinearForm", |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
48 "Problem"}; |
| 220 | 49 |
| 50 if (numel (varargin) < 1) | |
|
222
6671cb83a2dd
Make use of is_valid_file_id for fid checks
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
220
diff
changeset
|
51 if (! is_valid_file_id (fid)) |
| 220 | 52 print_usage (); |
| 53 endif | |
| 54 elseif (! all (cellfun ("ischar", varargin))) | |
| 55 error ("ufl: all arguments should be strings"); | |
| 56 elseif (strcmpi (varargin{1}, "start")) | |
|
222
6671cb83a2dd
Make use of is_valid_file_id for fid checks
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
220
diff
changeset
|
57 if (is_valid_file_id (fid)) |
| 220 | 58 error ("ufl: a file is already open"); |
| 59 else | |
| 60 if (numel (varargin) > 1) | |
| 61 filename = varargin{2}; | |
| 62 if (isempty (regexpi (filename, ".ufl$", "match"))) | |
| 63 filename = [filename, ".ufl"]; | |
|
224
a108a0bfe201
Ensure that ufl opens a file with .ufl extension
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
223
diff
changeset
|
64 else |
|
a108a0bfe201
Ensure that ufl opens a file with .ufl extension
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
223
diff
changeset
|
65 filename = [filename(1:end-4), ".ufl"]; |
| 220 | 66 endif |
| 67 endif | |
|
223
8b159045242f
ufl errors out if filename exists when start command is issued
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
222
diff
changeset
|
68 [~, err, ~] = stat (filename); |
|
8b159045242f
ufl errors out if filename exists when start command is issued
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
222
diff
changeset
|
69 if (err == 0) |
|
8b159045242f
ufl errors out if filename exists when start command is issued
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
222
diff
changeset
|
70 error (["ufl: a file named ", filename, " already exists"]); |
|
8b159045242f
ufl errors out if filename exists when start command is issued
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
222
diff
changeset
|
71 else |
|
8b159045242f
ufl errors out if filename exists when start command is issued
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
222
diff
changeset
|
72 fid = fopen (filename, "w"); |
|
8b159045242f
ufl errors out if filename exists when start command is issued
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
222
diff
changeset
|
73 endif |
| 220 | 74 endif |
|
222
6671cb83a2dd
Make use of is_valid_file_id for fid checks
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
220
diff
changeset
|
75 if (! is_valid_file_id (fid)) |
| 220 | 76 error (["ufl: could not open file ", filename]); |
| 77 endif | |
| 78 elseif (strcmpi (varargin{1}, "end")) | |
|
222
6671cb83a2dd
Make use of is_valid_file_id for fid checks
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
220
diff
changeset
|
79 if (! is_valid_file_id (fid)) |
| 220 | 80 error ("ufl: no open file"); |
| 81 else | |
| 82 fclose (fid); | |
| 83 fid = -1; | |
| 84 endif | |
|
225
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
85 compile = "Problem"; |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
86 if (numel (varargin) > 1) |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
87 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
|
88 endif |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
89 eval (["import_ufl_", compile, "(\"", filename(1:end-4), "\");"]); |
|
72af6354bf02
Compile right away after ufl end
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
224
diff
changeset
|
90 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
|
91 elseif (! is_valid_file_id (fid)) |
| 220 | 92 error ("ufl: no open file"); |
| 93 else | |
|
222
6671cb83a2dd
Make use of is_valid_file_id for fid checks
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
220
diff
changeset
|
94 if (is_valid_file_id (fid)) |
|
6671cb83a2dd
Make use of is_valid_file_id for fid checks
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
220
diff
changeset
|
95 fprintf (fid, "%s ", varargin{:}); |
|
6671cb83a2dd
Make use of is_valid_file_id for fid checks
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
220
diff
changeset
|
96 fprintf (fid, "\n"); |
|
6671cb83a2dd
Make use of is_valid_file_id for fid checks
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
220
diff
changeset
|
97 else |
|
6671cb83a2dd
Make use of is_valid_file_id for fid checks
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
220
diff
changeset
|
98 error ("ufl: no open file"); |
|
6671cb83a2dd
Make use of is_valid_file_id for fid checks
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
220
diff
changeset
|
99 endif |
| 220 | 100 endif |
| 101 | |
| 102 endfunction |