Mercurial > fem-fenics-eugenio
annotate src/save_mf.cc @ 268:61830a4f9ab9
Improve formatting
author | Eugenio Gianniti <eugenio.gianniti@mail.polimi.it> |
---|---|
date | Thu, 14 Aug 2014 12:26:55 +0200 |
parents | 0b9dc516ba36 |
children |
rev | line source |
---|---|
259
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
1 /* |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
2 Copyright (C) 2014 Eugenio Gianniti <eugenio.gianniti@mail.polimi.it> |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
3 |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
4 This program is free software; you can redistribute it and/or modify it under |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
5 the terms of the GNU General Public License as published by the Free Software |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
6 Foundation; either version 3 of the License, or (at your option) any later |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
7 version. |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
8 |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
9 This program is distributed in the hope that it will be useful, but WITHOUT |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
10 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
11 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
12 details. |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
13 |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
14 You should have received a copy of the GNU General Public License along with |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
15 this program; if not, see <http://www.gnu.org/licenses/>. |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
16 */ |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
17 |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
18 #include "meshfunction.h" |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
19 |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
20 DEFUN_DLD (save, args, nargout, "-*- texinfo -*-\n\ |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
21 @deftypefn {Function File} \ |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
22 save (@var{MeshFunction}, @var{Name})\n\ |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
23 Save a meshfunction in XDMF format.\n\ |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
24 The input parameters are\n\ |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
25 @itemize @bullet \n\ |
265
0b9dc516ba36
Correct doc-string for @meshfunction/save
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
259
diff
changeset
|
26 @item @var{MeshFunction} is the meshfunction you want to save\n\ |
259
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
27 @item @var{Name} is a string for the output name\n\ |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
28 @end itemize\n\ |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
29 The output is a .xdmf file.\n\ |
265
0b9dc516ba36
Correct doc-string for @meshfunction/save
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
259
diff
changeset
|
30 @seealso{MeshFunction}\n\ |
259
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
31 @end deftypefn") |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
32 { |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
33 int nargin = args.length (); |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
34 octave_value retval; |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
35 |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
36 if (nargin < 2 || nargin > 2 || nargout > 1) |
268
61830a4f9ab9
Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
265
diff
changeset
|
37 { print_usage (); } |
259
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
38 else |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
39 { |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
40 if (! meshfunction_type_loaded) |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
41 { |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
42 meshfunction::register_type (); |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
43 meshfunction_type_loaded = true; |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
44 mlock (); |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
45 } |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
46 |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
47 if (args(0).type_id () == meshfunction::static_type_id () && |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
48 args(1).is_string ()) |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
49 { |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
50 meshfunction const & mf_arg = |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
51 static_cast<meshfunction const &> (args(0).get_rep ()); |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
52 std::string str = args(1).string_value (); |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
53 |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
54 if (!error_state) |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
55 { |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
56 dolfin::MeshFunction <std::size_t> const & |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
57 mf = mf_arg.get_mf (); |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
58 str += ".xdmf"; |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
59 try |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
60 { |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
61 dolfin::File file (str); |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
62 file << mf; |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
63 } |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
64 catch (std::runtime_error &) |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
65 { error ("error saving meshfunction"); } |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
66 retval = 0; |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
67 } |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
68 } |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
69 else |
268
61830a4f9ab9
Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
265
diff
changeset
|
70 { error ("invalid input arguments"); } |
259
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
71 } |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
72 |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
73 return retval; |
598c5e9e0a9e
Add utilities for meshfunction
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
74 } |