Mercurial > fem-fenics-eugenio
annotate src/function.h @ 234:94895618701f
Add function evaluation with coordinates inside brackets
author | Eugenio Gianniti <eugenio.gianniti@mail.polimi.it> |
---|---|
date | Sat, 21 Jun 2014 16:07:42 +0200 |
parents | 9e944b0d0fc8 |
children | 568aff85dfa9 |
rev | line source |
---|---|
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
1 /* |
151 | 2 Copyright (C) 2013 Marco Vassallo <gedeone-octave@users.sourceforge.net> |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
3 |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
4 This program is free software; you can redistribute it and/or modify it under |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
5 the terms of the GNU General Public License as published by the Free Software |
151 | 6 Foundation; either version 3 of the License, or (at your option) any later |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
7 version. |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
8 |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
9 This program is distributed in the hope that it will be useful, but WITHOUT |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
10 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
11 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
12 details. |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
13 |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
14 You should have received a copy of the GNU General Public License along with |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
15 this program; if not, see <http://www.gnu.org/licenses/>. |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
16 */ |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
17 |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
18 #ifndef _FUNCTION_OCTAVE_ |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
19 #define _FUNCTION_OCTAVE_ |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
20 |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
21 #include <memory> |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
22 #include <vector> |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
23 #include <dolfin.h> |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
24 #include <octave/oct.h> |
234
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
25 #include <octave/parse.h> |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
26 |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
27 class function : public octave_base_value |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
28 { |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
29 |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
30 public: |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
31 |
173
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
32 function () |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
33 : octave_base_value (), fun () {} |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
34 |
173
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
35 function (std::string & _str, |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
36 boost::shared_ptr <const dolfin::Function> _fun) |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
37 : octave_base_value (), str(_str), fun (_fun) {} |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
38 |
234
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
39 function (function const& _func) |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
40 : octave_base_value () |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
41 { |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
42 str = _func.get_str (); |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
43 fun = _func.get_pfun (); |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
44 } |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
45 |
173
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
46 void |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
47 print (std::ostream& os, |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
48 bool pr_as_read_syntax = false) const |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
49 { os << "Function " << str << ": " << fun->str (true) << std::endl; } |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
50 |
173
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
51 ~function(void) {} |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
52 |
173
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
53 bool |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
54 is_defined (void) const |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
55 { return true; } |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
56 |
173
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
57 const dolfin::Function & |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
58 get_fun (void) const |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
59 { return (*fun); } |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
60 |
173
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
61 const boost::shared_ptr <const dolfin::Function> & |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
62 get_pfun (void) const |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
63 { return fun; } |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
64 |
173
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
65 void |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
66 set_fun (dolfin::Function & _fun) |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
67 { |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
68 dolfin::Function * p = new dolfin::Function (_fun); |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
69 fun = boost::shared_ptr<const dolfin::Function> (p); |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
70 } |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
71 |
173
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
72 const std::string & |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
73 get_str (void) const |
9e944b0d0fc8
Some Formatting improvements (?)
gedeone-octave <marcovass89@hotmail.it>
parents:
151
diff
changeset
|
74 { return str; } |
61
93e538063278
Add a string field to the class function
gedeone-octave <marco.vassallo@outlook.com>
parents:
46
diff
changeset
|
75 |
234
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
76 octave_value |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
77 subsref (const std::string& type, |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
78 const std::list<octave_value_list>& idx) |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
79 { |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
80 octave_value retval; |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
81 |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
82 switch (type[0]) |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
83 { |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
84 case '(': |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
85 { |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
86 std::list<octave_value_list> args; |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
87 args.push_back (octave_value (new function (*this))); |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
88 args.push_back (idx.front ()); |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
89 retval = feval ("feval", args); |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
90 } |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
91 break; |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
92 |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
93 case '{': |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
94 case '.': |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
95 { |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
96 std::string nm = type_name (); |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
97 error ("%s cannot be indexed with %c", nm.c_str (), type[0]); |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
98 } |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
99 break; |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
100 |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
101 default: |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
102 panic_impossible (); |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
103 } |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
104 |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
105 return retval; |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
106 } |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
107 |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
108 octave_value_list |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
109 subsref (const std::string& type, |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
110 const std::list<octave_value_list>& idx, |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
111 int) |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
112 { |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
113 octave_value_list retval; |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
114 retval = subsref (type, idx); |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
115 return retval; |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
116 } |
94895618701f
Add function evaluation with coordinates inside brackets
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
173
diff
changeset
|
117 |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
118 private: |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
119 |
61
93e538063278
Add a string field to the class function
gedeone-octave <marco.vassallo@outlook.com>
parents:
46
diff
changeset
|
120 std::string str; |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
121 boost::shared_ptr <const dolfin::Function> fun; |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
122 |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
123 DECLARE_OCTAVE_ALLOCATOR; |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
124 DECLARE_OV_TYPEID_FUNCTIONS_AND_DATA; |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
125 |
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
126 }; |
80
16ccfaf1476a
The function DEFINE_OCTAVE_ALLOCATOR and DEFINE_OV_TYPEID_FUNCTIONS_AND_DATA
gedeone-octave <marcovass89@hotmail.it>
parents:
61
diff
changeset
|
127 static bool function_type_loaded = false; |
16ccfaf1476a
The function DEFINE_OCTAVE_ALLOCATOR and DEFINE_OV_TYPEID_FUNCTIONS_AND_DATA
gedeone-octave <marcovass89@hotmail.it>
parents:
61
diff
changeset
|
128 |
16ccfaf1476a
The function DEFINE_OCTAVE_ALLOCATOR and DEFINE_OV_TYPEID_FUNCTIONS_AND_DATA
gedeone-octave <marcovass89@hotmail.it>
parents:
61
diff
changeset
|
129 DEFINE_OCTAVE_ALLOCATOR (function); |
16ccfaf1476a
The function DEFINE_OCTAVE_ALLOCATOR and DEFINE_OV_TYPEID_FUNCTIONS_AND_DATA
gedeone-octave <marcovass89@hotmail.it>
parents:
61
diff
changeset
|
130 DEFINE_OV_TYPEID_FUNCTIONS_AND_DATA (function, "function", "function"); |
46
a64e195d0611
Wrapper class for dolfin::Function
gedeone-octave <marco.vassallo@outlook.com>
parents:
diff
changeset
|
131 #endif |