Abstract:We observe a disconnect between the developers and the end users of linear algebra libraries. On the one hand, the numerical linear algebra and the high-performance communities invest significant effort in the development and optimization of highly sophisticated numerical kernels and libraries, aiming at the maximum exploitation of both the properties of the input matrices, and the architectural features of the target computing platform. On the other hand, end users are progressively less likely to go through the error-prone and time consuming process of directly using said libraries by writing their code in C or Fortran; instead, languages and libraries such as Matlab, Julia, Eigen and Armadillo, which offer a higher level of abstraction, are becoming more and more popular. Users are given the opportunity to code matrix computations with a syntax that closely resembles the mathematical description; it is then a compiler or an interpreter that internally maps the input program to lower level kernels, as provided by libraries such as BLAS and LAPACK. Unfortunately, our experience suggests that in terms of performance, this translation is typically vastly suboptimal. In this paper, we first introduce the Linear Algebra Mapping Problem, and then investigate how effectively a benchmark of test problems is solved by popular high-level programming languages. Specifically, we consider Matlab, Octave, Julia, R, Armadillo (C++), Eigen (C++), and NumPy (Python); the benchmark is meant to test both standard compiler optimizations such as common subexpression elimination and loop-invariant code motion, as well as linear algebra specific optimizations such as optimal parenthesization of a matrix product and kernel selection for matrices with properties. The aim of this study is to give concrete guidelines for the development of languages and libraries that support linear algebra computations.

The Challenges of Teaching Elementary Linear Algebra in a Modern Matrix Based Way

Teaching Linear Algebra in a Mechanized Mathematical Environment

Teaching linear algebra at university

Designing Problems for Improved Instruction and Learning -- Linear Algebra

Computer aided teaching in Linear Algebra

Research on Group Cooperative Learning in Linear Algebra Teaching

Moving Away from Lecture in Undergraduate Mathematics: Managing Tensions within a Coordinated Inquiry-Based Linear Algebra Course

Teaching Linear Algebra in Engineering Courses Using Critical Thinking

On the Teaching of the Linear Algebra for Engineering Students

Perspectives and reflections on teaching linear algebra

On Teaching Reform of Linear Algebra

Flipped classroom in Introductory Linear Algebra by utilizing Computer Algebra System {\sl SageMath} and a free electronic book

Pedagogical Strategies to Enhance Learning in a Linear Algebra Course

Assessing Confidence and Certainty of Students in an Undergraduate Linear Algebra Course

Engaging activities for teaching linear algebra

Applications of Matrix Theories in the Review Teaching of Linear Algebra Course

The Linear Algebra Mapping Problem. Current state of linear algebra languages and libraries

Mathematical analysis in high school : a fundamental dilemma

Traditional lectures versus active learning -- a false dichotomy?

The Discussion of Starting Experimental Courses in the Teaching of Linear Algrbra

Teaching Writing and Communication in a Mathematical Modeling Course