matrix.hpp

Go to the documentation of this file.

// vim: set tabstop=2:softtabstop=2:shiftwidth=2:expandtab

#ifndef LIN_CORE_TYPES_MATRIX_HPP_
#define LIN_CORE_TYPES_MATRIX_HPP_

#include "../config.hpp"
#include "../traits.hpp"
#include "tensor.hpp"

#include <type_traits>

namespace lin {

template <typename T, size_t R, size_t C, size_t MR = R, size_t MC = C>
class Matrix : public internal::Tensor<Matrix<T, R, C, MR, MC>> {
  static_assert(internal::is_matrix<Matrix<T, R, C, MR, MC>>::value,
      "Invalid Matrix<...> parameters");

 public:
  typedef internal::traits<Matrix<T, R, C, MR, MC>> Traits;

 protected:
  using internal::Tensor<Matrix<T, R, C, MR, MC>>::derived;

 public:
  using internal::Tensor<Matrix<T, R, C, MR, MC>>::Tensor;
  using internal::Tensor<Matrix<T, R, C, MR, MC>>::rows;
  using internal::Tensor<Matrix<T, R, C, MR, MC>>::cols;
  using internal::Tensor<Matrix<T, R, C, MR, MC>>::size;
  using internal::Tensor<Matrix<T, R, C, MR, MC>>::data;
  using internal::Tensor<Matrix<T, R, C, MR, MC>>::eval;
  using internal::Tensor<Matrix<T, R, C, MR, MC>>::resize;
  using internal::Tensor<Matrix<T, R, C, MR, MC>>::operator=;
  using internal::Tensor<Matrix<T, R, C, MR, MC>>::operator();

  constexpr Matrix() = default;
  constexpr Matrix(Matrix<T, R, C, MR, MC> const &) = default;
  constexpr Matrix(Matrix<T, R, C, MR, MC> &&) = default;
  constexpr Matrix<T, R, C, MR, MC> &operator=(Matrix<T, R, C, MR, MC> const &) = default;
  constexpr Matrix<T, R, C, MR, MC> &operator=(Matrix<T, R, C, MR, MC> &&) = default;
};

template <size_t R, size_t C, size_t MR = R, size_t MC = C>
using Matrixf = Matrix<float, R, C, MR, MC>;

typedef Matrixf<2, 2> Matrix2x2f; 
typedef Matrixf<3, 2> Matrix3x2f; 
typedef Matrixf<4, 2> Matrix4x2f; 
typedef Matrixf<2, 3> Matrix2x3f; 
typedef Matrixf<3, 3> Matrix3x3f; 
typedef Matrixf<4, 3> Matrix4x3f; 
typedef Matrixf<2, 4> Matrix2x4f; 
typedef Matrixf<3, 4> Matrix3x4f; 
typedef Matrixf<4, 4> Matrix4x4f; 

template <size_t R, size_t C, size_t MR = R, size_t MC = C>
using Matrixd = Matrix<double, R, C, MR, MC>;

typedef Matrixd<2, 2> Matrix2x2d; 
typedef Matrixd<3, 2> Matrix3x2d; 
typedef Matrixd<4, 2> Matrix4x2d; 
typedef Matrixd<2, 3> Matrix2x3d; 
typedef Matrixd<3, 3> Matrix3x3d; 
typedef Matrixd<4, 3> Matrix4x3d; 
typedef Matrixd<2, 4> Matrix2x4d; 
typedef Matrixd<3, 4> Matrix3x4d; 
typedef Matrixd<4, 4> Matrix4x4d; 

namespace internal {

template <typename T, size_t R, size_t C, size_t MR, size_t MC>
struct _elem<Matrix<T, R, C, MR, MC>> {
  typedef T type;
};

template <typename T, size_t R, size_t C, size_t MR, size_t MC>
struct _dims<Matrix<T, R, C, MR, MC>> {
  static constexpr size_t rows = R;
  static constexpr size_t cols = C;
  static constexpr size_t max_rows = MR;
  static constexpr size_t max_cols = MC;
};

template <class C>
struct _eval<C, std::enable_if_t<is_matrix<C>::value>> {
  typedef Matrix<
      _elem_t<C>,
      _dims<C>::rows,
      _dims<C>::cols,
      _dims<C>::max_rows,
      _dims<C>::max_cols
    > type;
};
}  // namespace internal
}  // namespace lin

#endif