/*************************************************************************** * blitz/mattoep.h Declarations for Toeplitz matrices * * $Id: mattoep.h,v 1.1.1.1 1999-04-09 17:59:00 ansari Exp $ * * Copyright (C) 1997,1998 Todd Veldhuizen * * 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 the Free Software Foundation; either version 2 * of the License, or (at your option) any later version. * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * GNU General Public License for more details. * * Suggestions: blitz-suggest@cybervision.com * Bugs: blitz-bugs@cybervision.com * * For more information, please see the Blitz++ Home Page: * http://seurat.uwaterloo.ca/blitz/ * *************************************************************************** * $Log: not supported by cvs2svn $ * Revision 1.3 1998/03/14 00:04:47 tveldhui * 0.2-alpha-05 * * Revision 1.2 1997/01/24 14:42:00 tveldhui * Periodic RCS update * * Revision 1.1 1997/01/13 22:19:58 tveldhui * Periodic RCS update * * */ #ifndef BZ_MATTOEP_H #define BZ_MATTOEP_H #ifndef BZ_MSTRUCT_H #error must be included via #endif BZ_NAMESPACE(blitz) // Toeplitz matrix // [ 0 1 2 3 ] // [ 1 2 3 4 ] // [ 2 3 4 5 ] // [ 3 4 5 6 ] class ToeplitzIterator { public: ToeplitzIterator(unsigned rows, unsigned cols) { rows_ = rows; cols_ = cols; i_ = 0; j_ = 0; good_ = true; offset_ = 0; } operator _bz_bool() const { return good_; } void operator++() { ++offset_; if (i_ < rows_ - 1) ++i_; else if (j_ < cols_ - 1) ++j_; else good_ = false; } unsigned row() const { return i_; } unsigned col() const { return j_; } unsigned offset() const { return offset_; } protected: unsigned offset_; unsigned i_, j_; unsigned rows_, cols_; _bz_bool good_; }; class Toeplitz : public GeneralMatrix { public: typedef ToeplitzIterator T_iterator; Toeplitz() : rows_(0), cols_(0) { } Toeplitz(unsigned rows, unsigned cols) : rows_(rows), cols_(cols) { } unsigned columns() const { return cols_; } unsigned coordToOffset(unsigned i, unsigned j) const { BZPRECONDITION(inRange(i,j)); return i + j; } unsigned firstInRow(unsigned i) const { return 0; } template T_numtype get(const T_numtype * _bz_restrict data, unsigned i, unsigned j) const { BZPRECONDITION(inRange(i,j)); return data[coordToOffset(i,j)]; } template T_numtype& get(T_numtype * _bz_restrict data, unsigned i, unsigned j) { BZPRECONDITION(inRange(i,j)); return data[coordToOffset(i,j)]; } unsigned lastInRow(unsigned i) const { return cols_ - 1; } unsigned firstInCol(unsigned j) const { return 0; } unsigned lastInCol(unsigned j) const { return rows_ - 1; } _bz_bool inRange(unsigned i, unsigned j) const { return (i < rows_) && (j < cols_); } unsigned numElements() const { return rows_ + cols_ - 1; } unsigned rows() const { return rows_; } void resize(unsigned rows, unsigned cols) { rows_ = rows; cols_ = cols; } private: unsigned rows_, cols_; }; BZ_NAMESPACE_END #endif // BZ_MATSYMM_H