3 * Wrapper template for making GiNaC classes out of C++ structures. */
6 * GiNaC Copyright (C) 1999-2007 Johannes Gutenberg University Mainz, Germany
8 * This program is free software; you can redistribute it and/or modify
9 * it under the terms of the GNU General Public License as published by
10 * the Free Software Foundation; either version 2 of the License, or
11 * (at your option) any later version.
13 * This program is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 * GNU General Public License for more details.
18 * You should have received a copy of the GNU General Public License
19 * along with this program; if not, write to the Free Software
20 * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
23 #ifndef __GINAC_STRUCTURE_H__
24 #define __GINAC_STRUCTURE_H__
31 #include "operators.h"
37 /** Comparison policy: all structures of one type are equal */
39 class compare_all_equal {
41 static bool struct_is_equal(const T * t1, const T * t2) { return true; }
42 static int struct_compare(const T * t1, const T * t2) { return 0; }
44 // disallow destruction of structure through a compare_all_equal*
46 ~compare_all_equal() {}
50 /** Comparison policy: use std::equal_to/std::less (defaults to operators
51 * == and <) to compare structures. */
53 class compare_std_less {
55 static bool struct_is_equal(const T * t1, const T * t2)
57 return std::equal_to<T>()(*t1, *t2);
60 static int struct_compare(const T * t1, const T * t2)
62 if (std::less<T>()(*t1, *t2))
64 else if (std::less<T>()(*t2, *t1))
70 // disallow destruction of structure through a compare_std_less*
72 ~compare_std_less() {}
76 /** Comparison policy: use bit-wise comparison to compare structures. */
78 class compare_bitwise {
80 static bool struct_is_equal(const T * t1, const T * t2)
82 const char * cp1 = reinterpret_cast<const char *>(t1);
83 const char * cp2 = reinterpret_cast<const char *>(t2);
85 return std::equal(cp1, cp1 + sizeof(T), cp2);
88 static int struct_compare(const T * t1, const T * t2)
90 const unsigned char * cp1 = reinterpret_cast<const unsigned char *>(t1);
91 const unsigned char * cp2 = reinterpret_cast<const unsigned char *>(t2);
92 typedef std::pair<const unsigned char *, const unsigned char *> cppair;
94 cppair res = std::mismatch(cp1, cp1 + sizeof(T), cp2);
96 if (res.first == cp1 + sizeof(T))
98 else if (*res.first < *res.second)
104 // disallow destruction of structure through a compare_bitwise*
106 ~compare_bitwise() {}
110 // Select default comparison policy
111 template <class T, template <class> class ComparisonPolicy = compare_all_equal> class structure;
114 /** Wrapper template for making GiNaC classes out of C++ structures. */
115 template <class T, template <class> class ComparisonPolicy>
116 class structure : public basic, public ComparisonPolicy<T> {
117 GINAC_DECLARE_REGISTERED_CLASS(structure, basic)
120 static tinfo_t get_tinfo() { return reg_info.options.get_id(); }
121 static const char *get_class_name() { return "structure"; }
122 static tinfo_t next_structure_tinfo_key() { return new tinfo_static_t; }
126 /** Construct structure as a copy of a given C++ structure. */
127 structure(const T & t) : inherited(get_tinfo()), obj(t) { }
129 // functions overriding virtual functions from base classes
130 // All these are just defaults that can be specialized by the user
133 ex eval(int level = 0) const { return hold(); }
134 ex evalf(int level = 0) const { return inherited::evalf(level); }
135 ex evalm() const { return inherited::evalm(); }
137 ex eval_ncmul(const exvector & v) const { return hold_ncmul(v); }
139 ex eval_indexed(const basic & i) const { return i.hold(); }
142 void print(const print_context & c, unsigned level = 0) const { inherited::print(c, level); }
143 unsigned precedence() const { return 70; }
146 bool info(unsigned inf) const { return false; }
149 size_t nops() const { return 0; }
150 ex op(size_t i) const { return inherited::op(i); }
151 ex operator[](const ex & index) const { return inherited::operator[](index); }
152 ex operator[](size_t i) const { return inherited::operator[](i); }
153 ex & let_op(size_t i) { return inherited::let_op(i); }
154 ex & operator[](const ex & index) { return inherited::operator[](index); }
155 ex & operator[](size_t i) { return inherited::operator[](i); }
158 bool has(const ex & other, unsigned options = 0) const { return inherited::has(other, options); }
159 bool match(const ex & pattern, lst & repl_lst) const { return inherited::match(pattern, repl_lst); }
161 bool match_same_type(const basic & other) const { return true; }
165 ex subs(const exmap & m, unsigned options = 0) const { return inherited::subs(m, options); }
168 ex map(map_function & f) const { return inherited::map(f); }
171 int degree(const ex & s) const { return inherited::degree(s); }
172 int ldegree(const ex & s) const { return inherited::ldegree(s); }
173 ex coeff(const ex & s, int n = 1) const { return inherited::coeff(s, n); }
176 ex expand(unsigned options = 0) const { return inherited::expand(options); }
177 ex collect(const ex & s, bool distributed = false) const { return inherited::collect(s, distributed); }
179 // differentiation and series expansion
181 ex derivative(const symbol & s) const { return inherited::derivative(s); }
183 ex series(const relational & r, int order, unsigned options = 0) const { return inherited::series(r, order, options); }
185 // rational functions
186 ex normal(exmap & repl, exmap & rev_lookup, int level = 0) const { return inherited::normal(repl, rev_lookup, level); }
187 ex to_rational(exmap & repl) const { return inherited::to_rational(repl); }
188 ex to_polynomial(exmap & repl) const { return inherited::to_polynomial(repl); }
190 // polynomial algorithms
191 numeric integer_content() const { return 1; }
192 ex smod(const numeric & xi) const { return *this; }
193 numeric max_coefficient() const { return 1; }
196 exvector get_free_indices() const { return exvector(); }
197 ex add_indexed(const ex & self, const ex & other) const { return self + other; }
198 ex scalar_mul_indexed(const ex & self, const numeric & other) const { return self * ex(other); }
199 bool contract_with(exvector::iterator self, exvector::iterator other, exvector & v) const { return false; }
202 unsigned return_type() const { return return_types::commutative; }
203 tinfo_t return_type_tinfo() const { return this; }
206 bool is_equal_same_type(const basic & other) const
208 GINAC_ASSERT(is_a<structure>(other));
209 const structure & o = static_cast<const structure &>(other);
211 return struct_is_equal(&obj, &o.obj);
214 unsigned calchash() const { return inherited::calchash(); }
216 // non-virtual functions in this class
218 // access to embedded structure
219 const T *operator->() const { return &obj; }
220 T &get_struct() { return obj; }
221 const T &get_struct() const { return obj; }
227 /** Default constructor */
228 template <class T, template <class> class CP>
229 structure<T, CP>::structure() : inherited(get_tinfo()) { }
231 /** Construct object from archive_node. */
232 template <class T, template <class> class CP>
233 structure<T, CP>::structure(const archive_node &n, lst &sym_lst) : inherited(n, sym_lst) {}
235 /** Unarchive the object. */
236 template <class T, template <class> class CP>
237 ex structure<T, CP>::unarchive(const archive_node &n, lst &sym_lst)
239 return (new structure(n, sym_lst))->setflag(status_flags::dynallocated);
242 /** Archive the object. */
243 template <class T, template <class> class CP>
244 void structure<T, CP>::archive(archive_node &n) const
246 inherited::archive(n);
249 /** Compare two structures of the same type. */
250 template <class T, template <class> class CP>
251 int structure<T, CP>::compare_same_type(const basic & other) const
253 GINAC_ASSERT(is_a<structure>(other));
254 const structure & o = static_cast<const structure &>(other);
256 return struct_compare(&obj, &o.obj);
259 template <class T, template <class> class CP>
260 registered_class_info structure<T, CP>::reg_info = registered_class_info(registered_class_options(structure::get_class_name(), "basic", structure::next_structure_tinfo_key(), &structure::unarchive));
265 #endif // ndef __GINAC_STRUCTURE_H__