Mercurial > fem-fenics-eugenio

annotate inst/private/generate_mf_save.m @ 260:1e2a9be8083a

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Support template dolfin::MeshFunction
author Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
date Wed, 06 Aug 2014 19:52:03 +0200
parents
children f22588ae37af
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
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
21 function output = generate_mf_save (typename)
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 ""meshfunction_@@TYPENAME@@.h""\n\
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\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
35 if (! meshfunction_@@TYPENAME@@_type_loaded)\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
36 {\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
37 meshfunction_@@TYPENAME@@::register_type ();\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
38 meshfunction_@@TYPENAME@@_type_loaded = true;\n\
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\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
42 if (args(0).type_id () == meshfunction_@@TYPENAME@@::static_type_id ()\n\
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\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
45 meshfunction_@@TYPENAME@@ const & mf_arg =\n\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
46 static_cast<meshfunction_@@TYPENAME@@ const &>\n\
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\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
52 dolfin::MeshFunction <@@TYPENAME@@> const &\n\
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\
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
61 { error (""error saving meshfunction@@TYPENAME@@""); }\n\
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
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
72 STRING = strrep (STRING, "@@TYPENAME@@", typename);
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
73
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
74 fid = fopen (["save_", typename, ".cc"], "w");
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
75 if (is_valid_file_id (fid))
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
76 fputs (fid, STRING);
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
77 output = fclose (fid);
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
78 else
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
79 error ("cannot open file");
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
80 output = 1;
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
81 endif
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
82
1e2a9be8083a Support template dolfin::MeshFunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
83 endfunction