]> www.ginac.de Git - ginac.git/blobdiff - ginac/expairseq.cpp
* (time_fateman_expand.cpp): new file.
[ginac.git] / ginac / expairseq.cpp
index c2dfc6a57cebd26363980b574739c89a7003d8f2..a9c3f12336b2d54a818fe3e9e076cb03bf629375 100644 (file)
@@ -3,7 +3,7 @@
  *  Implementation of sequences of expression pairs. */
 
 /*
- *  GiNaC Copyright (C) 1999-2001 Johannes Gutenberg University Mainz, Germany
+ *  GiNaC Copyright (C) 1999-2002 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
@@ -20,6 +20,7 @@
  *  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
  */
 
+#include <iostream>
 #include <algorithm>
 #include <string>
 #include <stdexcept>
@@ -32,7 +33,6 @@
 #include "wildcard.h"
 #include "print.h"
 #include "archive.h"
-#include "debugmsg.h"
 #include "utils.h"
 
 #if EXPAIRSEQ_USE_HASHTAB
@@ -58,20 +58,18 @@ public:
 };
 
 //////////
-// default ctor, dtor, copy ctor assignment operator and helpers
+// default ctor, dtor, copy ctor, assignment operator and helpers
 //////////
 
 // public
 
 expairseq::expairseq(const expairseq &other)
 {
-       debugmsg("expairseq copy ctor",LOGLEVEL_CONSTRUCT);
        copy(other);
 }
 
 const expairseq &expairseq::operator=(const expairseq &other)
 {
-       debugmsg("expairseq operator=",LOGLEVEL_ASSIGNMENT);
        if (this != &other) {
                destroy(true);
                copy(other);
@@ -115,14 +113,12 @@ DEFAULT_DESTROY(expairseq)
 
 expairseq::expairseq(const ex &lh, const ex &rh) : inherited(TINFO_expairseq)
 {
-       debugmsg("expairseq ctor from ex,ex",LOGLEVEL_CONSTRUCT);
        construct_from_2_ex(lh,rh);
        GINAC_ASSERT(is_canonical());
 }
 
 expairseq::expairseq(const exvector &v) : inherited(TINFO_expairseq)
 {
-       debugmsg("expairseq ctor from exvector",LOGLEVEL_CONSTRUCT);
        construct_from_exvector(v);
        GINAC_ASSERT(is_canonical());
 }
@@ -130,7 +126,7 @@ expairseq::expairseq(const exvector &v) : inherited(TINFO_expairseq)
 expairseq::expairseq(const epvector &v, const ex &oc)
   : inherited(TINFO_expairseq), overall_coeff(oc)
 {
-       debugmsg("expairseq ctor from epvector,ex",LOGLEVEL_CONSTRUCT);
+       GINAC_ASSERT(is_a<numeric>(oc));
        construct_from_epvector(v);
        GINAC_ASSERT(is_canonical());
 }
@@ -138,8 +134,8 @@ expairseq::expairseq(const epvector &v, const ex &oc)
 expairseq::expairseq(epvector *vp, const ex &oc)
   : inherited(TINFO_expairseq), overall_coeff(oc)
 {
-       debugmsg("expairseq ctor from epvector *,ex",LOGLEVEL_CONSTRUCT);
        GINAC_ASSERT(vp!=0);
+       GINAC_ASSERT(is_a<numeric>(oc));
        construct_from_epvector(*vp);
        delete vp;
        GINAC_ASSERT(is_canonical());
@@ -154,7 +150,6 @@ expairseq::expairseq(const archive_node &n, const lst &sym_lst) : inherited(n, s
        , hashtabsize(0)
 #endif
 {
-       debugmsg("expairseq ctor from archive_node", LOGLEVEL_CONSTRUCT);
        for (unsigned int i=0; true; i++) {
                ex rest;
                ex coeff;
@@ -188,15 +183,12 @@ DEFAULT_UNARCHIVE(expairseq)
 
 basic *expairseq::duplicate() const
 {
-       debugmsg("expairseq duplicate",LOGLEVEL_DUPLICATE);
        return new expairseq(*this);
 }
 
 void expairseq::print(const print_context &c, unsigned level) const
 {
-       debugmsg("expairseq print",LOGLEVEL_PRINT);
-
-       if (is_of_type(c, print_tree)) {
+       if (is_a<print_tree>(c)) {
 
                unsigned delta_indent = static_cast<const print_tree &>(c).delta_indent;
 
@@ -346,7 +338,7 @@ bool expairseq::match(const ex & pattern, lst & repl_lst) const
                bool has_global_wildcard = false;
                ex global_wildcard;
                for (unsigned int i=0; i<pattern.nops(); i++) {
-                       if (is_ex_exactly_of_type(pattern.op(i), wildcard)) {
+                       if (is_exactly_a<wildcard>(pattern.op(i))) {
                                has_global_wildcard = true;
                                global_wildcard = pattern.op(i);
                                break;
@@ -412,7 +404,7 @@ ex expairseq::subs(const lst &ls, const lst &lr, bool no_pattern) const
 {
        epvector *vp = subschildren(ls, lr, no_pattern);
        if (vp)
-               return ex_to<basic>(thisexpairseq(vp, overall_coeff)).basic::subs(ls, lr, no_pattern);
+               return ex_to<basic>(thisexpairseq(vp, overall_coeff));
        else
                return basic::subs(ls, lr, no_pattern);
 }
@@ -651,7 +643,7 @@ void expairseq::printseq(const print_context & c, char delim,
  *  @see expairseq::recombine_pair_to_ex() */
 expair expairseq::split_ex_to_pair(const ex &e) const
 {
-       return expair(e,_ex1());
+       return expair(e,_ex1);
 }
 
 
@@ -691,7 +683,7 @@ bool expairseq::expair_needs_further_processing(epp it)
 
 ex expairseq::default_overall_coeff(void) const
 {
-       return _ex0();
+       return _ex0;
 }
 
 void expairseq::combine_overall_coeff(const ex &c)
@@ -785,8 +777,8 @@ void expairseq::construct_from_2_ex(const ex &lh, const ex &rh)
        hashtabsize = 0;
 #endif // EXPAIRSEQ_USE_HASHTAB
        
-       if (is_ex_exactly_of_type(lh,numeric)) {
-               if (is_ex_exactly_of_type(rh,numeric)) {
+       if (is_exactly_a<numeric>(lh)) {
+               if (is_exactly_a<numeric>(rh)) {
                        combine_overall_coeff(lh);
                        combine_overall_coeff(rh);
                } else {
@@ -794,7 +786,7 @@ void expairseq::construct_from_2_ex(const ex &lh, const ex &rh)
                        seq.push_back(split_ex_to_pair(rh));
                }
        } else {
-               if (is_ex_exactly_of_type(rh,numeric)) {
+               if (is_exactly_a<numeric>(rh)) {
                        combine_overall_coeff(rh);
                        seq.push_back(split_ex_to_pair(lh));
                } else {
@@ -881,7 +873,7 @@ void expairseq::construct_from_expairseq_ex(const expairseq &s,
                                                                                        const ex &e)
 {
        combine_overall_coeff(s.overall_coeff);
-       if (is_ex_exactly_of_type(e,numeric)) {
+       if (is_exactly_a<numeric>(e)) {
                combine_overall_coeff(e);
                seq = s.seq;
                return;
@@ -1006,7 +998,7 @@ void expairseq::make_flat(const exvector &v)
                                ++cit_s;
                        }
                } else {
-                       if (is_ex_exactly_of_type(*cit,numeric))
+                       if (is_exactly_a<numeric>(*cit))
                                combine_overall_coeff(*cit);
                        else
                                seq.push_back(split_ex_to_pair(*cit));
@@ -1067,7 +1059,7 @@ void expairseq::make_flat(const epvector &v)
 /** Brings this expairseq into a sorted (canonical) form. */
 void expairseq::canonicalize(void)
 {
-       sort(seq.begin(), seq.end(), expair_is_less());
+       std::sort(seq.begin(), seq.end(), expair_rest_is_less());
 }
 
 
@@ -1257,7 +1249,7 @@ void expairseq::build_hashtab_and_combine(epvector::iterator &first_numeric,
        epp current = seq.begin();
 
        while (current!=first_numeric) {
-               if (is_ex_exactly_of_type(current->rest,numeric)) {
+               if (is_exactly_a<numeric>(current->rest)) {
                        --first_numeric;
                        iter_swap(current,first_numeric);
                } else {
@@ -1445,11 +1437,11 @@ bool expairseq::is_canonical() const
        epvector::const_iterator it_last = it;
        for (++it; it!=itend; it_last=it, ++it) {
                if (!(it_last->is_less(*it) || it_last->is_equal(*it))) {
-                       if (!is_ex_exactly_of_type(it_last->rest,numeric) ||
-                               !is_ex_exactly_of_type(it->rest,numeric)) {
+                       if (!is_exactly_a<numeric>(it_last->rest) ||
+                               !is_exactly_a<numeric>(it->rest)) {
                                // double test makes it easier to set a breakpoint...
-                               if (!is_ex_exactly_of_type(it_last->rest,numeric) ||
-                                       !is_ex_exactly_of_type(it->rest,numeric)) {
+                               if (!is_exactly_a<numeric>(it_last->rest) ||
+                                       !is_exactly_a<numeric>(it->rest)) {
                                        printpair(std::clog, *it_last, 0);
                                        std::clog << ">";
                                        printpair(std::clog, *it, 0);
@@ -1476,7 +1468,7 @@ bool expairseq::is_canonical() const
  *  if no members were changed. */
 epvector * expairseq::expandchildren(unsigned options) const
 {
-       epvector::const_iterator last = seq.end();
+       const epvector::const_iterator last = seq.end();
        epvector::const_iterator cit = seq.begin();
        while (cit!=last) {
                const ex &expanded_ex = cit->rest.expand(options);
@@ -1578,7 +1570,7 @@ epvector * expairseq::subschildren(const lst &ls, const lst &lr, bool no_pattern
        // because the numeric coefficients may be part of the search pattern.
        bool complex_subs = false;
        for (unsigned i=0; i<ls.nops(); ++i)
-               if (is_ex_exactly_of_type(ls.op(i), mul) || is_ex_exactly_of_type(ls.op(i), power)) {
+               if (is_exactly_a<mul>(ls.op(i)) || is_exactly_a<power>(ls.op(i))) {
                        complex_subs = true;
                        break;
                }