* Wrapper template for making GiNaC classes out of STL containers. */
/*
- * GiNaC Copyright (C) 1999-2015 Johannes Gutenberg University Mainz, Germany
+ * GiNaC Copyright (C) 1999-2020 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
#include "print.h"
#include "archive.h"
#include "assertion.h"
+#include "compiler.h"
#include <algorithm>
#include <iterator>
container_storage() {}
container_storage(size_t n, const ex & e) : seq(n, e) {}
+ container_storage(std::initializer_list<ex> il) : seq(il) {}
template <class In>
container_storage(In b, In e) : seq(b, e) {}
inline void container_storage<std::vector>::reserve(std::vector<ex> & v, size_t n) { v.reserve(n); }
-/** Helper template to allow initialization of containers via an overloaded
- * comma operator (idea stolen from Blitz++). */
-template <typename T, typename STLT>
-class container_init {
-public:
- container_init(STLT & s) : stlt(s) {}
-
- container_init<T, STLT> operator,(const T & x)
- {
- stlt.push_back(x);
- return container_init<T, STLT>(stlt);
- }
-
- // The following specializations produce much tighter code than the
- // general case above
-
- container_init<T, STLT> operator,(int x)
- {
- stlt.push_back(x);
- return container_init<T, STLT>(stlt);
- }
-
- container_init<T, STLT> operator,(unsigned int x)
- {
- stlt.push_back(x);
- return container_init<T, STLT>(stlt);
- }
-
- container_init<T, STLT> operator,(long x)
- {
- stlt.push_back(x);
- return container_init<T, STLT>(stlt);
- }
-
- container_init<T, STLT> operator,(unsigned long x)
- {
- stlt.push_back(x);
- return container_init<T, STLT>(stlt);
- }
-
- container_init<T, STLT> operator,(double x)
- {
- stlt.push_back(x);
- return container_init<T, STLT>(stlt);
- }
-
- container_init<T, STLT> operator,(const symbol & x)
- {
- stlt.push_back(T(x));
- return container_init<T, STLT>(stlt);
- }
-
-private:
- container_init();
- STLT & stlt;
-};
-
/** Wrapper template for making GiNaC classes out of STL containers. */
template <template <class T, class = std::allocator<T>> class C>
class container : public basic, public container_storage<C> {
setflag(get_default_flags());
}
- explicit container(const ex & p1)
- : container_storage<C>(1, p1)
- {
- setflag(get_default_flags());
- }
-
- container(const ex & p1, const ex & p2)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 2);
- this->seq.push_back(p1); this->seq.push_back(p2);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 3);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3,
- const ex & p4)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 4);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- this->seq.push_back(p4);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3,
- const ex & p4, const ex & p5)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 5);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- this->seq.push_back(p4); this->seq.push_back(p5);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3,
- const ex & p4, const ex & p5, const ex & p6)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 6);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- this->seq.push_back(p4); this->seq.push_back(p5); this->seq.push_back(p6);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3,
- const ex & p4, const ex & p5, const ex & p6,
- const ex & p7)
+ container(std::initializer_list<ex> il)
+ : container_storage<C>(il)
{
setflag(get_default_flags());
- this->reserve(this->seq, 7);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- this->seq.push_back(p4); this->seq.push_back(p5); this->seq.push_back(p6);
- this->seq.push_back(p7);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3,
- const ex & p4, const ex & p5, const ex & p6,
- const ex & p7, const ex & p8)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 8);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- this->seq.push_back(p4); this->seq.push_back(p5); this->seq.push_back(p6);
- this->seq.push_back(p7); this->seq.push_back(p8);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3,
- const ex & p4, const ex & p5, const ex & p6,
- const ex & p7, const ex & p8, const ex & p9)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 9);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- this->seq.push_back(p4); this->seq.push_back(p5); this->seq.push_back(p6);
- this->seq.push_back(p7); this->seq.push_back(p8); this->seq.push_back(p9);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3,
- const ex & p4, const ex & p5, const ex & p6,
- const ex & p7, const ex & p8, const ex & p9,
- const ex & p10)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 10);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- this->seq.push_back(p4); this->seq.push_back(p5); this->seq.push_back(p6);
- this->seq.push_back(p7); this->seq.push_back(p8); this->seq.push_back(p9);
- this->seq.push_back(p10);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3,
- const ex & p4, const ex & p5, const ex & p6,
- const ex & p7, const ex & p8, const ex & p9,
- const ex & p10, const ex & p11)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 11);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- this->seq.push_back(p4); this->seq.push_back(p5); this->seq.push_back(p6);
- this->seq.push_back(p7); this->seq.push_back(p8); this->seq.push_back(p9);
- this->seq.push_back(p10); this->seq.push_back(p11);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3,
- const ex & p4, const ex & p5, const ex & p6,
- const ex & p7, const ex & p8, const ex & p9,
- const ex & p10, const ex & p11, const ex & p12)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 12);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- this->seq.push_back(p4); this->seq.push_back(p5); this->seq.push_back(p6);
- this->seq.push_back(p7); this->seq.push_back(p8); this->seq.push_back(p9);
- this->seq.push_back(p10); this->seq.push_back(p11); this->seq.push_back(p12);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3,
- const ex & p4, const ex & p5, const ex & p6,
- const ex & p7, const ex & p8, const ex & p9,
- const ex & p10, const ex & p11, const ex & p12,
- const ex & p13)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 13);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- this->seq.push_back(p4); this->seq.push_back(p5); this->seq.push_back(p6);
- this->seq.push_back(p7); this->seq.push_back(p8); this->seq.push_back(p9);
- this->seq.push_back(p10); this->seq.push_back(p11); this->seq.push_back(p12);
- this->seq.push_back(p13);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3,
- const ex & p4, const ex & p5, const ex & p6,
- const ex & p7, const ex & p8, const ex & p9,
- const ex & p10, const ex & p11, const ex & p12,
- const ex & p13, const ex & p14)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 14);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- this->seq.push_back(p4); this->seq.push_back(p5); this->seq.push_back(p6);
- this->seq.push_back(p7); this->seq.push_back(p8); this->seq.push_back(p9);
- this->seq.push_back(p10); this->seq.push_back(p11); this->seq.push_back(p12);
- this->seq.push_back(p13); this->seq.push_back(p14);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3,
- const ex & p4, const ex & p5, const ex & p6,
- const ex & p7, const ex & p8, const ex & p9,
- const ex & p10, const ex & p11, const ex & p12,
- const ex & p13, const ex & p14, const ex & p15)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 15);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- this->seq.push_back(p4); this->seq.push_back(p5); this->seq.push_back(p6);
- this->seq.push_back(p7); this->seq.push_back(p8); this->seq.push_back(p9);
- this->seq.push_back(p10); this->seq.push_back(p11); this->seq.push_back(p12);
- this->seq.push_back(p13); this->seq.push_back(p14); this->seq.push_back(p15);
- }
-
- container(const ex & p1, const ex & p2, const ex & p3,
- const ex & p4, const ex & p5, const ex & p6,
- const ex & p7, const ex & p8, const ex & p9,
- const ex & p10, const ex & p11, const ex & p12,
- const ex & p13, const ex & p14, const ex & p15,
- const ex & p16)
- {
- setflag(get_default_flags());
- this->reserve(this->seq, 16);
- this->seq.push_back(p1); this->seq.push_back(p2); this->seq.push_back(p3);
- this->seq.push_back(p4); this->seq.push_back(p5); this->seq.push_back(p6);
- this->seq.push_back(p7); this->seq.push_back(p8); this->seq.push_back(p9);
- this->seq.push_back(p10); this->seq.push_back(p11); this->seq.push_back(p12);
- this->seq.push_back(p13); this->seq.push_back(p14); this->seq.push_back(p15);
- this->seq.push_back(p16);
- }
-
- // First step of initialization of container with a comma-separated
- // sequence of expressions. Subsequent steps are handled by
- // container_init<>::operator,().
- container_init<ex, STLT> operator=(const ex & x)
- {
- this->seq.push_back(x);
- return container_init<ex, STLT>(this->seq);
}
// functions overriding virtual functions from base classes
public:
- bool info(unsigned inf) const { return inherited::info(inf); }
- unsigned precedence() const { return 10; }
- size_t nops() const { return this->seq.size(); }
- ex op(size_t i) const;
- ex & let_op(size_t i);
- ex eval(int level = 0) const;
- ex subs(const exmap & m, unsigned options = 0) const;
-
- void read_archive(const archive_node &n, lst &sym_lst)
+ bool info(unsigned inf) const override { return inherited::info(inf); }
+ unsigned precedence() const override { return 10; }
+ size_t nops() const override { return this->seq.size(); }
+ ex op(size_t i) const override;
+ ex & let_op(size_t i) override;
+ ex subs(const exmap & m, unsigned options = 0) const override;
+
+ void read_archive(const archive_node &n, lst &sym_lst) override
{
inherited::read_archive(n, sym_lst);
setflag(get_default_flags());
- archive_node::archive_node_cit first = n.find_first("seq");
- archive_node::archive_node_cit last = n.find_last("seq");
- ++last;
- this->reserve(this->seq, last - first);
- for (archive_node::archive_node_cit i=first; i<last; ++i) {
+ auto range = n.find_property_range("seq", "seq");
+ this->reserve(this->seq, range.end - range.begin);
+ for (archive_node::archive_node_cit i=range.begin; i<range.end; ++i) {
ex e;
n.find_ex_by_loc(i, e, sym_lst);
- this->seq.push_back(e);
+ this->seq.emplace_back(e);
}
}
/** Archive the object. */
- void archive(archive_node &n) const
+ void archive(archive_node &n) const override
{
inherited::archive(n);
for (auto & i : this->seq) {
}
protected:
- ex conjugate() const
+ ex conjugate() const override
{
STLT *newcont = nullptr;
for (const_iterator i=this->seq.begin(); i!=this->seq.end(); ++i) {
return *this;
}
- ex real_part() const
+ ex real_part() const override
{
STLT cont;
this->reserve(cont, nops());
return thiscontainer(cont);
}
- ex imag_part() const
+ ex imag_part() const override
{
STLT cont;
this->reserve(cont, nops());
return thiscontainer(cont);
}
- bool is_equal_same_type(const basic & other) const;
+ bool is_equal_same_type(const basic & other) const override;
// new virtual functions which can be overridden by derived classes
protected:
void do_print_tree(const print_tree & c, unsigned level) const;
void do_print_python(const print_python & c, unsigned level) const;
void do_print_python_repr(const print_python_repr & c, unsigned level) const;
- STLT evalchildren(int level) const;
STLT subschildren(const exmap & m, unsigned options = 0) const;
};
setflag(get_default_flags());
}
-
template <template <class T, class = std::allocator<T>> class C>
void container<C>::do_print(const print_context & c, unsigned level) const
{
return *it;
}
-template <template <class T, class = std::allocator<T>> class C>
-ex container<C>::eval(int level) const
-{
- if (level == 1)
- return hold();
- else
- return thiscontainer(evalchildren(level));
-}
-
template <template <class T, class = std::allocator<T>> class C>
ex container<C>::subs(const exmap & m, unsigned options) const
{
c.s << closebracket;
}
-template <template <class T, class = std::allocator<T>> class C>
-typename container<C>::STLT container<C>::evalchildren(int level) const
-{
- if (level == 1)
- return this->seq;
- else if (level == -max_recursion_level)
- throw std::runtime_error("max recursion level reached");
-
- STLT s;
- this->reserve(s, this->seq.size());
-
- --level;
- const_iterator it = this->seq.begin(), itend = this->seq.end();
- while (it != itend) {
- s.push_back(it->eval(level));
- ++it;
- }
-
- return s;
-}
-
template <template <class T, class = std::allocator<T>> class C>
typename container<C>::STLT container<C>::subschildren(const exmap & m, unsigned options) const
{