dense_io2

Source: tests/dense/io2.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.      //
// ========================================================================= //

// Container types.
#include <myramath/utility/Number.h>
#include <myramath/dense/Matrix.h>
#include <myramath/dense/MatrixRange.h>
#include <myramath/dense/Vector.h>
#include <myramath/dense/VectorRange.h>
#include <myramath/dense/LowerMatrix.h>
#include <myramath/dense/LowerMatrixRange.h>

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

// Serialization backend.
#include <myramath/io/VectorStream.h>

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

using namespace myra;

namespace {

// Tests serialization roundtrip for MatrixRange -> file -> Matrix
template<class Number> void test1()
  {
  myra::out() << typestring<Number>() << std::endl;
  // Make random A.
  auto A = Matrix<Number>::random(5,5);  
  // Serialize a window of it, read back into B.
  VectorStream inout;
  MatrixRange<Number> A_range = A.bottom(3).right(2);
  A_range.write(inout);
  Matrix<Number> B(inout);
  // Verify A_range == B.
  typename ReflectPrecision<Number>::type error = frobenius(A_range-B);
  myra::out() << "  |A(2:4,3:4)-B| = " << error << std::endl;
  REQUIRE(error == 0);
  }

// Tests serialization roundtrip for VectorRange -> file -> Vector
template<class Number> void test2()
  {
  myra::out() << typestring<Number>() << std::endl;
  // Matrix random a.
  auto a = Vector<Number>::random(5);
  // Serialize a window of it, read back into B.
  VectorStream inout;
  VectorRange<Number> a_range = a.window(1,3);
  a_range.write(inout);
  Vector<Number> b(inout);
  // Verify a_range == b.
  typename ReflectPrecision<Number>::type error = frobenius(a_range-b);
  myra::out() << "  |a(1:3)-b| = " << error << std::endl;
  REQUIRE(error == 0);
  }

// Tests serialization roundtrip for LowerMatrixRange -> file -> LowerMatrix
template<class Number> void test3()
  {
  myra::out() << typestring<Number>() << std::endl;
  // Matrix random A.
  auto A = LowerMatrix<Number>::random(5);
  // Serialize a window of it, read back into B.
  VectorStream inout;
  LowerMatrixRange<Number> A_range = A.range();
  A_range.write(inout);
  LowerMatrix<Number> B(inout);
  // Verify a_range == B.
  typename ReflectPrecision<Number>::type error = frobenius(A_range-B);
  myra::out() << "  |a(:,:)-b| = " << error << std::endl;
  REQUIRE(error == 0);
  }

} // namespace

ADD_TEST("dense_io2","[dense]")
  {
  // Test MatrixRange/Matrix..
  test1<NumberS>();
  test1<NumberD>();
  test1<NumberC>();
  test1<NumberZ>();
  // Test VectorRange/Vector..
  test2<NumberS>();
  test2<NumberD>();
  test2<NumberC>();
  test2<NumberZ>();
  // Test LowerMatrixRange/LowerMatrix..
  test3<NumberS>();
  test3<NumberD>();
  test3<NumberC>();
  test3<NumberZ>();
  }

Results: [PASS]

float
  |A(2:4,3:4)-B| = 0
double
  |A(2:4,3:4)-B| = 0
std::complex<float>
  |A(2:4,3:4)-B| = 0
std::complex<double>
  |A(2:4,3:4)-B| = 0
float
  |a(1:3)-b| = 0
double
  |a(1:3)-b| = 0
std::complex<float>
  |a(1:3)-b| = 0
std::complex<double>
  |a(1:3)-b| = 0
float
  |a(:,:)-b| = 0
double
  |a(:,:)-b| = 0
std::complex<float>
  |a(:,:)-b| = 0
std::complex<double>
  |a(:,:)-b| = 0

Go back to Summary of /test programs.