frobenius

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

// Algorithms.
#include <myramath/sparse/frobenius.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)
  {
  myra::out() << typestring<Number>() << std::endl;
  // Make random sparse A operand.
  auto A =  SparseMatrix<Number>::random(I,J,N); 
  // Compare frobenius(A) and frobenius(A.make_Matrix())
  typedef typename ReflectPrecision<Number>::type Precision;
  Precision fro_sparse = frobenius(A);
  Precision fro_dense = frobenius(A.make_Matrix());
  Precision error = fro_sparse - fro_dense;
  myra::out() << "  |frobenius(A) [sparse-dense]| = " << error << std::endl;
  REQUIRE(error < tolerance);  
  }

} // namespace

ADD_TEST("frobenius","[sparse]")
  {
  test<NumberS>(15,10,40,1.0e-4f);
  test<NumberD>(15,10,40,1.0e-9);
  test<NumberC>(15,10,40,1.0e-4f);
  test<NumberZ>(15,10,40,1.0e-9);
  }

Results: [PASS]

float
  |frobenius(A) [sparse-dense]| = 0
double
  |frobenius(A) [sparse-dense]| = 0
std::complex<float>
  |frobenius(A) [sparse-dense]| = 0
std::complex<double>
  |frobenius(A) [sparse-dense]| = 0

Go back to Summary of /test programs.