X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Ffderivative.cpp;h=2ce60fb3efddf571546bd401a479cd80041e5e91;hp=17817e108671ffd4bc612ee343d78937455e2c49;hb=ac0ef893e48af99dc632adcdf21edd6a402327b3;hpb=6c946d4c762f5a0d6a3b742f03556dd018d63886 diff --git a/ginac/fderivative.cpp b/ginac/fderivative.cpp index 17817e10..2ce60fb3 100644 --- a/ginac/fderivative.cpp +++ b/ginac/fderivative.cpp @@ -3,7 +3,7 @@ * Implementation of abstract derivatives of functions. */ /* - * GiNaC Copyright (C) 1999-2015 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2016 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 @@ -81,7 +81,7 @@ GINAC_BIND_UNARCHIVER(fderivative); void fderivative::archive(archive_node &n) const { inherited::archive(n); - paramset::const_iterator i = parameter_set.begin(), end = parameter_set.end(); + auto i = parameter_set.begin(), end = parameter_set.end(); while (i != end) { n.add_unsigned("param", *i); ++i; @@ -102,7 +102,7 @@ void fderivative::print(const print_context & c, unsigned level) const void fderivative::do_print(const print_context & c, unsigned level) const { c.s << "D["; - paramset::const_iterator i = parameter_set.begin(), end = parameter_set.end(); + auto i = parameter_set.begin(), end = parameter_set.end(); --end; while (i != end) { c.s << *i++ << ","; @@ -114,7 +114,7 @@ void fderivative::do_print(const print_context & c, unsigned level) const void fderivative::do_print_csrc(const print_csrc & c, unsigned level) const { c.s << "D_"; - paramset::const_iterator i = parameter_set.begin(), end = parameter_set.end(); + auto i = parameter_set.begin(), end = parameter_set.end(); --end; while (i != end) c.s << *i++ << "_"; @@ -129,13 +129,13 @@ void fderivative::do_print_tree(const print_tree & c, unsigned level) const << std::hex << ", hash=0x" << hashvalue << ", flags=0x" << flags << std::dec << ", nops=" << nops() << ", params="; - paramset::const_iterator i = parameter_set.begin(), end = parameter_set.end(); + auto i = parameter_set.begin(), end = parameter_set.end(); --end; while (i != end) c.s << *i++ << ","; c.s << *i << std::endl; - for (size_t i=0; ihold(); } -/** Numeric evaluation falls back to evaluation of arguments. - * @see basic::evalf */ -ex fderivative::evalf(int level) const -{ - return basic::evalf(level); -} - /** The series expansion of derivatives falls back to Taylor expansion. * @see basic::series */ ex fderivative::series(const relational & r, int order, unsigned options) const @@ -222,4 +215,20 @@ bool fderivative::match_same_type(const basic & other) const return parameter_set == o.parameter_set && inherited::match_same_type(other); } +/** Expose this object's derivative structure. + * + * Parameter numbers occurring more than once stand for repeated + * differentiation with respect to that parameter. If a symbolic function + * f(x,y) is differentiated with respect to x, this method will return {0}. + * If f(x,y) is differentiated twice with respect to y, it will return {1,1}. + * (This corresponds to the way this object is printed.) + * + * @return multiset of function's parameter numbers that are abstractly + * differentiated. */ +const paramset& fderivative::derivatives() const +{ + return parameter_set; +} + + } // namespace GiNaC