Matrix.cpp
00001 /******************************************************************************
00002 
00003  Copyright 2008 Departamento de Realidad Virtual
00004  y Unidad de Cómputo Aplicado DGSGA, UNAM.
00005 
00006 
00007  This file is part of RBF++.
00008 
00009  RBF++ is free software: you can redistribute it and/or modify
00010  it under the terms of the GNU General Public License as published by
00011  the Free Software Foundation, either version 3 of the License, or
00012  (at your option) any later version.
00013 
00014  RBF++ is distributed in the hope that it will be useful,
00015  but WITHOUT ANY WARRANTY; without even the implied warranty of
00016  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
00017  GNU General Public License for more details.
00018 
00019  You should have received a copy of the GNU General Public License
00020  along with RBF++. If not, see .
00021 
00022 
00023 *******************************************************************************/
00024 
00025 /* 
00026  Author: Daniel Cervantes Cabrera.
00027  Project: RBF++
00028  Institution: DGSCA, UNAM.
00029  Date: 3/05/08
00030  Description: Template matrix class.
00031  */
00032 
00033 
00034 #include 
00035 #include 
00036 #include 
00037 
00038 
00039 
00040 
00041 namespace LA{
00042 
00043 #define SWAP(a,b) {temp=(a);(a)=(b);(b)=temp;}
00044 
00045 template
00046 Matrix::Matrix() : data(0), row(0), col(0) {
00047 }
00048 
00049 template
00050 Matrix::Matrix(int r, int c) :row(r), col(c) {
00051   
00052   data = new T* [row];
00053   for(int i = 0; i < row; ++i) 
00054     data[i] = new T [col];
00055 }
00056 
00057 template
00058 Matrix::Matrix(const Matrix& otra) { 
00059 
00060   if(this != &otra) {
00061     row = otra.row;
00062     col = otra.col;
00063     data = new T* [row];
00064     for(int i = 0; i < row; ++i)
00065       data[i] = new T [col];
00066     
00067     for(int i = 0; i < row; ++i)
00068       for(int j = 0; j < col; ++j)
00069         data[i][j] = otra.data[i][j];
00070   }
00071 }
00072 
00073 template
00074 Matrix::~Matrix() { 
00075   for(int i = 0; i < row; ++i) 
00076     delete [] data[i];
00077   delete [] data;
00078 }
00079 
00080 template
00081 void Matrix::resize(int r, int c) {
00082 
00083   if( row != 0 && col != 0 ) {
00084     for(int i = 0; i < row; ++i) 
00085       delete [] data[i];
00086     delete [] data;
00087   }
00088 
00089   row = r;
00090   col = c;
00091 
00092   data = new T* [row];
00093   for(int i = 0; i < row; ++i) 
00094     data[i] = new T [col];
00095 }
00096 
00097 template
00098 T& Matrix::operator()(int i, int j) const { 
00099   return data[i][j]; 
00100 }
00101 
00102 template
00103 Matrix Matrix::transpose() {
00104  
00105 
00106   Matrix TM(col, row);
00107   for(int i = 0; i < col; ++i) 
00108     for(int j = 0; j < row; ++j)
00109       TM(i,j) = data[j][i];
00110 
00111   return TM;
00112 }
00113 
00114 //Gauss con pivoteo.
00115 template
00116 vector Matrix::operator/(vector &b) {
00117   
00118   Matrix a;
00119   vector x;
00120   int n = row;
00121   a.resize(n,n+1);
00122   x.resize(n);
00123 
00124   T  pivot, temp,q;  
00125   int ipivot,ip1;
00126 
00127   for(int i=0; i
00182 ostream& operator<<(ostream& s, const Matrix& m) {
00183   for(int i = 0; i