00001 // $Id: NOX_LAPACK_Matrix.H,v 1.10 2006/08/22 00:01:31 rppawlo Exp $ 00002 // $Source: /space/CVS/Trilinos/packages/nox/src-lapack/NOX_LAPACK_Matrix.H,v $ 00003 00004 //@HEADER 00005 // ************************************************************************ 00006 // 00007 // NOX: An Object-Oriented Nonlinear Solver Package 00008 // Copyright (2002) Sandia Corporation 00009 // 00010 // LOCA: Library of Continuation Algorithms Package 00011 // Copyright (2005) Sandia Corporation 00012 // 00013 // Under terms of Contract DE-AC04-94AL85000, there is a non-exclusive 00014 // license for use of this work by or on behalf of the U.S. Government. 00015 // 00016 // This library is free software; you can redistribute it and/or modify 00017 // it under the terms of the GNU Lesser General Public License as 00018 // published by the Free Software Foundation; either version 2.1 of the 00019 // License, or (at your option) any later version. 00020 // 00021 // This library is distributed in the hope that it will be useful, but 00022 // WITHOUT ANY WARRANTY; without even the implied warranty of 00023 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00024 // Lesser General Public License for more details. 00025 // 00026 // You should have received a copy of the GNU Lesser General Public 00027 // License along with this library; if not, write to the Free Software 00028 // Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 00029 // USA 00030 // 00031 // Questions? Contact Roger Pawlowski (rppawlo@sandia.gov) or 00032 // Eric Phipps (etphipp@sandia.gov), Sandia National Laboratories. 00033 // ************************************************************************ 00034 // CVS Information 00035 // $Source: /space/CVS/Trilinos/packages/nox/src-lapack/NOX_LAPACK_Matrix.H,v $ 00036 // $Author: rppawlo $ 00037 // $Date: 2006/08/22 00:01:31 $ 00038 // $Revision: 1.10 $ 00039 // ************************************************************************ 00040 //@HEADER 00041 00042 #ifndef NOX_LAPACK_MATRIX_H 00043 #define NOX_LAPACK_MATRIX_H 00044 00045 #include "NOX_Common.H" 00046 00047 namespace NOX { 00048 00049 namespace LAPACK { 00050 00052 00057 template <typename T> 00058 class Matrix { 00059 00060 public: 00061 00063 Matrix() : p(0), q(0), entries() {} 00064 00066 Matrix(int m, int n) : p(m), q(n), entries(m*n) {} 00067 00069 Matrix(const Matrix& a) : 00070 p(a.p), q(a.q), entries(a.entries) {} 00071 00073 ~Matrix() {} 00074 00076 T& operator()(int i, int j) { return entries[i + (p*j)]; } 00077 00079 const T& operator()(int i, int j) const { return entries[i + (p*j)]; } 00080 00082 00085 void scale(T v ) { 00086 for (int i=0; i<p*q; i++) 00087 entries[i] *= v; 00088 } 00089 00091 bool print(std::ostream& stream) const { 00092 for (int i=0; i<p; i++) { 00093 stream << "[ "; 00094 for (int j=0; j<q; j++) 00095 stream << operator()(i,j) << " "; 00096 stream << "]" << std::endl; 00097 } 00098 return stream; 00099 } 00100 00102 int numRows() const { return p; } 00103 00105 int numCols() const { return q; } 00106 00107 private: 00108 00110 int p, q; 00111 00113 std::vector<T> entries; 00114 00115 }; 00116 00117 } // namespace LAPACK 00118 00119 } // namespace NOX 00120 00121 template <typename T> 00122 ostream& operator<<(ostream& stream, const NOX::LAPACK::Matrix<T>& v) { 00123 v.print(stream); 00124 return stream; 00125 } 00126 00127 00128 #endif