Wiki

Clone wiki

blaze / Addition


Vector/Vector Addition

The addition of vectors is as intuitive as the addition of scalar values. For the addition of any two vectors the addition operator (i.e. operator+()) can be used. It even enables the addition of dense and sparse vectors:

blaze::DynamicVector<int>      v1( 5UL ), v3;
blaze::CompressedVector<float> v2( 5UL );

// ... Initializing the vectors

v3 = v1 + v2;  // Addition of a dense and a sparse column vector of different data type

Note that it is necessary that both operands have exactly the same dimensions. Violating this precondition results in an exception. Also note that it is only possible to add vectors with the same transpose flag:

using blaze::columnVector;
using blaze::rowVector;

blaze::DynamicVector<int,columnVector>   v1( 5UL );
blaze::CompressedVector<float,rowVector> v2( 5UL );

v1 + v2;           // Compilation error: Cannot add a column vector and a row vector
v1 + trans( v2 );  // OK: Addition of two column vectors

Also note that the addition of two vectors with the same element type is favorable due to possible vectorization of the operation:

blaze::DynamicVector<double> v1( 100UL ), v2( 100UL ), v3;

// ... Initialization of the vectors

v3 = v1 + v2;  // Vectorized addition of two double precision vectors

Outer Sum

The addition between a column vector and a row vector results in the outer sum of the two vectors:

blaze::StaticVector<int,3UL,columnVector> v1{ 2, 5, -1 };
blaze::DynamicVector<int,rowVector> v2{ -1, 3, -2, 4 };

// Results in the matrix
//
//       (  1  5  0  6 )
//   A = (  4  8  3  9 )
//       ( -2  2 -3  3 )
//
blaze::StaticMatrix<int,3UL,4UL> M1 = v1 + v2;

The trans() function can be used to transpose a vector as necessary:

blaze::StaticVector<int,2UL,rowVector> v1{ 2, 5, -1 };
blaze::DynamicVector<int,rowVector> v2{ -1, 3, -2, 4 };

blaze::StaticMatrix<int,3UL,4UL> M1 = trans( v1 ) + v2;

Matrix/Matrix Addition

For the addition of any two matrices the addition operator (i.e. operator+()) can be used. It even enables the addition of dense and sparse matrices:

using blaze::rowMajor;
using blaze::columnMajor;

blaze::CompressedMatrix<size_t,columnMajor> M1( 7UL, 3UL );
blaze::DynamicMatrix<float,rowMajor>        M2( 7UL, 3UL ), M3;

// ... Initializing the matrices

M3 = M1 + M2;  // Addition of a sparse column-major and a dense row-major matrix of different data type

Note that it is necessary that both operands have exactly the same dimensions. Violating this precondition results in an exception. It is possible to add row-major and column-major matrices. Note however that in favor of performance the addition of two matrices with the same storage order is favorable. The same argument holds for the element type: In case two matrices with the same element type are added, the performance can be much higher due to vectorization of the operation.

!#c++
blaze::DynamicMatrix<float> M1( 50UL, 70UL ), M2( 50UL, 70UL ), M3;

// ... Initialization of the matrices

M3 = M1 + M2;  // Vectorized addition of two row-major, single precision dense matrices

Scalar Addition

For convenience it is also possible to add a scalar value to a dense vector or dense matrix, which has the same effect as adding a uniform vector or matrix. In Blaze it is possible to use all built-in/fundamental data types except bool as scalar values. Additionally, it is possible to use std::complex values with the same built-in data types as element type. Examples:

blaze::StaticVector<int,3UL> v1{ 3, 2, 5, -4, 1, 6 };

blaze::DynamicVector<int>    v2 = v1 + 2;  // Results in { 5, 4, 7, -2, 3, 8 }
blaze::CompressedVector<int> v3 = 3 + v1;  // Results in { 6, 5, 8, -1, 4, 9 }
blaze::StaticMatrix<int,2UL,3UL> M1{ {  3, 2, 5 },
                                     { -4, 1, 6 } };

blaze::DynamicMatrix<int>    M2 = M1 + 2;  // Results in { { 5, 4, 7 }, { -2, 3, 8 } }
blaze::CompressedMatrix<int> M3 = 3 + M1;  // Results in { { 6, 5, 8 }, { -1, 4, 9 } }

Previous: Arithmetic Operations ---- Next: Subtraction

Updated