expr_tensor

Source: tests/expression/arithmetic.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.      //
// ========================================================================= //

#include <myramath/expression/Expression.h>
#include <myramath/expression/expr.h>
#include <myramath/expression/arithmetic.h>
#include <myramath/expression/frobenius.h>

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

using namespace myra;

#ifdef MYRAMATH_ENABLE_CPP11

ADD_TEST("expr_arithmetic","[expression]")
  {
  double tol = 1.0e-12;
  auto E = expr({1.0,2.0,3.0});
  // operator+
  REQUIRE( frobenius( (E+E) - expr({2.0,4.0,6.0}) ) < tol );
  REQUIRE( frobenius( (E+2.0) - expr({3.0,4.0,5.0}) ) < tol );
  REQUIRE( frobenius( (3.0+E) - expr({4.0,5.0,6.0}) ) < tol );
  // operator-
  REQUIRE( frobenius( (E-E) - expr({0.0,0.0,0.0}) ) < tol );
  REQUIRE( frobenius( (E-2.0) - expr({-1.0,0.0,1.0}) ) < tol );
  REQUIRE( frobenius( (3.0-E) - expr({2.0,1.0,0.0}) ) < tol );
  // operator*
  REQUIRE( frobenius( (E*E) - expr({1.0,4.0,9.0}) ) < tol );
  REQUIRE( frobenius( (E*2.0) - expr({2.0,4.0,6.0}) ) < tol );
  REQUIRE( frobenius( (3.0*E) - expr({3.0,6.0,9.0}) ) < tol );
  // operator/
  REQUIRE( frobenius( (E/E) - expr({1.0,1.0,1.0}) ) < tol );
  REQUIRE( frobenius( (E/2.0) - expr({0.5,1.0,1.5}) ) < tol );
  REQUIRE( frobenius( (3.0/E) - expr({3.0,1.5,1.0}) ) < tol );
  // unary-
  REQUIRE( frobenius( (-E) - expr({-1.0,-2.0,-3.0}) ) < tol );
  }

ADD_TEST("expr_tensor","[expression]")
  {
  double tol = 1.0e-12;
  std::vector<double> v1({1.0,2.0,3.0});
  std::vector<double> v2({4.0,5.0,6.0});
  auto E = expr(v1) ^ expr(v2);
  for (int i = 0; i < 3; ++i)
    for (int j = 0; j < 3; ++j)
      REQUIRE( std::abs(E.evaluate({i,j})-v1[i]*v2[j]) < tol );
  }

#endif

Results: [PASS]

Go back to Summary of /test programs.