## BLAS

The *Basic Linear Algebra Subprograms* (BLAS) is a common library
including high-quality routines for performing basic operations on vectors and
matrices. The first FORTRAN version was released in 1979. As the routines are
efficient and portable, they are often used in linear algebra software, like
LAPACK. Most implementations are written in C, C++, or
FORTRAN 77.

- ATLAS
- High-performance implementation of the BLAS API, for C and FORTRAN 77 (BSD licence).
- BLAS
- The official reference implementation by Netlib (C and FORTRAN 77). Also
known as
*refblas*(public domain). - IMKL
- The Intel Math Kernel Library includes BLAS for Intel processors. Professional and community editions for C, C++, and Fortran are distributed (proprietary).

Please see Wikipedia for a more comprehensive overview of BLAS libraries.

### Operations

The functionality of BLAS is divided into three sets of routines, called
*levels*. Below, *A*, *B*, *C*, and *T* are
matrices, *x* and *y* are vectors, and *α* and *β*
are scalars.

- Level 1: Scalar-Vector and Vector-Vector Operations
The level provides low-level operations, like dot product of vectors and vector additions, such as

*y*←*αx*+*y*.- Level 2: Matrix-Vector Operations
The level contains basic matrix-vector operations, such as

*y*←*αAx*+*βy*.It also includes a solver for

*x*in the triangular equation*Tx*=*y*,with triangular matrix

*T*.- Level 3: Matrix-Matrix Operations
The level includes matrix-matrix operations, such as

*C*←*αAB*+*βC,*and routines for solving

*B*←*αT*^{−1}*B*,with triangular matrix

*T*.

### Installation

On FreeBSD, the LAPACK implementation of BLAS is available as a port:

`# pkg install math/blas`

However, BLAS may be built also from source as well. Download version 3.8.0
from Netlib, unpack the archive, and
compile the Fortran code manually. In order to create a shared library
`libblas.so`

, run:

`$ gfortran10 -O2 -shared -fPIC -o libblas.so *.f`

Or, if a static library is preferred, compile and archive
`libblas.a`

with:

```
$ gfortran10 -O2 -c *.f
$ ar cr libblas.a *.o
```

### Example

The following example program scales a vector by a constant using the BLAS
routine `sscal()`

. FORTRAN 77 procedures have to be be imported
individually using the `external`

statement if
`implicit none (type, external)`

is set. Otherwise, all linked
procedures are visible by default.

```
! example.f90
program main
implicit none (type, external)
external :: sscal
integer, parameter :: N = 3
real :: x(N) = [ 5., 6., 7. ]
real :: a = 5.
print '("a = ", f0.1)' a
print '("X = [ ", 3(f0.1, " "), "]")', x
call sscal(N, a, x, 1)
print '(/, "X = a * X")'
print '("X = [ ", 3(f0.1, " "), "]")', x
end program main
```

Compile and link the example with:

`$ gfortran10 -L/usr/local/lib/ -o example example.f90 -lblas`

If you prefer to link against the static library `libblas.a`

,
replace `-lblas`

with the actual path. The program outputs:

```
$ ./example
a = 5.0
X = [ 5.0 6.0 7.0 ]
X = a * X
X = [ 25.0 30.0 35.0 ]
```

The compilation can be simplified by writing an appropriate Makefile:

```
.POSIX:
FC = gfortran10
FFLAGS = -std=f2018
LDFLAGS =
LDLIBS = -lblas
TARGET = example
.PHONY: all clean
all:
$(FC) $(FFLAGS) $(LDFLAGS) -o $(TARGET) example.f90 $(LDLIBS)
clean:
rm $(TARGET)
```

We then just have to run `make`

to build the executable.

< POSIX Threads | [Index] | LAPACK > |