Mercurial > fem-fenics-eugenio

annotate inst/private/generate_mf_constructor.m @ 268:61830a4f9ab9

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Improve formatting
author Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
date Thu, 14 Aug 2014 12:26:55 +0200
parents 41b76530fe5e
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_constructor (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\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
24 #include <dolfin.h>\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
25 #include <fem-fenics/mesh.h>\n\
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
26 #include ""meshfunction_@@LABEL@@.h""\n\
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
27 #include <fem-fenics/dolfin_compat.h>\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
28 \n\
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 264
diff changeset
29 DEFUN_DLD (MeshFunction_@@LABEL@@, args, nargout,\
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 264
diff changeset
30 ""MF = MeshFunction_@@LABEL@@ (TYPE, MESH, FILENAME)"")\n\
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
31 {\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
32 \n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
33 int nargin = args.length ();\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
34 octave_value retval;\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
35 \n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
36 if (nargin < 2 || nargin > 2 || nargout > 1)\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
37 print_usage ();\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
38 else\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
39 {\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
40 if (! mesh_type_loaded)\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
41 {\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
42 mesh::register_type ();\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
43 mesh_type_loaded = true;\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
44 mlock ();\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
45 }\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
46 \n\
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
47 if (! meshfunction_@@LABEL@@_type_loaded)\n\
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
48 {\n\
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
49 meshfunction_@@LABEL@@::register_type ();\n\
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
50 meshfunction_@@LABEL@@_type_loaded = true;\n\
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
51 mlock ();\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
52 }\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
53 \n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
54 if (args(0).type_id () == mesh::static_type_id ()\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
55 && args(1).is_string ())\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 std::string filename = args(1).string_value ();\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
58 mesh const & msh_arg =\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
59 static_cast<mesh const &> (args(0).get_rep ());\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
60 \n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
61 if (!error_state)\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
62 {\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
63 SHARED_PTR <dolfin::Mesh const> const &\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
64 pmsh = msh_arg.get_pmsh ();\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
65 \n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
66 try\n\
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
67 { retval = new meshfunction_@@LABEL@@ (pmsh, filename); }\n\
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
68 catch (std::runtime_error &)\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
69 { error (""error reading file""); }\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
70 }\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
71 }\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
72 else\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
73 error (""invalid input arguments"");\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
74 }\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
75 \n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
76 return retval;\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
77 }";
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
78
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
79 STRING = strrep (STRING, "@@LABEL@@", label);
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
80
261
f22588ae37af Improve template meshfunction implementation
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 260
diff changeset
81 fid = fopen (["MeshFunction_", label, ".cc"], "w");
260
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
82 if (is_valid_file_id (fid))
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
83 fputs (fid, STRING);
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
84 output = fclose (fid);
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
85 else
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
86 error ("cannot open file");
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
87 output = 1;
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
88 endif
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
89
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
90 endfunction