Abstract:When writing high-performance code for numerical computation in a scripting language like MATLAB, it is crucial to have the operations in a large for-loop vectorized. If not, the code becomes too slow to use, even for a moderately large problem. However, in the process of vectorizing, the code often loses its original structure and becomes less readable. This is particularly true in the case of a finite element implementation, even though finite element methods are inherently structured. A basic remedy to this is the separation of the vectorization part from the mathematics part of the code, which is easily achieved through building the code on top of the basic linear algebra subprograms that are already vectorized codes, an idea that has been used in a series of papers over the last fifteen years, developing codes that are fast and still structured and readable. We discuss the vectorized basic linear algebra package and introduce a formalism using multi-linear algebra to explain and define formally the functions in the package, as well as MATLAB pagetime functions. We provide examples from computations of varying complexity, including the computation of normal vectors, volumes, and finite element methods. Benchmarking shows that we also get fast computations. Using the library, we can write codes that closely follow our mathematical thinking, making writing, following, reusing, and extending the code easier.

RandFile package for Mathematica for accessing file-based sources of randomness

Recent and Upcoming Developments in Randomized Numerical Linear Algebra for Machine Learning

An Introduction to Mrandom 3.0

OpenRAND: A Performance Portable, Reproducible Random Number Generation Library for Parallel Computations

Randentropy: a software to measure inequality in random systems

Randen - fast backtracking-resistant random generator with AES+Feistel+Reverie

Exhibiting Randomness in Arithmetic using Mathematica and C

Probability package for Macaulay2

A computer algebra system for R: Macaulay2 and the m2r package

Determinantal Point Processes in Randomized Numerical Linear Algebra

rlsm: R package for least squares Monte Carlo

Sketch 'n Solve: An Efficient Python Package for Large-Scale Least Squares Using Randomized Numerical Linear Algebra

Geometric Random Inner Products: A New Family of Tests for Random Number Generators

An R Implementation of the Polya-Aeppli Distribution

On a vectorized basic linear algebra package for prototyping codes in MATLAB

A Novel Design of Random Number Generators Using Chaos-Based Extremum Coding

MINT: a Computer Program for Adaptive Monte Carlo Integration and Generation of Unweighted Distributions

Nonparametric Estimation of the Random Coefficients Model in Python

Software/hardware framework for generating parallel Gaussian random numbers based on the Monty Python method

Hardware Random number Generator for cryptography

A Class of Random Number Generators Based on Weyl Sequence