diagcat_DiagonalMatrix

Source: tests/dense/diagcat.cpp

// ========================================================================= //
// This file is part of MyraMath, copyright (c) 2014-2019 by Ryan A Chilton  //
// and distributed by MyraCore, LLC. See LICENSE.txt for license terms.      //
// ========================================================================= //

// Containers.
#include <myramath/dense/Matrix.h>
#include <myramath/dense/MatrixRange.h>
#include <myramath/dense/DiagonalMatrix.h>
#include <myramath/dense/DiagonalMatrixRange.h>

// Algorithms.
#include <myramath/dense/diagcat.h>
#include <myramath/dense/frobenius.h>

// Reporting.
#include <tests/myratest.h>

using namespace myra;

namespace {

// Test diagcat(A,B)
template<class Number> void test1(typename ReflectPrecision<Number>::type tolerance)
  {
  typedef typename ReflectPrecision<Number>::type Precision;
  auto A = Matrix<Number>::random(6,6);
  A.bottom(4).left(2).zero();
  A.top(2).right(4).zero();
  auto B = diagcat( A.top(2).left(2), A.bottom(4).right(4) );
  Precision error = frobenius(A-B);
  myra::out() << "  |A - diagcat(a,a)| = " << error << std::endl;
  REQUIRE(error < tolerance);
  }

// Test diagcat(a,b)
template<class Number> void test2(typename ReflectPrecision<Number>::type tolerance)
  {
  typedef typename ReflectPrecision<Number>::type Precision;
  auto A = DiagonalMatrix<Number>::random(6);
  auto B = diagcat( A.first(2), A.last(4) );
  Precision error = frobenius(A-B);
  myra::out() << "  |A - diagcat(...)| = " << error << std::endl;
  REQUIRE(error < tolerance);
  }

} // namespace

ADD_TEST("diagcat_Matrix","[dense]")
  {
  test1<NumberS>(1.0e-5f);
  test1<NumberD>(1.0e-10);
  test1<NumberC>(1.0e-5f);  
  test1<NumberZ>(1.0e-10);
  }

ADD_TEST("diagcat_DiagonalMatrix","[dense]")
  {
  test2<NumberS>(1.0e-5f);
  test2<NumberD>(1.0e-10);
  test2<NumberC>(1.0e-5f);  
  test2<NumberZ>(1.0e-10);
  }

Results: [PASS]

  |A - diagcat(...)| = 0
  |A - diagcat(...)| = 0
  |A - diagcat(...)| = 0
  |A - diagcat(...)| = 0

Go back to Summary of /test programs.