sparse_io

Source: tests/sparse/io.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/sparse/SparseMatrix.h>
#include <myramath/sparse/SparseMatrixRange.h>
#include <myramath/sparse/SparseMatrixBuilder.h>
#include <myramath/sparse/Pattern.h>
#include <myramath/sparse/PatternRange.h>
#include <myramath/sparse/PatternBuilder.h>

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

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

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

using namespace myra;

namespace {

// Tests serialization roundtrip for SparseMatrixRange -> file -> SparseMatrix
template<class Number> void test1()
  {
  myra::out() << typestring<Number>() << std::endl;
  // Make random A.
  auto A = SparseMatrix<Number>::random(20,20,100);
  // Serialize a window of it, read back into B.
  VectorStream inout;
  SparseMatrixRange<Number> A_range = A.window(4,15,6,11);
  A_range.write(inout);
  SparseMatrix<Number> B(inout);
  // Verify A_range == B.
  typename ReflectPrecision<Number>::type error = frobenius(A_range-B.range().add_const());
  myra::out() << "  |A(4:15,6:11)-B| = " << error << std::endl;
  REQUIRE(error == 0);
  }

// Tests serialization roundtrip for SparseMatrixBuilder -> file -> SparseMatrix 
template<class Number> void test2()
  {
  myra::out() << typestring<Number>() << std::endl;
  // Save random SparseMatrixBuilder A.
  auto A = SparseMatrixBuilder<Number>::random(20,20,100);
  VectorStream inout;
  A.write(inout);
  SparseMatrix<Number> B(inout);
  // Verify A == B
  typename ReflectPrecision<Number>::type error = frobenius(A.make_SparseMatrix()-B);
  myra::out() << "  |A-B| = " << error << std::endl;
  REQUIRE(error == 0);
  }


// Tests serialization roundtrip for PatternBuilder -> file -> Pattern 
void test3()
  {
  // Save random PatternBuilder A.
  auto A = PatternBuilder::random(20,20,100);
  VectorStream inout;
  A.write(inout);
  // Verify A == B
  Pattern B(inout);
  REQUIRE(A.make_Pattern() == B);
  }

} // namespace

ADD_TEST("sparse_io","[sparse]")
  {
  // Test SparseMatrixRange/SparseMatrix..
  test1<NumberS>();
  test1<NumberD>();
  test1<NumberC>();
  test1<NumberZ>();
  // Test SparseMatrixBuilder/SparseMatrix..
  test2<NumberS>();
  test2<NumberD>();
  test2<NumberC>();
  test2<NumberZ>();
  // Test PatternBuilder/Pattern
  test3();
  }

Results: [PASS]

float
  |A(4:15,6:11)-B| = 0
double
  |A(4:15,6:11)-B| = 0
std::complex<float>
  |A(4:15,6:11)-B| = 0
std::complex<double>
  |A(4:15,6:11)-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.