Mercurial > fem-fenics-eugenio

annotate inst/private/generate_mf_save.m @ 261:f22588ae37af

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Improve template meshfunction implementation * inst/import_meshfunction_type.m: provide to the auxiliary functions a type name and a valid identifier * inst/private/generate_mf_*.m: add space before closing angle bracket to avoid parse errors
author Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
date Thu, 07 Aug 2014 11:13:54 +0200
parents 1e2a9be8083a
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
1 ## Copyright (C) 2014 Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
2 ##
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
3 ## This program is free software; you can redistribute it and/or modify it under
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
4 ## the terms of the GNU General Public License as published by the Free Software
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
5 ## Foundation; either version 3 of the License, or (at your option) any later
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
6 ## version.
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
7 ##
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
8 ## This program is distributed in the hope that it will be useful, but WITHOUT
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
9 ## ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
10 ## FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
11 ## details.
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
12 ##
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
13 ## You should have received a copy of the GNU General Public License along with
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
14 ## this program; if not, see <http://www.gnu.org/licenses/>.
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
15
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
16
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
17 ## -*- texinfo -*-
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
18 ## function for internal usage only
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
19 ## @end deftypefn
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
20
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
21 function output = generate_mf_save (typename, label)
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
22
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
23 STRING ="\n\
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
24 #include ""meshfunction_@@LABEL@@.h""\n\
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
25 \n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
26 DEFUN_DLD (save, args, nargout, ""save (MESHFUNCTION, FILENAME)"")\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
27 {\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
28 int nargin = args.length ();\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
29 octave_value retval;\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
30 \n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
31 if (nargin < 2 || nargin > 2 || nargout > 1)\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
32 print_usage ();\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
33 else\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
34 {\n\
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
35 if (! meshfunction_@@LABEL@@_type_loaded)\n\
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
36 {\n\
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
37 meshfunction_@@LABEL@@::register_type ();\n\
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
38 meshfunction_@@LABEL@@_type_loaded = true;\n\
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
39 mlock ();\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
40 }\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
41 \n\
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
42 if (args(0).type_id () == meshfunction_@@LABEL@@::static_type_id ()\n\
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
43 && args(1).is_string ())\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
44 {\n\
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
45 meshfunction_@@LABEL@@ const & mf_arg =\n\
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
46 static_cast<meshfunction_@@LABEL@@ const &>\n\
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
47 (args(0).get_rep ());\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
48 std::string str = args(1).string_value ();\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
49 \n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
50 if (!error_state)\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
51 {\n\
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
52 dolfin::MeshFunction <@@TYPENAME@@ > const &\n\
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
53 mf = mf_arg.get_mf ();\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
54 str += "".xdmf"";\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
55 try\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
56 {\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
57 dolfin::File file (str);\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
58 file << mf;\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
59 }\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
60 catch (std::runtime_error &)\n\
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
61 { error (""error saving meshfunction@@LABEL@@""); }\n\
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
62 retval = 0;\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
63 }\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
64 }\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
65 else\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
66 error (""invalid input arguments"");\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
67 }\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
68 \n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
69 return retval;\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
70 }";
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
71
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
72 STRING = strrep (STRING, "@@TYPENAME@@", typename);
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
73 STRING = strrep (STRING, "@@LABEL@@", label);
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
74
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
75 fid = fopen (["save_", label, ".cc"], "w");
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
76 if (is_valid_file_id (fid))
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
77 fputs (fid, STRING);
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
78 output = fclose (fid);
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
79 else
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
80 error ("cannot open file");
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
81 output = 1;
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
82 endif
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
83
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
84 endfunction