GramMatrix.h

Go to the documentation of this file.

/* This file is part of gbLAB.
 *
 * gbLAB is distributed without any warranty under the MIT License.
 */
 
 
#ifndef gbLAB_GramMatrix_h_
#define gbLAB_GramMatrix_h_
 
#include <Eigen/Core>
#include <array>
 
namespace oILAB {
template <typename T, int N> class GramMatrix : public Eigen::Matrix<T, N, N> {
 
  template <int dim>
  static Eigen::Matrix<T, N, N>
  getMatrix(const std::array<Eigen::Matrix<T, dim, 1>, N> &a) {
    Eigen::Matrix<T, dim, N> X(Eigen::Matrix<T, dim, N>::Zero());
    for (int k = 0; k < N; ++k) {
      X.col(k) = a[k];
    }
    return X.transpose() * X;
  }
 
public:
  template <int dim>
  GramMatrix(const std::array<Eigen::Matrix<T, dim, 1>, N> &a)
      : /* init */ Eigen::Matrix<T, N, N>(getMatrix(a)) {}
        
        
    };
    
 
    
//    Rational rat(const long int& n,const long int& d)
//    {
//        return Rational(n,d);
//    }
 
    } // namespace oILAB
 
#endif