* Archiving of GiNaC expressions. */
/*
- * GiNaC Copyright (C) 1999 Johannes Gutenberg University Mainz, Germany
+ * GiNaC Copyright (C) 1999-2001 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 "config.h"
#include "utils.h"
-#ifndef NO_GINAC_NAMESPACE
namespace GiNaC {
-#endif // ndef NO_GINAC_NAMESPACE
-/** Archive an expression.
- * @param ex the expression to be archived
- * @param name name under which the expression is stored */
void archive::archive_ex(const ex &e, const char *name)
{
// Create root node (which recursively archives the whole expression tree)
archive_node_id archive::add_node(const archive_node &n)
{
// Search for node in nodes vector
- vector<archive_node>::const_iterator i = nodes.begin(), iend = nodes.end();
+ std::vector<archive_node>::const_iterator i = nodes.begin(), iend = nodes.end();
archive_node_id id = 0;
while (i != iend) {
if (i->has_same_ex_as(n))
}
-/** Retrieve expression from archive by name.
- * @param sym_lst list of pre-defined symbols */
ex archive::unarchive_ex(const lst &sym_lst, const char *name) const
{
// Find root node
- string name_string = name;
+ std::string name_string = name;
archive_atom id = atomize(name_string);
- vector<archived_ex>::const_iterator i = exprs.begin(), iend = exprs.end();
+ std::vector<archived_ex>::const_iterator i = exprs.begin(), iend = exprs.end();
while (i != iend) {
if (i->name == id)
goto found;
i++;
}
- throw (std::logic_error("expression with name '" + name_string + "' not found in archive"));
+ throw (std::runtime_error("expression with name '" + name_string + "' not found in archive"));
found:
// Recursively unarchive all nodes, starting at the root node
return nodes[i->root].unarchive(sym_lst);
}
-/** Retrieve expression from archive by index.
- * @param sym_lst list of pre-defined symbols */
ex archive::unarchive_ex(const lst &sym_lst, unsigned int index) const
{
if (index >= exprs.size())
return nodes[exprs[index].root].unarchive(sym_lst);
}
-/** Retrieve expression and its name from archive by index.
- * @param sym_lst list of pre-defined symbols */
-ex archive::unarchive_ex(const lst &sym_lst, string &name, unsigned int index) const
+ex archive::unarchive_ex(const lst &sym_lst, std::string &name, unsigned int index) const
{
if (index >= exprs.size())
throw (std::range_error("index of archived expression out of range"));
return nodes[exprs[index].root].unarchive(sym_lst);
}
-
-/** Return number of archived expressions. */
unsigned int archive::num_expressions(void) const
{
return exprs.size();
}
+const archive_node &archive::get_top_node(unsigned int index) const
+{
+ if (index >= exprs.size())
+ throw (std::range_error("index of archived expression out of range"));
+
+ return nodes[exprs[index].root];
+}
+
/*
* Archive file format
* - unsigned root node ID
* - unsigned number of nodes
* - unsigned number of properties
- * - unsigned type (PTYPE_*)
- * - unsigned name atom
+ * - unsigned containing type (PTYPE_*) in its lower 3 bits and
+ * name atom in the upper bits
* - unsigned property value
*
* Unsigned quantities are stored in a compressed format:
*/
/** Write unsigned integer quantity to stream. */
-static void write_unsigned(ostream &os, unsigned int val)
+static void write_unsigned(std::ostream &os, unsigned int val)
{
- while (val > 0x80) {
+ while (val >= 0x80) {
os.put((val & 0x7f) | 0x80);
val >>= 7;
}
}
/** Read unsigned integer quantity from stream. */
-static unsigned int read_unsigned(istream &is)
+static unsigned int read_unsigned(std::istream &is)
{
unsigned char b;
unsigned int ret = 0;
unsigned int shift = 0;
do {
- is.get(b);
+ char b2;
+ is.get(b2);
+ b = b2;
ret |= (b & 0x7f) << shift;
shift += 7;
} while (b & 0x80);
}
/** Write archive_node to binary data stream. */
-ostream &operator<<(ostream &os, const archive_node &n)
+std::ostream &operator<<(std::ostream &os, const archive_node &n)
{
// Write properties
unsigned int num_props = n.props.size();
write_unsigned(os, num_props);
for (unsigned int i=0; i<num_props; i++) {
- write_unsigned(os, n.props[i].type);
- write_unsigned(os, n.props[i].name);
+ write_unsigned(os, n.props[i].type | (n.props[i].name << 3));
write_unsigned(os, n.props[i].value);
}
+ return os;
}
/** Write archive to binary data stream. */
-ostream &operator<<(ostream &os, const archive &ar)
+std::ostream &operator<<(std::ostream &os, const archive &ar)
{
// Write header
os.put('G'); // Signature
unsigned int num_atoms = ar.atoms.size();
write_unsigned(os, num_atoms);
for (unsigned int i=0; i<num_atoms; i++)
- os << ar.atoms[i] << ends;
+ os << ar.atoms[i] << std::ends;
// Write expressions
unsigned int num_exprs = ar.exprs.size();
write_unsigned(os, num_nodes);
for (unsigned int i=0; i<num_nodes; i++)
os << ar.nodes[i];
+ return os;
}
/** Read archive_node from binary data stream. */
-istream &operator>>(istream &is, archive_node &n)
+std::istream &operator>>(std::istream &is, archive_node &n)
{
// Read properties
unsigned int num_props = read_unsigned(is);
n.props.resize(num_props);
for (unsigned int i=0; i<num_props; i++) {
- n.props[i].type = (archive_node::property_type)read_unsigned(is);
- n.props[i].name = read_unsigned(is);
+ unsigned int name_type = read_unsigned(is);
+ n.props[i].type = (archive_node::property_type)(name_type & 7);
+ n.props[i].name = name_type >> 3;
n.props[i].value = read_unsigned(is);
}
+ return is;
}
/** Read archive from binary data stream. */
-istream &operator>>(istream &is, archive &ar)
+std::istream &operator>>(std::istream &is, archive &ar)
{
// Read header
char c1, c2, c3, c4;
ar.nodes.resize(num_nodes, ar);
for (unsigned int i=0; i<num_nodes; i++)
is >> ar.nodes[i];
+ return is;
}
/** Atomize a string (i.e. convert it into an ID number that uniquely
* represents the string). */
-archive_atom archive::atomize(const string &s) const
+archive_atom archive::atomize(const std::string &s) const
{
// Search for string in atoms vector
- vector<string>::const_iterator i = atoms.begin(), iend = atoms.end();
+ std::vector<std::string>::const_iterator i = atoms.begin(), iend = atoms.end();
archive_atom id = 0;
while (i != iend) {
if (*i == s)
}
/** Unatomize a string (i.e. convert the ID number back to the string). */
-const string &archive::unatomize(archive_atom id) const
+const std::string &archive::unatomize(archive_atom id) const
{
if (id >= atoms.size())
throw (std::range_error("archive::unatomizee(): atom ID out of range"));
/** Copy constructor of archive_node. */
archive_node::archive_node(const archive_node &other)
- : a(other.a), props(other.props), has_expression(other.has_expression), e(other.e)
+ : a(other.a), props(other.props), has_expression(other.has_expression), e(other.e)
{
}
/** Recursively construct archive node from expression. */
archive_node::archive_node(archive &ar, const ex &expr)
- : a(ar), has_expression(true), e(expr)
+ : a(ar), has_expression(true), e(expr)
{
expr.bp->archive(*this);
}
}
-/** Add property of type "bool" to node. */
-void archive_node::add_bool(const string &name, bool value)
+void archive_node::add_bool(const std::string &name, bool value)
{
props.push_back(property(a.atomize(name), PTYPE_BOOL, value));
}
-/** Add property of type "unsigned int" to node. */
-void archive_node::add_unsigned(const string &name, unsigned int value)
+void archive_node::add_unsigned(const std::string &name, unsigned int value)
{
props.push_back(property(a.atomize(name), PTYPE_UNSIGNED, value));
}
-/** Add property of type "string" to node. */
-void archive_node::add_string(const string &name, const string &value)
+void archive_node::add_string(const std::string &name, const std::string &value)
{
props.push_back(property(a.atomize(name), PTYPE_STRING, a.atomize(value)));
}
-/** Add property of type "ex" to node. */
-void archive_node::add_ex(const string &name, const ex &value)
+void archive_node::add_ex(const std::string &name, const ex &value)
{
// Recursively create an archive_node and add its ID to the properties of this node
archive_node_id id = a.add_node(archive_node(a, value));
}
-/** Retrieve property of type "bool" from node.
- * @return "true" if property was found, "false" otherwise */
-bool archive_node::find_bool(const string &name, bool &ret) const
+bool archive_node::find_bool(const std::string &name, bool &ret) const
{
archive_atom name_atom = a.atomize(name);
- vector<property>::const_iterator i = props.begin(), iend = props.end();
+ std::vector<property>::const_iterator i = props.begin(), iend = props.end();
while (i != iend) {
if (i->type == PTYPE_BOOL && i->name == name_atom) {
ret = i->value;
return false;
}
-/** Retrieve property of type "unsigned" from node.
- * @return "true" if property was found, "false" otherwise */
-bool archive_node::find_unsigned(const string &name, unsigned int &ret) const
+bool archive_node::find_unsigned(const std::string &name, unsigned int &ret) const
{
archive_atom name_atom = a.atomize(name);
- vector<property>::const_iterator i = props.begin(), iend = props.end();
+ std::vector<property>::const_iterator i = props.begin(), iend = props.end();
while (i != iend) {
if (i->type == PTYPE_UNSIGNED && i->name == name_atom) {
ret = i->value;
return false;
}
-/** Retrieve property of type "string" from node.
- * @return "true" if property was found, "false" otherwise */
-bool archive_node::find_string(const string &name, string &ret) const
+bool archive_node::find_string(const std::string &name, std::string &ret) const
{
archive_atom name_atom = a.atomize(name);
- vector<property>::const_iterator i = props.begin(), iend = props.end();
+ std::vector<property>::const_iterator i = props.begin(), iend = props.end();
while (i != iend) {
if (i->type == PTYPE_STRING && i->name == name_atom) {
ret = a.unatomize(i->value);
return false;
}
-/** Retrieve property of type "ex" from node.
- * @return "true" if property was found, "false" otherwise */
-bool archive_node::find_ex(const string &name, ex &ret, const lst &sym_lst, unsigned int index) const
+bool archive_node::find_ex(const std::string &name, ex &ret, const lst &sym_lst, unsigned int index) const
{
archive_atom name_atom = a.atomize(name);
- vector<property>::const_iterator i = props.begin(), iend = props.end();
+ std::vector<property>::const_iterator i = props.begin(), iend = props.end();
unsigned int found_index = 0;
while (i != iend) {
if (i->type == PTYPE_NODE && i->name == name_atom) {
return true;
}
+const archive_node &archive_node::find_ex_node(const std::string &name, unsigned int index) const
+{
+ archive_atom name_atom = a.atomize(name);
+ std::vector<property>::const_iterator i = props.begin(), iend = props.end();
+ unsigned int found_index = 0;
+ while (i != iend) {
+ if (i->type == PTYPE_NODE && i->name == name_atom) {
+ if (found_index == index)
+ return a.get_node(i->value);
+ found_index++;
+ }
+ i++;
+ }
+ throw (std::runtime_error("property with name '" + name + "' not found in archive node"));
+}
+
+
+void archive_node::get_properties(propinfovector &v) const
+{
+ v.clear();
+ std::vector<property>::const_iterator i = props.begin(), iend = props.end();
+ while (i != iend) {
+ property_type type = i->type;
+ std::string name = a.unatomize(i->name);
+
+ propinfovector::iterator a = v.begin(), aend = v.end();
+ bool found = false;
+ while (a != aend) {
+ if (a->type == type && a->name == name) {
+ a->count++;
+ found = true;
+ break;
+ }
+ a++;
+ }
+ if (!found)
+ v.push_back(property_info(type, name));
+ i++;
+ }
+}
+
/** Convert archive node to GiNaC expression. */
ex archive_node::unarchive(const lst &sym_lst) const
return e;
// Find instantiation function for class specified in node
- string class_name;
+ std::string class_name;
if (!find_string("class", class_name))
throw (std::runtime_error("archive node contains no class name"));
unarch_func f = find_unarch_func(class_name);
}
+/** Assignment operator of property_info. */
+const archive_node::property_info &archive_node::property_info::operator=(const property_info &other)
+{
+ if (this != &other) {
+ type = other.type;
+ name = other.name;
+ count = other.count;
+ }
+ return *this;
+}
+
/** Assignment operator of property. */
const archive_node::property &archive_node::property::operator=(const property &other)
{
}
-/** Clear all archived expressions. */
void archive::clear(void)
{
atoms.clear();
/** Delete cached unarchived expressions in all archive_nodes (mainly for debugging). */
void archive::forget(void)
{
- vector<archive_node>::iterator i = nodes.begin(), iend = nodes.end();
+ std::vector<archive_node>::iterator i = nodes.begin(), iend = nodes.end();
while (i != iend) {
i->forget();
i++;
/** Print archive to stream in ugly raw format (for debugging). */
-void archive::printraw(ostream &os) const
+void archive::printraw(std::ostream &os) const
{
// Dump atoms
os << "Atoms:\n";
{
- vector<string>::const_iterator i = atoms.begin(), iend = atoms.end();
+ std::vector<std::string>::const_iterator i = atoms.begin(), iend = atoms.end();
archive_atom id = 0;
while (i != iend) {
- os << " " << id << " " << *i << endl;
+ os << " " << id << " " << *i << std::endl;
i++; id++;
}
}
- os << endl;
+ os << std::endl;
// Dump expressions
os << "Expressions:\n";
{
- vector<archived_ex>::const_iterator i = exprs.begin(), iend = exprs.end();
+ std::vector<archived_ex>::const_iterator i = exprs.begin(), iend = exprs.end();
unsigned int index = 0;
while (i != iend) {
- os << " " << index << " \"" << unatomize(i->name) << "\" root node " << i->root << endl;
+ os << " " << index << " \"" << unatomize(i->name) << "\" root node " << i->root << std::endl;
i++; index++;
}
}
- os << endl;
+ os << std::endl;
// Dump nodes
os << "Nodes:\n";
{
- vector<archive_node>::const_iterator i = nodes.begin(), iend = nodes.end();
+ std::vector<archive_node>::const_iterator i = nodes.begin(), iend = nodes.end();
archive_node_id id = 0;
while (i != iend) {
os << " " << id << " ";
}
/** Output archive_node to stream in ugly raw format (for debugging). */
-void archive_node::printraw(ostream &os) const
+void archive_node::printraw(std::ostream &os) const
{
// Dump cached unarchived expression
if (has_expression)
os << "(basic * " << e.bp << " = " << e << ")\n";
else
- os << "(no expression)\n";
+ os << "\n";
// Dump properties
- vector<property>::const_iterator i = props.begin(), iend = props.end();
+ std::vector<property>::const_iterator i = props.begin(), iend = props.end();
while (i != iend) {
os << " ";
switch (i->type) {
case PTYPE_NODE: os << "node"; break;
default: os << "<unknown>"; break;
}
- os << " \"" << a.unatomize(i->name) << "\" " << i->value << endl;
+ os << " \"" << a.unatomize(i->name) << "\" " << i->value << std::endl;
i++;
}
}
+/** Create a dummy archive. The intention is to fill archive_node's default
+ * ctor, which is currently a Cint-requirement. */
+archive* archive_node::dummy_ar_creator(void)
+{
+ static archive* some_ar = new archive;
+ return some_ar;
+}
+
-#ifndef NO_GINAC_NAMESPACE
} // namespace GiNaC
-#endif // ndef NO_GINAC_NAMESPACE