]> www.ginac.de Git - ginac.git/blobdiff - ginac/basic.cpp
fixed typo
[ginac.git] / ginac / basic.cpp
index 725f28a2e6b0ea40965ebbc45d50e6ca9d5e7a98..7a0634760354727c77aa06166e1922764b6f8338 100644 (file)
@@ -3,7 +3,7 @@
  *  Implementation of GiNaC's ABC. */
 
 /*
- *  GiNaC Copyright (C) 1999-2003 Johannes Gutenberg University Mainz, Germany
+ *  GiNaC Copyright (C) 1999-2004 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
@@ -55,7 +55,7 @@ GINAC_IMPLEMENT_REGISTERED_CLASS_OPT(basic, void,
 /** basic copy constructor: implicitly assumes that the other class is of
  *  the exact same type (as it's used by duplicate()), so it can copy the
  *  tinfo_key and the hash value. */
-basic::basic(const basic & other) : tinfo_key(other.tinfo_key), flags(other.flags & ~status_flags::dynallocated), hashvalue(other.hashvalue), refcount(0)
+basic::basic(const basic & other) : tinfo_key(other.tinfo_key), flags(other.flags & ~status_flags::dynallocated), hashvalue(other.hashvalue)
 {
        GINAC_ASSERT(typeid(*this) == typeid(other));
 }
@@ -68,13 +68,13 @@ const basic & basic::operator=(const basic & other)
                // The other object is of a derived class, so clear the flags as they
                // might no longer apply (especially hash_calculated). Oh, and don't
                // copy the tinfo_key: it is already set correctly for this object.
-               flags = 0;
+               fl &= ~(status_flags::evaluated | status_flags::expanded | status_flags::hash_calculated);
        } else {
                // The objects are of the exact same class, so copy the hash value.
                hashvalue = other.hashvalue;
        }
        flags = fl;
-       refcount = 0;
+       set_refcount(0);
        return *this;
 }
 
@@ -93,7 +93,7 @@ const basic & basic::operator=(const basic & other)
 //////////
 
 /** Construct object from archive_node. */
-basic::basic(const archive_node &n, lst &sym_lst) : flags(0), refcount(0)
+basic::basic(const archive_node &n, lst &sym_lst) : flags(0)
 {
        // Reconstruct tinfo_key from class name
        std::string class_name;
@@ -140,11 +140,9 @@ void basic::print_dispatch(const registered_class_info & ri, const print_context
        const print_context_class_info * pc_info = &c.get_class_info();
 
 next_class:
-std::clog << "searching class " << reg_info->options.get_name() << std::endl;
        const std::vector<print_functor> & pdt = reg_info->options.get_print_dispatch_table();
 
 next_context:
-std::clog << "searching context " << pc_info->options.get_name() << ", ID " << pc_info->options.get_id() << std::endl;
        unsigned id = pc_info->options.get_id();
        if (id >= pdt.size() || !(pdt[id].is_valid())) {
 
@@ -173,8 +171,6 @@ std::clog << "searching context " << pc_info->options.get_name() << ", ID " << p
        } else {
 
                // Call method
-std::clog << "method found, calling" << std::endl;
-std::clog << " this = " << class_name() << ", context = " << c.class_name() << std::endl;
                pdt[id](*this, c, level);
        }
 }
@@ -188,10 +184,11 @@ void basic::do_print(const print_context & c, unsigned level) const
 /** Tree output to stream. */
 void basic::do_print_tree(const print_tree & c, unsigned level) const
 {
-       c.s << std::string(level, ' ') << class_name()
-           << std::hex << ", hash=0x" << hashvalue << ", flags=0x" << flags << std::dec
-           << ", nops=" << nops()
-           << std::endl;
+       c.s << std::string(level, ' ') << class_name() << " @" << this
+           << std::hex << ", hash=0x" << hashvalue << ", flags=0x" << flags << std::dec;
+       if (nops())
+               c.s << ", nops=" << nops();
+       c.s << std::endl;
        for (size_t i=0; i<nops(); ++i)
                op(i).print(c, level + c.delta_indent);
 }
@@ -552,8 +549,8 @@ bool basic::match(const ex & pattern, lst & repl_lst) const
        if (is_exactly_a<wildcard>(pattern)) {
 
                // Wildcard matches anything, but check whether we already have found
-               // a match for that wildcard first (if so, it the earlier match must
-               // be the same expression)
+               // a match for that wildcard first (if so, the earlier match must be
+               // the same expression)
                for (lst::const_iterator it = repl_lst.begin(); it != repl_lst.end(); ++it) {
                        if (it->op(0).is_equal(pattern))
                                return is_equal(ex_to<basic>(it->op(1)));
@@ -677,6 +674,11 @@ exvector basic::get_free_indices() const
        return exvector(); // return an empty exvector
 }
 
+ex basic::conjugate() const
+{
+       return *this;
+}
+
 ex basic::eval_ncmul(const exvector & v) const
 {
        return hold_ncmul(v);
@@ -806,10 +808,16 @@ ex basic::expand(unsigned options) const
  *  1 greater. */
 int basic::compare(const basic & other) const
 {
+#ifdef GINAC_COMPARE_STATISTICS
+       compare_statistics.total_basic_compares++;
+#endif
        const unsigned hash_this = gethash();
        const unsigned hash_other = other.gethash();
        if (hash_this<hash_other) return -1;
        if (hash_this>hash_other) return 1;
+#ifdef GINAC_COMPARE_STATISTICS
+       compare_statistics.compare_same_hashvalue++;
+#endif
 
        const unsigned typeid_this = tinfo();
        const unsigned typeid_other = other.tinfo();
@@ -825,6 +833,9 @@ int basic::compare(const basic & other) const
 //                     std::cout << std::endl;
 //             }
 //             return cmpval;
+#ifdef GINAC_COMPARE_STATISTICS
+               compare_statistics.compare_same_type++;
+#endif
                return compare_same_type(other);
        } else {
 //             std::cout << "hash collision, different types: " 
@@ -845,13 +856,22 @@ int basic::compare(const basic & other) const
  *  @see is_equal_same_type */
 bool basic::is_equal(const basic & other) const
 {
+#ifdef GINAC_COMPARE_STATISTICS
+       compare_statistics.total_basic_is_equals++;
+#endif
        if (this->gethash()!=other.gethash())
                return false;
+#ifdef GINAC_COMPARE_STATISTICS
+       compare_statistics.is_equal_same_hashvalue++;
+#endif
        if (this->tinfo()!=other.tinfo())
                return false;
        
        GINAC_ASSERT(typeid(*this)==typeid(other));
        
+#ifdef GINAC_COMPARE_STATISTICS
+       compare_statistics.is_equal_same_type++;
+#endif
        return is_equal_same_type(other);
 }
 
@@ -869,7 +889,7 @@ const basic & basic::hold() const
  *  is not the case. */
 void basic::ensure_if_modifiable() const
 {
-       if (refcount > 1)
+       if (get_refcount() > 1)
                throw(std::runtime_error("cannot modify multiply referenced object"));
        clearflag(status_flags::hash_calculated | status_flags::evaluated);
 }
@@ -880,4 +900,27 @@ void basic::ensure_if_modifiable() const
 
 int max_recursion_level = 1024;
 
+
+#ifdef GINAC_COMPARE_STATISTICS
+compare_statistics_t::~compare_statistics_t()
+{
+       std::clog << "ex::compare() called " << total_compares << " times" << std::endl;
+       std::clog << "nontrivial compares: " << nontrivial_compares << " times" << std::endl;
+       std::clog << "basic::compare() called " << total_basic_compares << " times" << std::endl;
+       std::clog << "same hashvalue in compare(): " << compare_same_hashvalue << " times" << std::endl;
+       std::clog << "compare_same_type() called " << compare_same_type << " times" << std::endl;
+       std::clog << std::endl;
+       std::clog << "ex::is_equal() called " << total_is_equals << " times" << std::endl;
+       std::clog << "nontrivial is_equals: " << nontrivial_is_equals << " times" << std::endl;
+       std::clog << "basic::is_equal() called " << total_basic_is_equals << " times" << std::endl;
+       std::clog << "same hashvalue in is_equal(): " << is_equal_same_hashvalue << " times" << std::endl;
+       std::clog << "is_equal_same_type() called " << is_equal_same_type << " times" << std::endl;
+       std::clog << std::endl;
+       std::clog << "basic::gethash() called " << total_gethash << " times" << std::endl;
+       std::clog << "used cached hashvalue " << gethash_cached << " times" << std::endl;
+}
+
+compare_statistics_t compare_statistics;
+#endif
+
 } // namespace GiNaC