SparseMatrixRange

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

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

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

using namespace myra;

namespace {

template<class Number> void test(int I, int J, int N, typename ReflectPrecision<Number>::type tolerance)
  {
  typedef typename ReflectPrecision<Number>::type Precision;
  // Generate random A (sparse) and B (dense).
  auto A = SparseMatrix<Number>::random(I,J,N);
  auto B = Matrix<Number>::random(J,10);
  // Multiply only using the lower right corner, DenseMatrix (reference)
  //  [ ~  ~ ] * [ ~  ]
  //  [ ~ A22]   [ B2 ]
  Matrix<Number> C1 = gemm(A.make_Matrix().window(I/2,I,J/2,J), 'N', B.rows(J/2,J), 'N');
  // Repeat the multiply with a CSparseMatrixRange.
  Matrix<Number> C2 = gemm(A.window(I/2,I,J/2,J), 'N', B.rows(J/2,J), 'N');
  // Repeat the multiply with a SparseMatrix, constructed from the range.
  SparseMatrix<Number> A2( A.window(I/2,I,J/2,J) );  
  Matrix<Number> C3 = gemm(A2, 'N', B.rows(J/2,J), 'N');
  // Check error C2 - C1.
  Precision error2 = frobenius(C2-C1);
  myra::out() << "  |A.window()*B [sparse-dense]| = " << error2 << std::endl;
  REQUIRE(error2 < tolerance);
  // Check error C3 - C1.
  Precision error3 = frobenius(C3-C1);  
  myra::out() << "  |A*B.window() [sparse-dense]| = " << error3 << std::endl;
  REQUIRE(error3 < tolerance);    
  }

} // namespace

ADD_TEST("SparseMatrixRange","[sparse]")
  {
  myra::out() << typestring<NumberS>() << std::endl;
  test<NumberS> (100,100,1000,1.0e-4f);
  test<NumberS > (10,1,3,1.0e-4f);  
  test<NumberS > (1,10,3,1.0e-4f);  
  test<NumberS > (0,0,0,1.0e-4f);
  myra::out() << typestring<NumberD>() << std::endl;
  test<NumberD> (100,100,1000,1.0e-9);
  test<NumberD> (10,1,3,1.0e-9);  
  test<NumberD> (1,10,3,1.0e-9);  
  test<NumberD> (0,0,0,1.0e-9);
  myra::out() << typestring<NumberC>() << std::endl;
  test<NumberC> (100,100,1000,1.0e-4f);
  test<NumberC> (10,1,3,1.0e-4f);
  test<NumberC> (1,10,3,1.0e-4f);  
  test<NumberC> (0,0,0,1.0e-4f);
  myra::out() << typestring<NumberZ>() << std::endl;
  test<NumberZ> (100,100,1000,1.0e-9);
  test<NumberZ> (10,1,3,1.0e-9);  
  test<NumberZ> (1,10,3,1.0e-9);  
  test<NumberZ> (0,0,0,1.0e-9);
  }

Results: [PASS]

float
  |A.window()*B [sparse-dense]| = 8.40485e-07
  |A*B.window() [sparse-dense]| = 8.40485e-07
  |A.window()*B [sparse-dense]| = 0
  |A*B.window() [sparse-dense]| = 0
  |A.window()*B [sparse-dense]| = 0
  |A*B.window() [sparse-dense]| = 0
  |A.window()*B [sparse-dense]| = 0
  |A*B.window() [sparse-dense]| = 0
double
  |A.window()*B [sparse-dense]| = 1.86056e-15
  |A*B.window() [sparse-dense]| = 1.86056e-15
  |A.window()*B [sparse-dense]| = 0
  |A*B.window() [sparse-dense]| = 0
  |A.window()*B [sparse-dense]| = 3.46945e-17
  |A*B.window() [sparse-dense]| = 3.46945e-17
  |A.window()*B [sparse-dense]| = 0
  |A*B.window() [sparse-dense]| = 0
std::complex<float>
  |A.window()*B [sparse-dense]| = 2.57639e-06
  |A*B.window() [sparse-dense]| = 2.57639e-06
  |A.window()*B [sparse-dense]| = 0
  |A*B.window() [sparse-dense]| = 0
  |A.window()*B [sparse-dense]| = 0
  |A*B.window() [sparse-dense]| = 0
  |A.window()*B [sparse-dense]| = 0
  |A*B.window() [sparse-dense]| = 0
std::complex<double>
  |A.window()*B [sparse-dense]| = 4.8191e-15
  |A*B.window() [sparse-dense]| = 4.8191e-15
  |A.window()*B [sparse-dense]| = 0
  |A*B.window() [sparse-dense]| = 0
  |A.window()*B [sparse-dense]| = 2.77556e-16
  |A*B.window() [sparse-dense]| = 2.77556e-16
  |A.window()*B [sparse-dense]| = 0
  |A*B.window() [sparse-dense]| = 0

Go back to Summary of /test programs.