X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Foperators.cpp;h=eb91c7c19b8cfd83dadb3d5064a47b82b594bde0;hp=345d8f59eb1b552e22a41c37431762770962dbf5;hb=3d09388a8ea144a19949c216fae43ccba92e47aa;hpb=c8feefe95a6c219195aea22050f17e2294656f32 diff --git a/ginac/operators.cpp b/ginac/operators.cpp index 345d8f59..eb91c7c1 100644 --- a/ginac/operators.cpp +++ b/ginac/operators.cpp @@ -3,7 +3,7 @@ * Implementation of GiNaC's overloaded operators. */ /* - * GiNaC Copyright (C) 1999-2003 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2010 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 @@ -17,12 +17,9 @@ * * You should have received a copy of the GNU General Public License * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA */ -#include -#include - #include "operators.h" #include "numeric.h" #include "add.h" @@ -33,6 +30,9 @@ #include "print.h" #include "utils.h" +#include +#include + namespace GiNaC { /** Used internally by operator+() to add two ex objects together. */ @@ -48,10 +48,11 @@ static inline const ex exmul(const ex & lh, const ex & rh) // ncmul::eval()'s rule to pull out commutative elements we need to check // only one of the elements. if (rh.return_type()==return_types::commutative || - lh.return_type()==return_types::commutative) + lh.return_type()==return_types::commutative) { return (new mul(lh,rh))->setflag(status_flags::dynallocated); - else + } else { return (new ncmul(lh,rh))->setflag(status_flags::dynallocated); + } } /** Used internally by operator-() and friends to change the sign of an argument. */ @@ -175,7 +176,7 @@ const numeric operator+(const numeric & lh) const numeric operator-(const numeric & lh) { - return _num_1.mul(lh); + return _num_1_p->mul(lh); } @@ -214,14 +215,14 @@ const ex operator--(ex & lh, int) /** Numeric prefix increment. Adds 1 and returns incremented number. */ numeric& operator++(numeric & rh) { - rh = rh.add(_num1); + rh = rh.add(*_num1_p); return rh; } /** Numeric prefix decrement. Subtracts 1 and returns decremented number. */ numeric& operator--(numeric & rh) { - rh = rh.add(_num_1); + rh = rh.add(*_num_1_p); return rh; } @@ -230,7 +231,7 @@ numeric& operator--(numeric & rh) const numeric operator++(numeric & lh, int) { numeric tmp(lh); - lh = lh.add(_num1); + lh = lh.add(*_num1_p); return tmp; } @@ -239,7 +240,7 @@ const numeric operator++(numeric & lh, int) const numeric operator--(numeric & lh, int) { numeric tmp(lh); - lh = lh.add(_num_1); + lh = lh.add(*_num_1_p); return tmp; } @@ -349,6 +350,92 @@ std::ostream & operator<<(std::ostream & os, const ex & e) return os; } +std::ostream & operator<<(std::ostream & os, const exvector & e) +{ + print_context *p = get_print_context(os); + exvector::const_iterator i = e.begin(); + exvector::const_iterator vend = e.end(); + + if (i==vend) { + os << "[]"; + return os; + } + + os << "["; + while (true) { + if (p == 0) + i -> print(print_dflt(os)); + else + i -> print(*p); + ++i; + if (i==vend) + break; + os << ","; + } + os << "]"; + + return os; +} + +std::ostream & operator<<(std::ostream & os, const exset & e) +{ + print_context *p = get_print_context(os); + exset::const_iterator i = e.begin(); + exset::const_iterator send = e.end(); + + if (i==send) { + os << "<>"; + return os; + } + + os << "<"; + while (true) { + if (p == 0) + i->print(print_dflt(os)); + else + i->print(*p); + ++i; + if (i == send) + break; + os << ","; + } + os << ">"; + + return os; +} + +std::ostream & operator<<(std::ostream & os, const exmap & e) +{ + print_context *p = get_print_context(os); + exmap::const_iterator i = e.begin(); + exmap::const_iterator mend = e.end(); + + if (i==mend) { + os << "{}"; + return os; + } + + os << "{"; + while (true) { + if (p == 0) + i->first.print(print_dflt(os)); + else + i->first.print(*p); + os << "=="; + if (p == 0) + i->second.print(print_dflt(os)); + else + i->second.print(*p); + ++i; + if( i==mend ) + break; + os << ","; + } + os << "}"; + + return os; +} + std::istream & operator>>(std::istream & is, ex & e) { throw (std::logic_error("expression input from streams not implemented"));