X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Fmatrix.h;h=6081e900257131db359b0f40fc89e0b9a9e6b19e;hp=ab8d7997a6a93a192f16316ff165968eedc2ce37;hb=8f283de519668b70b2e675a7055c7f1bf7ba197c;hpb=92203097b593322eda950cf8c831177a0c204129 diff --git a/ginac/matrix.h b/ginac/matrix.h index ab8d7997..6081e900 100644 --- a/ginac/matrix.h +++ b/ginac/matrix.h @@ -3,7 +3,7 @@ * Interface to symbolic matrices */ /* - * GiNaC Copyright (C) 1999-2006 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2015 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,16 +20,17 @@ * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA */ -#ifndef __GINAC_MATRIX_H__ -#define __GINAC_MATRIX_H__ +#ifndef GINAC_MATRIX_H +#define GINAC_MATRIX_H -#include -#include #include "basic.h" #include "ex.h" +#include "archive.h" -namespace GiNaC { +#include +#include +namespace GiNaC { /** Helper template to allow initialization of matrices via an overloaded * comma operator (idea stolen from Blitz++). */ @@ -97,36 +98,42 @@ class matrix : public basic // other constructors public: matrix(unsigned r, unsigned c); - matrix(unsigned r, unsigned c, const exvector & m2); matrix(unsigned r, unsigned c, const lst & l); - // First step of initialization of matrix with a comma-separated seqeuence + // First step of initialization of matrix with a comma-separated sequence // of expressions. Subsequent steps are handled by matrix_init<>::operator,(). matrix_init operator=(const ex & x) { m[0] = x; return matrix_init(++m.begin()); } - + +protected: + matrix(unsigned r, unsigned c, const exvector & m2); + matrix(unsigned r, unsigned c, exvector && m2); // functions overriding virtual functions from base classes public: - size_t nops() const; - ex op(size_t i) const; - ex & let_op(size_t i); - ex eval(int level=0) const; - ex evalm() const {return *this;} - ex subs(const exmap & m, unsigned options = 0) const; - ex eval_indexed(const basic & i) const; - ex add_indexed(const ex & self, const ex & other) const; - ex scalar_mul_indexed(const ex & self, const numeric & other) const; - bool contract_with(exvector::iterator self, exvector::iterator other, exvector & v) const; - ex conjugate() const; - ex real_part() const; - ex imag_part() const; - + size_t nops() const override; + ex op(size_t i) const override; + ex & let_op(size_t i) override; + ex eval(int level=0) const override; + ex evalm() const override {return *this;} + ex subs(const exmap & m, unsigned options = 0) const override; + ex eval_indexed(const basic & i) const override; + ex add_indexed(const ex & self, const ex & other) const override; + ex scalar_mul_indexed(const ex & self, const numeric & other) const override; + bool contract_with(exvector::iterator self, exvector::iterator other, exvector & v) const override; + ex conjugate() const override; + ex real_part() const override; + ex imag_part() const override; + + /** Save (a.k.a. serialize) object into archive. */ + void archive(archive_node& n) const override; + /** Read (a.k.a. deserialize) object from archive. */ + void read_archive(const archive_node& n, lst& syms) override; protected: - bool match_same_type(const basic & other) const; - unsigned return_type() const { return return_types::noncommutative; }; + bool match_same_type(const basic & other) const override; + unsigned return_type() const override { return return_types::noncommutative; }; // non-virtual functions in this class public: @@ -170,6 +177,7 @@ protected: unsigned col; ///< number of columns exvector m; ///< representation (cols indexed first) }; +GINAC_DECLARE_UNARCHIVER(matrix); // wrapper functions around member functions @@ -244,4 +252,4 @@ inline ex symbolic_matrix(unsigned r, unsigned c, const std::string & base_name) } // namespace GiNaC -#endif // ndef __GINAC_MATRIX_H__ +#endif // ndef GINAC_MATRIX_H