num_dimensions.cpp File Reference

#include "main.h"
#include <Eigen/SparseCore>

Typedefs
template<typename Scalar , int Rows, int Cols>
using	TArray = Array< Scalar, Rows, Cols >
template<typename Scalar , int Rows, int Cols>
using	TMatrix = Matrix< Scalar, Rows, Cols >

Functions
template<int ExpectedDim, typename Xpr >
void	check_dim (const Xpr &)
template<template< typename, int, int > class Object>
void	map_num_dimensions ()
	EIGEN_DECLARE_TEST (num_dimensions)

Typedef Documentation

TArray

template<typename Scalar , int Rows, int Cols>

using TArray = Array<Scalar, Rows, Cols>

TMatrix

template<typename Scalar , int Rows, int Cols>

using TMatrix = Matrix<Scalar, Rows, Cols>

Function Documentation

check_dim()

template<int ExpectedDim, typename Xpr >

void check_dim

(

const Xpr &

)

                           {
  STATIC_CHECK(Xpr::NumDimensions == ExpectedDim);
}

References STATIC_CHECK.

EIGEN_DECLARE_TEST()

EIGEN_DECLARE_TEST

(

num_dimensions

)

                                   {
  int n = 10;
  ArrayXXd A(n, n);
  CALL_SUBTEST(check_dim<2>(A));
  CALL_SUBTEST(check_dim<2>(A.block(1, 1, 2, 2)));
  CALL_SUBTEST(check_dim<1>(A.col(1)));
  CALL_SUBTEST(check_dim<1>(A.row(1)));
 
  MatrixXd M(n, n);
  CALL_SUBTEST(check_dim<0>(M.row(1) * M.col(1)));
 
  SparseMatrix<double> S(n, n);
  CALL_SUBTEST(check_dim<2>(S));
  CALL_SUBTEST(check_dim<2>(S.block(1, 1, 2, 2)));
  CALL_SUBTEST(check_dim<1>(S.col(1)));
  CALL_SUBTEST(check_dim<1>(S.row(1)));
 
  SparseVector<double> s(n);
  CALL_SUBTEST(check_dim<1>(s));
  CALL_SUBTEST(check_dim<1>(s.head(2)));
 
  CALL_SUBTEST(map_num_dimensions<TArray>());
  CALL_SUBTEST(map_num_dimensions<TMatrix>());
}

References CALL_SUBTEST, n, s, and oomph::QuadTreeNames::S.

map_num_dimensions()

template<template< typename, int, int > class Object>

void map_num_dimensions

(

)

                          {
  typedef Object<double, 1, 1> ArrayScalarType;
  typedef Object<double, 2, 1> ArrayVectorType;
  typedef Object<double, 1, 2> TransposeArrayVectorType;
  typedef Object<double, 2, 2> ArrayType;
  typedef Object<double, Eigen::Dynamic, 1> DynamicArrayVectorType;
  typedef Object<double, 1, Eigen::Dynamic> DynamicTransposeArrayVectorType;
  typedef Object<double, Eigen::Dynamic, Eigen::Dynamic> DynamicArrayType;
 
  STATIC_CHECK(ArrayScalarType::NumDimensions == 0);
  STATIC_CHECK(ArrayVectorType::NumDimensions == 1);
  STATIC_CHECK(TransposeArrayVectorType::NumDimensions == 1);
  STATIC_CHECK(ArrayType::NumDimensions == 2);
  STATIC_CHECK(DynamicArrayVectorType::NumDimensions == 1);
  STATIC_CHECK(DynamicTransposeArrayVectorType::NumDimensions == 1);
  STATIC_CHECK(DynamicArrayType::NumDimensions == 2);
 
  typedef Eigen::Map<ArrayScalarType> ArrayScalarMap;
  typedef Eigen::Map<ArrayVectorType> ArrayVectorMap;
  typedef Eigen::Map<TransposeArrayVectorType> TransposeArrayVectorMap;
  typedef Eigen::Map<ArrayType> ArrayMap;
  typedef Eigen::Map<DynamicArrayVectorType> DynamicArrayVectorMap;
  typedef Eigen::Map<DynamicTransposeArrayVectorType> DynamicTransposeArrayVectorMap;
  typedef Eigen::Map<DynamicArrayType> DynamicArrayMap;
 
  STATIC_CHECK(ArrayScalarMap::NumDimensions == 0);
  STATIC_CHECK(ArrayVectorMap::NumDimensions == 1);
  STATIC_CHECK(TransposeArrayVectorMap::NumDimensions == 1);
  STATIC_CHECK(ArrayMap::NumDimensions == 2);
  STATIC_CHECK(DynamicArrayVectorMap::NumDimensions == 1);
  STATIC_CHECK(DynamicTransposeArrayVectorMap::NumDimensions == 1);
  STATIC_CHECK(DynamicArrayMap::NumDimensions == 2);
}

References STATIC_CHECK.