Mercurial > fem-fenics-eugenio
annotate src/mark.cc @ 268:61830a4f9ab9
Improve formatting
author | Eugenio Gianniti <eugenio.gianniti@mail.polimi.it> |
---|---|
date | Thu, 14 Aug 2014 12:26:55 +0200 |
parents | 53039ac90368 |
children |
rev | line source |
---|---|
267
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
1 /* |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
2 Copyright (C) 2014 Eugenio Gianniti <eugenio.gianniti@mail.polimi.it> |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
3 |
53039ac90368
Mark meshfunction using subdomain
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 |
53039ac90368
Mark meshfunction using subdomain
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 |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
6 Foundation; either version 3 of the License, or (at your option) any later |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
7 version. |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
8 |
53039ac90368
Mark meshfunction using subdomain
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 |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
10 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
11 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
12 details. |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
13 |
53039ac90368
Mark meshfunction using subdomain
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 |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
15 this program; if not, see <http://www.gnu.org/licenses/>. |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
16 */ |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
17 |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
18 #include <dolfin.h> |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
19 #include "meshfunction.h" |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
20 #include "subdomain.h" |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
21 #include "dolfin_compat.h" |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
22 |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
23 DEFUN_DLD (mark, args, nargout, |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
24 "-*- texinfo -*-\n\ |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
25 @deftypefn {Function File} {[@var{markers}]} = \ |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
26 mark (@var{subdomain}, @var{meshfunction}, @var{number})\n\ |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
27 Mark @var{meshfunction} with @var{number} on the entities contained in \ |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
28 @var{subdomain}. The output, @var{markers}, is a copy of @var{meshfunction} \ |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
29 where entities inside @var{subdomain} are marked.\n\ |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
30 @seealso{MeshFunction, SubDomain}\n\ |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
31 @end deftypefn") |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
32 { |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
33 int nargin = args.length (); |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
34 octave_value retval; |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
35 |
268
61830a4f9ab9
Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
267
diff
changeset
|
36 if (nargin < 3 || nargin > 3 || nargout > 1) |
61830a4f9ab9
Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
267
diff
changeset
|
37 { print_usage (); } |
267
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
38 else |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
39 { |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
40 if (! meshfunction_type_loaded) |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
41 { |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
42 meshfunction::register_type (); |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
43 meshfunction_type_loaded = true; |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
44 mlock (); |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
45 } |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
46 if (! subdomain_type_loaded) |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
47 { |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
48 subdomain::register_type (); |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
49 subdomain_type_loaded = true; |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
50 mlock (); |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
51 } |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
52 |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
53 if (args(0).type_id () == subdomain::static_type_id () && |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
54 args(1).type_id () == meshfunction::static_type_id () && |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
55 args(2).is_real_scalar ()) |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
56 { |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
57 std::size_t number = args(2).ulong_value (); |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
58 subdomain const & sd = static_cast <subdomain const &> |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
59 (args(0).get_rep ()); |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
60 meshfunction const & mf = static_cast <meshfunction const &> |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
61 (args(1).get_rep ()); |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
62 |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
63 if (! error_state) |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
64 { |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
65 dolfin::MeshFunction <std::size_t> out = mf.get_mf (); |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
66 SHARED_PTR <subdomain_rep const> const & sdrep = sd.get_psd (); |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
67 |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
68 sdrep->mark (out, number); |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
69 retval = new meshfunction (mf.get_str (), out); |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
70 } |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
71 } |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
72 else |
268
61830a4f9ab9
Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
267
diff
changeset
|
73 { error ("invalid input arguments"); } |
267
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
74 } |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
75 |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
76 return retval; |
53039ac90368
Mark meshfunction using subdomain
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff
changeset
|
77 } |