[e51dba]: / src / function.h  Maximize  Restore  History

Download this file

135 lines (105 with data), 3.3 kB

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