bothcat_Pattern

Source: tests/sparse/bothcat.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/Array2.h>
#include <myramath/dense/intRange.h>
#include <myramath/sparse/SparseMatrix.h>
#include <myramath/sparse/SparseMatrixRange.h>
#include <myramath/sparse/Pattern.h>
#include <myramath/sparse/PatternRange.h>

// Algorithms.
#include <myramath/sparse/bothcat.h>
#include <myramath/sparse/frobenius.h>

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

using namespace myra;

namespace {

// SparseMatrix/SparseMatrixRange

// Test bothcat(A,B,C,D)
template<class Number> void test1(typename ReflectPrecision<Number>::type tolerance)
  {
  typedef typename ReflectPrecision<Number>::type Precision;
  auto A = SparseMatrix<Number>::random(20,20,50);
  auto B = bothcat( A.top(12).left(6), A.top(12).right(14), A.bottom(8).left(6), A.bottom(8).right(14) );
  Precision error = frobenius(A-B);
  myra::out() << "  |A - bothcat(a,a,a,a)| = " << error << std::endl;
  REQUIRE(error < tolerance);
  }

#ifdef MYRAMATH_ENABLE_CPP11

// Test bothcat({A,B,C,D,E,F})
template<class Number> void test2(typename ReflectPrecision<Number>::type tolerance)
  {
  typedef typename ReflectPrecision<Number>::type Precision;
  auto A = SparseMatrix<Number>::random(20,20,50);
  auto B = bothcat( A.add_const().windows({10,5,5},{6,8,6}) );
  Precision error = frobenius(A-B);
  myra::out() << "  |A - bothcat(...)| = " << error << std::endl;
  REQUIRE(error < tolerance);
  }

#endif

// Pattern/PatternRange

// Test bothcat(A,B,C,D)
void test3()
  {
  auto A = Pattern::random(20,20,50);
  auto B = bothcat( A.top(12).left(6), A.top(12).right(14), A.bottom(8).left(6), A.bottom(8).right(14) );
  bool equal = (A==B);
  myra::out() << "  A == bothcat(a,a,a,a) : " << (equal?"true":"false") << std::endl;
  REQUIRE(equal);
  }

#ifdef MYRAMATH_ENABLE_CPP11

// Test bothcat({A,B,C,D,E,F})
void test4()
  {
  auto A = Pattern::random(20,20,50);
  auto B = bothcat( A.add_const().windows({10,5,5},{6,8,6}) );
  bool equal = (A==B);
  myra::out() << "  A == bothcat(...) : " << (equal?"true":"false") << std::endl;
  REQUIRE(equal);
  }

#endif

} // namespace

ADD_TEST("bothcat_SparseMatrix","[sparse]")
  {
  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("bothcat_SparseMatrix_cpp11","[sparse]")
  {
  test2<NumberS>(1.0e-5f);
  test2<NumberD>(1.0e-10);
  test2<NumberC>(1.0e-5f);  
  test2<NumberZ>(1.0e-10);
  }

#endif

ADD_TEST("bothcat_Pattern","[sparse]")
  { test3(); }

#ifdef MYRAMATH_ENABLE_CPP11

ADD_TEST("bothcat_Pattern_cpp11","[sparse]")
  { test4(); }

#endif

Results: [PASS]

  A == bothcat(a,a,a,a) : true

Go back to Summary of /test programs.