Skip to content

Fastor V0.6.1
Compare
Choose a tag to compare
@romeric romeric released this 07 May 21:26
· 223 commits to master since this release
V0.6.1
5654474

Fastor V0.6.1 is an incremental change over V0.6 release that introduced a significant overhaul in Fastor's internal design and exposed API. This release includes

  1. lu function introduced for LU decomposition of 2D tensors. Multiple variants of LU decomposition is available including no pivoting, partial pivoting with a permutation vector and partial pivoting with a permutation matrix. This is perhaps the most performant implementation of the LU decomposition available today for small matrices of up to 64x64. If no pivoting is used it the performance is unbeaten for all sizes up to the stack limit however given that the implementation is based on compile time loop recursion for sizes up to 32x32 and further it uses block recursion which in turn uses block-triangular-inversion compilation would be quite time consuming for bigger sizes
  2. ut_inverse and lut_inverse for fast triangular inversion of upper and unit lower matrices using block-wise inversion
  3. tmatmul function equivalent to BLAS's TRMM function for triangular matrix-matrix (or vector) multiplication which allows either or both operand to be upper/lower triangular. The function can be used to specifiy which matrix is lower/upper at compile time like tmatmul<matrix_type::lower_tri,matrix_type::general>(A,B). Proper 2X speed up over matmul for when one operand is triangular and 4X when both are triangular can be achieved for bigger sizes
  4. det/determinant can now be computed for all sizes using the LU decomposition [default for matrix sizes bigger than 4x4].
  5. inv/inverse and solve can be performed with any variant of the LU decomposition
  6. There is now a unified interface for choosing the computation type of linear algebra functions for instance det<DetCompType::BlockLU>(A) or inv<InvCompType::SimpleLUPiv>(A) or solve<SolveCompType::BlockLUPiv> etc
  7. tril/triu functions added for getting the lower/upper part of a 2D tensor
  8. Comprehensive unit tests and benchmarks are added and are available for these newly added (and some old) routines
Assets 2
Loading