Mercurial > fem-fenics-eugenio
annotate src/DirichletBC.cc @ 143:9486cbdc0a2e
Maint: update the documentation
author | gedeone-octave <marcovass89@hotmail.it> |
---|---|
date | Mon, 09 Sep 2013 22:40:46 +0200 |
parents | 3a0f952a79ec |
children | 5fe2a157f4eb |
rev | line source |
---|---|
31
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
1 /* |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
2 Copyright (C) 2013 Marco Vassallo |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
3 |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
4 This program is free software; you can redistribute it and/or modify it under |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
5 the terms of the GNU General Public License as published by the Free Software |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
6 Foundation; either version 2 of the License, or (at your option) any later |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
7 version. |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
8 |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
9 This program is distributed in the hope that it will be useful, but WITHOUT |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
10 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
11 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
12 details. |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
13 |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
14 You should have received a copy of the GNU General Public License along with |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
15 this program; if not, see <http://www.gnu.org/licenses/>. |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
16 */ |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
17 |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
18 #include "boundarycondition.h" |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
19 #include "functionspace.h" |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
20 #include "expression.h" |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
21 |
78
670a5d91c397
Rename the function visible from the user accordingly to dolfin interface.
gedeone-octave <marcovass89@hotmail.it>
parents:
75
diff
changeset
|
22 DEFUN_DLD (DirichletBC, args, , "-*- texinfo -*-\n\ |
31
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
23 @deftypefn {Function File} {[@var{bc}]} = \ |
135 | 24 DirichletBC (@var{Functional Space}, @var{Boundary Label}, @var{Function handle}) \n\ |
31
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
25 The input parameters are\n\ |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
26 @itemize @bullet \n\ |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
27 @item @var{Functional Space} is a fem-fenics functional space where\ |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
28 we want to apply the BC\n\ |
117
70322135b25f
Maint: Add copyright notice
gedeone-octave <marcovass89@hotmail.it>
parents:
78
diff
changeset
|
29 @item @var{Function handle} is a function handle which contains the expression\ |
70322135b25f
Maint: Add copyright notice
gedeone-octave <marcovass89@hotmail.it>
parents:
78
diff
changeset
|
30 that we want to apply as BC. If we have a Vector field, we can just use a\n\ |
125
a80ac536c78a
Fix bug in the texinfo code
gedeone-octave <marcovass89@hotmail.it>
parents:
117
diff
changeset
|
31 vector of function handle: @var{Function handle} = ['@'(x, y) f1, '@'(x, y) f2, ...]\n\ |
117
70322135b25f
Maint: Add copyright notice
gedeone-octave <marcovass89@hotmail.it>
parents:
78
diff
changeset
|
32 @item @var{Boundary} is an Array which contains the labels of the \ |
31
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
33 sides where the BC is applied\n\ |
50
fcfecdd3a9b5
Maint: improve the documentation
gedeone-octave <marco.vassallo@outlook.com>
parents:
44
diff
changeset
|
34 @end itemize\n\ |
31
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
35 The output @var{BC} is an object which contains the boundary conditions\n\ |
143
9486cbdc0a2e
Maint: update the documentation
gedeone-octave <marcovass89@hotmail.it>
parents:
137
diff
changeset
|
36 @seealso{Mesh, FunctionSpace}\n\ |
31
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
37 @end deftypefn") |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
38 { |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
39 int nargin = args.length (); |
75
d6df5cc8ef53
Improve the expression class to deal also with vector fields.
gedeone-octave <marcovass89@hotmail.it>
parents:
50
diff
changeset
|
40 octave_value retval; |
31
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
41 |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
42 if (nargin < 3 || nargin > 3) |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
43 print_usage (); |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
44 else |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
45 { |
44
fca8c3d75036
register_type is called before every type is used.
gedeone-octave <marco.vassallo@outlook.com>
parents:
31
diff
changeset
|
46 |
fca8c3d75036
register_type is called before every type is used.
gedeone-octave <marco.vassallo@outlook.com>
parents:
31
diff
changeset
|
47 if (! functionspace_type_loaded) |
fca8c3d75036
register_type is called before every type is used.
gedeone-octave <marco.vassallo@outlook.com>
parents:
31
diff
changeset
|
48 { |
fca8c3d75036
register_type is called before every type is used.
gedeone-octave <marco.vassallo@outlook.com>
parents:
31
diff
changeset
|
49 functionspace::register_type (); |
fca8c3d75036
register_type is called before every type is used.
gedeone-octave <marco.vassallo@outlook.com>
parents:
31
diff
changeset
|
50 functionspace_type_loaded = true; |
fca8c3d75036
register_type is called before every type is used.
gedeone-octave <marco.vassallo@outlook.com>
parents:
31
diff
changeset
|
51 mlock (); |
fca8c3d75036
register_type is called before every type is used.
gedeone-octave <marco.vassallo@outlook.com>
parents:
31
diff
changeset
|
52 } |
fca8c3d75036
register_type is called before every type is used.
gedeone-octave <marco.vassallo@outlook.com>
parents:
31
diff
changeset
|
53 |
117
70322135b25f
Maint: Add copyright notice
gedeone-octave <marcovass89@hotmail.it>
parents:
78
diff
changeset
|
54 if (! boundarycondition_type_loaded) |
70322135b25f
Maint: Add copyright notice
gedeone-octave <marcovass89@hotmail.it>
parents:
78
diff
changeset
|
55 { |
70322135b25f
Maint: Add copyright notice
gedeone-octave <marcovass89@hotmail.it>
parents:
78
diff
changeset
|
56 boundarycondition::register_type (); |
70322135b25f
Maint: Add copyright notice
gedeone-octave <marcovass89@hotmail.it>
parents:
78
diff
changeset
|
57 boundarycondition_type_loaded = true; |
70322135b25f
Maint: Add copyright notice
gedeone-octave <marcovass89@hotmail.it>
parents:
78
diff
changeset
|
58 mlock (); |
70322135b25f
Maint: Add copyright notice
gedeone-octave <marcovass89@hotmail.it>
parents:
78
diff
changeset
|
59 } |
70322135b25f
Maint: Add copyright notice
gedeone-octave <marcovass89@hotmail.it>
parents:
78
diff
changeset
|
60 |
31
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
61 if (args(0).type_id () == functionspace::static_type_id ()) |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
62 { |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
63 const functionspace & fspo = static_cast<const functionspace&> (args(0).get_rep ()); |
137
3a0f952a79ec
Avois useless output.
gedeone-octave <marcovass89@hotmail.it>
parents:
135
diff
changeset
|
64 octave_fcn_handle * fh = args(1).fcn_handle_value (); |
3a0f952a79ec
Avois useless output.
gedeone-octave <marcovass89@hotmail.it>
parents:
135
diff
changeset
|
65 Array<octave_idx_type> side = args(2).array_value (); |
3a0f952a79ec
Avois useless output.
gedeone-octave <marcovass89@hotmail.it>
parents:
135
diff
changeset
|
66 |
31
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
67 |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
68 if (!error_state) |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
69 { |
117
70322135b25f
Maint: Add copyright notice
gedeone-octave <marcovass89@hotmail.it>
parents:
78
diff
changeset
|
70 std::cout <<"Building essential boundary condition..."<< std::endl; |
31
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
71 const boost::shared_ptr <const dolfin::FunctionSpace> & V (fspo.get_pfsp ()); |
75
d6df5cc8ef53
Improve the expression class to deal also with vector fields.
gedeone-octave <marcovass89@hotmail.it>
parents:
50
diff
changeset
|
72 |
d6df5cc8ef53
Improve the expression class to deal also with vector fields.
gedeone-octave <marcovass89@hotmail.it>
parents:
50
diff
changeset
|
73 octave_value_list b (3, 1); |
d6df5cc8ef53
Improve the expression class to deal also with vector fields.
gedeone-octave <marcovass89@hotmail.it>
parents:
50
diff
changeset
|
74 octave_value_list tmp = feval (fh->function_value (), b); |
d6df5cc8ef53
Improve the expression class to deal also with vector fields.
gedeone-octave <marcovass89@hotmail.it>
parents:
50
diff
changeset
|
75 Array<double> res = tmp(0).array_value (); |
d6df5cc8ef53
Improve the expression class to deal also with vector fields.
gedeone-octave <marcovass89@hotmail.it>
parents:
50
diff
changeset
|
76 std::size_t l = res.length (); |
d6df5cc8ef53
Improve the expression class to deal also with vector fields.
gedeone-octave <marcovass89@hotmail.it>
parents:
50
diff
changeset
|
77 |
d6df5cc8ef53
Improve the expression class to deal also with vector fields.
gedeone-octave <marcovass89@hotmail.it>
parents:
50
diff
changeset
|
78 expression * pf; |
d6df5cc8ef53
Improve the expression class to deal also with vector fields.
gedeone-octave <marcovass89@hotmail.it>
parents:
50
diff
changeset
|
79 if (l > 1) |
d6df5cc8ef53
Improve the expression class to deal also with vector fields.
gedeone-octave <marcovass89@hotmail.it>
parents:
50
diff
changeset
|
80 pf = new expression (*fh, l); |
d6df5cc8ef53
Improve the expression class to deal also with vector fields.
gedeone-octave <marcovass89@hotmail.it>
parents:
50
diff
changeset
|
81 else |
d6df5cc8ef53
Improve the expression class to deal also with vector fields.
gedeone-octave <marcovass89@hotmail.it>
parents:
50
diff
changeset
|
82 pf = new expression (*fh); |
d6df5cc8ef53
Improve the expression class to deal also with vector fields.
gedeone-octave <marcovass89@hotmail.it>
parents:
50
diff
changeset
|
83 |
31
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
84 boost::shared_ptr<const expression> f (pf); |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
85 boundarycondition * pbc = new boundarycondition (); |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
86 |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
87 for (octave_idx_type i = 0; i < side.length (); ++i) |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
88 pbc->add_bc (V, f, side(i)); |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
89 retval = octave_value (pbc); |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
90 } |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
91 } |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
92 } |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
93 return retval; |
74cbe97f55af
New class for treatment of BC
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
94 } |