#include "power.h"
#include "archive.h"
#include "inifcns.h"
-#include "tostring.h"
#include "utils.h"
#include "hash_seed.h"
#include "remember.h"
set_name("unnamed_function", "\\\\mbox{unnamed}");
nparams = 0;
eval_f = evalf_f = real_part_f = imag_part_f = conjugate_f = expand_f
- = derivative_f = expl_derivative_f = power_f = series_f = 0;
+ = derivative_f = expl_derivative_f = power_f = series_f = nullptr;
info_f = 0;
evalf_params_first = true;
use_return_type = false;
// the following lines have been generated for max. @maxargs@ parameters
+++ for N in range(1, maxargs + 1):
function::function(unsigned ser, @seq('const ex & param%(n)d', N)@)
- : exprseq(@seq('param%(n)d', N)@), serial(ser)
+ : exprseq{@seq('param%(n)d', N)@}, serial(ser)
{
}
---
clearflag(status_flags::evaluated);
}
-function::function(unsigned ser, const exvector & v, bool discardable)
- : exprseq(v,discardable), serial(ser)
+function::function(unsigned ser, const exvector & v)
+ : exprseq(v), serial(ser)
{
}
ex function::eval(int level) const
{
+ if ((level == 1) && (flags & status_flags::evaluated)) {
+ return *this;
+ }
+
if (level>1) {
// first evaluate children, then we will end up here again
return function(serial,evalchildren(level));
}
}
- if (opt.eval_f==0) {
+ if (opt.eval_f==nullptr) {
return this->hold();
}
GINAC_ASSERT(serial<registered_functions().size());
const function_options &opt = registered_functions()[serial];
- // No derivative defined? Then return abstract derivative object
- if (opt.derivative_f == nullptr)
- return fderivative(serial, diff_param, seq);
-
- current_serial = serial;
- if (opt.derivative_use_exvector_args)
- return ((derivative_funcp_exvector)(opt.derivative_f))(seq, diff_param);
- switch (opt.nparams) {
- // the following lines have been generated for max. @maxargs@ parameters
+ if (opt.derivative_f) {
+ // Invoke the defined derivative function.
+ current_serial = serial;
+ if (opt.derivative_use_exvector_args)
+ return ((derivative_funcp_exvector)(opt.derivative_f))(seq, diff_param);
+ switch (opt.nparams) {
+ // the following lines have been generated for max. @maxargs@ parameters
+++ for N in range(1, maxargs + 1):
- case @N@:
- return ((derivative_funcp_@N@)(opt.derivative_f))(@seq('seq[%(n)d]', N, 0)@, diff_param);
+ case @N@:
+ return ((derivative_funcp_@N@)(opt.derivative_f))(@seq('seq[%(n)d]', N, 0)@, diff_param);
---
- // end of generated lines
+ // end of generated lines
+ }
}
- throw(std::logic_error("function::pderivative(): no diff function defined"));
+ // No derivative defined? Fall back to abstract derivative object.
+ return fderivative(serial, diff_param, seq);
}
ex function::expl_derivative(const symbol & s) const // explicit differentiation
GINAC_ASSERT(serial<registered_functions().size());
const function_options &opt = registered_functions()[serial];
- // No explicit derivative defined? Then this function shall not be called!
- if (opt.expl_derivative_f == nullptr)
- throw(std::logic_error("function::expl_derivative(): explicit derivation is called, but no such function defined"));
-
- current_serial = serial;
- if (opt.expl_derivative_use_exvector_args)
- return ((expl_derivative_funcp_exvector)(opt.expl_derivative_f))(seq, s);
- switch (opt.nparams) {
- // the following lines have been generated for max. @maxargs@ parameters
+ if (opt.expl_derivative_f) {
+ // Invoke the defined explicit derivative function.
+ current_serial = serial;
+ if (opt.expl_derivative_use_exvector_args)
+ return ((expl_derivative_funcp_exvector)(opt.expl_derivative_f))(seq, s);
+ switch (opt.nparams) {
+ // the following lines have been generated for max. @maxargs@ parameters
+++ for N in range(1, maxargs + 1):
- case @N@:
- return ((expl_derivative_funcp_@N@)(opt.expl_derivative_f))(@seq('seq[%(n)d]', N, 0)@, s);
+ case @N@:
+ return ((expl_derivative_funcp_@N@)(opt.expl_derivative_f))(@seq('seq[%(n)d]', N, 0)@, s);
---
- // end of generated lines
+ // end of generated lines
+ }
}
+ // There is no fallback for explicit deriviative.
+ throw(std::logic_error("function::expl_derivative(): explicit derivation is called, but no such 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];
- if (opt.power_f == nullptr)
- return (new GiNaC::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
+ if (opt.power_f) {
+ // Invoke the defined power function.
+ 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
+++ for N in range(1, maxargs + 1):
- case @N@:
- return ((power_funcp_@N@)(opt.power_f))(@seq('seq[%(n)d]', N, 0)@, power_param);
+ case @N@:
+ return ((power_funcp_@N@)(opt.power_f))(@seq('seq[%(n)d]', N, 0)@, power_param);
---
- // end of generated lines
+ // end of generated lines
+ }
}
- throw(std::logic_error("function::power(): no power function defined"));
+ // No power function defined? Fall back to returning a power object.
+ return dynallocate<GiNaC::power>(*this, power_param).setflag(status_flags::evaluated);
}
ex function::expand(unsigned options) const
GINAC_ASSERT(serial<registered_functions().size());
const function_options &opt = registered_functions()[serial];
- // No expand defined? Then return the same function with expanded arguments (if required)
- if (opt.expand_f == nullptr) {
- // Only expand arguments when asked to do so
- if (options & expand_options::expand_function_args)
- return inherited::expand(options);
- else
- return (options == 0) ? setflag(status_flags::expanded) : *this;
- }
-
- current_serial = serial;
- if (opt.expand_use_exvector_args)
- return ((expand_funcp_exvector)(opt.expand_f))(seq, options);
- switch (opt.nparams) {
- // the following lines have been generated for max. @maxargs@ parameters
+ if (opt.expand_f) {
+ // Invoke the defined expand function.
+ current_serial = serial;
+ if (opt.expand_use_exvector_args)
+ return ((expand_funcp_exvector)(opt.expand_f))(seq, options);
+ switch (opt.nparams) {
+ // the following lines have been generated for max. @maxargs@ parameters
+++ for N in range(1, maxargs + 1):
- case @N@:
- return ((expand_funcp_@N@)(opt.expand_f))(@seq('seq[%(n)d]', N, 0)@, options);
+ case @N@:
+ return ((expand_funcp_@N@)(opt.expand_f))(@seq('seq[%(n)d]', N, 0)@, options);
---
- // end of generated lines
+ // end of generated lines
+ }
}
- throw(std::logic_error("function::expand(): no expand of function defined"));
+ // No expand function defined? Return the same function with expanded arguments (if required)
+ if (options & expand_options::expand_function_args)
+ return inherited::expand(options);
+ else
+ return (options == 0) ? setflag(status_flags::expanded) : *this;
}
std::vector<function_options> & function::registered_functions()
return serial;
++serial;
}
- throw (std::runtime_error("no function '" + name + "' with " + ToString(nparams) + " parameters defined"));
+ throw (std::runtime_error("no function '" + name + "' with " + std::to_string(nparams) + " parameters defined"));
}
/** Return the print name of the function. */
git://www.ginac.de / ginac.git / blobdiff