Sparse Matrices and Vectors
Enhancements have been made in Maple 17 for the performance of stacking and augmenting sparse floating-point, hardware datatype Matrices and Vectors. Computing linear combinations of pairs of sparse floating-point, hardware datatype Matrices using the LinearAlgebra:-MatrixAdd command has also been improved.
In the examples below, the Matrix and Vector constructors are used for stacking or augmentation (concatenation).
Benchmarks computed in 64-bit Maple 17 for Linux on an Intel Core i5 CPU 760 @ 2.80GHz.
Sparse Vector concatenation
Sparse Matrix concatenation
Simple linear combinations of sparse Matrices
See Also
In Maple 16, the concatenation of Vector V with itself three times would take approximately 28 seconds and would result in an output Vector with full rectangular storage. The allocated memory would increase by 1.27GiB.
In Maple 17, it takes less than 0.2 seconds and produces a Vector with sparse storage. Memory allocation increases by 30.5MiB.
Example
restart:
withLinearAlgebra:
V:=RandomVector107,generator=0.0..1.0,density=0.05,storage=sparse,datatype=float8
rtable_num_elemsV,NonZeroStored
500366
biggerV:=CodeTools:-UsageVectorV,V,V,V
memory used=76.37MiB, alloc change=30.54MiB, cpu time=100.00ms, real time=121.00ms
rtable_num_elemsbiggerV,NonZeroStored
2001464
In Maple 16, the concatenation of Matrix M with itself as performed below would take approximately 20 seconds and would result in an output Matrix with full rectangular storage.
In Maple 17, it takes less than 0.1 seconds and produces a Matrix with sparse storage. Memory allocation increases by 48MiB.
restart:withLinearAlgebra:M:=RandomMatrix103,103,generator=0.0..1.0,density=0.01,storage=sparse,datatype=float8
rtable_num_elemsM,NonZeroStored
10053
biggerM:=CodeTools:-UsageMatrixM,M
memory used=0.76MiB, alloc change=48.00MiB, cpu time=30.00ms, real time=55.00ms
rtable_num_elemsbiggerM,NonZeroStored
20106
In Maple 16, the following command would take 3 minutes and allocated memory would increase by 400MiB.
In Maple 17, it takes less than 1 second and allocated memory increases by 76MiB.
restart:withLinearAlgebra:M:=RandomMatrix104,104,generator=0.0..1.0,density=0.01,storage=sparse,datatype=float8;CodeTools:-UsageMatrixAddM,M,2,−3.3
memory used=76.67MiB, alloc change=76.35MiB, cpu time=360.00ms, real time=417.00ms
Matrix, Vector
Download Help Document
