/* * $Id: matrix.cc,v 1.1.1.1 1999-04-09 17:59:00 ansari Exp $ * * Copyright (C) 1997 Todd Veldhuizen * All rights reserved. Please see for terms and * conditions of use. * * $Log: not supported by cvs2svn $ * Revision 1.4 1998/03/14 00:04:47 tveldhui * 0.2-alpha-05 * * Revision 1.3 1997/07/16 14:51:20 tveldhui * Update: Alpha release 0.2 (Arrays) * * Revision 1.2 1997/01/24 14:42:00 tveldhui * Periodic RCS update * */ #ifndef BZ_MATRIX_CC #define BZ_MATRIX_CC #ifndef BZ_MATRIX_H #include #endif BZ_NAMESPACE(blitz) // Matrix expression operand template template Matrix& Matrix::operator=(_bz_MatExpr expr) { // Check for compatible structures. // Fast evaluation (compatible structures) // (not implemented) // Slow evaluation _bz_typename P_structure::T_iterator iter(rows(), cols()); while (iter) { data_[iter.offset()] = expr(iter.row(), iter.col()); ++iter; } return *this; } template ostream& operator<<(ostream& os, const Matrix& matrix) { os << "[ "; for (int i=0; i < matrix.rows(); ++i) { for (int j=0; j < matrix.columns(); ++j) { os << setw(10) << matrix(i,j); if ((!((j+1)%7)) && (j < matrix.cols()-1)) os << endl << " ..."; } if (i != matrix.rows() - 1) os << endl << " "; } os << " ]"; return os; } BZ_NAMESPACE_END #endif // BZ_MATRIX_CC