doc/html/test_ProductAction.html

 // ========================================================================= //
 // 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/gemm.h>
 #include <myramath/dense/frobenius.h>

 // Iterative solve/action stuff.
 #include <myramath/iterative/Action.h>
 #include <myramath/iterative/ProductAction.h>
 #include <myramath/iterative/GemmAction.h>

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

 using namespace myra;

 namespace {

 template<class Number> void test(typename ReflectPrecision<Number>::type tolerance)
   {
   typedef typename ReflectPrecision<Number>::type Precision;
   auto A = Matrix<Number>::random(3,10);
   auto B = Matrix<Number>::random(10,3);
   auto C1 = A*B;
   auto C2 = make_GemmAction(A)*make_GemmAction(B);
   Precision error = frobenius(C1-C2.make_Matrix());
   myra::out() << "|(A*B)[dense] - (A*B)[action]| = " << error << std::endl;
   REQUIRE(error < tolerance);
   }

 } // namespace

 ADD_TEST("ProductAction","[iterative]")
   {
   test<NumberS>(1.0e-4f);
   test<NumberD>(1.0e-12);
   test<NumberC>(1.0e-4f);
   test<NumberZ>(1.0e-12);
   }
MatrixRange.h
Interface class for representing subranges of dense Matrix&#39;s.

Action.h
Applies the "Action" of a linear operator, b := A*x, used in iterative solution algorithms.

frobenius.h
Routines for computing Frobenius norms of various algebraic containers.

myra::Matrix::random
static Matrix< Number > random(int I, int J)
Generates a random Matrix of specified size.
Definition: Matrix.cpp:353

myra
Definition: syntax.dox:1

Number.h
Various utility functions/classes related to scalar Number types.

ProductAction.h
Composes two Action&#39;s A and B, yielding an Action that cascades A*(B*X)

Matrix.h
General purpose dense matrix container, O(i*j) storage.

myra::ReflectPrecision
Reflects Precision trait for a Number, scalar Number types should specialize it.
Definition: Number.h:33

GemmAction.h
An Action for multiplying by a dense Matrix or SparseMatrix using gemm().

gemm.h
Variety of routines all for dense Matrix*Matrix multiplies. Delegates to the BLAS.