Mercurial > fem-fenics-eugenio

annotate src/interpolate.cc @ 221:470586565dc7

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Changed order in interpolate signature
author Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
date Mon, 19 May 2014 23:18:33 +0200
parents a13b7d744b86
children 958a0e0e8102
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
1 /*
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
2 Copyright (C) 2014 Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
3
8a3361bfa434 interpolate function added
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
8a3361bfa434 interpolate function added
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
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
6 Foundation; either version 3 of the License, or (at your option) any later
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
7 version.
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
8
8a3361bfa434 interpolate function added
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
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
10 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
11 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
12 details.
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
13
8a3361bfa434 interpolate function added
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
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
15 this program; if not, see <http://www.gnu.org/licenses/>.
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
16 */
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
17
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
18 #include "coefficient.h"
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
19 #include "function.h"
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
20 #include "functionspace.h"
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
21
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
22 DEFUN_DLD (interpolate, args, nargout, "-*- texinfo -*-\n\
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
23 @deftypefn {Function File} @var{interp} = \
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
24 interpolate (@var{name}, @var{u}, @var{v})\n\
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
25 Interpolate a function on a FunctionSpace.\n\
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
26 @var{u} is the Function, Expression or Constant \
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
27 to interpolate. @var{v} may be a FunctionSpace \
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
28 or a Function. In the latter case @var{u} is interpolated \
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
29 on the same FunctionSpace where @var{v} is defined.\n\
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
30 Note that @var{name} is optional: if not provided, \
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
31 it will default to the one assigned to @var{u}.\n\
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
32 @seealso{Function, Expression, Constant, FunctionSpace}\n\
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
33 @end deftypefn")
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
34 {
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
35
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
36 int nargin = args.length ();
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
37 octave_value retval;
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
38
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
39 if (nargin < 2 || nargin > 3 || nargout > 1)
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
40 print_usage ();
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
41 else
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
42 {
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
43 if (! function_type_loaded)
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
44 {
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
45 function::register_type ();
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
46 function_type_loaded = true;
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
47 mlock ();
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
48 }
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
49
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
50 if (! coefficient_type_loaded)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
51 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
52 coefficient::register_type ();
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
53 coefficient_type_loaded = true;
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
54 mlock ();
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
55 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
56
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
57 if (! functionspace_type_loaded)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
58 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
59 functionspace::register_type ();
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
60 functionspace_type_loaded = true;
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
61 mlock ();
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
62 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
63
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
64 octave_idx_type offset = 0;
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
65 std::string name;
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
66 if (args(0).is_string ())
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
67 {
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
68 name = args(0).string_value ();
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
69 offset = 1;
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
70 }
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
71
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
72 if (args(1+offset).type_id () == functionspace::static_type_id ())
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
73 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
74 const functionspace & u1 =
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
75 static_cast<const functionspace&> (args(1+offset).get_rep ());
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
76
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
77 if (! error_state)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
78 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
79 boost::shared_ptr<dolfin::Function>
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
80 output (new dolfin::Function (u1.get_pfsp ()));
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
81
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
82 if (args(0+offset).type_id () == function::static_type_id ())
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
83 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
84 const function & u0 =
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
85 static_cast<const function&> (args(0+offset).get_rep ());
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
86
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
87 if (! error_state)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
88 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
89 output->interpolate (u0.get_fun ());
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
90 if (name.empty ())
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
91 name = u0.get_str ();
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
92 retval = new function (name, output);
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
93 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
94 }
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
95 else if (args(0+offset).type_id () ==
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
96 coefficient::static_type_id ())
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
97 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
98 const coefficient & u0 =
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
99 static_cast<const coefficient&> (args(0+offset).get_rep ());
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
100
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
101 if (! error_state)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
102 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
103 output->interpolate (* u0.get_expr ());
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
104 if (name.empty ())
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
105 name = u0.get_str ();
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
106 retval = new function (name, output);
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
107 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
108 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
109 else
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
110 error ("interpolate: invalid arguments");
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
111 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
112 }
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
113 else if (args(1+offset).type_id () == function::static_type_id ())
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
114 {
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
115 const function & u0 =
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
116 static_cast<const function&> (args(1+offset).get_rep ());
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
117
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
118 if (! error_state)
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
119 {
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
120 boost::shared_ptr<dolfin::Function>
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
121 output (new dolfin::Function (u0.get_fun ()));
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
122
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
123 if (args(0+offset).type_id () == function::static_type_id ())
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
124 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
125 const function & u1 =
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
126 static_cast<const function&> (args(0+offset).get_rep ());
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
127
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
128 if (! error_state)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
129 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
130 output->interpolate (u1.get_fun ());
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
131 if (name.empty ())
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
132 name = u1.get_str ();
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
133 retval = new function (name, output);
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
134 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
135 }
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
136 else if (args(0+offset).type_id () ==
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
137 coefficient::static_type_id ())
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
138 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
139 const coefficient & u1 =
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
140 static_cast<const coefficient&> (args(0+offset).get_rep ());
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
141
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
142 if (! error_state)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
143 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
144 output->interpolate (* u1.get_expr ());
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
145 if (name.empty ())
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
146 name = u1.get_str ();
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
147 retval = new function (name, output);
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
148 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
149 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
150 else
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
151 error ("interpolate: invalid arguments");
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
152 }
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
153 }
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
154 else
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
155 error ("interpolate: invalid arguments");
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
156 }
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
157 return retval;
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
158 }