GramMatrix.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: Gramm matrix class
00031 */
00032
00033
00034 #include
00035 #include
00036
00037 namespace RBF{
00038
00039 template
00040 GrammMatrix::GrammMatrix(LA::Matrix &M, LA::Matrix &P)
00041 {
00042 dimension = M.getCol() + P.getCol();
00043 this -> resize(dimension, dimension);
00044
00045 int N = M.getRow();
00046
00047 for(int i = 0; i < N; ++i)
00048 for(int j = 0; j < N; ++j)
00049 this->data[i][j] = M(i, j);
00050
00051 for(int i = 0; i < N; ++i)
00052 for(int j = N; j < dimension; ++j)
00053 this->data[i][j] = P(i, j - N);
00054
00055 LA::Matrix PT = P.transpose();
00056
00057 for(int i = N; i < dimension; ++i)
00058 for(int j = 0; j < N; ++j)
00059 this->data[i][j] = PT(i - N, j);
00060 }
00061
00063
00076 template
00077 GrammMatrix::GrammMatrix(LA::Matrix &WL, LA::Matrix &PL, LA::Matrix &WB, LA::Matrix &PB, LA::Matrix &PT) {
00078
00079 int NI = WL.getRow(), NF = WB.getRow();
00080 int N = NI + NF;
00081 int NC = PL.getCol();
00082
00083 dimension = N + NC;
00084
00085 this->resize(dimension, dimension);
00086
00087 for(int i = 0; i < NI; ++i)
00088 for(int j = 0; j < N; ++j)
00089 this->data[i][j] = WL(i, j);
00090
00091 for(int i = 0; i < NI; ++i)
00092 for(int j = N; j < dimension; ++j)
00093 this->data[i][j] = PL(i, j - N);
00094
00095 for(int i = NI; i < N ; ++i)
00096 for(int j = 0; j < N; ++j)
00097 this->data[i][j] = WB(i - NI, j);
00098
00099 for(int i = NI; i < N ; ++i)
00100 for(int j = N; j < dimension; ++j)
00101 this->data[i][j] = PB(i - NI, j - N);
00102
00103 for(int i = N; i < dimension ; ++i)
00104 for(int j = 0; j < N; ++j)
00105 this->data[i][j] = PT(i-N,j);
00106
00107 for(int i=0; i< NC; ++i)
00108 for(int j = 0; j < NC; ++j)
00109 this->data[N+i][N+j] = 0.0;
00110
00111 }
00112 };