# function.pl options: \$maxargs=${maxargs}
#
-# GiNaC Copyright (C) 1999-2005 Johannes Gutenberg University Mainz, Germany
+# GiNaC Copyright (C) 1999-2006 Johannes Gutenberg University Mainz, Germany
#
# 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
'typedef ex (* derivative_funcp_${N})(${SEQ1}, unsigned);'."\n",
'const ex &','','');
+$typedef_power_funcp=generate(
+'typedef ex (* power_funcp_${N})(${SEQ1}, const ex &);'."\n",
+'const ex &','','');
+
$typedef_series_funcp=generate(
'typedef ex (* series_funcp_${N})(${SEQ1}, const relational &, int, unsigned);'."\n",
'const ex &','','');
$derivative_func_interface=generate(' function_options & derivative_func(derivative_funcp_${N} d);'."\n",'','','');
+$power_func_interface=generate(' function_options & power_func(power_funcp_${N} d);'."\n",'','','');
+
$series_func_interface=generate(' function_options & series_func(series_funcp_${N} s);'."\n",'','','');
$print_func_interface=generate(
function::function(unsigned ser, ${SEQ1})
: exprseq(${SEQ2}), serial(ser)
{
- tinfo_key = TINFO_function;
+ tinfo_key = &function::tinfo_static;
}
END_OF_CONSTRUCTORS_IMPLEMENTATION
return ((derivative_funcp_${N})(opt.derivative_f))(${SEQ1},diff_param);
END_OF_DIFF_SWITCH_STATEMENT
+$power_switch_statement=generate(
+ <<'END_OF_POWER_SWITCH_STATEMENT','seq[${N}-1]','','');
+ case ${N}:
+ return ((power_funcp_${N})(opt.power_f))(${SEQ1},power_param);
+END_OF_POWER_SWITCH_STATEMENT
+
$series_switch_statement=generate(
<<'END_OF_SERIES_SWITCH_STATEMENT','seq[${N}-1]','','');
case ${N}:
}
END_OF_DERIVATIVE_FUNC_IMPLEMENTATION
+$power_func_implementation=generate(
+ <<'END_OF_POWER_FUNC_IMPLEMENTATION','','','');
+function_options & function_options::power_func(power_funcp_${N} d)
+{
+ test_and_set_nparams(${N});
+ power_f = power_funcp(d);
+ return *this;
+}
+END_OF_POWER_FUNC_IMPLEMENTATION
+
$series_func_implementation=generate(
<<'END_OF_SERIES_FUNC_IMPLEMENTATION','','','');
function_options & function_options::series_func(series_funcp_${N} s)
typedef ex (* evalf_funcp)();
typedef ex (* conjugate_funcp)();
typedef ex (* derivative_funcp)();
+typedef ex (* power_funcp)();
typedef ex (* series_funcp)();
typedef void (* print_funcp)();
$typedef_evalf_funcp
$typedef_conjugate_funcp
$typedef_derivative_funcp
+$typedef_power_funcp
$typedef_series_funcp
$typedef_print_funcp
// end of generated lines
typedef ex (* evalf_funcp_exvector)(const exvector &);
typedef ex (* conjugate_funcp_exvector)(const exvector &);
typedef ex (* derivative_funcp_exvector)(const exvector &, unsigned);
+typedef ex (* power_funcp_exvector)(const exvector &, const ex &);
typedef ex (* series_funcp_exvector)(const exvector &, const relational &, int, unsigned);
typedef void (* print_funcp_exvector)(const exvector &, const print_context &);
$evalf_func_interface
$conjugate_func_interface
$derivative_func_interface
+$power_func_interface
$series_func_interface
$print_func_interface
// end of generated lines
function_options & evalf_func(evalf_funcp_exvector ef);
function_options & conjugate_func(conjugate_funcp_exvector d);
function_options & derivative_func(derivative_funcp_exvector d);
+ function_options & power_func(power_funcp_exvector d);
function_options & series_func(series_funcp_exvector s);
template <class Ctx> function_options & print_func(print_funcp_exvector p)
return *this;
}
- function_options & set_return_type(unsigned rt, unsigned rtt=0);
+ function_options & set_return_type(unsigned rt, tinfo_t rtt=NULL);
function_options & do_not_evalf_params();
function_options & remember(unsigned size, unsigned assoc_size=0,
unsigned strategy=remember_strategies::delete_never);
protected:
bool has_derivative() const { return derivative_f != NULL; }
+ bool has_power() const { return power_f != NULL; }
void test_and_set_nparams(unsigned n);
void set_print_func(unsigned id, print_funcp f);
evalf_funcp evalf_f;
conjugate_funcp conjugate_f;
derivative_funcp derivative_f;
+ power_funcp power_f;
series_funcp series_f;
std::vector<print_funcp> print_dispatch_table;
bool use_return_type;
unsigned return_type;
- unsigned return_type_tinfo;
+ tinfo_t return_type_tinfo;
bool use_remember;
unsigned remember_size;
bool evalf_use_exvector_args;
bool conjugate_use_exvector_args;
bool derivative_use_exvector_args;
+ bool power_use_exvector_args;
bool series_use_exvector_args;
bool print_use_exvector_args;
bool is_equal_same_type(const basic & other) const;
bool match_same_type(const basic & other) const;
unsigned return_type() const;
- unsigned return_type_tinfo() const;
+ tinfo_t return_type_tinfo() const;
// new virtual functions which can be overridden by derived classes
// none
bool lookup_remember_table(ex & result) const;
void store_remember_table(ex const & result) const;
public:
+ ex power(const ex & exp) const;
static unsigned register_new(function_options const & opt);
static unsigned current_serial;
static unsigned find_function(const std::string &name, unsigned nparams);
// utility functions/macros
-/** Specialization of is_exactly_a<function>(obj) for objects of type function. */
-template<> inline bool is_exactly_a<function>(const basic & obj)
-{
- return obj.tinfo()==TINFO_function;
-}
-
template <typename T>
inline bool is_the_function(const ex & x)
{
#include "lst.h"
#include "symmetry.h"
#include "print.h"
+#include "power.h"
#include "archive.h"
#include "inifcns.h"
#include "tostring.h"
{
set_name("unnamed_function", "\\\\mbox{unnamed}");
nparams = 0;
- eval_f = evalf_f = conjugate_f = derivative_f = series_f = 0;
+ eval_f = evalf_f = conjugate_f = derivative_f = power_f = series_f = 0;
evalf_params_first = true;
use_return_type = false;
eval_use_exvector_args = false;
evalf_use_exvector_args = false;
conjugate_use_exvector_args = false;
derivative_use_exvector_args = false;
+ power_use_exvector_args = false;
series_use_exvector_args = false;
print_use_exvector_args = false;
use_remember = false;
$evalf_func_implementation
$conjugate_func_implementation
$derivative_func_implementation
+$power_func_implementation
$series_func_implementation
// end of generated lines
derivative_f = derivative_funcp(d);
return *this;
}
+function_options& function_options::power_func(power_funcp_exvector d)
+{
+ power_use_exvector_args = true;
+ power_f = power_funcp(d);
+ return *this;
+}
function_options& function_options::series_func(series_funcp_exvector s)
{
series_use_exvector_args = true;
return *this;
}
-function_options & function_options::set_return_type(unsigned rt, unsigned rtt)
+function_options & function_options::set_return_type(unsigned rt, tinfo_t rtt)
{
use_return_type = true;
return_type = rt;
function::function() : serial(0)
{
- tinfo_key = TINFO_function;
+ tinfo_key = &function::tinfo_static;
}
//////////
function::function(unsigned ser) : serial(ser)
{
- tinfo_key = TINFO_function;
+ tinfo_key = &function::tinfo_static;
}
// the following lines have been generated for max. ${maxargs} parameters
function::function(unsigned ser, const exprseq & es) : exprseq(es), serial(ser)
{
- tinfo_key = TINFO_function;
+ tinfo_key = &function::tinfo_static;
// Force re-evaluation even if the exprseq was already evaluated
// (the exprseq copy constructor copies the flags)
function::function(unsigned ser, const exvector & v, bool discardable)
: exprseq(v,discardable), serial(ser)
{
- tinfo_key = TINFO_function;
+ tinfo_key = &function::tinfo_static;
}
function::function(unsigned ser, std::auto_ptr<exvector> vp)
: exprseq(vp), serial(ser)
{
- tinfo_key = TINFO_function;
+ tinfo_key = &function::tinfo_static;
}
//////////
unsigned function::calchash() const
{
- unsigned v = golden_ratio_hash(golden_ratio_hash(tinfo()) ^ serial);
+ unsigned v = golden_ratio_hash(golden_ratio_hash((p_int)tinfo()) ^ serial);
for (size_t i=0; i<nops(); i++) {
v = rotate_left(v);
v ^= this->op(i).gethash();
}
}
-unsigned function::return_type_tinfo() const
+tinfo_t function::return_type_tinfo() const
{
GINAC_ASSERT(serial<registered_functions().size());
const function_options &opt = registered_functions()[serial];
// Default behavior is to use the return type of the first
// argument. Thus, exp() of a matrix behaves like a matrix, etc.
if (seq.empty())
- return tinfo_key;
+ return this;
else
return seq.begin()->return_type_tinfo();
}
throw(std::logic_error("function::pderivative(): no diff function defined"));
}
+ex function::power(const ex & power_param) const // power of function
+{
+ GINAC_ASSERT(serial<registered_functions().size());
+ const function_options &opt = registered_functions()[serial];
+
+ // No derivative defined? Then return abstract derivative object
+ if (opt.power_f == NULL)
+ return (new power::power(*this, power_param))->setflag(status_flags::dynallocated |
+ status_flags::evaluated);
+
+ current_serial = serial;
+ if (opt.power_use_exvector_args)
+ return ((power_funcp_exvector)(opt.power_f))(seq, power_param);
+ switch (opt.nparams) {
+ // the following lines have been generated for max. ${maxargs} parameters
+${power_switch_statement}
+ // end of generated lines
+ }
+ throw(std::logic_error("function::power(): no power function defined"));
+}
+
std::vector<function_options> & function::registered_functions()
{
static std::vector<function_options> * rf = new std::vector<function_options>;
git://www.ginac.de / ginac.git / blobdiff