00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034 #include <RBF/GrammMatrix.h>
00035 #include <LA/Matrix.h>
00036
00037 namespace RBF{
00038
00039 template<typename T>
00040 GrammMatrix<T>::GrammMatrix(LA::Matrix<T> &M, LA::Matrix<T> &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<T> 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<typename T>
00077 GrammMatrix<T>::GrammMatrix(LA::Matrix<T> &WL, LA::Matrix<T> &PL, LA::Matrix<T> &WB, LA::Matrix<T> &PB, LA::Matrix<T> &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 };