X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Fmatrix.cpp;h=37c4927786299f1ad84453897258c6b0401fe6c3;hp=d108bd4c04be6bb315a3966c3f4a716d77a45504;hb=3627c2c0f2fe8fe04816703c5fc6dad9c96fb0a6;hpb=a377cee53b71348235ec36f83afeced7e10288a8 diff --git a/ginac/matrix.cpp b/ginac/matrix.cpp index d108bd4c..37c49277 100644 --- a/ginac/matrix.cpp +++ b/ginac/matrix.cpp @@ -53,7 +53,7 @@ GINAC_IMPLEMENT_REGISTERED_CLASS_OPT(matrix, basic, ////////// /** Default ctor. Initializes to 1 x 1-dimensional zero-matrix. */ -matrix::matrix() : inherited(&matrix::tinfo_static), row(1), col(1), m(1, _ex0) +matrix::matrix() : row(1), col(1), m(1, _ex0) { setflag(status_flags::not_shareable); } @@ -68,8 +68,7 @@ matrix::matrix() : inherited(&matrix::tinfo_static), row(1), col(1), m(1, _ex0) * * @param r number of rows * @param c number of cols */ -matrix::matrix(unsigned r, unsigned c) - : inherited(&matrix::tinfo_static), row(r), col(c), m(r*c, _ex0) +matrix::matrix(unsigned r, unsigned c) : row(r), col(c), m(r*c, _ex0) { setflag(status_flags::not_shareable); } @@ -78,7 +77,7 @@ matrix::matrix(unsigned r, unsigned c) /** Ctor from representation, for internal use only. */ matrix::matrix(unsigned r, unsigned c, const exvector & m2) - : inherited(&matrix::tinfo_static), row(r), col(c), m(m2) + : row(r), col(c), m(m2) { setflag(status_flags::not_shareable); } @@ -88,7 +87,7 @@ matrix::matrix(unsigned r, unsigned c, const exvector & m2) * If the list has more elements than the matrix, the excessive elements are * thrown away. */ matrix::matrix(unsigned r, unsigned c, const lst & l) - : inherited(&matrix::tinfo_static), row(r), col(c), m(r*c, _ex0) + : row(r), col(c), m(r*c, _ex0) { setflag(status_flags::not_shareable); @@ -106,22 +105,25 @@ matrix::matrix(unsigned r, unsigned c, const lst & l) // archiving ////////// -matrix::matrix(const archive_node &n, lst &sym_lst) : inherited(n, sym_lst) +void matrix::read_archive(const archive_node &n, lst &sym_lst) { - setflag(status_flags::not_shareable); + inherited::read_archive(n, sym_lst); if (!(n.find_unsigned("row", row)) || !(n.find_unsigned("col", col))) throw (std::runtime_error("unknown matrix dimensions in archive")); m.reserve(row * col); + // XXX: default ctor inserts a zero element, we need to erase it here. + m.pop_back(); archive_node::archive_node_cit first = n.find_first("m"); archive_node::archive_node_cit last = n.find_last("m"); ++last; - for (archive_node::archive_node_cit i=first; i