diagcat_Action_cpp11

Source: tests/iterative/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/utility/Number.h>
#include <myramath/utility/Array1.h>
#include <myramath/dense/Matrix.h>
#include <myramath/dense/MatrixRange.h>

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

// Iterative solve/action stuff.
#include <myramath/iterative/Action.h>
#include <myramath/iterative/GemmAction.h>
#include <myramath/iterative/diagcat.h>

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

using namespace myra;

namespace {

// Tests diagcat(A,B)
template<class Number> void test1(typename ReflectPrecision<Number>::type tolerance)
  {
  typedef typename ReflectPrecision<Number>::type Precision;
  auto A0 = Matrix<Number>::random(5,3);
  auto B0 = Matrix<Number>::random(5,7);
  auto C0 = diagcat(A0,B0);
  auto A = make_GemmAction(A0);
  auto B = make_GemmAction(B0);
  auto C = diagcat(A,B);
  Precision error = frobenius(C0-C.make_Matrix());
  myra::out() << "  |C0-C| = " << error << std::endl;
  REQUIRE(error < tolerance);
  }

#ifdef MYRAMATH_ENABLE_CPP11

// Tests diagcat(A,B,C)
template<class Number> void test2(typename ReflectPrecision<Number>::type tolerance)
  {
  typedef typename ReflectPrecision<Number>::type Precision;
  auto A0 = Matrix<Number>::random(5,3);
  auto B0 = Matrix<Number>::random(5,7);
  auto C0 = Matrix<Number>::random(5,6);
  auto D0 = diagcat({A0,B0,C0});
  auto A = make_GemmAction(A0);
  auto B = make_GemmAction(B0);
  auto C = make_GemmAction(C0);
  auto D = diagcat({A,B,C});
  Precision error = frobenius(D0-D.make_Matrix());
  myra::out() << "  |D0-D| = " << error << std::endl;
  REQUIRE(error < tolerance);
  }

#endif

}

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

#ifdef MYRAMATH_ENABLE_CPP11

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

#endif

Results: [PASS]

  |D0-D| = 0
  |D0-D| = 0
  |D0-D| = 0
  |D0-D| = 0

Go back to Summary of /test programs.