convert_sds

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

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

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

using namespace myra;

ADD_TEST("convert_sds","[dense]")
  {
  // Verify that you can convert real float->double->float, without loss.
  int I = 10;
  int J = 12;
  auto S1 = Matrix<NumberS>::random(I,J);
  auto S2 = lower_precision(raise_precision(S1));
  NumberS S_error = frobenius(S1-S2);
  myra::out() << "  |S - float(double(S))| = " << S_error << std::endl;  
  REQUIRE(S_error == 0.0f);
  }

ADD_TEST("convert_czc","[dense]")
  {  
  // Verify that you can convert complex float->double->float, without loss.
  int I = 10;
  int J = 12;
  auto C1 = Matrix<NumberC>::random(I,J);
  auto C2 = lower_precision(raise_precision(C1));
  NumberS C_error = frobenius(C1-C2);
  myra::out() << "  |C - float(double(C))| = " << C_error << std::endl;  
  REQUIRE(C_error == 0.0f);
  }

ADD_TEST("convert_dsd","[dense]")
  {
  // Converting a real double->float->double implies some loss.
  int I = 10;
  int J = 12;
  auto D1 = Matrix<NumberD>::random(I,J);
  auto D2 = raise_precision(lower_precision(D1));
  NumberD D_error = frobenius(D1-D2);
  myra::out() << "  |D - double(float(D))| = " << D_error << std::endl;
  REQUIRE(D_error < 1.0e-6);
  }

ADD_TEST("convert_zcz","[dense]")
  {  
  // Converting a complex double->float->double implies some loss.
  int I = 10;
  int J = 12;
  auto Z1 = Matrix<NumberZ>::random(I,J);
  auto Z2 = raise_precision(lower_precision(Z1));
  NumberD Z_error = frobenius(Z1-Z2);
  myra::out() << "  |Z - double(float(Z))| = " << Z_error << std::endl;
  REQUIRE(Z_error < 1.0e-6);
  }

Results: [PASS]

  |S - float(double(S))| = 0

Go back to Summary of /test programs.