Mercurial > fem-fenics-eugenio

annotate src/interpolate.cc @ 268:61830a4f9ab9

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Improve formatting
author Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
date Thu, 14 Aug 2014 12:26:55 +0200
parents 5e9b5bbdc56b
children
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"
242
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
21 #include <stdexcept>
253
5e9b5bbdc56b Support both DOLFIN 1.3.0 and 1.4.0
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 242
diff changeset
22 #include "dolfin_compat.h"
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
23
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
24 DEFUN_DLD (interpolate, args, nargout, "-*- texinfo -*-\n\
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
25 @deftypefn {Function File} @var{interp} = \
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
26 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
27 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
28 @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
29 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
30 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
31 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
32 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
33 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
34 @seealso{Function, Expression, Constant, FunctionSpace}\n\
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
35 @end deftypefn")
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
36 {
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
37
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
38 int nargin = args.length ();
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
39 octave_value retval;
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
40
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
41 if (nargin < 2 || nargin > 3 || nargout > 1)
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
42 { print_usage (); }
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
43 else
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 if (! function_type_loaded)
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
46 {
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
47 function::register_type ();
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
48 function_type_loaded = true;
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
49 mlock ();
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
50 }
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
51
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
52 if (! coefficient_type_loaded)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
53 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
54 coefficient::register_type ();
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
55 coefficient_type_loaded = true;
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
56 mlock ();
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
57 }
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 if (! functionspace_type_loaded)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
60 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
61 functionspace::register_type ();
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
62 functionspace_type_loaded = true;
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
63 mlock ();
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
64 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
65
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
66 octave_idx_type offset = 0;
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
67 std::string name;
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
68 if (args(0).is_string ())
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
69 {
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
70 name = args(0).string_value ();
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
71 offset = 1;
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
72 }
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
73
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
74 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
75 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
76 const functionspace & u1 =
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
77 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
78
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
79 if (! error_state)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
80 {
253
5e9b5bbdc56b Support both DOLFIN 1.3.0 and 1.4.0
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 242
diff changeset
81 SHARED_PTR <dolfin::Function>
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
82 output (new dolfin::Function (u1.get_pfsp ()));
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
83
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
84 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
85 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
86 const function & u0 =
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
87 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
88
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
89 if (! error_state)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
90 {
242
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
91 try { output->interpolate (u0.get_fun ()); }
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
92 catch (std::runtime_error &)
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
93 {
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
94 error ("unable to interpolate on this function space");
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
95 }
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
96 if (! error_state)
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
97 {
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
98 if (name.empty ())
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
99 { name = u0.get_str (); }
242
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
100 retval = new function (name, output);
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
101 }
219
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 }
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
104 else if (args(0 + offset).type_id () ==
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
105 coefficient::static_type_id ())
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
106 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
107 const coefficient & u0 =
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
108 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
109
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
110 if (! error_state)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
111 {
242
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
112 try { output->interpolate (* u0.get_expr ()); }
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
113 catch (std::runtime_error &)
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
114 {
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
115 error ("unable to interpolate on this function space");
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
116 }
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
117 if (! error_state)
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
118 {
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
119 if (name.empty ())
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
120 { name = u0.get_str (); }
242
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
121 retval = new function (name, output);
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
122 }
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
123 }
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 else
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
126 { error ("interpolate: invalid arguments"); }
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
127 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
128 }
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
129 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
130 {
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
131 const function & u0 =
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
132 static_cast<const function &> (args(1 + offset).get_rep ());
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
133
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
134 if (! error_state)
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
135 {
253
5e9b5bbdc56b Support both DOLFIN 1.3.0 and 1.4.0
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 242
diff changeset
136 SHARED_PTR <dolfin::Function>
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
137 output (new dolfin::Function (u0.get_fun ()));
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
138
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
139 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
140 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
141 const function & u1 =
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
142 static_cast<const function &> (args(0 + offset).get_rep ());
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
143
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
144 if (! error_state)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
145 {
242
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
146 try { output->interpolate (u1.get_fun ()); }
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
147 catch (std::runtime_error &)
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
148 {
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
149 error ("unable to interpolate on this function space");
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
150 }
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
151 if (! error_state)
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
152 {
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
153 if (name.empty ())
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
154 { name = u1.get_str (); }
242
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
155 retval = new function (name, output);
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
156 }
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
157 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
158 }
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
159 else if (args(0 + offset).type_id () ==
221
470586565dc7 Changed order in interpolate signature
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 219
diff changeset
160 coefficient::static_type_id ())
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
161 {
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
162 const coefficient & u1 =
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
163 static_cast<const coefficient &> (args(0 + offset).get_rep ());
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
164
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
165 if (! error_state)
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
166 {
242
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
167 try { output->interpolate (* u1.get_expr ()); }
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
168 catch (std::runtime_error &)
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
169 {
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
170 error ("unable to interpolate on this function space");
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
171 }
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
172 if (! error_state)
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
173 {
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
174 if (name.empty ())
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
175 { name = u1.get_str (); }
242
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
176 retval = new function (name, output);
958a0e0e8102 interpolate.cc: catch exception if interpolation fails
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 221
diff changeset
177 }
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
178 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
179 }
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
180 else
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
181 { error ("interpolate: invalid arguments"); }
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
182 }
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
183 }
219
a13b7d744b86 Added Python-like interpolate
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 218
diff changeset
184 else
268
61830a4f9ab9 Improve formatting
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents: 253
diff changeset
185 { error ("interpolate: invalid arguments"); }
218
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
186 }
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
187 return retval;
8a3361bfa434 interpolate function added
Eugenio Gianniti <eugenio.gianniti@mail.polimi.it>
parents:
diff changeset
188 }